Problems in four-dimensional data assimilation

With the introduction of new observing systems based on asynoptic observations, the analysis problem has changed in character. In the near future we may expect that a considerable part of meteorological observations will be unevenly distributed in four dimensions, i.e. three dimensions in space and one in time. The term analysis, or objective analysis in meteorology, means the process of interpolating observed meteorological observations from unevenly distributed locations to a network of regularly spaced grid points. Necessitated by the requirement of numerical weather prediction models to solve the governing finite difference equations on such a grid lattice, the objective analysis is a three-dimensional (or mostly two-dimensional) interpolation technique. As a consequence of the structure of the conventional synoptic network with separated data-sparse and data-dense areas, four-dimensional analysis has in fact been intensively used for many years. Weather services have thus based their analysis not only on synoptic data at the time of the analysis and climatology, but also on the fields predicted from the previous observation hour and valid at the time of the analysis. The inclusion of the time dimension in objective analysis will be called four-dimensional data assimilation. From one point of view it seems possible to apply the conventional technique on the new data sources by simply reducing the time interval in the analysis-forecasting cycle. This could in fact be justified also for the conventional observations. We have a fairly good coverage of surface observations 8 times a day and several upper air stations are making radiosonde and radiowind observations 4 times a day. If we have a 3-hour step in the analysis-forecasting cycle instead of 12 hours, which is applied most often, we may without any difficulties treat all observations as synoptic. No observation would thus be more than 90 minutes off time and the observations even during strong transient motion would fall within a horizontal mesh of 500 km * 500 km.

Multivariate analysis of the dielectric response of materials modeled using networks of resistors and capacitors

We discuss the modeling of dielectric responses of electromagnetically excited networks which are composed of a mixture of capacitors and resistors. Such networks can be employed as lumped-parameter circuits to model the response of composite materials containing conductive and insulating grains. The dynamics of the excited network systems are studied using a state space model derived from a randomized incidence matrix. Time and frequency domain responses from synthetic data sets generated from state space models are analyzed for the purpose of estimating the fraction of capacitors in the network. Good results were obtained by using either the time-domain response to a pulse excitation or impedance data at selected frequencies. A chemometric framework based on a Successive Projections Algorithm (SPA) enables the construction of multiple linear regression (MLR) models which can efficiently determine the ratio of conductive to insulating components in composite material samples. The proposed method avoids restrictions commonly associated with Archie’s law, the application of percolation theory or Kohlrausch-Williams-Watts models and is applicable to experimental results generated by either time domain transient spectrometers or continuous-wave instruments. Furthermore, it is quite generic and applicable to tomography, acoustics as well as other spectroscopies such as nuclear magnetic resonance, electron paramagnetic resonance and, therefore, should be of general interest across the dielectrics community.

Numerical prediction of three-dimensional time-dependent viscoelastic extrudate swell using differential and algebraic models

This study investigates the numerical simulation of three-dimensional time-dependent viscoelastic free surface flows using the Upper-Convected Maxwell (UCM) constitutive equation and an algebraic explicit model. This investigation was carried out to develop a simplified approach that can be applied to the extrudate swell problem. The relevant physics of this flow phenomenon is discussed in the paper and an algebraic model to predict the extrudate swell problem is presented. It is based on an explicit algebraic representation of the non-Newtonian extra-stress through a kinematic tensor formed with the scaled dyadic product of the velocity field. The elasticity of the fluid is governed by a single transport equation for a scalar quantity which has dimension of strain rate. Mass and momentum conservations, and the constitutive equation (UCM and algebraic model) were solved by a three-dimensional time-dependent finite difference method. The free surface of the fluid was modeled using a marker-and-cell approach. The algebraic model was validated by comparing the numerical predictions with analytic solutions for pipe flow. In comparison with the classical UCM model, one advantage of this approach is that computational workload is substantially reduced: the UCM model employs six differential equations while the algebraic model uses only one. The results showed stable flows with very large extrudate growths beyond those usually obtained with standard differential viscoelastic models. (C) 2010 Elsevier Ltd. All rights reserved.

Stability of numerical schemes on staggered grids

