Numerical simulation of combined heat transfer phenomena (conduction, convection, and radiation). Application to he optimisation of thermal systems

Thermal systems interchanging heat and mass by conduction, convection, radiation (solar and thermal ) occur in many engineering applications like energy storage by solar collectors, window glazing in buildings, refrigeration of plastic moulds, air handling units etc. Often these thermal systems are composed of various elements for example a building with wall, windows, rooms, etc. It would be of particular interest to have a modular thermal system which is formed by connecting different modules for the elements, flexibility to use and change models for individual elements, add or remove elements without changing the entire code. A numerical approach to handle the heat transfer and fluid flow in such systems helps in saving the full scale experiment time, cost and also aids optimisation of parameters of the system. In subsequent sections are presented a short summary of the work done until now on the orientation of the thesis in the field of numerical methods for heat transfer and fluid flow applications, the work in process and the future work.

Long-term outcome after cognitive and behavioral regression in nonlesional epilepsy with continuous spike-waves during slow-wave sleep.

PURPOSE: To present the long-term follow-up of 10 adolescents and young adults with documented cognitive and behavioral regression as children due to nonlesional focal, mainly frontal, epilepsy with continuous spike-waves during slow wave sleep (CSWS). METHODS: Past medical and electroencephalography (EEG) data were reviewed and neuropsychological tests exploring main cognitive functions were administered. KEY FINDINGS: After a mean duration of follow-up of 15.6 years (range, 8-23 years), none of the 10 patients had recovered fully, but four regained borderline to normal intelligence and were almost independent. Patients with prolonged global intellectual regression had the worst outcome, whereas those with more specific and short-lived deficits recovered best. The marked behavioral disorders resolved in all but one patient. Executive functions were neither severely nor homogenously affected. Three patients with a frontal syndrome during the active phase (AP) disclosed only mild residual executive and social cognition deficits. The main cognitive gains occurred shortly after the AP, but qualitative improvements continued to occur. Long-term outcome correlated best with duration of CSWS. SIGNIFICANCE: Our findings emphasize that cognitive recovery after cessation of CSWS depends on the severity and duration of the initial regression. None of our patients had major executive and social cognition deficits with preserved intelligence, as reported in adults with early destructive lesions of the frontal lobes. Early recognition of epilepsy with CSWS and rapid introduction of effective therapy are crucial for a best possible outcome.

Fast numerical algorithms for the computation of invariant tori in Hamiltonian systems

In this paper, we develop numerical algorithms that use small requirements of storage and operations for the computation of invariant tori in Hamiltonian systems (exact symplectic maps and Hamiltonian vector fields). The algorithms are based on the parameterization method and follow closely the proof of the KAM theorem given in [LGJV05] and [FLS07]. They essentially consist in solving a functional equation satisfied by the invariant tori by using a Newton method. Using some geometric identities, it is possible to perform a Newton step using little storage and few operations. In this paper we focus on the numerical issues of the algorithms (speed, storage and stability) and we refer to the mentioned papers for the rigorous results. We show how to compute efficiently both maximal invariant tori and whiskered tori, together with the associated invariant stable and unstable manifolds of whiskered tori. Moreover, we present fast algorithms for the iteration of the quasi-periodic cocycles and the computation of the invariant bundles, which is a preliminary step for the computation of invariant whiskered tori. Since quasi-periodic cocycles appear in other contexts, this section may be of independent interest. The numerical methods presented here allow to compute in a unified way primary and secondary invariant KAM tori. Secondary tori are invariant tori which can be contracted to a periodic orbit. We present some preliminary results that ensure that the methods are indeed implementable and fast. We postpone to a future paper optimized implementations and results on the breakdown of invariant tori.

Bilateral symmetry of radial pulse in high-level tennis players: implications for the validity of central aortic pulse wave analysis.

