A generalized alternating-direction implicit scheme for incompressible magnetohydrodynamic viscous flows at low magnetic Reynolds number

This paper presents numerical simulations of incompressible fluid flows in the presence of a magnetic field at low magnetic Reynolds number. The equations governing the flow are the Navier-Stokes equations of fluid motion coupled with Maxwell's equations of electromagnetics. The study of fluid flows under the influence of a magnetic field and with no free electric charges or electric fields is known as magnetohydrodynamics. The magnetohydrodynamics approximation is considered for the formulation of the non-dimensional problem and for the characterization of similarity parameters. A finite-difference technique is used to discretize the equations. In particular, an extension of the generalized Peaceman and Rachford alternating-direction implicit (ADI) scheme for simulating two-dimensional fluid flows is presented. The discretized conservation equations are solved in stream function-vorticity formulation. We compare the ADI and generalized ADI schemes, and show that the latter is more efficient in simulating low Reynolds number and magnetic Reynolds number problems. Numerical results demonstrating the applicability of this technique are also presented. The simulation of incompressible magneto hydrodynamic fluid flows is illustrated by numerical solution for two-dimensional cases. (c) 2007 Elsevier B.V. All rights reserved.

Modelagem matemática e simulação numérica do resfriamento rápido de morango com ar forçado

Este trabalho abordou o resfriamento rápido com ar forçado de morango via simulação numérica. Para tanto, foi empregado o modelo matemático que descreve o processo de transferência de calor, com base na lei de Fourier, escrito em coordenadas esféricas e simplificado para descrever o processo unidimensional. A resolução da equação expressa pelo modelo matemático deu-se por meio da implementação de um algoritmo, fundamentado no esquema explícito do método numérico das diferenças finitas, executado no ambiente de computação científica MATLAB 6.1. A validação do modelo matemático foi realizada a partir da comparação de dados teóricos com dados obtidos num experimento, no qual morangos foram resfriados com ar forçado. Os resultados mostraram que esse tipo de investigação para a determinação do coeficiente de transferência de calor por convecção é promissora como ferramenta no suporte à decisão do uso ou desenvolvimento de equipamentos na área de resfriamento rápido de frutos esféricos com ar forçado.

Influence of drying conditions on ascorbic acid during convective drying of whole persimmons

Degradation kinetics of food constituents may be related to the matrix molecular mobility by glass transition temperature. Our objective was to test this approach to describe ascorbic acid degradation during drying of persimmons in an automatically controlled tray dryer with temperatures (40 to 70 degrees C) and air velocities (0.8 to 2.0 m/s) varying according a second order central composite design. The Williams-Landel-Ferry model was satisfactorily adjusted to degradation curves for both control strategies adopted-constant air temperature and temperature fixed inside the fruit. Degradation rates were higher at higher drying temperatures, independent of the necessary time to attain the desired moisture content.

A critical investigation of the Tanford-Kirkwood scheme by means of Monte Carlo simulations

Monte Carlo simulations are used to assess the adequacy of the Tanford-Kirkwood prescription for electrostatic interactions in macromolecules. Within a continuum dielectric framework, the approach accurately describes salt screening of electrostatic interactions for moderately charged systems consistent with common proteins at physiological conditions. The limitations of the Debye-Huckel theory, which forms the statistical mechanical basis for the Tanford-Kirkwood result, become apparent for highly charged systems. It is shown, both by an analysis of the Debye-Huckel theory and by numerical simulations, that the difference in dielectric permittivity between macromolecule and surrounding solvent does not play a significant role for salt effects if the macromolecule is highly charged. By comparison to experimental data, the continuum dielectric model (combined with either an approximate effective Hamiltonian as in the Tanford-Kirkwood treatment or with exact Monte Carlo simulations) satisfactorily predicts the effects of charge mutation on metal ion binding constants, but only if the macromolecule and solvent are assigned the same or similar permittivities.

A case study of convective organization into precipitating lines in the Southwest Amazon during the WETAMC and TRMM-LBA

A case study of convective development in the Southwest Amazon region during the Wet Season Atmospheric Mesoscale Campaign (WETAMC) and Tropical Rainfall Measuring Mission (TRMM)/Large-Scale Biosphere-Atmosphere (LBA) Experiment in Amazonia is presented. The convective development during 7 February 1999 is shown to occur during a period of very weak large-scale forcing in the presence of topography and deforestation. The available data include dual Doppler radar analysis, radiosonde launches, and surface and boundary layer observations. The observational analysis is complemented with a series of model simulations using the RAMS with 2-km resolution over a 300 km 300 km area forced by a morning radiosonde profile. A comparison of the observed and simulated thermodynamic transformation of the boundary layer and of the formation of convective lines, and of their kinematic and microphysical properties is presented. It is shown that only a few very deep and intense convective cells are necessary to explain the overall precipitating line formation and that discrete propagation and coupling with upper atmosphere circulations may explain the appearance of several lines. The numerical simulation indicates that topography may be the cause of initial convective development, although later on the convective line is parallel to the midlevel shear. There are indications that small-scale deforestation may have an effect on increasing rainfall in the wet season when the large-scale forcing is very weak.

