Criteria for mixed grids in computational fluid dynamics

The finite volume method is used as a numerical method for solving the fluid flow equations. This method is appropriate to employ under structured and unstructured meshes. Mixed grids, combining both types of grids, are investigated. The coupling of different grids is done by overlapping strategy. The computational effort for the mixed grid is evaluated by the CPU-time, with different percentage of covering area of the unstructured mesh. The present scheme is tested for the driven cavity problem, where the incompressible fluid is integrated by calculating the velocity fields and computing the pressure field in each time step. Several schemes for unstructured grid are examined, and the compatibility condition is applied to check their consistency. A scheme to verify the compatibility condition for the unstructured grids is presented. (c) 2006 IMACS. Published by Elsevier B.V. All rights reserved.

Transient stability indices from a hybrid approach

This paper deals with hybrid method for transient stability analysis combining time domain simulation and a direct method. Nowadays, the step-by-step simulation is the best available tool for allowing the uses of detailed models and for providing reliable results. The main limitation of this approach involves the large time of computational simulations and the absence of stability margin. On the other hand, direct methods, that demand less CPU time, did not show ample reliability and applicability yet. The best way seems to be using hybrid solutions, in which a direct method is incorporated in a time domain simulation tool. This work has studied a direct method using the transient potential and kinetic energy of the critical machine only. In this paper the critical machine is identified by a fast and efficient method, and the proposal is new for using to get stability margins from hybrid approaches. Results from systems, like 16-machine, show stability indices to dynamic security assessment. © 2001 IEEE.

Simulação numérica do processo de secagem de madeiras

Pós-graduação em Engenharia Mecânica - FEG

Spurious transients of projection methods in microflow simulations

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

Stability analysis of Crank-Nicolson and Euler schemes for time-dependent diffusion equations

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

Utilização da modelagem sistêmica para analisar o efeito da variação climática no parâmetro produtividade da água na laranja Natal (Citrus sinensis L. Osbeck) no Estado de São Paulo

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

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 SIMPLE NUMERICAL-METHOD FOR HYDROSTATIC INCOMPRESSIBLE MODELS WITH RIGID LIDS

A simple and easily implemented method is developed to keep the vertical velocity equal to zero at the bottom and top of hydrostatic incompressible numerical models. The pressure is computed at the top by correcting its value given in the previous time step so that the vertical integral of the horizontal divergence is zero at each column. Numerical experiments that exhibit small time variations of pressure at the top are able to simplify the algorithm and save computer time. Numerical simulations illustrate the method effectiveness for a horizontal deformation-induced frontogenesis.

Modelling study of the impact of deep convection on the utls air composition - Part I: Analysis of ozone precursors

The aim of this work is to study the local impact on the upper troposphere/lower stratosphere air composition of an extreme deep convective system. For this purpose, we performed a simulation of a convective cluster composed of many individual deep convective cells that occurred near Bauru (Brazil). The simulation is performed using the 3-D mesoscale model RAMS coupled on-line with a chemistry model. The comparisons with meteorological measurements show that the model produces meteorological fields generally consistent with the observations. The present paper (part I) is devoted to the analysis of the ozone precursors (CO, NO x and non-methane volatile organic compounds) and HO x in the UTLS. The simulation results show that the distribution of CO with altitude is closely related to the upward convective motions and consecutive outflow at the top of the convective cells leading to a bulge of CO between 7 km altitude and the tropopause (around 17km altitude). The model results for CO are consistent with satellite-borne measurements at 700 hPa. The simulation also indicates enhanced amounts of NO x up to 2 ppbv in the 7-17 km altitude layer mainly produced by the lightning associated with the intense convective activity. For insoluble non-methane volatile organic compounds, the convective activity tends to significantly increase their amount in the 7-17km layer by dynamical effects. During daytime in the presence of lightning NO x, this bulge is largely reduced in the upper part of the layer for reactive species (e.g. isoprene, ethene) because of their reactions with OH that is increased on average during daytime. Lightning NO x also impacts on the oxydizing capacity of the upper troposphere by reducing on average HO x, HO 2, H 2O 2 and organic hydroperoxides. During the simulation time, the impact of convection on the air composition of the lower stratosphere is negligible for all ozone precursors although several of the simulated convective cells nearly reach the tropopause. There is no significant transport from the upper troposphere to the lower stratosphere, the isentropic barrier not being crossed by convection. The impact of the increase of ozone precursors and HO x in the upper troposphere on the ozone budget in the LS is discussed in part II of this series of papers.

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.

Modelling fungus dispersal scenarios using cellular automata

This study focused on representing spatio-temporal patterns of fungal dispersal using cellular automata. Square lattices were used, with each site representing a host for a hypothetical fungus population. Four possible host states were allowed: resistant, permissive, latent or infectious. In this model, the probability of infection for each of the healthy states (permissive or resistant) in a time step was determined as a function of the host's susceptibility, seasonality, and the number of infectious sites and the distance between them. It was also assumed that infected sites become infectious after a pre-specified latency period, and that recovery is not possible. Several scenarios were simulated to understand the contribution of the model's parameters and the spatial structure on the dynamic behaviour of the modelling system. The model showed good capability for representing the spatio-temporal pattern of fungus dispersal over planar surfaces. With a specific problem in mind, the model can be easily modified and used to describe field behaviour, which can contribute to the conservation and development of management strategies for both natural and agricultural systems. © 2012 Elsevier B.V.

Artificial neural network model of discharge lamps in the power quality context

This paper presents a methodology for modeling high intensity discharge lamps based on artificial neural networks. The methodology provides a model which is able to represent the device operating in the frequency of distribution systems, facing events related to power quality. With the aid of a data acquisition system to monitor the laboratory experiment, and using $$\text{ MATLAB }^{\textregistered }$$ software, data was obtained for the training of two neural networks. These neural networks, working together, were able to represent with high fidelity the behavior of a discharge lamp. The excellent performance obtained by these models allowed the simulation of a group of lamps in a distribution system with shorter simulation time when compared to mathematical models. This fact justified the application of this family of loads in electric power systems. The representation of the device facing power quality disturbances also proved to be a useful tool for more complex studies in distribution systems. © 2013 Brazilian Society for Automatics - SBA.

Fenômenos não-lineares, incluindo-se os não-ideias, em captura de energia utilizando-se dispositivos piezoelétricos

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

Simulação do enovelamento de proteínas desnaturadas utilizando o método de crescimento de cadeias em rede tetraédrica

Pós-graduação em Biofísica Molecular - IBILCE

Substâncias húmicas: interações com nutrientes e contaminantes

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