Evaluation of a bounded high order upwind scheme for 3D incompressible free surface flow computations

fit the context of normalized variable formulation (NVF) of Leonard and total variation diminishing (TVD) constraints of Harten. this paper presents an extension of it previous work by the authors for solving unsteady incompressible flow problems. The main contributions of the paper are threefold. First, it presents the results of the development and implementation of a bounded high order upwind adaptative QUICKEST scheme in the 3D robust code (Freeflow), for the numerical solution of the full incompressible Navier-Stokes equations. Second, it reports numerical simulation results for 1D hock tube problem, 2D impinging jet and 2D/3D broken clam flows. Furthermore, these results are compared with existing analytical and experimental data. and third, it presents the application of the numerical method for solving 3D free surface flow problems. (C) 2007 IMACS. Published by Elsevier B.V. All rights reserved,

An application of the Element Free Galerkin Method to the analysis of quantum well structures

In order to obtain the quantum-mechanical properties of layered semicondutor structures (quantum well and superlattice structures, for instance), solutions of the Schrodinger equation should be obtained for arbitrary potential profiles. In this paper, it is shown that such problems may be also studied by the Element Free Galerkin Method.

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 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.

Optic Flow Contribution to Locomotion Adjustments in Obstacle Avoidance

Locomotion generates a visual movement pattern characterized as optic flow. To explore how the locomotor adjustments are affected by this pattern, an experimental paradigm was developed to eliminate optic flow during obstacle avoidance. The aim was to investigate the contribution of optic flow in obstacle avoidance by using a stroboscopic lamp. Ten young adults walked on an 8m pathway and stepped over obstacles at two heights. Visual sampling was determined by a stroboscopic lamp (static and dynamic visual sampling). Three-dimensional kinematics data showed that the visual information about self-motion provided by the optic flow was crucial for estimating the distance from and the height of the obstacle. Participants presented conservative behavior for obstacle avoidance under experimental visual sampling conditions, which suggests that optic flow favors the coupling of vision to adaptive behavior for obstacle avoidance.

Simulation of free surface flows in a distributed memory environment

A parallel technique, for a distributed memory machine, based on domain decomposition for solving the Navier-Stokes equations in cartesian and cylindrical coordinates in two dimensions with free surfaces is described. It is based on the code by Tome and McKee (J. Comp. Phys. 110 (1994) 171-186) and Tome (Ph.D. Thesis, University of Strathclyde, Glasgow, 1993) which in turn is based on the SMAC method by Amsden and Harlow (Report LA-4370, Los Alamos Scientific Laboratory, 1971), which solves the Navier-Stokes equations in three steps: the momentum and Poisson equations and particle movement, These equations are discretized by explicit and 5-point finite differences. The parallelization is performed by splitting the computation domain into vertical panels and assigning each of these panels to a processor. All the computation can then be performed using nearest neighbour communication. Test runs comparing the performance of the parallel with the serial code, and a discussion of the load balancing question are presented. PVM is used for communication between processes. (C) 1999 Elsevier B.V. B.V. All rights reserved.

Numerical assessment of mass conservation on a MAC-type method for viscoelastic free surface flows

A numerical study of mass conservation of MAC-type methods is presented, for viscoelastic free-surface flows. We use an implicit formulation which allows for greater time steps, and therefore time marching schemes for advecting the free surface marker particles have to be accurate in order to preserve the good mass conservation properties of this methodology. We then present an improvement by using a Runge-Kutta scheme coupled with a local linear extrapolation on the free surface. A thorough study of the viscoelastic impacting drop problem, for both Oldroyd-B and XPP fluid models, is presented, investigating the influence of timestep, grid spacing and other model parameters to the overall mass conservation of the method. Furthermore, an unsteady fountain flow is also simulated to illustrate the low mass conservation error obtained.

Application of a bounded upwinding scheme to complex fluid dynamics problems

