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.

An algorithm for arranging response surface designs in small blocks

We consider the problem of blocking response surface designs when the block sizes are prespecified to control variation efficiently and the treatment set is chosen independently of the block structure. We show how the loss of information due to blocking is related to scores defined by Mead and present an interchange algorithm based on scores to improve a given blocked design. Examples illustrating the performance of the algorithm are given and some comparisons with other designs are made. (C) 2000 Elsevier B.V. B.V. All rights reserved.

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.

d(-)-Lactic Acid Production by Leuconostoc mesenteroides B512 Using Different Carbon and Nitrogen Sources

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

Optimization of xylanase production by Bacillus circulans D1 in submerged fermentation using response surface methodology

In this work, a 3(3) factorial design was performed with the aim of optimizing the culture conditions for xylanase production by an alkalophilic thermophilic strain of Bacillus circulans, using response surface methodology. The variables involved in this study were xylan concentration (X-1), pH (X-2) and cultivation time (X-3). The optimal response region was approached without using paths of steepest ascent. Statistical analysis of results showed that, in the range studied, only pH did not have a significant effect on xylanase production. A second-order model was proposed to represent the enzymic activity as a function of xylan concentration (X-1) and cultivation time (X-3). The optimum xylan concentration and cultivation time were 5 g/l and 48 h, respectively. Under these conditions, the model predicted a xylanase activity of 19.1 U/ml. (C) 2002 Elsevier B.V. Ltd. All rights reserved.

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.

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.

Ag-rich phase formation in the Cu-6 wt% Al alloy with Ag additions

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

PHASE-TRANSITION IN A Q-DEFORMED LIPKIN MODEL

Assuming q-deformed commutation relations for the fermions, an extension of the standard Lipkin Hamiltonian is presented. The usual quasi-spin representation of the standard Lipkin model is also obtained in this q-deformed framework. A variationally obtained energy functional is used to analyse the phase transition associated with the spherical symmetry breaking. The only phase transitions in this q-deformed model are of second order. As an outcome of this analysis a critical parameter is obtained which is dependent on the deformation of the algebra and on the number of particles.

Nonsmooth continuous-time optimization problems: Sufficient conditions

We discuss sufficient conditions of optimality for nonsmooth continuous-time nonlinear optimization problems under generalized convexity assumptions. These include both first-order and second-order criteria. (C) 1998 Academic Press.

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.

Automation and cybernetics: Control of a flexible one-link manipulator

States that control is of the essence in cybernetics. Summarizes the dynamic equations for a flexible one-link manipulator moving in the horizontal plane. Employs the finite element method, based on elementary beam theory, during the process of formulation. Develops and instruments a one-link flexible manipulator in order to control its vibration modes. Uses a simple second-order vibration model which permits vibrations on the rod to be estimated using the hub angle. The validation of the dynamic model and the structural analysis of the flexible manipulator is reached using proper infrared cameras and active light sources for determining actual positions of objects in space. Shows that the performance of the control is satisfactory, even under perturbation action.

An insight into the theoretical investigation of possible piezoelectric effect in samarium titanate (SmTiO3)

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Disulfide Biochemistry in 2-Cys Peroxiredoxin: Requirement of Glu50 and Arg146 for the Reduction of Yeast Tsa1 by Thioredoxin

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

Quantum chemical modeling of perovskite: An investigation of piezoelectricity in ferrite of yttrium

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)