Using state-space techniques to represent frequency dependent single-phase lines directly in time domain

This paper proposes to use a state-space technique to represent a frequency dependent line for simulating electromagnetic transients directly in time domain. The distributed nature of the line is represented by a multiple 1t section network made up of the lumped parameters and the frequency dependence of the per unit longitudinal parameters is matched by using a rational function. The rational function is represented by its equivalent circuit with passive elements. This passive circuit is then inserted in each 1t circuit of the cascade that represents the line. Because the system is very sparse, it is possible to use a sparsity technique to store only nonzero elements of this matrix for saving space and running time. The model was used to simulate the energization process of a 10 km length single-phase line. ©2008 IEEE.

Alternate treatments of jacobian singularities in polar coordinates within finite-difference schemes

Jacobian singularities of differential operators in curvilinear coordinates occur when the Jacobian determinant of the curvilinear-to-Cartesian mapping vanishes, thus leading to unbounded coefficients in partial differential equations. Within a finite-difference scheme, we treat the singularity at the pole of polar coordinates by setting up complementary equations. Such equations are obtained by either integral or smoothness conditions. They are assessed by application to analytically solvable steady-state heat-conduction problems.

A Time-domain Harmonic Power-Flow Analysis in Electrical Energy Distribution Networks, Using Norton Models for Non-linear Loading

The purpose of this paper is to present the application of a three-phase harmonic propagation analysis time-domain tool, using the Norton model to approach the modeling of non-linear loads, making the harmonics currents flow more appropriate to the operation analysis and to the influence of mitigation elements analysis. This software makes it possible to obtain results closer to the real distribution network, considering voltages unbalances, currents imbalances and the application of mitigation elements for harmonic distortions. In this scenario, a real case study with network data and equipments connected to the network will be presented, as well as the modeling of non-linear loads based on real data obtained from some PCCs (Points of Common Coupling) of interests for a distribution company.

Comparison of frequency and time domain electromyography parameters in women with patellofemoral pain

Despite its high incidence, patellofemoral pain etiology remains unclear. No prior study has compared surface electromyography frequency domain parameters and surface electromyography time domain variables, which have been used as a classic analysis of patellofemoral pain. Thirty one women with patellofemoral pain and twenty eight pain-free women were recruited. Each participant was asked to descend a seven step staircase and data from five successful trials were collected. During the task, the vastus medialis and vastus lateralis muscle activities were monitored by surface electromyography. The data were processed and analyzed in four variables of the frequency domain (median frequency, low, medium and high frequency bands) and three time domain variables (Automatic, Cross-correlation and Visual Onset between the vastus medialis and vastus lateralis muscles). Reliability, Receiver Operating Characteristic curves and regression models were performed. The medium frequency band was the most reliable variable and different between the groups for both muscles, also demonstrated the best values of sensitivity and sensibility, 72% and 69% for the vastus medialis and 68% and 62% for the vastus lateralis, respectively. The frequency variables predicted the pain of individuals with patellofemoral pain, 26% for the vastus medialis and 20% for the vastus lateralis, being better than the time variables, which achieved only 7%. The frequency domain parameters presented greater reliability, diagnostic accuracy and capacity to predict pain than the time domain variables during stair descent and might be a useful tool to diagnose individuals with patellofemoral pain.

Avaliação do consumo de água no sistema de irrigação por inundação com distribuição de múltiplas entradas na cultura de arroz (Oryza sativa L.)

Pós-graduação em Agronomia (Irrigação e Drenagem) - FCA

Numerical assessment of mass conservation on a MAC-type method for viscoelastic free surface flows

A numerical study of mass conservation of MAC-type methods is presented, for viscoelastic free-surface flows. We use an implicit formulation which allows for greater time steps, and therefore time marching schemes for advecting the free surface marker particles have to be accurate in order to preserve the good mass conservation properties of this methodology. We then present an improvement by using a Runge-Kutta scheme coupled with a local linear extrapolation on the free surface. A thorough study of the viscoelastic impacting drop problem, for both Oldroyd-B and XPP fluid models, is presented, investigating the influence of timestep, grid spacing and other model parameters to the overall mass conservation of the method. Furthermore, an unsteady fountain flow is also simulated to illustrate the low mass conservation error obtained.

On determining spectral parameters, tracking jitter, and GPS positioning improvement by scintillation mitigation

