A General Solution of the Fokker-Planck Equation

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

On an approximate analytical solution to a nonlinear vibrating problem, excited by a nonideal motor

This paper analyzes through Multiple Scales Method a response of a simplified nonideal and nonlinear vibrating system. Here, one verifies the interactions between the dynamics of the DC motor (excitation) and the dynamics of the foundation (spring, damper, and mass). We remarked that we consider cubic nonlinearity (spring) and quadratic nonlinearity (DC motor) of the same order of magnitude according to experimental results. Both analytical and numerical results that we have obtained had good agreement.

Localization of the optimal solution and a posteriori bounds for aggregation

After an aggregated problem has been solved, it is often desirable to estimate the accuracy loss due to the fact that a simpler problem than the original one has been solved. One way of measuring this loss in accuracy is the difference in objective function values. To get the bounds for this difference, Zipkin (Operations Research 1980;28:406) has assumed, that a simple (knapsack-type) localization of an original optimal solution is known. Since then various extensions of Zipkin's bound have been proposed, but under the same assumption. A method to compute the bounds for variable aggregation for convex problems, based on general localization of the original solution is proposed. For some classes of the original problem it is shown how to construct the localization. Examples are given to illustrate the main constructions and a small numerical study is presented.

Spatial distribution of solutes and water in sucrose solution dehydrated apples

Half-fresh apples were immersed in sucrose solution (50% w/w, 27 degrees C) during different times of exposition (2, 4, and 8 h). Then each fruit was sliced from the transversal exposed surface. Density, water, and sugar content were determined for each slice. A mathematical model was fitted to experimental data of water and sucrose content considering the global flux and the tissue shrinkage. By numerical analysis, the binary effective diffusion coefficients as a function of concentration were calculated, using material coordinates and integrating simultaneously two differential equations (for water and sucrose). The coefficients obtained are one or even two orders of magnitude lower than the ones for pure solutions and present an unusual concentration dependence. This comparison shows the influence of the tissue resistance to the diffusion.

Numerical study of the spherically symmetric Gross-Pitaevskii equation in two space dimensions

The Gross-Pitaevskii equation for Bose-Einstein condensation (BEC) in two space dimensions under the action of a harmonic oscillator trap potential for bosonic atoms with attractive and repulsive interparticle interactions was numerically studied by using time-dependent and time-independent approaches. In both cases, numerical difficulty appeared for large nonlinearity. Nonetheless, the solution of the time-dependent approach exhibited intrinsic oscillation with time iteration which is independent of space and time steps used in discretization.

A new solution to the optimal power flow problem

A new approach to solving the Optimal Power Flow problem is described, making use of some recent findings, especially in the area of primal-dual methods for complex programming. In this approach, equality constraints are handled by Newton's method inequality constraints for voltage and transformer taps by the logarithmic barrier method and the other inequality constraints by the augmented Lagrangian method. Numerical test results are presented, showing the effective performance of this algorithm. © 2001 IEEE.

On the solution of mathematical programming problems with equilibrium constraints

Mathematical programming problems with equilibrium constraints (MPEC) are nonlinear programming problems where the constraints have a form that is analogous to first-order optimality conditions of constrained optimization. We prove that, under reasonable sufficient conditions, stationary points of the sum of squares of the constraints are feasible points of the MPEC. In usual formulations of MPEC all the feasible points are nonregular in the sense that they do not satisfy the Mangasarian-Fromovitz constraint qualification of nonlinear programming. Therefore, all the feasible points satisfy the classical Fritz-John necessary optimality conditions. In principle, this can cause serious difficulties for nonlinear programming algorithms applied to MPEC. However, we show that most feasible points do not satisfy a recently introduced stronger optimality condition for nonlinear programming. This is the reason why, in general, nonlinear programming algorithms are successful when applied to MPEC.

Electroosmotic pumping in rectangular microchannels: A numerical treatment by the finite element method

In this work, the analysis of electroosmotic pumping mechanisms in microchannels is performed through the solution of Poisson-Boltzmann and Navier Stokes equations by the Finite Element Method. This approach is combined with a Newton-Raphson iterative scheme, allowing a full treatment of the non-linear Poisson-Boltzmann source term which is normally approximated by linearizations in other methods.

