Some experiments with high order compact methods using a computer algebra software - Part I

We consider a procedure for obtaining a compact fourth order method to the steady 2D Navier-Stokes equations in the streamfunction formulation using the computer algebra system Maple. The resulting code is short and from it we obtain the Fortran program for the method. To test the procedure we have solved many cavity-type problems which include one with an analytical solution and the results are compared with results obtained by second order central differences to moderate Reynolds numbers. (c) 2005 Elsevier B.V. All rights reserved.

Some practical advice on polynomial regression analysis from blocked response surface designs

It is often necessary to run response surface designs in blocks. In this paper the analysis of data from such experiments, using polynomial regression models, is discussed. The definition and estimation of pure error in blocked designs are considered. It is recommended that pure error is estimated by assuming additive block and treatment effects, as this is more consistent with designs without blocking. The recovery of inter-block information using REML analysis is discussed, although it is shown that it has very little impact if thc design is nearly orthogonally blocked. Finally prediction from blocked designs is considered and it is shown that prediction of many quantities of interest is much simpler than prediction of the response itself.

A study of a numerical solution of the steady two dimensions Navier-Stokes equations in a constricted channel problem by a compact fourth order method

We present a numerical solution for the steady 2D Navier-Stokes equations using a fourth order compact-type method. The geometry of the problem is a constricted symmetric channel, where the boundary can be varied, via a parameter, from a smooth constriction to one possessing a very sharp but smooth corner allowing us to analyse the behaviour of the errors when the solution is smooth or near singular. The set of non-linear equations is solved by the Newton method. Results have been obtained for Reynolds number up to 500. Estimates of the errors incurred have shown that the results are accurate and better than those of the corresponding second order method. (C) 2002 Elsevier B.V. All rights reserved.

Simultaneous recording of shadowgrams and schlieren images for the study of plasma assisted aerodynamics

This work describes an optical device for the simultaneous recording of shadowgrams and schlieren images, and some results are presented concerning its applications to the study of plasma assisted flow control in airfoil models. This approach offers many advantages in comparison to other methods, specially because the use of tracer particles (like smoke in wind tunnels) is not required for the experiments, thus avoiding contaminations in the electric discharges or air flows. Besides, while schlieren images reveal the refractive index gradients in the area of study, shadowgrams detect the second order spatial derivatives of the refractive indexes. Therefore, the simultaneous recording of these different images may give interesting information about the phenomena under study. In this paper, these images were used to confirm the existence of vortex structures in the flow induced by corona discharges on airfoil models. These structures are a possible explanation for the effects of drag reduction and lift force increasing, which have been reported in experiments of plasma assisted Aerodynamics.

Simultaneous Recording of Shadowgrams and Schlieren Images for the Study of Plasma Assisted Aerodynamics

This work describes an optical device for the simultaneous recording of shadowgrams and schlieren images, and some results are presented concerning its application to the study of plasma assisted flow control in airfoil models. This approach offers many advantages in comparison to other methods, specially because the use of tracer particles (like smoke in wind tunnels) is not required for the experiments, thus avoiding contaminations in the electric discharges or air flows. Besides, while schlieren images reveal the refractive index gradients in the area of study, shadowgrams detect the second order spatial derivatives of the refractive indexes. Therefore, the simultaneous recording of these different images may give interesting information about the phenomena under study.

Critical temperature for first-order phase transitions in confined systems

We consider the Euclidean D-dimensional -lambda vertical bar phi vertical bar(4)+eta vertical bar rho vertical bar(6) (lambda,eta > 0) model with d (d <= D) compactified dimensions. Introducing temperature by means of the Ginzburg-Landau prescription in the mass term of the Hamiltonian, this model can be interpreted as describing a first-order phase transition for a system in a region of the D-dimensional space, limited by d pairs of parallel planes, orthogonal to the coordinates axis x(1), x(2),..., x(d). The planes in each pair are separated by distances L-1, L-2, ... , L-d. We obtain an expression for the transition temperature as a function of the size of the system, T-c({L-i}), i = 1, 2, ..., d. For D = 3 we particularize this formula, taking L-1 = L-2 = ... = L-d = L for the physically interesting cases d = 1 (a film), d = 2 (an infinitely long wire having a square cross-section), and for d = 3 (a cube). For completeness, the corresponding formulas for second-order transitions are also presented. Comparison with experimental data for superconducting films and wires shows qualitative agreement with our theoretical expressions.

Thermal photon production in Au plus Au collisions: Viscous corrections in two different hydrodynamic formalisms

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

Shear viscosity from thermal fluctuations in relativistic conformal fluid dynamics

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

Gauge formulation for higher order gravity

This work is an application of the second order gauge theory for the Lorentz group, where a description of the gravitational interaction is obtained that includes derivatives of the curvature. We analyze the form of the second field strength, G=partial derivative F+fAF, in terms of geometrical variables. All possible independent Lagrangians constructed with quadratic contractions of F and quadratic contractions of G are analyzed. The equations of motion for a particular Lagrangian, which is analogous to Podolsky's term of his generalized electrodynamics, are calculated. The static isotropic solution in the linear approximation was found, exhibiting the regular Newtonian behavior at short distances as well as a meso-large distance modification.