This paper considers the stability of explicit, implicit and Crank-Nicolson schemes for the one-dimensional heat equation on a staggered grid. Furthemore, we consider the cases when both explicit and implicit approximations of the boundary conditions arc employed. Why we choose to do this is clearly motivated and arises front solving fluid flow equations with free surfaces when the Reynolds number can be very small. in at least parts of the spatial domain. A comprehensive stability analysis is supplied: a novel result is the precise stability restriction on the Crank-Nicolson method when the boundary conditions are approximated explicitly, that is, at t =n delta t rather than t = (n + 1)delta t. The two-dimensional Navier-Stokes equations were then solved by a marker and cell approach for two simple problems that had analytic solutions. It was found that the stability results provided in this paper were qualitatively very similar. thereby providing insight as to why a Crank-Nicolson approximation of the momentum equations is only conditionally, stable. Copyright (C) 2008 John Wiley & Sons, Ltd.

Design and analysis of a low actuation voltage electrowetting-on-dielectric device

This paper presents an Electrowetting-on-Dielectric (EWOD) device with optimized insulating layers operated by low actuation voltage. The device consists of an electrode array on a silicon substrate, covered by a dielectric layer and a hydrophobic layer. To characterize the performance of the device, simulations are performed for the dielectric layer of Sio₂ and the hydrophobic layer of Sio₂, Su-8 and Parylene C at different voltages. The volume finite difference approach of the Coventorware software was used to carry out the simulations. Two different molar of di-ionized water droplet were considered in the simulations. It was observed that the device having the Sio₂ dielectric layer and the Parylene C hydrophobic layer moved the 1M KCL (potassium chloride) droplet at the actuation voltage of 25V.

A comparative study of different power system stabilizers for dynamic stability analysis

Power system stabilizers (PSSs) are extensively used to ensure the dynamic stability of power systems through the modulation of excitation signals supplied to synchronous generators. This paper presents a comparative study of two different PSSs: STAB1 and IEEEST. The stabilizers are designed for the linearized model of a single machine infinite bus (SMIB) system with different loads. Both time-and frequency-domain simulations are carried out to investigate the performance of these stabilizers. For all PSSs, the time-domain simulations are performed by applying a three-phase short-circuit fault at the terminal of the synchronous generator. These simulation results are compared against the open-loop characteristics of the SMIB system where no PSS is implemented. Simulation results demonstrate that the speed-fed PSS provides more damping as compared to frequency- and power-fed stabilizers.

Disfunção autonômica e doença vascular periférica em pacientes com diabete melito tipo 2

A disfunção autonômica está associada com aumento da mortalidade em pacientes diabéticos, especialmente naqueles com doença cardiovascular. Neuropatia periférica, mau controle glicêmico, dislipidemia e hipertensão são alguns dos fatores de risco para o desenvolvimento de doença vascular periférica (DVP) nestes pacientes. O objetivo deste estudo foi avaliar os fatores de risco associados com a presença de DVP em pacientes com DM tipo 2. Um estudo transversal foi realizado em 84 pacientes com DM tipo 2 ( 39 homens, idade média de 64,9 ± 7,5 anos). Os pacientes foram submetidos a uma avaliação clínica e laboratorial. A presença de DVP foi definida, utilizando-se um um aparelho manual de ultrasom com doppler (índice perna-braço < 0,9). A atividade autonômica foi avaliada através da análise da variabilidade da freqüência cardíaca (HRV) por métodos no domínio do tempo e da freqüência (análise espectral), e pelo mapa de retorno tridimensional durante o período do dia e da noite. Para a análise da HRV, um eletrocardiograma de 24 horas foi gravado e as fitas analisadas em um analisador de Holter Mars 8000 (Marquete). A potência espectral foi quantificada pela área em duas bandas de freqüência: 0,04-0,15 Hz – baixa freqüência (BF), 0,015-0,5 Hz – alta freqüência (AF). A razão BF/AF foi calculada em cada paciente. O mapa de retorno tridimensional foi construído através de um modelo matemático onde foram analisados os intervalos RR versus a diferença entre os intervalos RR adjacentes versus o número de contagens verificadas, e quantificado por três índices refletindo a modulação simpática (P1) e vagal (P2 e P3). DVP estava presente em 30 (36%) pacientes. Na análise univariada, pacientes com DVP apresentaram índices que refletem a modulação autonômica (análise espectral) diminuídos quando comparados aos pacientes sem DVP, respectivamente: BF = 0,19 ± 0,07 m/s2 vs. 0,29 ± 0,11 m/s2 P = 0,0001; BF/AF = 1,98 ± 0,9 m/s2 vs. 3,35 ± 1,83 m/s2 p = 0,001. Além disso, o índice que reflete a atividade simpática no mapa de retorno tridimensional (P1), foi mais baixo em pacientes com DVP (61,7 ± 9,4 vs. 66,8 ± 9,7 unidades arbitrárias, P = 0,04) durante a noite, refletindo maior ativação simpática neste período. Estes pacientes também apresentavam uma maior duração do diabetes (20 ± 8,1 vs. 15,3 ± 6,7 anos, P = 0,006), níveis de pressão arterial sistólica (154 ± 20 vs. 145 ± 20 mmHg, P = 0,04), razão cintura-quadril ( 0,98 ± 0,09 vs.0,92 ± 0,08, P = 0,01), e níveis de HbA1c mais elevados (7,7 ± 1,6 vs. 6,9 ± 1,7 %, P = 0,04), bem como valores de triglicerídeos ( 259 ± 94 vs. 230 ± 196 mg/dl, P= 0,03) e de excreção urinária de albumina ( 685,5 ± 1359,9 vs. 188,2 ± 591,1 μ/min, P = 0,02) superiores aos dos pacientes sem DVP.. Nos pacientes com DVP observou-se uma presença aumentada de nefropatia diabética (73,3% vs. 29,6% P = 0,0001), de retinopatia (73,3% vs. 44,4% P = 0,02) e neuropatia periférica (705 vs. 35,1% P = 0,006). Os grupos não diferiram quanto à idade, índice de massa corporal, tabagismo e presença de doença arterial coronariana. Na análise logística multivariada, a DVP permaneceu associada com a disfunção autonômica, mesmo após ter sido controlada pela pressão arterial sistólica, duração do DM, HbA1c, triglicerídeos e excreção urinária de albumina. Concluindo, pacientes com DVP e DM tipo 2 apresentam índices que refletem a modulação autonômica diminuídos, o que pode representar um fator de risco adicional para o aumento da mortalidade nestes pacientes.

