Three-dimensional transient complex free surface flows: Numerical simulation of XPP fluid

In this paper we present a finite difference MAC-type approach for solving three-dimensional viscoelastic incompressible free surface flows governed by the eXtended Pom-Pom (XPP) model, considering a wide range of parameters. The numerical formulation presented in this work is an extension to three-dimensions of our implicit technique [Journal of Non-Newtonian Fluid Mechanics 166 (2011) 165-179] for solving two-dimensional viscoelastic free surface flows. To enhance the stability of the numerical method, we employ a combination of the projection method with an implicit technique for treating the pressure on the free surfaces. The differential constitutive equation of the fluid is solved using a second-order Runge-Kutta scheme. The numerical technique is validated by performing a mesh refinement study on a pipe flow, and the numerical results presented include the simulation of two complex viscoelastic free surface flows: extrudate-swell problem and jet buckling phenomenon. © 2013 Elsevier B.V.

Optimisation of High Hydrostatic Pressure Processing of Pêra Rio Orange Juice

The influence of high hydrostatic pressure (HHP) on Pêra Rio orange juice was investigated using response surface methodology. A central composite design was used to evaluate the effects of three processing conditions (independent variables), namely pressure (100-600 MPa), temperature (30-60 °C) and time (30-360 s), on the native microflora and pectin methylesterase (PME) activity of orange juice. Analysis of variance showed that second-order polynomial models fitted well with the experimental data for PME residual activity (R2 = 0.9586, p < 0.001) and aerobic microorganism count (R2 = 0.9879, p < 0.001). The optimum HHP processing conditions to produce orange juice with PME residual activity of less than 20 % and low microorganism count (<2 log cycles CFU/mL) were 550 to 600 MPa, 55 to 60 °C and 330 to 360 s. © 2013 Springer Science+Business Media New York.

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

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

Ethanolysis optimisation of Jupati (Raphia taedigera Mart.) oil to biodiesel using response surface methodology

In this work, the transesterification of jupati (Raphia taedigera Mart.) oil using ethanol and acid catalyst was examined. The production of biodiesel was performed using a central composite design (CCD). A range of values for catalyst concentration (1 to 4.21%), temperature (70-80 °C), and the molar ratio of alcohol to oil (6:1-13.83:1) were tested, and ester content, viscosity, and yield were the response variables. The synthesis process was optimised using response surface methodology (RSM), resulting in the following optimal conditions for the production of jupati ethyl esters: a catalyst concentration of 3.85% at 80 °C and an alcohol-to-oil molar ratio of 10:1.

Numerical modeling of flow through an industrial burner orifice

This paper presents numerical modeling of a turbulent natural gas flow through a non-premixed industrial burner of a slab reheating furnace. The furnace is equipped with diffusion side swirl burners capable of utilizing natural gas or coke oven gas alternatively through the same nozzles. The study is focused on one of the burners of the preheating zone. Computational Fluid Dynamics simulation has been used to predict the burner orifice turbulent flow. Flow rate and pressure at burner upstream were validated by experimental measurements. The outcomes of the numerical modeling are analyzed for the different turbulence models in terms of pressure drop, velocity profiles, and orifice discharge coefficient. The standard, RNG, and Realizable k-epsilon models and Reynolds Stress Model (RSM) have been used. The main purpose of the numerical investigation is to determine the turbulence model that more consistently reproduces the experimental results of the flow through an industrial non-premixed burner orifice. The comparisons between simulations indicate that all the models tested satisfactorily and represent the experimental conditions. However, the Realizable k-epsilon model seems to be the most appropriate turbulence model, since it provides results that are quite similar to the RSM and RNG k-epsilon models, requiring only slightly more computational power than the standard k-epsilon model. (C) 2014 Elsevier Ltd. All rights reserved.

A real-time thermal field theoretical analysis of Kubo-type shear viscosity: Numerical understanding with simple examples

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

