A combination of implicit and adaptative upwind tools for the numerical solution of incompressible free surface flows

This paper is concerned with the numerical solutions of time dependent two-dimensional incompressible flows. By using the primitive variables of velocity and pressure, the Navier-Stokes and mass conservation equations are solved by a semi-implicit finite difference projection method. A new bounded higher order upwind convection scheme is employed to deal with the non-linear (advective) terms. The procedure is an adaptation of the GENSMAC (J. Comput. Phys. 1994; 110: 171-186) methodology for calculating confined and free surface fluid flows at both low and high Reynolds numbers. The calculations were performed by using the 2D version of the Freeflow simulation system (J. Comp. Visual. Science 2000; 2:199-210). In order to demonstrate the capabilities of the numerical method, various test cases are presented. These are the fully developed flow in a channel, the flow over a backward facing step, the die-swell problem, the broken dam flow, and an impinging jet onto a flat plate. The numerical results compare favourably with the experimental data and the analytical solutions. Copyright (c) 2006 John Wiley & Sons, Ltd.

Asymptotic soliton train solutions of Kaup-Boussinesq equations

Asymptotic soliton trains arising from a 'large and smooth' enough initial pulse are investigated by the use of the quasiclassical quantization method for the case of Kaup-Boussinesq shallow water equations. The parameter varying along the soliton train is determined by the Bohr-Sommerfeld quantization rule which generalizes the usual rule to the case of 'two potentials' h(0)(x) and u(0)(x) representing initial distributions of height and velocity, respectively. The influence of the initial velocity u(0)(x) on the asymptotic stage of the evolution is determined. Excellent agreement of numerical solutions of the Kaup-Boussinesq equations with predictions of the asymptotic theory is found. (C) 2003 Elsevier B.V. All rights reserved.

Numerical and variational solutions of the dipolar Gross-Pitaevskii equation in reduced dimensions

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

Asymptotic soliton train solutions of the defocusing nonlinear Schrodinger equation

Asymptotic behavior of initially large and smooth pulses is investigated at two typical stages of their evolution governed by the defocusing nonlinear Schrodinger equation. At first, wave breaking phenomenon is studied in the limit of small dispersion. A solution of the Whitham modulational equations is found for the case of dissipationless shock wave arising after the wave breaking point. Then, asymptotic soliton trains arising eventually from a large and smooth initial pulse are studied by means of a semiclassical method. The parameter varying along the soliton train is calculated from the generalized Bohr-Sommerfeld quantization rule, so that the distribution of eigenvalues depends on two functions-intensity rho(0)(x) of the initial pulse and its initial chirp v(0)(x). The influence of the initial chirp on the asymptotic state is investigated. Excellent agreement of the numerical solution of the defocusing NLS equation with predictions of the asymptotic theory is found.

A note about the appearance of non-hyperbolic solutions in a mechanical pendulum system

The investigation of the behavior of a nonlinear system consists in the analysis of different stages of its motion, where the complexity varies with the proximity of a resonance region. Near this region the stability domain of the system undergoes sudden changes due basically to competition and interaction between periodic and saddle solutions inside the phase portrait, leading to the occurrence of the most different phenomena. Depending of the domain of the chosen control parameter, these events can reveal interesting geometric features of the system so that the phase portrait is not capable to express all them, since the projection of these solutions on the two-dimensional surface can hide some aspects of these events. In this work we will investigate the numerical solutions of a particular pendulum system close to a secondary resonance region, where we vary the control parameter in a restrict domain in order to draw a preliminary identification about what happens with this system. This domain includes the appearance of non-hyperbolic solutions where the basin of attraction in the center of the phase portrait diminishes considerably, almost disappearing, and afterwards its size increases with the direction of motion inverted. This phenomenon delimits a boundary between low and high frequency of the external excitation.

An analytical and numerical study of plane change maneuvers using aerodynamic force

The goal of the present work is to analyze space missions that use the terrestrial atmosphere to accomplish orbital maneuvers that involve a plane change. A set of analytical solutions is presented for the variation of the orbital elements due to a single passage through the atmosphere, assuming that the interval the spacecraft travels through the atmosphere is not too large. The study considers both the lift influence on the spacecraft orbit as well as drag. The final equations are tested with numerical integration and can be considered in accordance with the numerical results whenever the perigee height is larger than a critical value. Next, a numerical study of the ratio between the velocity increment required to correct the semimajor axis decay due to the atmospheric passage and the velocity variation required to obtain the change in the inclination is also presented. This analysis can be used to decide if a maneuver passing through the atmosphere can decrease the fuel consumption of the mission and, in the cases where this technique can be used, if a multiple passage is more efficient than a single passage.

A study about the solutions of Schrödinger equation with an asymptotically linear potential

We determine the solutions of the Schrödinger equation for an asymptotically linear potential. Analytical solutions are obtained by superalgebra in quantum mechanics and we establish when these solutions are possible. Numerical solutions for the spectra are obtained by the shifted 1/N expansion method.

On certain new exact solutions of a diffusive predator-prey system

We construct exact solutions for a system of two coupled nonlinear partial differential equations describing the spatio-temporal dynamics of a predator-prey system where the prey per capita growth rate is subject to the Allee effect. Using the G'/G expansion method, we derive exact solutions to this model for two different wave speeds. For each wave velocity we report three different forms of solutions. We also discuss the biological relevance of the solutions obtained. © 2012 Elsevier B.V.

A Numerical Study of the FENE-CR Model Applied to a Jet Flow Problem

The FENE-CR model is investigated through a numerical algorithm to simulate the time-dependent moving free surface flow produced by a jet impinging on a flat surface. The objective is to demonstrate that by increasing the extensibility parameter L, the numerical solutions converge to the solutions obtained with the Oldroyd-B model. The governing equations are solved by an established free surface flow solver based on the finite difference and marker-and-cell methods. Numerical predictions of the extensional viscosity obtained with several values of the parameter L are presented. The results show that if the extensibility parameter L is sufficiently large then the extensional viscosities obtained with the FENE-CR model approximate the corresponding Oldroyd-B viscosity. Moreover, the flow from a jet impinging on a flat surface is simulated with various values of the extensibility parameter L and the fluid flow visualizations display convergence to the Oldroyd-B jet flow results.

A numerical study of powered Swing-Bys around the Moon

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