Covalidation of Integral Transforms and Method of Lines in Nonlinear Convection-Diffusion with Mathematica

The Mathematica system (version 4.0) is employed in the solution of nonlinear difusion and convection-difusion problems, formulated as transient one-dimensional partial diferential equations with potential dependent equation coefficients. The Generalized Integral Transform Technique (GITT) is first implemented for the hybrid numerical-analytical solution of such classes of problems, through the symbolic integral transformation and elimination of the space variable, followed by the utilization of the built-in Mathematica function NDSolve for handling the resulting transformed ODE system. This approach ofers an error-controlled final numerical solution, through the simultaneous control of local errors in this reliable ODE's solver and of the proposed eigenfunction expansion truncation order. For covalidation purposes, the same built-in function NDSolve is employed in the direct solution of these partial diferential equations, as made possible by the algorithms implemented in Mathematica (versions 3.0 and up), based on application of the method of lines. Various numerical experiments are performed and relative merits of each approach are critically pointed out.

Estimating Attractor Dimension on the Nonlinear Pendulum Time Series

Chaotic dynamical systems exhibit trajectories in their phase space that converges to a strange attractor. The strangeness of the chaotic attractor is associated with its dimension in which instance it is described by a noninteger dimension. This contribution presents an overview of the main definitions of dimension discussing their evaluation from time series employing the correlation and the generalized dimension. The investigation is applied to the nonlinear pendulum where signals are generated by numerical integration of the mathematical model, selecting a single variable of the system as a time series. In order to simulate experimental data sets, a random noise is introduced in the time series. State space reconstruction and the determination of attractor dimensions are carried out regarding periodic and chaotic signals. Results obtained from time series analyses are compared with a reference value obtained from the analysis of mathematical model, estimating noise sensitivity. This procedure allows one to identify the best techniques to be applied in the analysis of experimental data.

Simulation of Steady-State Nonlinear Heat Transfer Problems Through the Minimization of Quadratic Functionals

In this work it is presented a systematic procedure for constructing the solution of a large class of nonlinear conduction heat transfer problems through the minimization of quadratic functionals like the ones usually employed for linear descriptions. The proposed procedure gives rise to an efficient and easy way for carrying out numerical simulations of nonlinear heat transfer problems by means of finite elements. To illustrate the procedure a particular problem is simulated by means of a finite element approximation.

An Investigation into Mechanisms of Loss of Safe Basins in a 2 D.O.F. Nonlinear Oscillator

In this work we consider the transient stability of coupled motions of a 2 D.O.F. nonlinear oscillator that can represent, for example, the motions of a sea vessel under the action of trains of regular lateral waves. Instability is studied as the escape of the system from a safe potential well. The set of initial conditions in phase space that lead to acceptable motions constitutes its safe basin. We investigate the evolution of these safe basins under variation of parameters such as frequency and amplitude of waves, and an internal tuning parameter. Complex nonlinear phenomena are known to play an important role in determining the loss of safe basins as, say, wave amplitude is increased. We therefore investigate those processes, and attempt to classify them in terms of their speed relative to changes in parameter values. "Mechanism basins" are produced depicting regions of parameter space in which rapid or slow losses of safe basin are observed. We propose that a comprehensive understanding of mechanisms of loss of safe basins can be a valuable tool in assessing stability properties of these systems, and we give a conceptual view of how such information could be used.

Experimental Vibration Characteristics of a Beam with Nonlinear Support Using Receptance Coupling Analysis

This paper applies the Multi-Harmonic Nonlinear Receptance Coupling Approach (MUHANORCA) (Ferreira 1998) to evaluate the frequency response characteristics of a beam which is clamped at one end and supported at the other end by a nonlinear cubic stiffness joint. In order to apply the substructure coupling technique, the problem was characterised by coupling a clamped linear beam with a nonlinear cubic stiffness joint. The experimental results were obtained by a sinusoidal excitation with a special force control algorithm where the level of the fundamental force is kept constant and the level of the harmonics is kept zero for all the frequencies measured.

Application of the center manifold theory to the study of slewing flexible Non-ideal structures with nonlinear curvature: a case study

In this paper is Analyzed the local dynamical behavior of a slewing flexible structure considering nonlinear curvature. The dynamics of the original (nonlinear) governing equations of motion are reduced to the center manifold in the neighborhood of an equilibrium solution with the purpose of locally study the stability of the system. In this critical point, a Hopf bifurcation occurs. In this region, one can find values for the control parameter (structural damping coefficient) where the system is unstable and values where the system stability is assured (periodic motion). This local analysis of the system reduced to the center manifold assures the stable / unstable behavior of the original system around a known solution.

Lung tissue mechanics in the early stages of induced paracoccidioidomycosis in rats

