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.

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.

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 experiments with high order compact methods using a computer algebra software - Part II (non-uniform grid)

We generalize a procedure proposed by Mancera and Hunt [P.F.A. Mancera, R. Hunt, Some experiments with high order compact methods using a computer algebra software-Part 1, Appl. Math. Comput., in press, doi: 10.1016/j.amc.2005.05.015] for obtaining a compact fourth-order method to the steady 2D Navier-Stokes equations in the streamfunction formulation-vorticity using the computer algebra system Maple, which includes conformal mappings and non-uniform grids. To analyse the procedure we have solved a constricted stepped channel problem, where a fine grid is placed near the re-entrant corner by transformation of the independent variables. (c) 2006 Elsevier B.V. All rights reserved.

NUMERICAL DETERMINATION OF THE ORDER OF PHASE-TRANSITIONS IN ISING SYSTEMS WITH MULTISPIN INTERACTIONS

The method of the fourth-order cumulant of Challa, Landau, and Binder is used together with the Monte Carlo histogram technique of Ferrenberg and Swendsen to study the order of the phase transitions of two-dimensional Ising systems with multispin interactions in the horizontal direction and two-body interactions in the vertical direction.

Chawla-Numerov method revisited

The well-known two-step fourth-order Numerov method was shown to have better interval of periodicity when made explicit, see Chawla (1984). It is readily verifiable that the improved method still has phase-lag of order 4. We suggest a slight modification from which linear problems could benefit. Phase-lag of any order can be achieved, but only order 6 is derived. © 1991.

Persistência na lactação para vacas da raça Holandesa criadas no Estado do Rio Grande do Sul via modelos de regressão aleatória

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

Random regression models using different functions to model test-day milk yield of Brazilian Holstein cows

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

Study of Stability of Rotational Motion of Spacecraft with Canonical Variables

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

Bose-Einstein condensation dynamics in three dimensions by the pseudospectral and finite-difference methods

We suggest a pseudospectral method for solving the three-dimensional time-dependent Gross-Pitaevskii (GP) equation, and use it to study the resonance dynamics of a trapped Bose-Einstein condensate induced by a periodic variation in the atomic scattering length. When the frequency of oscillation of the scattering length is an even multiple of one of the trapping frequencies along the x, y or z direction, the corresponding size of the condensate executes resonant oscillation. Using the concept of the differentiation matrix, the partial-differential GP equation is reduced to a set of coupled ordinary differential equations, which is solved by a fourth-order adaptive step-size control Runge-Kutta method. The pseudospectral method is contrasted with the finite-difference method for the same problem, where the time evolution is performed by the Crank-Nicholson algorithm. The latter method is illustrated to be more suitable for a three-dimensional standing-wave optical-lattice trapping potential.

Tm and Tm-Tb-doped germanate glasses for S-band amplifiers

The mechanism involved in the Tm(3+)((3)F(4)) -> Tb(3+)((7)F(0,1,2)) energy transfer as a function of the Tb concentration was investigated in Tm:Tb-doped germanate (GLKZ) glass. The experimental transfer rate was determined from the best fit of the (3)F(4) luminescence decay due to the Tm -> Tb energy transfer using the Burshtein model. The result showed that the 1700 nm emission from (3)F(4) can be completely quenched by 0.8 mol% of Tb(3+). As a consequence, the (7)F(3) state of Tb(3+) interacts with the (3)H(4) upper excited state of TM(3+) slighting decreasing its population. The effective amplification coefficient beta(cm(-1)) that depends on the population density difference Delta n = n((3)H(4))-n((3)F(4)) involved in the optical transition of Tm(3+) (S-band) was calculated by solving the rate equations of the system for continuous pumping with laser at 792 nm, using the Runge-Kutta numerical method including terms of fourth order. The population density inversion An as a function of Tb(3+) concentration was calculated by computational simulation for three pumping intensities, 0.2, 2.2 and 4.4 kWcm(-2). These calculations were performed using the experimental Tm -> Tb transfer rates and the optical constants of the Tm (0.1 mol%) system. It was demonstrated that 0.2 mol% of Tb(3+) propitiates best population density inversion of Tin(3+) maximizing the amplification coefficient of Tm-doped (0.1 mol%) GLKZ glass when operating as laser intensity amplification at 1.47 mu m. (C) 2007 Elsevier B.V. All rights reserved.

Solução numérica das equações de Navier-Stokes em um canal-tipo estenose usando métodos compactos e não compacos de alta ordem

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

Estudo da estabilidade do movimento rotacional de satélites artificiais com variáveis canônicas

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

Modelagem matemática de tumores e estequiometria biológica

Cancer biology is a complex and expanding field of science study. Due its complexity, there is a strong motivation to integrate many fields of knowledge to study cancer biology, and biological stoichiometry can make this. Biological stoichiometry is the study of the balance of multiple chemical elements in biological systems. A key idea in biological stoichiometry is the growth rate hypothesis, which states that variation in the carbon:nitrogen:phosphorus stoichiometry of living things is associated with growth rate because of the elevated demands for phosphorusrich ribosomal RNA and other elements necessary to protein synthesis. As tumor cells has high rate proliferation, the growth rate hypothesis can be used in cancer study. In this work the dynamic of two tumors (primary and secondary) and the chemical elements carbon and nitrogen are simulate and analyzed through mathematical models that utilize as central idea biological stoichiometry. Differential equations from mathematical model are solved by numerical method Runge-Kutta fourth order

Morfometria de microbacias do Córrego Rico, afluente do Rio Mogi-Guaçu, Estado de São Paulo, Brasil

Este trabalho teve como objetivo avaliar as características morfométricas das microbacias (2ª, 3ª, 4ª e 5ª ordens de magnitude) da bacia hidrográfica do córrego Rico, sub-bacia do Rio Mogi-Guaçu, localizada na região administrativa de Ribeirão Preto, Estado de São Paulo, Brasil. Para tanto, foram determinados os parâmetros físicos e a configuração topográfica natural do sistema de drenagem. Os procedimentos para a obtenção dos dados foram fundamentados em técnicas de sensoriamento remoto e geoprocessamento. A partir da vetorização das cartas topográficas correspondentes à área de estudo, realizou-se a análise morfométrica quanto às características dimensionais, do padrão de drenagem e do relevo no sistema de informação geográfica ArcView. A microbacia é considerada de sexta ordem de magnitude, com área estimada de 542 km², com 85 microbacias de segunda ordem, 22 de terceira, sete de quarta ordem e duas de quinta. Utilizando o critério geométrico, na disposição fluvial das sub-bacias de cabeceiras observou-se a predominância dos modelos dendríticos e subdendríticos, enquanto a jusante predominava o modelo subparalelo, respectivamente, nas áreas de ocorrências dos arenitos Bauru e rochas efusivas básicas.