An exact series solution for the vibration analysis of cylindrical shells with arbitrary boundary conditions

In this paper, an exact series solution for the vibration analysis of circular cylindrical shells with arbitrary boundary conditions is obtained, using the elastic equations based on Flügge's theory. Each of the three displacements is represented by a Fourier series and auxiliary functions and sought in a strong form by letting the solution exactly satisfy both the governing differential equations and the boundary conditions on a point-wise basis. Since the series solution has to be truncated for numerical implementation, the term exactly satisfying should be understood as a satisfaction with arbitrary precision. One of the important advantages of this approach is that it can be universally applied to shells with a variety of different boundary conditions, without the need of making any corresponding modifications to the solution algorithms and implementation procedures as typically required in other techniques. Furthermore, the current method can be easily used to deal with more complicated boundary conditions such as point supports, partial supports, and non-uniform elastic restraints. Numerical examples are presented regarding the modal parameters of shells with various boundary conditions. The capacity and reliability of this solution method are demonstrated through these examples. © 2012 Elsevier Ltd. All rights reserved.

Numerical study of the Ginzburg-Landau-Langevin equation: coherent structures and noise perturbation theory

Pós-graduação em Física - IFT

Generalized integral transform solution for hydrodynamically developing non-Newtonian flows in circular tubes

ABSTRACT: The Generalized Integral Transform Technique (GITT) is applied to the solution of the momentum equations in a hydrodynamically developing laminar flow of a non-Newtonian power-law fluid inside a circular duct. A primitive variables formulation is adopted in order to avoid the singularity of the auxiliary eigenvalue problem in terms of Bessel functions at the centerline of the duct when the GITT approach is applied. Results for the velocity field and friction factor-Reynolds number product are computed for different power-law indices, which are tabulated and graphically presented as functions of the dimensionless coordinates. Critical comparisons with previous results in the literature are also performed, in order to validate the numerical codes developed in the present work and to demonstrate the consistency of the final results.

Integral transform solution for the forced convection of herschel-bulkley fluids in circular tubes and parallel-plates ducts

ABSTRACT: The thermal entry region in laminar forced convection of Herschel-Bulkley fluids is solved analytically through the integral transform technique, for both circular and parallel-plates ducts, which are maintained at a prescribed wall temperature or at a prescribed wall heat flux. The local Nusselt numbers are obtained with high accuracy in both developing and fully-developed thermal regions, and critical comparisons with previously reported numerical results are performed.

Exact Solution of the Elastica with Transverse Shear Effects

A new derivation of Euler's Elastica with transverse shear effects included is presented. The elastic potential energy of bending and transverse shear is set up. The work of the axial compression force is determined. The equation of equilibrium is derived using the variation of the total potential. Using substitution of variables an exact solution is derived. The equation is transcendental and does not have a closed form solution. It is evaluated in a dimensionless form by using a numerical procedure. Finally, numerical examples of laminates made of composite material (fiber reinforced) and sandwich panels are provided.

A numerical study of powered Swing-Bys around the Moon

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

URANS calculations for smooth circular cylinder flow in a wide range of Reynolds Numbers: solution verification and validation

The flow around circular smooth fixed cylinder in a large range of Reynolds numbers is considered in this paper. In order to investigate this canonical case, we perform CFD calculations and apply verification & validation (V&V) procedures to draw conclusions regarding numerical error and, afterwards, assess the modeling errors and capabilities of this (U)RANS method to solve the problem. Eight Reynolds numbers between Re = 10 and Re 5 x 10(5) will be presented with, at least, four geometrically similar grids and five discretization in time for each case (when unsteady), together with strict control of iterative and round-off errors, allowing a consistent verification analysis with uncertainty estimation. Two-dimensional RANS, steady or unsteady, laminar or turbulent calculations are performed. The original 1994 k - omega SST turbulence model by Menter is used to model turbulence. The validation procedure is performed by comparing the numerical results with an extensive set of experimental results compiled from the literature. [DOI: 10.1115/1.4007571]