BACKGROUND: Reconstruction of the central aortic pressure wave from the noninvasive recording of the radial pulse with applanation tonometry has become a standard tool in the field of hypertension. It is not presently known whether recording the radial pulse on the dominant or the nondominant side has any effect on such reconstruction. METHOD: We carried out radial applanation tonometry on both forearms in young, healthy, male volunteers, who were either sedentary (n = 11) or high-level tennis players (n = 10). The purpose of including tennis players was to investigate individuals with extreme asymmetry between the dominant and nondominant upper limb. RESULTS: In the sedentary individuals, forearm circumference and handgrip strength were slightly larger on the dominant (mean +/- SD respectively 27.9 +/- 1.5 cm and 53.8 +/- 10 kg) than on nondominant side (27.3 +/- 1.6 cm, P < 0.001 vs. dominant, and 52.1 +/- 11 kg, P = NS). In the tennis players, differences between sides were more conspicuous (forearm circumference: dominant 28.0 +/- 1.7 cm nondominant 26.4 +/- 1.5 cm, P < 0.001; handgrip strength 61.4 +/- 10.8 vs. 53.4 +/- 9.7 kg, P < 0.001). We found that in both sedentary individuals and tennis players, the radial pulse had identical shape on both sides and, consequently, the reconstructed central aortic pressure waveforms, as well as derived indices of central pulsatility, were not dependent on the side where applanation tonometry was carried out. CONCLUSION: Evidence from individuals with maximal asymmetry of dominant vs. nondominant upper limb indicates that laterality of measurement is not a methodological issue for central pulse wave analysis carried out with radial applanation tonometry.

Quantifying rates of landscape evolution and tectonic processes by thermochronology and numerical modeling of crustal heat transport using PECUBE

PECUBE is a three-dimensional thermal-kinematic code capable of solving the heat production-diffusion-advection equation under a temporally varying surface boundary condition. It was initially developed to assess the effects of time-varying surface topography (relief) on low-temperature thermochronological datasets. Thermochronometric ages are predicted by tracking the time-temperature histories of rock-particles ending up at the surface and by combining these with various age-prediction models. In the decade since its inception, the PECUBE code has been under continuous development as its use became wider and addressed different tectonic-geomorphic problems. This paper describes several major recent improvements in the code, including its integration with an inverse-modeling package based on the Neighborhood Algorithm, the incorporation of fault-controlled kinematics, several different ways to address topographic and drainage change through time, the ability to predict subsurface (tunnel or borehole) data, prediction of detrital thermochronology data and a method to compare these with observations, and the coupling with landscape-evolution (or surface-process) models. Each new development is described together with one or several applications, so that the reader and potential user can clearly assess and make use of the capabilities of PECUBE. We end with describing some developments that are currently underway or should take place in the foreseeable future. (C) 2012 Elsevier B.V. All rights reserved.

A numerical algorithm for finding solutions of a generalized Nash equilibrium problem

A family of nonempty closed convex sets is built by using the data of the Generalized Nash equilibrium problem (GNEP). The sets are selected iteratively such that the intersection of the selected sets contains solutions of the GNEP. The algorithm introduced by Iusem-Sosa (2003) is adapted to obtain solutions of the GNEP. Finally some numerical experiments are given to illustrate the numerical behavior of the algorithm.

A numerical solver for a nonlinear Fokker-Planck equation representation of neuronal network dynamics

To describe the collective behavior of large ensembles of neurons in neuronal network, a kinetic theory description was developed in [13, 12], where a macroscopic representation of the network dynamics was directly derived from the microscopic dynamics of individual neurons, which are modeled by conductance-based, linear, integrate-and-fire point neurons. A diffusion approximation then led to a nonlinear Fokker-Planck equation for the probability density function of neuronal membrane potentials and synaptic conductances. In this work, we propose a deterministic numerical scheme for a Fokker-Planck model of an excitatory-only network. Our numerical solver allows us to obtain the time evolution of probability distribution functions, and thus, the evolution of all possible macroscopic quantities that are given by suitable moments of the probability density function. We show that this deterministic scheme is capable of capturing the bistability of stationary states observed in Monte Carlo simulations. Moreover, the transient behavior of the firing rates computed from the Fokker-Planck equation is analyzed in this bistable situation, where a bifurcation scenario, of asynchronous convergence towards stationary states, periodic synchronous solutions or damped oscillatory convergence towards stationary states, can be uncovered by increasing the strength of the excitatory coupling. Finally, the computation of moments of the probability distribution allows us to validate the applicability of a moment closure assumption used in [13] to further simplify the kinetic theory.