Remote sensing analysis of depositional landforms in alluvial settings: Method development and application to the Taquari megafan, Pantanal (Brazil)

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

Non-Hermitian multichannel theory of ultracold collisions modified by intense light fields tuned to the red of the trapping transition

We develop a systematic scheme to treat binary collisions between ultracold atoms in the presence of a strong laser field, tuned to the red of the trapping transition. We assume that the Rabi frequency is much less than the spacing between adjacent bound-state resonances, In this approach we neglect fine and hyperfine structures, but consider fully the three-dimensional aspects of the scattering process, up to the partial d wave. We apply the scheme to calculate the S matrix elements up to the second order in the ratio between the Rabi frequency and the laser detuning, We also obtain, fur this simplified multichannel model, the asymmetric line shapes of photoassociation spectroscopy, and the modification of the scattering length due to the light field at low, but finite, entrance kinetic energy. We emphasize that the present calculations can be generalized to treat more realistic models, and suggest how to carry out a thorough numerical comparison to this semianalytic theory. [S1050-2947(98)04902-6].

Viscoelastic behavior of persimmons dried at constant air temperature

The viscoelastic behavior of dried persimmons at different air-drying temperatures and velocities was evaluated. Air temperatures and velocities were varied according to a second-order central composite design, with temperature ranging from 40degreesC to 70degreesC and air velocity from 0.8 to 2.0 m/s. After drying, persimmons were equilibrated at four different water activities: 0.432, 0.576, 0.625 and 0.751. The rheological behavior of dried and conditioned persimmons was studied under uniaxial compression-relaxation tests. Three different rheological models were fitted to the experimental relaxation curves: Maxwell, Generalized Maxwell and Peleg and Normand. Based on the root mean square of residuals, the Generalized Maxwell model showed the best fit and a regression analysis was applied to obtain response surfaces for the model parameters. The dependence of the rheological properties on water activity was also analysed. Results showed that only the linear effect of air temperature was significant at a 5% level on the equilibrium stress and relaxation times. In a general way, these parameters increased with increasing air temperature and decreasing water activity. (C) 2004 Swiss Society of Food Science and Technology. Published by Elsevier Ltd. All rights reserved.

Non-stationary frequency analysis of extreme daily rainfall in São Paulo, Brazil

This work is an assessment of frequency of extreme values (EVs) of daily rainfall in the city of São Paulo. Brazil, over the period 1933-2005, based on the peaks-over-threshold (POT) and Generalized Pareto Distribution (GPD) approach. Usually. a GPD model is fitted to a sample of POT Values Selected With a constant threshold. However. in this work we use time-dependent thresholds, composed of relatively large p quantities (for example p of 0.97) of daily rainfall amounts computed from all available data. Samples of POT values were extracted with several Values of p. Four different GPD models (GPD-1, GPD-2, GPD-3. and GDP-4) were fitted to each one of these samples by the maximum likelihood (ML) method. The shape parameter was assumed constant for the four models, but time-varying covariates were incorporated into scale parameter of GPD-2. GPD-3, and GPD-4, describing annual cycle in GPD-2. linear trend in GPD-3, and both annual cycle and linear trend in GPD-4. The GPD-1 with constant scale and shape parameters is the simplest model. For identification of the best model among the four models WC used rescaled Akaike Information Criterion (AIC) with second-order bias correction. This criterion isolates GPD-3 as the best model, i.e. the one with positive linear trend in the scale parameter. The slope of this trend is significant compared to the null hypothesis of no trend, for about 98% confidence level. The non-parametric Mann-Kendall test also showed presence of positive trend in the annual frequency of excess over high thresholds. with p-value being virtually zero. Therefore. there is strong evidence that high quantiles of daily rainfall in the city of São Paulo have been increasing in magnitude and frequency over time. For example. 0.99 quantiles of daily rainfall amount have increased by about 40 mm between 1933 and 2005. Copyright (C) 2008 Royal Meteorological Society

Use of Biopolymeric Membranes for Adsorption of Paraquat Herbicide from Water

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

Fourth-order method for solving the Navier-Stokes equations in a constricting channel

A fourth-order numerical method for solving the Navier-Stokes equations in streamfunction/vorticity formulation on a two-dimensional non-uniform orthogonal grid has been tested on the fluid flow in a constricted symmetric channel. The family of grids is generated algebraically using a conformal transformation followed by a non-uniform stretching of the mesh cells in which the shape of the channel boundary can vary from a smooth constriction to one which one possesses a very sharp but smooth corner. The generality of the grids allows the use of long channels upstream and downstream as well as having a refined grid near the sharp corner. Derivatives in the governing equations are replaced by fourth-order central differences and the vorticity is eliminated, either before or after the discretization, to form a wide difference molecule for the streamfunction. Extra boundary conditions, necessary for wide-molecule methods, are supplied by a procedure proposed by Henshaw et al. The ensuing set of non-linear equations is solved using Newton iteration. Results have been obtained for Reynolds numbers up to 250 for three constrictions, the first being smooth, the second having a moderately sharp corner and the third with a very sharp corner. Estimates of the error incurred show that the results are very accurate and substantially better than those of the corresponding second-order method. The observed order of the method has been shown to be close to four, demonstrating that the method is genuinely fourth-order. © 1977 John Wiley & Sons, Ltd.