INFLUENCE of THE EVAPORATOR POSITION IN STATIC DOMESTIC REFRIGERATORS

The temperature and velocity distributions of the air inside the cabinet of domestic refrigerators affect the quality of food products. If the consumer knows the location of warm and cold zones in the refrigerator, the products can be placed in the right zone. In addition, the knowledge of the thickness of thermal and hydrodynamic boundary layers near the evaporator and the other walls is also important. If the product is too close to the evaporator wall, freezing can occur, and if it is too close to warm walls, the products can be deteriorated. The aim of the present work is to develop a steady state computational fluid dynamics (CFD) model for domestic refrigerators working on natural convection regime. The Finite Volume Methodology is chosen as numerical procedure for discretizing the governing equations. The SIMPLE-Semi-Implicit Method for Pressure-Linked Equations algorithm applied to a staggered mesh was used for solving the pressure-velocity coupling problem. The Power-Law scheme is employed as interpolation function for the convective-diffusive terms, and the TDMA-Tri-Diagonal Matrix Algorithm is used to solve the systems of algebraic equations. The model is applied to a commercial static refrigerator, where the cabinet is considered an empty three-dimensional rectangular cavity with one drawer at the bottom of the cabinet, but without shelves. In order to analyze the velocity and temperature fields of the air flow inside the cabinet the evaporator temperature, Te, was varied from -20 degrees C to 0 degrees C, and nine different evaporator positions are evaluated for evaporator temperature of -15 degrees C. The cooling capacity of the evaporator for the steady state regime is also computed for each case. One can conclude that the vertical positioning of the evaporator inside the cabinet plays an important role on the temperature distribution inside the cabinet.

Convective storm initiation in a moist tropical environment

Radar and satellite data from the Tropical Rainfall Measuring Mission-Large-Scale Biosphere-Atmosphere (TRMM-LBA) project have been examined to determine causes for convective storm initiation in the southwest Amazon region. The locations and times of storm initiation were based on the National Center for Atmospheric Research (NCAR) S-band dual-polarization Doppler radar (S-Pol). Both the radar and the Geostationary Operational Environmental Satellite-8 (GOES-8) visible data were used to identify cold pools produced by convective precipitation. These data along with high-resolution topographic data were used to determine possible convective storm triggering mechanisms. The terrain elevation varied from 100 to 600 m. Tropical forests cover the area with numerous clear-cut areas used for cattle grazing and farming. This paper presents the results from 5 February 1999. A total of 315 storms were initiated within 130 km of the S-Pol radar. This day was classified as a weak monsoon regime where convection developed in response to the diurnal cycle of solar heating. Scattered shallow cumulus during the morning developed into deep convection by early afternoon. Storm initiation began about 1100 LST and peaked around 1500-1600 LST. The causes of storm initiation were classified into four categories. The most common initiation mechanism was caused by forced lifting by a gust front (GF; 36%). Forcing by terrain (>300 m) without any other triggering mechanism accounted for 21% of the initiations and colliding GFs accounted for 16%. For the remaining 27% a triggering mechanism was not identified. Examination of all days during TRMM-LBA showed that this one detailed study day was representative of many days. A conceptual model of storm initiation and evolution is presented. The results of this study should have implications for other locations when synoptic-scale forcing mechanisms are at a minimum. These results should also have implications for very short-period forecasting techniques in any location where terrain, GFs, and colliding boundaries influence storm evolution.

Exact dynamics of single-qubit-gate fidelities under the measurement-based quantum computation scheme

Measurement-based quantum computation is an efficient model to perform universal computation. Nevertheless, theoretical questions have been raised, mainly with respect to realistic noise conditions. In order to shed some light on this issue, we evaluate the exact dynamics of some single-qubit-gate fidelities using the measurement-based quantum computation scheme when the qubits which are used as a resource interact with a common dephasing environment. We report a necessary condition for the fidelity dynamics of a general pure N-qubit state, interacting with this type of error channel, to present an oscillatory behavior, and we show that for the initial canonical cluster state, the fidelity oscillates as a function of time. This state fidelity oscillatory behavior brings significant variations to the values of the computational results of a generic gate acting on that state depending on the instants we choose to apply our set of projective measurements. As we shall see, considering some specific gates that are frequently found in the literature, the fast application of the set of projective measurements does not necessarily imply high gate fidelity, and likewise the slow application thereof does not necessarily imply low gate fidelity. Our condition for the occurrence of the fidelity oscillatory behavior shows that the oscillation presented by the cluster state is due exclusively to its initial geometry. Other states that can be used as resources for measurement-based quantum computation can present the same initial geometrical condition. Therefore, it is very important for the present scheme to know when the fidelity of a particular resource state will oscillate in time and, if this is the case, what are the best times to perform the measurements.

Transient non-Darcy forced convective heat transfer from a flat plate embedded in a fluid-saturated porous medium