This paper is concerned with an overview of upwinding schemes, and further nonlinear applications of a recently introduced high resolution upwind differencing scheme, namely the ADBQUICKEST [V.G. Ferreira, F.A. Kurokawa, R.A.B. Queiroz, M.K. Kaibara, C.M. Oishi, J.A.Cuminato, A.F. Castelo, M.F. Tomé, S. McKee, assessment of a high-order finite difference upwind scheme for the simulation of convection-diffusion problems, International Journal for Numerical Methods in Fluids 60 (2009) 1-26]. The ADBQUICKEST scheme is a new TVD version of the QUICKEST [B.P. Leonard, A stable and accurate convective modeling procedure based on quadratic upstream interpolation, Computer Methods in Applied Mechanics and Engineering 19 (1979) 59-98] for solving nonlinear balance laws. The scheme is based on the concept of NV and TVD formalisms and satisfies a convective boundedness criterion. The accuracy of the scheme is compared with other popularly used convective upwinding schemes (see, for example, Roe (1985) [19], Van Leer (1974) [18] and Arora & Roe (1997) [17]) for solving nonlinear conservation laws (for example, Buckley-Leverett, shallow water and Euler equations). The ADBQUICKEST scheme is then used to solve six types of fluid flow problems of increasing complexity: namely, 2D aerosol filtration by fibrous filters; axisymmetric flow in a tubular membrane; 2D two-phase flow in a fluidized bed; 2D compressible Orszag-Tang MHD vortex; axisymmetric jet onto a flat surface at low Reynolds number and full 3D incompressible flows involving moving free surfaces. The numerical simulations indicate that this convective upwinding scheme is a good generic alternative for solving complex fluid dynamics problems. © 2012.

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.

Anteroventral third ventricle lesions impair cardiovascular responses to intravenous hypertonic saline infusion

The anteroventral third ventricle (AV3V) region is a critical area of the forebrain, acting on fluid and electrolyte balance and maintaining cardiovascular homeostasis. The purpose of this study was to determine the effects of lesions to the anteroventral third ventricle region on cardiovascular responses to intravenous hypertonic saline (HS) infusion, Male Wistar rats were anesthetized with urethane. The femoral artery and jugular vein were cannulated to record mean arterial pressure (MAP) and infuse hypertonic saline (3M NaCl, 0.18 mL/100 g bw, over 1 min), respectively. Renal blood flow (RBF) was recorded by ultrasonic transit-time flow probes. Renal vascular conductance (RVC) was calculated as renal blood flow to mean arterial pressure ratio and expressed as percentage of baseline. After hypertonic saline infusion in sham animals, renal blood flow and renal vascular conductance increased to 137+10% and 125+7% (10 min), and 141 +/- 10% and 133 +/- 10% (60 min), respectively. Increases in mean arterial pressure (20-min peak: 12 +/- 3 mm Hg) were also observed. An acute lesion in the AV3V region (DC, 2 mA 25s) 30 min before infusion abrogated the effects of hypertonic saline. Mean arterial pressure was unchanged and renal blood flow and renal vascular conductance were 107 +/- 7% and 103 +/- 6% (10 min), and 107 +/- 4 and 106 +/- 4% (60 min), respectively. Marked tachycardia was observed immediately after lesion. Responses of chronic sham or lesioned rats were similar to those of acute animals. However, in chronic lesioned rats, hypertonic saline induced sustained hypertension. These results demonstrate that integrity of the AV3V region is essential for the renal vasodilation that follows acute changes in extracellular fluid compartment composition. (C) 2004 Elsevier B.V. All rights reserved.

An interactive environment for solid modeling based on free software

This paper describes an interactive environment built entirely upon public domain or free software, intended to be used as the preprocessor of a finite element package for the simulation of three-dimensional electromagnetic problems.

Influence from the free volume on the photoinduced birefringence in azocompound-containing polymers