Pulmonary dysfunction represents the most important cause of death in patients with paracoccidioidomycosis (PBM). In order to investigate the functional changes of the lungs in the early stages of PBM, a model of benign disease was developed by intratracheal challenge of 12-week old isogenic Wistar rats with 1 x 106 yeast forms of Paracoccidioides brasiliensis. Animals were studied 30 and 60 days after infection, when fully developed granulomas were demonstrable in the lungs. Measurements of airway resistance, lung elastance and tissue hysteresis were made during sinusoidal deformations (100 breaths/min, tidal volume = 2 ml) with direct measurement of alveolar pressure using the alveolar capsule technique. Infection caused a significant increase in hysteresis (infected: 1.69, N = 13; control: 1.13, N = 12, P = 0.024, ANOVA), with no alterations in airway resistance or lung elastance. Histopathological analysis revealed the presence of fully developed granulomas located in the axial compartment of the lung interstitial space. These results suggest that alterations of tissue mechanics represent an early event in experimental PBM

The effects of positive end-expiratory pressure on respiratory system mechanics and hemodynamics in postoperative cardiac surgery patients

We prospectively evaluated the effects of positive end-expiratory pressure (PEEP) on the respiratory mechanical properties and hemodynamics of 10 postoperative adult cardiac patients undergoing mechanical ventilation while still anesthetized and paralyzed. The respiratory mechanics was evaluated by the inflation inspiratory occlusion method and hemodynamics by conventional methods. Each patient was randomized to a different level of PEEP (5, 10 and 15 cmH2O), while zero end-expiratory pressure (ZEEP) was established as control. PEEP of 15-min duration was applied at 20-min intervals. The frequency dependence of resistance and the viscoelastic properties and elastance of the respiratory system were evaluated together with hemodynamic and respiratory indexes. We observed a significant decrease in total airway resistance (13.12 ± 0.79 cmH2O l-1 s-1 at ZEEP, 11.94 ± 0.55 cmH2O l-1 s-1 (P<0.0197) at 5 cmH2O of PEEP, 11.42 ± 0.71 cmH2O l-1 s-1 (P<0.0255) at 10 cmH2O of PEEP, and 10.32 ± 0.57 cmH2O l-1 s-1 (P<0.0002) at 15 cmH2O of PEEP). The elastance (Ers; cmH2O/l) was not significantly modified by PEEP from zero (23.49 ± 1.21) to 5 cmH2O (21.89 ± 0.70). However, a significant decrease (P<0.0003) at 10 cmH2O PEEP (18.86 ± 1.13), as well as (P<0.0001) at 15 cmH2O (18.41 ± 0.82) was observed after PEEP application. Volume dependence of viscoelastic properties showed a slight but not significant tendency to increase with PEEP. The significant decreases in cardiac index (l min-1 m-2) due to PEEP increments (3.90 ± 0.22 at ZEEP, 3.43 ± 0.17 (P<0.0260) at 5 cmH2O of PEEP, 3.31 ± 0.22 (P<0.0260) at 10 cmH2O of PEEP, and 3.10 ± 0.22 (P<0.0113) at 15 cmH2O of PEEP) were compensated for by an increase in arterial oxygen content owing to shunt fraction reduction (%) from 22.26 ± 2.28 at ZEEP to 11.66 ± 1.24 at PEEP of 15 cmH2O (P<0.0007). We conclude that increments in PEEP resulted in a reduction of both airway resistance and respiratory elastance. These results could reflect improvement in respiratory mechanics. However, due to possible hemodynamic instability, PEEP should be carefully applied to postoperative cardiac patients.

Efficient Optimization Algorithms for Nonlinear Data Analysis

Identification of low-dimensional structures and main sources of variation from multivariate data are fundamental tasks in data analysis. Many methods aimed at these tasks involve solution of an optimization problem. Thus, the objective of this thesis is to develop computationally efficient and theoretically justified methods for solving such problems. Most of the thesis is based on a statistical model, where ridges of the density estimated from the data are considered as relevant features. Finding ridges, that are generalized maxima, necessitates development of advanced optimization methods. An efficient and convergent trust region Newton method for projecting a point onto a ridge of the underlying density is developed for this purpose. The method is utilized in a differential equation-based approach for tracing ridges and computing projection coordinates along them. The density estimation is done nonparametrically by using Gaussian kernels. This allows application of ridge-based methods with only mild assumptions on the underlying structure of the data. The statistical model and the ridge finding methods are adapted to two different applications. The first one is extraction of curvilinear structures from noisy data mixed with background clutter. The second one is a novel nonlinear generalization of principal component analysis (PCA) and its extension to time series data. The methods have a wide range of potential applications, where most of the earlier approaches are inadequate. Examples include identification of faults from seismic data and identification of filaments from cosmological data. Applicability of the nonlinear PCA to climate analysis and reconstruction of periodic patterns from noisy time series data are also demonstrated. Other contributions of the thesis include development of an efficient semidefinite optimization method for embedding graphs into the Euclidean space. The method produces structure-preserving embeddings that maximize interpoint distances. It is primarily developed for dimensionality reduction, but has also potential applications in graph theory and various areas of physics, chemistry and engineering. Asymptotic behaviour of ridges and maxima of Gaussian kernel densities is also investigated when the kernel bandwidth approaches infinity. The results are applied to the nonlinear PCA and to finding significant maxima of such densities, which is a typical problem in visual object tracking.