High-resolution numerical modelling of flow-vegetation interactions

In this paper, we present and apply a new three-dimensional model for the prediction of canopy-flow and turbulence dynamics in open-channel flow. The approach uses a dynamic immersed boundary technique that is coupled in a sequentially staggered manner to a large eddy simulation. Two different biomechanical models are developed depending on whether the vegetation is dominated by bending or tensile forces. For bending plants, a model structured on the Euler-Bernoulli beam equation has been developed, whilst for tensile plants, an N-pendula model has been developed. Validation against flume data shows good agreement and demonstrates that for a given stem density, the models are able to simulate the extraction of energy from the mean flow at the stem-scale which leads to the drag discontinuity and associated mixing layer.

Diseño de filtros Bulk Acoustic Wave con topología en escalera de altas prestaciones

Degut als avenços als dispositius de telecomunicacions durant l’última dècada, els filtres integrats en aquests dispositius requereixen de millors prestacions, baix cost i, per sobre de tot, requereixen unes dimensions el més reduïdes possibles. Tot i que avui dia encara s’utilitzen el filtres SAW en aquests dispositius, cada cop més s’estan substituint pels filtres amb tecnologia BAW, ja que tenen millors prestacions. En l’actualitat la topologia BAW més extensa i utilitzada és la topologia en escala. En aquest projecte s’ha portat a terme un estudi en profunditat de les limitacions dels filtres en escala. A partir de les limitacions detectades s’ha presentat una nova estructura de disseny per aquest tipus de filtres que redueix les dimensions d’aquests i millora considerablement algunes de les limitacions de l’estructura convencional. Paral·lelament s’ha desenvolupat una metodologia sistemàtica pròpia pel disseny de la nova estructura.

Numerical Taxonomy of Old World Phlebotominae (Diptera: Psychodidae): 1. Considerations of Morphological Characters in the Genus Phlebotomus Rondani & Berté 1840

Numerical analyses (correspondence analysis, ascending hierarchical classification, cladistic approach) were applied to the morphological characters of the adults of the genus Phlebotomus Rondani & Berté 1840. They confirm the reliability of the classic classifications, and also redefine the taxonomic and phylogenetic position of certain taxa. Thus, Spelaeophlebotomus Theodor 1948, Idiophlebotomus Quate & Fairchild 1961 and Australophlebotomus Theodor 1948 deserve generic rank. Among the vectors of leishmaniasis, the subgenus Phlebotomus Rondani & Berté 1840 is probably ancient. The results attribute an intermediate taxonomic and phylogenetic position to the taxa Euphlebotomus Theodor 1948 and Anaphlebotomus Theodor 1948, and reveal the probable artificial nature of the latter. The comparatively large numbers of species of subgenera Paraphlebotomus Theodor 1948, Synphlebotomus Theodor 1948 and, above all, Larroussius Nitzulescu 1931 and Adlerius Nitzulescu 1931, suggest that they are relatively recent. The development of adult morphological characters, the validity of their use in taxonomy and proposals for further studies are discussed.

Numerical Taxonomy of Old World Phlebotominae (Diptera: Psychodidae): 2. Restatement of Classification upon Subgeneric Morphological Characters

Numerical analyses (correspondence analysis, ascending hierarchical classification, and cladistics) were done with morphological characters of adult phlebotomine sand flies. The resulting classification largely confirms that of classical taxonomy for supra-specific groups from the Old World, though the positions of some groups are adjusted. The taxa Spelaeophlebotomus Theodor 1948, Idiophlebotomus Quate & Fairchild 1961, Australophlebotomus Theodor 1948 and Chinius Leng 1987 are notably distinct from other Old World groups, particularly from the genus Phlebotomus Rondani & Berté 1840. Spelaeomyia Theodor 1948 and, in particular, Parvidens Theodor & Mesghali 1964 are clearly separate from Sergentomyia França & Parrot 1920.