A fully adaptive multiresolution scheme for shock computations

The scheme is based on Ami Harten's ideas (Harten, 1994), the main tools coming from wavelet theory, in the framework of multiresolution analysis for cell averages. But instead of evolving cell averages on the finest uniform level, we propose to evolve just the cell averages on the grid determined by the significant wavelet coefficients. Typically, there are few cells in each time step, big cells on smooth regions, and smaller ones close to irregularities of the solution. For the numerical flux, we use a simple uniform central finite difference scheme, adapted to the size of each cell. If any of the required neighboring cell averages is not present, it is interpolated from coarser scales. But we switch to ENO scheme in the finest part of the grids. To show the feasibility and efficiency of the method, it is applied to a system arising in polymer-flooding of an oil reservoir. In terms of CPU time and memory requirements, it outperforms Harten's multiresolution algorithm.The proposed method applies to systems of conservation laws in 1Dpartial derivative(t)u(x, t) + partial derivative(x)f(u(x, t)) = 0, u(x, t) is an element of R-m. (1)In the spirit of finite volume methods, we shall consider the explicit schemeupsilon(mu)(n+1) = upsilon(mu)(n) - Deltat/hmu ((f) over bar (mu) - (f) over bar (mu)-) = [Dupsilon(n)](mu), (2)where mu is a point of an irregular grid Gamma, mu(-) is the left neighbor of A in Gamma, upsilon(mu)(n) approximate to 1/mu-mu(-) integral(mu-)(mu) u(x, t(n))dx are approximated cell averages of the solution, (f) over bar (mu) = (f) over bar (mu)(upsilon(n)) are the numerical fluxes, and D is the numerical evolution operator of the scheme.According to the definition of (f) over bar (mu), several schemes of this type have been proposed and successfully applied (LeVeque, 1990). Godunov, Lax-Wendroff, and ENO are some of the popular names. Godunov scheme resolves well the shocks, but accuracy (of first order) is poor in smooth regions. Lax-Wendroff is of second order, but produces dangerous oscillations close to shocks. ENO schemes are good alternatives, with high order and without serious oscillations. But the price is high computational cost.Ami Harten proposed in (Harten, 1994) a simple strategy to save expensive ENO flux calculations. The basic tools come from multiresolution analysis for cell averages on uniform grids, and the principle is that wavelet coefficients can be used for the characterization of local smoothness.. Typically, only few wavelet coefficients are significant. At the finest level, they indicate discontinuity points, where ENO numerical fluxes are computed exactly. Elsewhere, cheaper fluxes can be safely used, or just interpolated from coarser scales. Different applications of this principle have been explored by several authors, see for example (G-Muller and Muller, 1998).Our scheme also uses Ami Harten's ideas. But instead of evolving the cell averages on the finest uniform level, we propose to evolve the cell averages on sparse grids associated with the significant wavelet coefficients. This means that the total number of cells is small, with big cells in smooth regions and smaller ones close to irregularities. This task requires improved new tools, which are described next.