The photoinduced birefringence is analyzed in a guest-host azobenzene-containing polymer in the temperature range from 20 to 330 K. An anomalous behavior arises in the low-temperature range, suggesting strong influence from the free volume for the chromophores in the polymer. This influence is so strong that quenched samples have a photoinduced signal ca. 5 times greater than the annealed ones at room temperature. An extended free volume model is presented based on two assumptions about thermal fluctuations in the cavities and their size distribution. This model, which is an extension of the model by Mita et al., can explain the main features of the photoinduced birefringence as a function of time, temperature, and initial free volume state. To account for the influence of free volume on the photoorientation, the detailed reorientation model by Sekkat's was used. We show that Sekkat's model leads to an exponential behavior at small orientation regimes, which simplifies the mathematical treatment and allows the mean free volume to be obtained from the data fitting.

Interaction between wave and coastal structure: Validation of two lagrangian numerical models with experimental results marine 2011

Numerical modeling of the interaction among waves and coastal structures is a challenge due to the many nonlinear phenomena involved, such as, wave propagation, wave transformation with water depth, interaction among incident and reflected waves, run-up / run-down and wave overtopping. Numerical models based on Lagrangian formulation, like SPH (Smoothed Particle Hydrodynamics), allow simulating complex free surface flows. The validation of these numerical models is essential, but comparing numerical results with experimental data is not an easy task. In the present paper, two SPH numerical models, SPHysics LNEC and SPH UNESP, are validated comparing the numerical results of waves interacting with a vertical breakwater, with data obtained in physical model tests made in one of the LNEC's flume. To achieve this validation, the experimental set-up is determined to be compatible with the Characteristics of the numerical models. Therefore, the flume dimensions are exactly the same for numerical and physical model and incident wave characteristics are identical, which allows determining the accuracy of the numerical models, particularly regarding two complex phenomena: wave-breaking and impact loads on the breakwater. It is shown that partial renormalization, i.e. renormalization applied only for particles near the structure, seems to be a promising compromise and an original method that allows simultaneously propagating waves, without diffusion, and modeling accurately the pressure field near the structure.

Ablation of the ability to control the right-to-left cardiac shunt does not affect oxygen uptake, specific dynamic action or growth in the rattlesnake Crotalus durissus

The morphologically undivided ventricle of the heart in non-crocodilian reptiles permits the mixing of oxygen-rich blood returning from the lungs and oxygen-poor blood from the systemic circulation. A possible functional significance for this intra-cardiac shunt has been debated for almost a century. Unilateral left vagotomy rendered the single effective pulmonary artery of the South American rattlesnake, Crotalus durissus, unable to adjust the magnitude of blood flow to the lung. The higher constant perfusion of the lung circulation and the incapability of adjusting the right-left shunt in left-denervated snakes persisted over time, providing a unique model for investigation of the long-term consequences of cardiac shunting in a squamate. Oxygen uptake recorded at rest and during spontaneous and forced activity was not affected by removing control of the cardiac shunt. Furthermore, metabolic rate and energetic balance during the post-prandial metabolic increment, plus the food conversion efficiency and growth rate, were all similarly unaffected. These results show that control of cardiac shunting is not associated with a clear functional advantage in adjusting metabolic rate, effectiveness of digestion or growth rates. © 2013. Published by The Company of Biologists Ltd.

Assessing the impact of menopause on salivary flow and xerostomia

Background: The aim of this study was to evaluate the symptoms of dry mouth and salivary flow in menarche and menopausal women. Methods: Objective and subjective assessment of salivary function were analysed by Xerostomia Inventory and Visual Analogue Scale questionnaire in menopausal and menarche women (control group). Salivary flow was evaluated by a chemical absorption stimulation test. Each subject provided three saliva samples: S1, non-stimulated saliva; S2, saliva initially stimulated with two drops of citric acid 2.5%; and S3, saliva super-stimulated with two drops of citric acid 2.5% every 30 seconds for two minutes. Results: No intergroup association was observed between Xerostomia Inventory and Visual Analogue Scale questionnaire. In both groups, the salivary flow was greatest at S3, followed by S2 and finally S1. Salivary flow was lower in the menopausal group compared to the control group only in S2 and S3. Conclusions: In the menopausal group, the salivary flow showed reduction but without clinical symptoms of dry mouth. It is important to normalize salivary flow to prevent oral abnormalities and maintain oral health. © 2013 Australian Dental Association.