Confounding effects of heart rate on pulse wave velocity in paced patients with a low degree of atherosclerosis.

BACKGROUND: Pulse wave velocity (PWV), an index of arterial wall stiffness, is modulated by blood pressure (BP). Whether heart rate (HR) is also a modulator of PWV is controversial. Recent research involving mainly patients with high aortic PWV have found either no change or a positive correlation between the two. Given that PWV is increasingly being measured in cardiovascular studies, the relationship between HR and PWV should be known in patients with preserved arterial wall elasticity. OBJECTIVE: The aim of this study was to evaluate the importance of HR as a determinant of the variability in PWV in patients with a low degree of atherosclerosis. DESIGN AND METHODS: Fourteen patients (five female, nine male; aged 68 +/- 8 years) were evaluated post pacemaker implantation due to sick sinus or carotid hypersensitivity syndromes. Carotid-femoral PWV was measured at rest and during atrial pacing at 80, 90 and 100 bpm (paced HR). Arterial femoral blood flow (AFBF) was measured by echodoppler. RESULTS: PWV increased from 6.2 +/- 1.5 m/s (mean +/- SD) during resting sinus rhythm (HR 62 +/- 8 bpm; mean +/- SD) to 6.8 +/- 1.0, 7.0 +/- 0.9, and 7.6 +/- 1.1 m/s at pacing rates of 80, 90 and 100 bpm, respectively (P < 0.0001). Systolic (SBP) and mean blood pressure (MBP) remained constant at all HR levels, whereas AFBF increased in a linear fashion. CONCLUSIONS: These results demonstrate that even in patients with a low degree of atherosclerosis, HR is a potential modulator of carotid-femoral PWV.

Analyzing complex rock slope deformation at Stampa, western Norway, by integrating geomorphology, kinematics and numerical modeling

The unstable rock slope, Stampa, above the village of Flåm, Norway, shows signs of both active and postglacial gravitational deformation over an area of 11 km2. Detailed structural field mapping, annual differential Global Navigation Satellite System (GNSS) surveys, as well as geomorphic analysis of high-resolution digital elevation models based on airborne and terrestrial laser scanning indicate that slope deformation is complex and spatially variable. Numerical modeling was used to investigate the influence of former rockslide activity and to better understand the failure mechanism. Field observations, kinematic analysis and numerical modeling indicate a strong structural control of the unstable area. Based on the integration of the above analyses, we propose that the failure mechanism is dominated by (1) a toppling component, (2) subsiding bilinear wedge failure and (3) planar sliding along the foliation at the toe of the unstable slope. Using differential GNSS, 18 points were measured annually over a period of up to 6 years. Two of these points have an average yearly movement of around 10 mm/year. They are located at the frontal cliff on almost completely detached blocks with volumes smaller than 300,000 m3. Large fractures indicate deep-seated gravitational deformation of volumes reaching several 100 million m3, but the movement rates in these areas are below 2 mm/year. Two different lobes of prehistoric rock slope failures were dated with terrestrial cosmogenic nuclides. While the northern lobe gave an average age of 4,300 years BP, the southern one resulted in two different ages (2,400 and 12,000 years BP), which represent most likely multiple rockfall events. This reflects the currently observable deformation style with unstable blocks in the northern part in between Joasete and Furekamben and no distinct blocks but a high rockfall activity around Ramnanosi in the south. With a relative susceptibility analysis it is concluded that small collapses of blocks along the frontal cliff will be more frequent. Larger collapses of free-standing blocks along the cliff with volumes > 100,000 m3, thus large enough to reach the fjord, cannot be ruled out. A larger collapse involving several million m3 is presently considered of very low likelihood.