Analysis of the Influence of Sparkout Time on Grinding Using Several Lubrication/Cooling Methods

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Analysis of the Electrical Characteristics of an Alternative Solution for the Brazilian-Amazon Transmission System

This article analyzes the electrical parameters of a 3-phase transmission line using a 280-m-high steel tower that has been proposed for the Amazon transmission system in Brazil. The height of the line conductors and the distance between them are intrinsically related to the longitudinal and transverse parameters of the line. Hence, an accurate study is carried out in order to show the electrical variations between a transmission line using the new technology and a conventional 3-phase 440-kV line, considering a wide range of frequencies and variable soil resistivity. First, a brief review of the fundamental theory of line parameters is presented. In addition, by using a digital line model, simulations are carried out in the time domain to analyze possible and critical over-voltage transients on the proposed line representation.

Numerical approximation of the Ginzburg-Landau equation with memory effects in the dynamics of phase transitions

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Heart rate variability in healthy adolescents at rest

Objective: to describe the cardiac autonomic function in healthy adolescents between 13 and 18 years old. Methods: data from 93 adolescents, of both sexes, were analyzed; they were divided into three groups according to age. Group A 13-14 years old (n=26), Group B 15-16 years old (n = 30) and Group C 17-18 years old (n = 17). The protocol consisted of raising the heart rate, beat by beat for twenty minutes of rest in supine position with spontaneous breathing. The cardiac autonomic behavior analysis was performed by the method of heart rate variability (HRV) through the indices of the Poincaré plot (SD1, SD2 and ratio SD1/SD2) and its qualitative analysis and indices in linear time domain (RMSSD and SDNN) and frequency (LFnu, HFnu and LF/HF). One-way ANOVA test or Kruskal-Wallis test was used for comparison between groups, with a significance level of 5%. Results: there was no significant difference in none of the indices compared the three groups. In addition,visual analysis of Poincaré plot has observed for all age groups large dispersion of RR intervals, indicating that these groups have good HRV. Conclusion: there were no changes in HRV in the different age groups analyzed, however the result allowed to establish a standard for each group that was studied, becomes an important tool for comparison of cardiac autonomic function between healthy and patients subjects or in different areas.

Proposal of a simplified process to correct the phase decoupling using modal analysis

Modal analysis is widely approached in the classic theory of transmission line modeling. This technique is applied to model the three-phase representation of conventional electric systems taking into account their self and mutual electrical parameters. However the methodology has some particularities and inaccuracies for specific applications which are not clearly described in the basic references of this topic. This paper provides a thorough review of modal analysis theory applied to line models followed by an original and simple procedure to overcome the possible errors embedded in the modal decoupling through the three-phase system modeling. © 2012 IEEE.

Development of a simplified transmission line model directly in the phase domain

A transmission line digital model is developed direct in the phase and time domains. The successive modal transformations considered in the three-phase representation are simplified and then the proposed model can be easily applied to several operation condition based only on the previous knowing of the line parameters, without a thorough theoretical knowledge of modal analysis. The proposed model is also developed based on lumped elements, providing a complete current and voltage profile at any point of the transmission system. This model makes possible the modeling of non-linear power devices and electromagnetic phenomena along the transmission line using simple electric circuit components, representing a great advantage when compared to several models based on distributed parameters and inverse transforms. In addition, an efficient integration method is proposed to solve the system of differential equations resulted from the line modeling by lumped elements, thereby making possible simulations of transient and steady state using a wide and constant integration step. © 2012 IEEE.

The use of Gramian matrices for aeroelastic stability analysis

Most of the established procedures for analysis of aeroelastic flutter in the development of aircraft are based on frequency domain methods. Proposing new methodologies in this field is always a challenge, because the new methods need to be validated by many experimental procedures. With the interest for new flight control systems and nonlinear behavior of aeroelastic structures, other strategies may be necessary to complete the analysis of such systems. If the aeroelastic model can be written in time domain, using state-space formulation, for instance, then many of the tools used in stability analysis of dynamic systems may be used to help providing an insight into the aeroelastic phenomenon. In this respect, this paper presents a discussion on the use of Gramian matrices to determine conditions of aeroelastic flutter. The main goal of this work is to introduce how observability gramian matrix can be used to identify the system instability. To explain the approach, the theory is outlined and simulations are carried out on two benchmark problems. Results are compared with classical methods to validate the approach and a reduction of computational time is obtained for the second example. © 2013 Douglas Domingues Bueno et al.