Transient non-Darcy forced convection on a flat plate embedded in a porous medium is investigated using the Forchheimer-extended Darcy law. A sudden uniform pressure gradient is applied along the flat plate, and at the same time, its wall temperature is suddenly raised to a high temperature. Both the momentum and energy equations are solved by retaining the unsteady terms. An exact velocity solution is obtained and substituted into the energy equation, which then is solved by means of a quasi-similarity transformation. The temperature field can be divided into the one-dimensional transient (downstream) region and the quasi-steady-state (upstream) region. Thus the transient local heat transfer coefficient can be described by connecting the quasi-steady-state solution and the one-dimensional transient solution. The non-Darcy porous inertia works to decrease the velocity level and the time required for reaching the steady-state velocity level. The porous-medium inertia delays covering of the plate by the steady-state thermal boundary layer. © 1990.

Statistics on tropical convective storms observed by radar

Convective storm rainfall is of special importance to urban hydrological studies due to its temporal and spatial variability. Although dense networks of recording rain gauges can be employed to characterize such rainfall, very few investigations of this type have been undertaken due to their prohibitive cost. This paper reports some data on characteristics of tropical convective storms obtained from radar at Bauru in the State of São Paulo, Brazil. Periods of convective precipitation were identified by exclusion of those related to frontal activity with the help of synoptic maps and the radar screen record. The occurrence and evolution of convective storms were observed in two 28 km × 28 km windows obtaining information on the life history of convective cells and the magnitude of rainfall. Frequency distributions of the time of occurrence of convective rainfall, cell size, area covered, life duration and maximum and average rainfall observed in the experimental areas are presented and discussed.

Feynman integrals with tensorial structure in the negative dimensional integration scheme

The negative-dimensional integration method (NDIM) is revealing itself as a very useful technique for computing massless and/or massive Feynman integrals, covariant and noncovanant alike. Up until now however, the illustrative calculations done using such method have been mostly covariant scalar integrals/without numerator factors. We show here how those integrals with tensorial structures also can be handled straightforwardly and easily. However, contrary to the absence of significant features in the usual approach, here the NDIM also allows us to come across surprising unsuspected bonuses. Toward this end, we present two alternative ways of working out the integrals and illustrate them by taking the easiest Feynman integrals in this category that emerge in the computation of a standard one-loop self-energy diagram. One of the novel and heretofore unsuspected bonuses is that there are degeneracies in the way one can express the final result for the referred Feynman integral.

An active leakage-injection scheme applied to low-voltage SRAMs

An active leakage-injection scheme (ALIS) for low-voltage (LV) high-density (HD) SRAMs is presented. By means of a feedback loop comprising a servo-amplifier and a common-drain MOSFET, a current matching the respective bit-line leakage is injected onto the line during precharge and sensing, preventing the respective capacitances from erroneous discharges. The technique is able to handle leakages up to hundreds of μA at high operating temperatures. Since no additional timing is required, read-out operations are performed at no speed penalty. A simplified 256×1bit array was designed in accordance with a 0.35 CMOS process and 1.2V-supply. A range of PSPICE simulation attests the efficacy of ALIS. With an extra power consumption of 242 μW, a 200 μA-leakage @125°C, corresponding to 13.6 times the cell current, is compensated.

Analysis of a sub-optimal scheme of drug dosage in the AIDS treatment

Here the results for CD4+T cells count and the viral load obtained from HIV sero-positive patients are compared with results from numerical simulations by computer. Also, the standard scheme of administration of drugs anti HIV (HAART schemes) which uses constant doses is compared with an alternative sub-optimal teatment scheme which uses variable drug dosage according to the evolution of a quantitative measure of the side effects. The quantitative analysis done here shows that it is possible to obtain, using the alternative scheme, the same performance of actual data but using variable dosage and having fewer side effects. Optimal control theory is used to solve and also to provide a prognosis related to the strategies for control of viraemia.

Gluon mass generation in the PT-BFM scheme

In this article we study the general structure and special properties of the Schwinger-Dyson equation for the gluon propagator constructed with the pinch technique, together with the question of how to obtain infrared finite solutions, associated with the generation of an effective gluon mass. Exploiting the known all-order correspondence between the pinch technique and the background field method, we demonstrate that, contrary to the standard formulation, the non-perturbative gluon self-energy is transverse order-by-order in the dressed loop expansion, and separately for gluonic and ghost contributions. We next present a comprehensive review of several subtle issues relevant to the search of infrared finite solutions, paying particular attention to the role of the seagull graph in enforcing transversality, the necessity of introducing massless poles in the three-gluon vertex, and the incorporation of the correct renormalization group properties. In addition, we present a method for regulating the seagull-type contributions based on dimensional regularization; its applicability depends crucially on the asymptotic behavior of the solutions in the deep ultraviolet, and in particular on the anomalous dimension of the dynamically generated gluon mass. A linearized version of the truncated Schwinger-Dyson equation is derived, using a vertex that satisfies the required Ward identity and contains massless poles belonging to different Lorentz structures. The resulting integral equation is then solved numerically, the infrared and ultraviolet properties of the obtained solutions are examined in detail, and the allowed range for the effective gluon mass is determined. Various open questions and possible connections with different approaches in the literature are discussed. © SISSA 2006.