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)
Aims. Sunspot distribution in the northern and southern solar hemispheres exibit striking synchronous behaviour on the scale of a Schwabe cycle. However, sometimes the bilateral symmetry of the Butterfly diagram relative to the solar equatorial plane breaks down. The investigation of this phenomenon is important to explaining the almost-periodic behaviour of solar cycles. Methods. We use cross-recurrence plots for the study of the time-varying phase asymmetry of the northern and southern hemisphere and compare our results with the latitudinal distribution of the sunspots. Results. We observe a long-term persistence of phase leading in one of the hemispheres, which lasts almost 4 solar cycles and probably corresponds to the Gleissberg cycle. Long-term variations in the hemispheric-leading do not demonstrate clear periodicity but are strongly anti-correlated with the long-term variations in the magnetic equator.
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.
Sensory information entering the nervous system follows independent paths of processing such that specific features are individually detected. However, sensory perception, awareness, and cognition emerge from the combination of information. Here we have analyzed the corticocortical network of the cat, looking for the anatomical substrate which permits the simultaneous segregation and integration of information in the brain. We find that cortical communications are mainly governed by three topological factors of the underlying network: (i) a large density of connections, (ii) segregation of cortical areas into clusters, and (iii) the presence of highly connected hubs aiding the multisensory processing and integration. Statistical analysis of the shortest paths reveals that, while information is highly accessible to all cortical areas, the complexity of cortical information processing may arise from the rich and intricate alternative paths in which areas can influence each other.
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.
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.
The response of scale-free networks with community structure to external stimuli is studied. By disturbing some nodes with different strategies, it is shown that the robustness of this kind of network can be enhanced due to the existence of communities in the networks. Some of the response patterns are found to coincide with topological communities. We show that such phenomena also occur in the cat brain network which is an example of a scale-free like network with community structure. Our results provide insights into the relationship between network topology and the functional organization in complex networks from another viewpoint.
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.
We have used techniques of nonlinear dynamics to compare a special model for the reversals of the Earth's magnetic field with the observational data. Although this model is rather simple, there is no essential difference to the data by means of well-known characteristics, such as correlation function and probability distribution. Applying methods of symbolic dynamics we have found that the considered model is not able to describe the dynamical properties of the observed process. These significant differences are expressed by algorithmic complexity and Renyi information.
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.
The incidence of cardiovascular diseases increases with the growth of the human population and an aging society, leading to very high expenses in the public health system. Therefore, it is challenging to develop sophisticated methods in order to improve medical diagnostics. The question whether the normal heart rate is chaotic or not is an attempt to elucidate the underlying mechanisms of cardiovascular dynamics and therefore a highly controversial topical challenge. In this contribution we demonstrate that linear and nonlinear parameters allow us to separate completely the data sets of the three groups provided for this controversial topic in nonlinear dynamics. The question whether these time series are chaotic or not cannot be answered satisfactorily without investigating the underlying mechanisms leading to them. We give an example of the dominant influence of respiration on heart beat dynamics, which shows that observed fluctuations can be mostly explained by respiratory modulations of heart rate and blood pressure (coefficient of determination: 96%). Therefore, we recommend reformulating the following initial question: "Is the normal heart rate chaotic?" We rather ask the following: " Is the normal heart rate 'chaotic' due to respiration?"
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
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.
Using a special technique of data analysis, we have found out 34 grand minima of solar activity in a 7,700 years long C14 record. The method used rests on a proper filtering of the C14 record and the extrapolation of verifiable results for the later history back in time. Additionally, we have applied a method of nonlinear dynamics, the recurrence rate, to back up the results. Our findings are not contradictory to the record of grand minima by Eddy, but constitute a considerable extension. Hence, it has become possible to look closer at the validity of models. This way, we have tested esp. the model of Barnes et al. There are hints for that the grand minima might solely be driven by the 209--year period found in the C14 record.
Using a special technique of data analysis, we have found out 34 grand minima of solar activity obtained from a 7,700 years long Δ14C record. The method used rests on a proper filtering of the Δ14C record and the extrapolation of verifiable results for the later history back in time. Additionally, we use a method of nonlinear dynamics, the recurrence rate, to back up the results. Our findings are not contradictory to the record of solar maxima resp. minima by Eddy [5], but constitute a considerable extension. Hence, it has become possible to look closer at the validity of models. This way, we have tested several models for solar activity, esp. the model of Barnes et al. [1]. There are hints for that the grand minima might solely be driven by the 209 year period found in the Δ14C record.
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.
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.
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.
The Voyager 2 Photopolarimeter experiment has yielded the highest resolved data of Saturn's rings, exhibiting a wide variety of features. The B-ring region between 105000 km and 110000 km distance from Saturn has been investigated. It has a high matter density and contains no significance features visible by eye. Analysis with statistical methods has let us to the detection of two significant events. These features are correlated with the inner 3:2 resonances of the F-ring shepherd satellites Pandora and Prometheus, and may be evidence of large ring paricles caught in the corotation resonances.
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.
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.
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 have numerically studied the bifurcation properties of a sheet pinch with impenetrable stress-free boundaries. An incompressible, electrically conducting fluid with spatially and temporally uniform kinematic viscosity and magnetic diffusivity is confined between planes at x1=0 and 1. Periodic boundary conditions are assumed in the x2 and x3 directions and the magnetofluid is driven by an electric field in the x3 direction, prescribed on the boundary planes. There is a stationary basic state with the fluid at rest and a uniform current J=(0,0,J3). Surprisingly, this basic state proves to be stable and apparently to be the only time-asymptotic state, no matter how strong the applied electric field and irrespective of the other control parameters of the system, namely, the magnetic Prandtl number, the spatial periods L2 and L3 in the x2 and x3 directions, and the mean values B¯2 and B¯3 of the magnetic-field components in these directions.
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 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.
It is shown that the ff effect of mean-field magnetohydrodynamics, which consists in the generation of a mean electromotive force along the mean magnetic field by turbulently fluctuating parts of velocity and magnetic field, is equivalent to the simultaneous generation of both turbulent and mean-field magnetic helicities, the generation rates being equal in magnitude and opposite in sign. In the particular case of statistically stationary and homogeneous fluctuations this implies that the ff effect can increase the energy in the mean magnetic field only under the condition that also magnetic helicity is accumulated there.
Using quantities of symbolic dynamics, such as mutual information, Shannon information and algorithmic complexity, we have searched for interrelations of spikes emitted simultaneously at different frequencies during the impulsive phase of a flare event. As the spikes are related to the flare energy release and are interpreted as emissions originating at different sites having different magnetic field strengths, any relation in frequency is interpretated as a relation in space. This approach is appropriate to characterize such spatio-temporal patterns, whereas the popular estimate of fractal dimensions can be applied to low-dimensional systems only. Depending on the energy release and emission processes, two types of fragmentation are possible: a scenario of global organization (spikes are emitted in a succession of similar events by the same system) or a scenario of local organization (many systems triggered by an initial event). Mutual information which is a generalization of correlation indicates a relation in frequency beyond the bandwidth of individual spikes. The scans in the spectrograms with large mutual information also show a low level of Shannon information and algorithmic complexity, indicating that the simultaneous appearance of spikes at other frequencies is not a completely stochastic phenomenon (white noise). It may be caused by a nonlinear deterministic system or by a Markov process. By means of mutual information we find a memory over frequency intervals up to 60 MHz. Shannon information and algorithmic complexity concern the mbox{whole} frequency region, i.e. the global source region. A global organization is also apparent in quasi-periodic changes of the Shannon information and algorithmic complexity in the range of 2 - 8 seconds. The finding is compatible with a scenario of local organization in which the information of one event spreads spatially and triggers further events at different places. The region is not an ensemble of independently flashing sources, each representing a system that cascades in energy after an initial trigger. On the contrary, there is a causal connection between the sources at any time. The analysis of the four spike events suggests that the structure in frequency is not stochastic but a process in which spikes at nearby locations are simultaneously triggered by a common exciter.
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
Im vorletzten Absatz des o.g. Kurzberichtes befindet sich eine falsche Aussage zur C14-Produktion waehrend des Maunder-Minimums. Wie aus der in meiner Abbildung gezeigten Delta C14-Haeufigkeit fuer den Zeitraum des Maunder-Minimums hervorgeht, war die C14-Produktion zu dieser Zeit erhoeht statt, wie von Herrn Buehrke und anderen Autoren in der Literatur behauptet, erniedrigt. Die allgemein akzeptierte Begruendung fuer die erhoehte C14-Produktion lautet: Der geringere Sonnenwind schirmt die Erde weniger stark von der kosmischen Strahlung ab.
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.
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.
The present paper is related to the problem of approximating the exact solution to the magnetohydrodynamic equations (MHD). The behaviour of a viscous, incompressible and resistive fluid is exemined for a long period of time. Contents: 1 The magnetohydrodynamic equations 2 Notations and precise functional setting of the problem 3 Existence, uniqueness and regularity results 4 Statement and Proof of the main theorem 5 The approximate inertial manifold 6 Summary
The EEG is one of the most commonly used tools in brain research. Though of high relevance in research, the data obtained is very noisy and nonstationary. In the present article we investigate the applicability of a nonlinear data analysis method, the recurrence quantification analysis (RQA), to Such data. The method solely rests on the natural property of recurrence which is a phenomenon inherent to complex systems, such as the brain. We show that this method is indeed suitable for the analysis of EEG data and that it might improve contemporary EEG analysis.
In the recent past, recurrence quantification analysis (RQA) has gained an increasing interest in various research areas. The complexity measures the RQA provides have been useful in describing and analysing a broad range of data. It is known to be rather robust to noise and nonstationarities. Yet, one key question in empirical research concerns the confidence bounds of measured data. In the present Letter we suggest a method for estimating the confidence bounds of recurrence-based complexity measures. We study the applicability of the suggested method with model and real- life data.
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 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.
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.
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.
The investigation of foetal reaction to internal and external conditions and stimuli is an important tool in the characterization of the developing neural integration of the foetus. An interesting example of this is the study of the interrelationship between the foetal and the maternal heart rate. Recent studies have shown a certain likelihood of occasional heart rate synchronization between mother and foetus. In the case of respiratory-induced heart rate changes, the comparison with maternal surrogates suggests that the evidence for detected synchronization is largely statistical and does not result from physiological interaction. Rather, they simply reflect a stochastic, temporary stability of two independent oscillators with time-variant frequencies. We reanalysed three datasets from that study for a more local consideration. Epochs of assumed synchronization associated with short-term regulation of the foetal heart rate were selected and compared with synchronization resulting from white noise instead of the foetal signal. Using data-driven modelling analysis, it was possible to identify the consistent influence of the heartbeat duration of maternal beats preceding the foetal beats during epochs of synchronization. These maternal beats occurred approximately one maternal respiratory cycle prior to the affected foetal beat. A similar effect could not be found in the epochs without synchronization. Simulations based on the fitted models led to a higher likelihood of synchronization in the data segments with assumed foetal-maternal interaction than in the segment without such assumed interaction. We conclude that the data-driven model-based analysis can be a useful tool for the identification of synchronization.
The analysis of baroreflex sensitivity (BRS) and heart rate variability (HRV) leads to additional insights into patients' prognosis after cardiovascular events. The following study was performed to assess the differences in the post-operative recovery of autonomic regulation after mitral valve (MV) and aortic valve (AV) surgery with a heart lung machine. Among the 43 consecutive male patients enrolled in a prospective study, 26 underwent isolated AV surgery and 17 isolated MV surgery. Blood pressure as well as ECG signals were recorded the day before, 24 hours after and one week after surgery. BRS was calculated according to the dual sequence method, and HRV was calculated using standard linear as well as nonlinear parameters. There were no major differences between the two groups in the pre-operative values. At 24 hours a comparable depression of HRV and BRS in both groups was observed, while at 7 days there was partial recovery in AV patients, which was absent in MV patients: p(AV versus MV) < 0.001. While the response of the autonomic system to surgery is similar in AV and MV patients, there is obviously a decreased ability to recover in MV patients, probably attributed to traumatic lesions of the autonomic nervous system by opening the atria. Ongoing research is required for further clarification of the pathophysiology of this phenomenon and to establish strategies to restore autonomic function.
We demonstrate the occurrence of regimes with singular continuous (fractal) Fourier spectra in autonomous dissipative dynamical systems. The particular example in an ODE system at the accumulation points of bifurcation sequences associated to the creation of complicated homoclinic orbits. Two different machanisms responsible for the appearance of such spectra are proposed. In the first case when the geometry of the attractor is symbolically represented by the Thue-Morse sequence, both the continuous-time process and its descrete Poincaré map have singular power spectra. The other mechanism owes to the logarithmic divergence of the first return times near the saddle point; here the Poincaré map possesses the discrete spectrum, while the continuous-time process displays the singular one. A method is presented for computing the multifractal characteristics of the singular continuous spectra with the help of the usual Fourier analysis technique.
Strange nonchaotic attractors typically appear in quasiperiodically driven nonlinear systems. Two methods of their characterization are proposed. The first one is based on the bifurcation analysis of the systems, resulting from periodic approximations of the quasiperiodic forcing. Secondly, we propose th characterize their strangeness by calculating a phase sensitivity exponent, that measures the sensitivity with respect to changes of the phase of the external force. It is shown, that phase sensitivity appears if there is a non-zero probability for positive local Lyapunov exponents to occur.
Analysis of blood pressure dynamics in male and female rats using the continuous wavelet transform
(2009)
We study gender-related particularities in cardiovascular responses to stress and nitric oxide (NO) deficiency in rats using HR, mean arterial pressure (MAP) and a proposed wavelet-based approach. Blood pressure dynamics is analyzed: (1) under control conditions, (2) during immobilization stress and recovery and (3) during nitric oxide blockade by N-G-nitro-L-arginine-methyl ester (L-NAME). We show that cardiovascular sensitivity to stress and NO deficiency depends upon gender. Actually, in females the chronotropic effect of stress is more pronounced, while the pressor effect is weakened compared with males. We conclude that females demonstrate more favorable patterns of cardiovascular responses to stress and more effective NO control of cardiovascular activity than males.
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.
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.
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.
In the recent article "Stochastic analysis of recurrence plots with applications to the detection of deterministic signals" (Physica D 237 (2008) 619-629), Rohde et al. stated that the performance of RQA in order to detect deterministic signals would be below traditional and well-known detectors. However, we have concerns about such a general statement. Based on our own studies we cannot confirm their conclusions. Our findings suggest that the measures of complexity provided by RQA are useful detectors outperforming well-known traditional detectors, in particular for the detection of signals of complex systems, with phase differences or signals modified due to the measurement process.
Contents: 1 Introduction 1.1 Tikhanov-Phillips Regularization of Ill-Posed Problems 1.2 A Compact Course to Wavelets 2 A Multilevel Iteration for Tikhonov-Phillips Regularization 2.1 Multilevel Splitting 2.2 The Multilevel Iteration 2.3 Multilevel Approach to Cone Beam Reconstuction 3 The use of approximating operators 3.1 Computing approximating families {Ah}
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.
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.
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
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.
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.
In the modern industrialized countries every year several hundred thousands of people die due to the sudden cardiac death. The individual risk for this sudden cardiac death cannot be defined precisely by common available, non-invasive diagnostic tools like Holter-monitoring, highly amplified ECG and traditional linear analysis of heart rate variability (HRV). Therefore, we apply some rather unconventional methods of nonlinear dynamics to analyse the HRV. Especially, some complexity measures that are basing on symbolic dynamics as well as a new measure, the renormalized entropy, detect some abnormalities in the HRV of several patients who have been classified in the low risk group by traditional methods. A combination of these complexity measures with the parameters in the frequency domain seems to be a promising way to get a more precise definition of the individual risk. These findings have to be validated by a representative number of patients.
We have discussed some tools from nonlinear dynamics which may help to analyze transient phenomena, such as solar bursts. The structure function known from turbulence theory is an appropriate method to find out some scaling behavior of fluctuations in time. More generally, the wavelet analysis, which is some generalization of the power spectrum, exhibits information on the location as well as the size of hidden characteristic features. Applying both techniques to microwave bursts, we have found some scaling properties that refer to the existence of hierarchic time structures. This is in good accordance with the electric circuit model for describing the flare-particle energization process.
The radiocarbon record that has been extended from 7199 BC to 1891 AD is of fundamental importance to understand century-scale variations of solar activity. We have, therefore, studied how to extract information from dynamic reconstructions of this observational record. Using some rather unusual methods of nonlinear dynamics, we have found that the data are significantly different from linear colored noise and that there is some evidence of nonlinear behavior. The method of recurrence plots exhibits that the grand minima of solar activity are quite different in their recurrence. Most remarkably, it suggests that the recent epoch seems to be similar to the Medieval maximum.
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.
We look for structural properties in the light curve of the dwarf nova SS Cyg by means of techniques from nonlinear dynamics. Applying the popular Grassberger-Procaccia procedure, Cannizzo and Goddings (1988) showed that there is no evidence for a low-dimensional attractor underlying this record. Because there are some hints for order in the light curve, we search for other signatures of deterministic systems. Therefore, we use other methods recently developed in this theory, such as local linear prediction and recurrence maps. Our main findings are: i] the prediction error grows exponentially during outburst phases, but via a power law in the quiescent states, ii] there are some rather regular patterns in this light curve which sometimes recur, but the recurrence is not regular. This leads to the following conclusions: i] The outburst dynamics shows a higher degree of order than the quiescent one. There are some hints for deterministic chaos in the outburst behavior. ii] The light curve is a complex mixture of deterministic and stochastic structures. The analysis presented in this paper shows that methods of nonlinear dynamics can be an efficient tool for the study of complex processes, even if there is no evidence for a low-dimensional attractor.
The application of chaos theory has become popular to understand the nature of various features of solar activity because most of them are far from regular. The usual approach, however, that is basing on finding low- dimensional structures of the underlying processes seems to be successful only in a few exceptional cases, such as in rather coherent phenomena as coronal pulsations. It is important to note that most phenomena in solar radio emission are more complex. We present two kinds of techniques from nonlinear dynamics which can be useful to analyse such phenomena: i] Fragmentation processes observed in solar spike events are studied by means of symbolic dynamics methods. Different measures of complexity calculated from such observations reveal that there is some order in this fragmentation. ii] Bursts are a typical transient phenomenon. To study energization processes causing impulsive microwave bursts, the wavelet analysis is applied. It exhibits structural differences of the pre- and post-impulsive phase in cases where the power spectra of both are not distinct.
Two deterministic processes leading to roughening interfaces are considered. It is shown that the dynamics of linear perturbations of turbulent regimes in coupled map lattices is governed by a discrete version of the Kardar-Parisi-Zhang equation. The asymptotic scaling behavior of the perturbation field is investigated in the case of large lattices. Secondly, the dynamics of an order-disorder interface is modelled with a simple two-dimensional coupled map lattice, possesing a turbulent and a laminar state. It is demonstrated, that in some range of parameters the spreading of the turbulent state is accompanied by kinetic roughening of the interface.
The paper presents a method that determines, by standard numerical means, the type of mutual relations of fold and flip bifurcations (configured as a so-called communication area) of a map. Equation systems are developed for the computation of points where a transition between areas of different types occurs. Furthermore, it is shown that saddle area<->spring area transitions can exist which have not yet been considered in the literature. Analytical conditions of that transition are derived.
Contents: I. Algorithms 1. Theoretical Backround 2. Numerical Procedures 3. Graph Representation of the Solutions 4. Applications and Example II. Users' Manual 5. About the Program 6. The Course of a Qualitative Analysis 7. The Model Module 8. Input description 9. Output Description 10. Example 11. Graphics
In one of the data mining techniques, change-point detection is of importance in evaluating time series measured in real world. For decades this technique has been developed as a nonlinear dynamics. We apply the method for detecting the change points, Singular Spectrum Transformation (SST), to the climate time series. To know where the structures of climate data sets change can reveal a climate background. In this paper we discuss the structures of precipitation data in Kenya and Wrangel Island (Arctic land) by using the SST.
Objectives: Scoring sleep visually based on polysomnography is an important but time-consuming element of sleep medicine. Where-as computer software assists human experts in the assignment of sleep stages to polysomnogram epochs, their performance is usually insufficient. This study evaluates the possibility to fully automatize sleep staging considering the reliability of the sleep stages available from human expert sleep scorers. Methods: We obtain features from EEG, ECG and respiratory signals of polysomnograms from ten healthy subjects. Using the sleep stages provided by three human experts, we evaluate the performance of linear discriminant analysis on the entire polysomnogram and:only on epochs where the three experts agree in their-sleep stage scoring. Results: We show that in polysomnogram intervals, to which all three scorers assign the same sleep stage, our algorithm achieves 90% accuracy. This high rate of agreement with the human experts is accomplished with only a small set of three frequency features from the EEG. We increase-the performance to 93% by including ECG and respiration features. In contrast, on intervals of ambiguous sleep stage, the sleep stage classification obtained from our algorithm, agrees with the human consensus scorer in approximately 61%. Conclusions: These findings suggest that machine classification is highly consistent with human sleep staging and that error in the algorithm's assignments is rather a problem of lack of well-defined criteria for human experts to judge certain polysomnogram epochs than an insufficiency of computational procedures
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
Recent studies have pointed out the importance of transient synchronization between widely distributed neural assemblies to understand conscious perception. These neural assemblies form intricate networks of neurons and synapses whose detailed map for mammals is still unknown and far from our experimental capabilities. Only in a few cases, for example the C. elegans, we know the complete mapping of the neuronal tissue or its mesoscopic level of description provided by cortical areas. Here we study the process of transient and global synchronization using a simple model of phase-coupled oscillators assigned to cortical areas in the cerebral cat cortex. Our results highlight the impact of the topological connectivity in the developing of synchronization, revealing a transition in the synchronization organization that goes from a modular decentralized coherence to a centralized synchronized regime controlled by a few cortical areas forming a Rich-Club connectivity pattern.
We employ a spectral decomposition method to analyze synchronization of a non-identical oscillator network. We study the case that a small parameter mismatch of oscillators is characterized by one parameter and phase synchronization is observed. We derive a linearized equation for each eigenmode of the coupling matrix. The parameter mismatch is reflected on inhomogeneous term in the linearized equation. We find that the oscillation of each mode is essentially characterized only by the eigenvalue of the coupling matrix with a suitable normalization. We refer to this property as spectral universality, because it is observed irrespective of network topology. Numerical results in various network topologies show good agreement with those based on linearized equation. This universality is also observed in a system driven by additive independent Gaussian noise.
Experimental evidences point Out the participation of nonsynaptic mechanisms (e.g., fluctuations in extracellular tons) in epileptiform bursting and spreading depression (SD). During these abnormal oscillatory patterns, it is observed an increase of extracellular potassium concentration [K+](o) and a decrease of extracellular calcium concentration [Ca2+](o) which raises the neuronal excitability. However, whether the high [K+](o) triggers and propagates these abnormal neuronal activities or plays a secondary role into this process is unclear. To better understand the influence of extracellular potassium dynamics in these oscillatory patterns, the experimental conditions of high [K+](o) and zero [Ca2+](o) were replicated in an extended Golomb model where we added important regulatory mechanisms of ion concentration as Na+-K+ pump, ion diffusion and glial buffering. Within these Conditions, simulations of the cell model exhibit seizure-like discharges (ictal bursting). The SD was elicited by the interruption of the Na+- K+ pump activity, mimicking the effect of cellular hypoxia (an experimental protocol to elicit SD, the hypoxia-induced SD). We used the bifurcation theory and the fast-slow method to analyze the interference of K+ dynamics in the cellular excitability. This analysis indicates that the system loses its stability at a high [K+](o), transiting to an elevated state of neuronal excitability. Effects of high [K+](o), are observed in different stages of ictal bursting and SD. In the initial stage, the increase of [K+](o) creates favorable conditions to trigger both oscillatory patterns. During the neuronal activity, a continuous growth of [K+](o) by outward K+ flow depresses K+ Currents in a positive feedback way. At the last stage, due to the depression of K+ currents, the Na+-K+ pump is the main mechanism in the end of neuronal activity. Thus, this work suggests that [K+](o) dynamics may play a fundamental role in these abnormal oscillatory patterns.
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.
We report on bifurcation studies for the incompressible magnetohydrodynamic equations in three space dimensions with periodic boundary conditions and a temporally constant external forcing. Fourier reprsentations of velocity, pressure and magnetic field have been used to transform the original partial differential equations into systems of ordinary differential equations (ODE), to which then special numerical methods for the qualitative analysis of systems of ODE have been applied, supplemented by the simulative calculation of solutions for selected initial conditions. In a part of the calculations, in order to reduce the number of modes to be retained, the concept of approximate inertial manifolds has been applied. For varying (incereasing from zero) strength of the imposed forcing, or varying Reynolds number, respectively, time-asymptotic states, notably stable stationary solutions, have been traced. A primary non-magnetic steady state loses, in a Hopf bifurcation, stability to a periodic state with a non-vanishing magnetic field, showing the appearance of a generic dynamo effect. From now on the magnetic field is present for all values of the forcing. The Hopf bifurcation is followed by furhter, symmetry-breaking, bifurcations, leading finally to chaos. We pay particular attention to kinetic and magnetic helicities. The dynamo effect is observed only if the forcing is chosen such that a mean kinetic helicity is generated; otherwise the magnetic field diffuses away, and the time-asymptotic states are non-magnetic, in accordance with traditional kinematic dynamo theory.
The bifurcation behaviour of the 3D magnetohydrodynamic equations has been studied for external forcings of varying degree of helicity. With increasing strength of the forcing a primary non-magnetic steady state loses stability to a magnetic periodic state if the helicity exceeds a threshold value and to different non-magnetic states otherwise.
We have studied the bifurcations in a three-dimensional incompressible magnetofluid with periodic boundary conditions and an external forcing of the Arnold-Beltrami-Childress (ABC) type. Bifurcation-analysis techniques have been applied to explore the qualitative behavior of solution branches. Due to the symmetry of the forcing, the equations are equivariant with respect to a group of transformations isomorphic to the octahedral group, and we have paid special attention to symmetry-breaking effects. As the Reynolds number is increased, the primary nonmagnetic steady state, the ABC flow, loses its stability to a periodic magnetic state, showing the appearance of a generic dynamo effect; the critical value of the Reynolds number for the instability of the ABC flow is decreased compared to the purely hydrodynamic case. The bifurcating magnetic branch in turn is subject to secondary, symmetry-breaking bifurcations. We have traced periodic and quasi- periodic branches until they end up in chaotic states. In particular detail we have analyzed the subgroup symmetries of the bifurcating periodic branches, which are closely related to the spatial structure of the magnetic field.
We report on bifurcation studies for the incompressible Navier-Stokes equations in two space dimensions with periodic boundary conditions and an external forcing of the Kolmogorov type. Fourier representations of velocity and pressure have been used to approximate the original partial differential equations by a finite-dimensional system of ordinary differential equations, which then has been studied by means of bifurcation-analysis techniques. A special route into chaos observed for increasing Reynolds number or strength of the imposed forcing is described. It includes several steady states, traveling waves, modulated traveling waves, periodic and torus solutions, as well as a period-doubling cascade for a torus solution. Lyapunov exponents and Kaplan-Yorke dimensions have been calculated to characterize the chaotic branch. While studying the dynamics of the system in Fourier space, we also have transformed solutions to real space and examined the relation between the different bifurcations in Fourier space and toplogical changes of the streamline portrait. In particular, the time-dependent solutions, such as, e.g., traveling waves, torus, and chaotic solutions, have been characterized by the associated fluid-particle motion (Lagrangian dynamics).
We have studied bifurcation phenomena for the incompressable Navier-Stokes equations in two space dimensions with periodic boundary conditions. Fourier representations of velocity and pressure have been used to transform the original partial differential equations into systems of ordinary differential equations (ODE), to which then numerical methods for the qualitative analysis of systems of ODE have been applied, supplemented by the simulative calculation of solutions for selected initial conditions. Invariant sets, notably steady states, have been traced for varying Reynolds number or strength of the imposed forcing, respectively. A complete bifurcation sequence leading to chaos is described in detail, including the calculation of the Lyapunov exponents that characterize the resulting chaotic branch in the bifurcation diagram.
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 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.