Respiratory mechanics and morphometric changes during anesthesia with ketamine in normal rats

Ketamine is believed to reduce airway and pulmonary tissue resistance. The aim of the present study was to determine the effects of ketamine on the resistive, elastic and viscoelastic/inhomogeneous mechanical properties of the respiratory system, lungs and chest wall, and to relate the mechanical data to findings from histological lung analysis in normal animals. Fifteen adult male Wistar rats were assigned randomly to two groups: control (N = 7) and ketamine (N = 8). All animals were sedated (diazepam, 5 mg, ip) and anesthetized with pentobarbital sodium (20 mg/kg, ip) or ketamine (30 mg/kg, ip). The rats were paralyzed and ventilated mechanically. Ketamine increased lung viscoelastic/inhomogeneous pressure (26%) compared to the control group. Dynamic and static elastances were similar in both groups, but the difference was greater in the ketamine than in the control group. Lung morphometry demonstrated dilation of alveolar ducts and increased areas of alveolar collapse in the ketamine group. In conclusion, ketamine did not act at the airway level but acted at the lung periphery increasing mechanical inhomogeneities possibly resulting from dilation of distal airways and alveolar collapse.

Initialization of Continuous Nonlinear Models Using Extended Kalman Filter

The two main objectives of Bayesian inference are to estimate parameters and states. In this thesis, we are interested in how this can be done in the framework of state-space models when there is a complete or partial lack of knowledge of the initial state of a continuous nonlinear dynamical system. In literature, similar problems have been referred to as diffuse initialization problems. This is achieved first by extending the previously developed diffuse initialization Kalman filtering techniques for discrete systems to continuous systems. The second objective is to estimate parameters using MCMC methods with a likelihood function obtained from the diffuse filtering. These methods are tried on the data collected from the 1995 Ebola outbreak in Kikwit, DRC in order to estimate the parameters of the system.

Linear and nonlinear analysis of heart rate variability in healthy subjects and after acute myocardial infarction in patients

The objectives of this study were to evaluate and compare the use of linear and nonlinear methods for analysis of heart rate variability (HRV) in healthy subjects and in patients after acute myocardial infarction (AMI). Heart rate (HR) was recorded for 15 min in the supine position in 10 patients with AMI taking β-blockers (aged 57 ± 9 years) and in 11 healthy subjects (aged 53 ± 4 years). HRV was analyzed in the time domain (RMSSD and RMSM), the frequency domain using low- and high-frequency bands in normalized units (nu; LFnu and HFnu) and the LF/HF ratio and approximate entropy (ApEn) were determined. There was a correlation (P < 0.05) of RMSSD, RMSM, LFnu, HFnu, and the LF/HF ratio index with the ApEn of the AMI group on the 2nd (r = 0.87, 0.65, 0.72, 0.72, and 0.64) and 7th day (r = 0.88, 0.70, 0.69, 0.69, and 0.87) and of the healthy group (r = 0.63, 0.71, 0.63, 0.63, and 0.74), respectively. The median HRV indexes of the AMI group on the 2nd and 7th day differed from the healthy group (P < 0.05): RMSSD = 10.37, 19.95, 24.81; RMSM = 23.47, 31.96, 43.79; LFnu = 0.79, 0.79, 0.62; HFnu = 0.20, 0.20, 0.37; LF/HF ratio = 3.87, 3.94, 1.65; ApEn = 1.01, 1.24, 1.31, respectively. There was agreement between the methods, suggesting that these have the same power to evaluate autonomic modulation of HR in both AMI patients and healthy subjects. AMI contributed to a reduction in cardiac signal irregularity, higher sympathetic modulation and lower vagal modulation.

Molecular mechanics calculations of tin-bis (triphenylphosphine oxide) complexes

The capability of molecular mechanics for modeling the wide distribution of bond angles and bond lengths characteristic of coordination complexes was investigatecl. This was the preliminary step for future modeling of solvent extraction. Several tin-phosphine oxide COrnI)le:){es were selected as the test groUl) for t.he d,esired range of geometry they eX!libi ted as \-vell as the ligands they cOD.tained r Wllich were c\f interest in connection with solvation. A variety of adjustments were made to Allinger's M:M2 force·-field ill order to inl.prove its performance in the treatment of these systems. A set of u,nique force constants was introduced for' those terms representing the metal ligand bond lengths, bond angles, and, torsion angles. These were significantly smaller than trad.itionallY used. with organic compounds. The ~1orse poteIlt.ial energ'Y function was incorporated for the M-X l')ond lE~ngths and the cosine harmonic potential erlerg-y function was invoked for the MOP bond angle. These functions were found to accomodate the wide distribution of observed values better than the traditional harmonic approximations~ Crystal packing influences on the MOP angle were explored thr"ollgh ttle inclusion of the isolated molecule withil1 a shell cc)ntaini11g tl1e nearest neigl1'bors duri.rlg energy rninimization experiments~ This was found to further improve the fit of the MOP angle.