A method of determining spectral parameters p (slope of the phase PSD) and T (phase PSD at 1 Hz) and hence tracking error variance in a GPS receiver PLL from just amplitude and phase scintillation indices and an estimated value of the Fresnel frequency has been previously presented. Here this method is validated using 50 Hz GPS phase and amplitude data from high latitude receivers in northern Norway and Svalbard. This has been done both using (1) a Fresnel frequency estimated using the amplitude PSD (in order to check the accuracy of the method) and (2) a constant assumed value of Fresnel frequency for the data set, convenient for the situation when contemporaneous phase PSDs are not available. Both of the spectral parameters (p, T) calculated using this method are in quite good agreement with those obtained by direct measurements of the phase spectrum as are tracking jitter variances determined for GPS receiver PLLs using these values. For the Svalbard data set, a significant difference in the scintillation level observed on the paths from different satellites received simultaneously was noted. Then, it is shown that the accuracy of relative GPS positioning can be improved by use of the tracking jitter variance in weighting the measurements from each satellite used in the positioning estimation. This has significant advantages for scintillation mitigation, particularly since the method can be accomplished utilizing only time domain measurements thus obviating the need for the phase PSDs in order to extract the spectral parameters required for tracking jitter determination.

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.

A combination of implicit and adaptative upwind tools for the numerical solution of incompressible free surface flows

This paper is concerned with the numerical solutions of time dependent two-dimensional incompressible flows. By using the primitive variables of velocity and pressure, the Navier-Stokes and mass conservation equations are solved by a semi-implicit finite difference projection method. A new bounded higher order upwind convection scheme is employed to deal with the non-linear (advective) terms. The procedure is an adaptation of the GENSMAC (J. Comput. Phys. 1994; 110: 171-186) methodology for calculating confined and free surface fluid flows at both low and high Reynolds numbers. The calculations were performed by using the 2D version of the Freeflow simulation system (J. Comp. Visual. Science 2000; 2:199-210). In order to demonstrate the capabilities of the numerical method, various test cases are presented. These are the fully developed flow in a channel, the flow over a backward facing step, the die-swell problem, the broken dam flow, and an impinging jet onto a flat plate. The numerical results compare favourably with the experimental data and the analytical solutions. Copyright (c) 2006 John Wiley & Sons, Ltd.

A high dynamic range method for the direct readout of a dynamic phase change in homodyne interferometers

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

Mechanistic approach to predict real machining time for milling free-form geometries applying high feed rate

An accurate estimate of machining time is very important for predicting delivery time, manufacturing costs, and also to help production process planning. Most commercial CAM software systems estimate the machining time in milling operations simply by dividing the entire tool path length by the programmed feed rate. This time estimate differs drastically from the real process time because the feed rate is not always constant due to machine and computer numerical controlled (CNC) limitations. This study presents a practical mechanistic method for milling time estimation when machining free-form geometries. The method considers a variable called machine response time (MRT) which characterizes the real CNC machine's capacity to move in high feed rates in free-form geometries. MRT is a global performance feature which can be obtained for any type of CNC machine configuration by carrying out a simple test. For validating the methodology, a workpiece was used to generate NC programs for five different types of CNC machines. A practical industrial case study was also carried out to validate the method. The results indicated that MRT, and consequently, the real machining time, depends on the CNC machine's potential: furthermore, the greater MRT, the larger the difference between predicted milling time and real milling time. The proposed method achieved an error range from 0.3% to 12% of the real machining time, whereas the CAM estimation achieved from 211% to 1244% error. The MRT-based process is also suggested as an instrument for helping in machine tool benchmarking.

Energy method to calculate the density of liquids using ultrasonic reflection techniques

In this work it is introduced a new approach to calculate the density of liquids in terms of the energies of the acoustic signals. This method is compared to other methods in the time domain (peak-to-peak amplitudes) and frequency domain magnitudes at a single frequency. It is used a measurement cell based on a multiple reflection technique, and it is developed an acoustic model for the cell. Simulations and experiments using several liquids are presented, showing that the energy method a less sensitive to noise than the other techniques. The relative errors in the density are smaller than 0.2% when compared to the values measured with a pycnometer.

Harmonic content identification based on neural method for single phase power systems

An alternative method is presented in this paper to identify the harmonic components of non-linear loads in single phase power systems based on artificial neural networks. The components are identified by analyzing the single phase current waveform in time domain in half-cycle of the ac voltage source. The proposed method is compared to the fast Fourier transform. Simulation and experimental results are presented to validate the proposed approach.

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.