The NINJA-2 project: detecting and characterizing gravitational waveforms modelled using numerical binary black hole simulations

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave (GW) astrophysics communities. The purpose of NINJA is to study the ability to detect GWs emitted from merging binary black holes (BBH) and recover their parameters with next-generation GW observatories. We report here on the results of the second NINJA project, NINJA-2, which employs 60 complete BBH hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. In a 'blind injection challenge' similar to that conducted in recent Laser Interferometer Gravitational Wave Observatory (LIGO) and Virgo science runs, we added seven hybrid waveforms to two months of data recoloured to predictions of Advanced LIGO (aLIGO) and Advanced Virgo (AdV) sensitivity curves during their first observing runs. The resulting data was analysed by GW detection algorithms and 6 of the waveforms were recovered with false alarm rates smaller than 1 in a thousand years. Parameter-estimation algorithms were run on each of these waveforms to explore the ability to constrain the masses, component angular momenta and sky position of these waveforms. We find that the strong degeneracy between the mass ratio and the BHs' angular momenta will make it difficult to precisely estimate these parameters with aLIGO and AdV. We also perform a large-scale Monte Carlo study to assess the ability to recover each of the 60 hybrid waveforms with early aLIGO and AdV sensitivity curves. Our results predict that early aLIGO and AdV will have a volume-weighted average sensitive distance of 300 Mpc (1 Gpc) for 10M circle dot + 10M circle dot (50M circle dot + 50M circle dot) BBH coalescences. We demonstrate that neglecting the component angular momenta in the waveform models used in matched-filtering will result in a reduction in sensitivity for systems with large component angular momenta. This reduction is estimated to be up to similar to 15% for 50M circle dot + 50M circle dot BBH coalescences with almost maximal angular momenta aligned with the orbit when using early aLIGO and AdV sensitivity curves.

Relaxation to Fixed Points in the Logistic and Cubic Maps: Analytical and Numerical Investigation

Convergence to a period one fixed point is investigated for both logistic and cubic maps. For the logistic map the relaxation to the fixed point is considered near a transcritical bifurcation while for the cubic map it is near a pitchfork bifurcation. We confirmed that the convergence to the fixed point in both logistic and cubic maps for a region close to the fixed point goes exponentially fast to the fixed point and with a relaxation time described by a power law of exponent -1. At the bifurcation point, the exponent is not universal and depends on the type of the bifurcation as well as on the nonlinearity of the map.

Numerical and experimental investigations of three-dimensional container filling with Newtonian viscous fluids

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

Numerical solution of the FENE-CR model in complex flows

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

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.

Numerical solution of the time dependent Ginzburg-Landau equations for mixed (d plus s)-wave superconductors

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

Discovering promising regions to help global numerical optimization algorithms

We have developed an algorithm using a Design of Experiments technique for reduction of search-space in global optimization problems. Our approach is called Domain Optimization Algorithm. This approach can efficiently eliminate search-space regions with low probability of containing a global optimum. The Domain Optimization Algorithm approach is based on eliminating non-promising search-space regions, which are identifyed using simple models (linear) fitted to the data. Then, we run a global optimization algorithm starting its population inside the promising region. The proposed approach with this heuristic criterion of population initialization has shown relevant results for tests using hard benchmark functions.

Optimisation of a sample preparation method for the determination of fumonisin B-1 in rice

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

The dynamic behavior of a cantilever beam coupled to a non-ideal unbalanced motor through numerical and experimental analysis

An excitation force that is not influenced by the system state is said to be an ideal energy source. In real situations, a direct and feedback coupling between the excitation source and the system must always exist at a certain level. This manifestation of the law of conservation of energy is known as the Sommerfeld effect. In the case of obtaining a mathematical model for such a system, additional equations are usually necessary to describe the vibration sources with limited power and its coupling with the mechanical system. In this work, a cantilever beam and a non-ideal DC motor fixed to its free end are analyzed. The motor has an unbalanced mass that provides excitation to the system which is proportional to the current applied to the motor. During the coast up operation of the motor, if the drive power is increased slowly, making the excitation frequency pass through the first natural frequency of the beam, the DC motor speed will remain the same until it suddenly jumps to a much higher value (simultaneously its amplitude jumps to a much lower value) upon exceeding a critical input power. It was found that the Sommerfeld effect depends on some system parameters and the motor operational procedures. These parameters are explored to avoid the resonance capture in the Sommerfeld effect. Numerical simulations and experimental tests are used to help gather insight of this dynamic behavior. (C) 2014 Elsevier Ltd. All rights reserved.