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,

Time reduction in a forced-convection thermal process using simplified transient analysis

This study investigates the utilisation of a simplified model in the transient analysis of a thermal cooling process. In such process the external thermal resistance between the surface and the surroundings is high compared to the system internal thermal resistance, so that the first controls the heat transfer process. In this case the Biot number is lower than 0.1. Aluminium reels were utilised, which, with proper internal instrumentation, furnished experimental results for the thermal cooling process. Based on experimental data, a simplified model for the determination of the process film coefficient was used. Subsequently, experimental and theoretical results were compared. The change of the airflow direction was also investigated for the cooling process, aiming at process time optimisation. (C) 2001 Elsevier B.V. Ltd.

NUMERICAL SIMULATION of CONVECTION-DIFFUSION PROBLEMS BY THE CONTROL-VOLUME-BASED FINITE-ELEMENT METHOD

This work presents a numerical study of the tri-dimensional convection-diffusion equation by the control-volume-based on finite-element method using quadratic hexahedral elements. Considering that the equation governing this problem in its main variable may represent several properties, including temperature, turbulent kinetic energy, viscous dissipation rate of the turbulent kinetic energy, specific dissipation rate of the turbulent kinetic energy, or even the concentration of a contaminant in a given medium, among others, the wide applicability of this problem is thus evidenced. Three cases of temperature distributions will be studied specifically in this work, in addition to one case of pollutant dispersion upon analysis of the concentration of a contaminant in a fixed flow point. Some comparisons will be carried out against works found in the open literature, while others will be done according to each phenomenon characteristics.

Comments on nonlinear dynamics of a non-ideal Duffing-Rayleigh oscillator: Numerical and analytical approaches

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

Forced convection to laminar flow of liquid egg yolk in circular and annular ducts

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

BOUSSINESQ-TYPE SYSTEM OF EQUATIONS IN THE BENARD-MARANGONI SYSTEM

A system of coupled evolution equations for the bulk velocity and the surface displacement is found to govern the long-wavelength perturbations in a Benard-Marangoni system. This system of equations, involving nonlinearity, dispersion, and dissipation, is a generalization of the usual Boussinesq system.

Numerical and experimental analysis of natural convection in a cavity heated from below

In this work, an analysis of the natural convection flow caused by heat sources dissipating energy at a constant rate simulating electronic components mounted at the bottom surface of a cavity symmetrically cooled from the sides and insulated at the top is performed. This problem was studied numerically and experimentally for several aspect ratios (height/width), for different levels of dissipation in the sources, and for different side wall temperatures. Temperature and velocity fields were determined as well as the temperature variation along the surface where the sources are mounted and the average Nusselt number in the source surfaces. Numerical and experimental results were found to agree.

DISSIPATIVE BOUSSINESQ SYSTEM OF EQUATIONS IN THE BENARD-MARANGONI PHENOMENON

By using the long-wavelength approximation, a system of coupled evolution equations for the bulk velocity and the surface perturbations of a Benard-Marangoni system is obtained. It includes nonlinearity, dispersion, and dissipation, and it can be interpreted as a dissipative generalization of the usual Boussinesq system of equations. As a particular case, a strictly dissipative version of the Boussinesq system is obtained.

Laminar forced convection to a pseudoplastic fluid food in circular and annular ducts

In this manuscript we investigated experimentally the steady-state heat transfer to an important pseudoplastic fluid food, the soursop juice, flowing in laminar regime through circular and concentric annular ducts. The mean convection heat transfer coefficients, determined by measuring the bulk temperatures before and after the heating sections with constant temperatures of the tube walls, were used to correlate simple new empiric expressions to estimate the average Nusselt number in the thermal entrance of the considered geometries. In addition, the thermophysical properties of the tested fluid food, as well as the rheological behavior, being essential for the heat transfer analyses, were experimentally determined. (c) 2006 Elsevier Ltd. All rights reserved.

Nonlinear parameter estimation in laminar forced convection within a circular sector tube

In this article we examine an inverse heat convection problem of estimating unknown parameters of a parameterized variable boundary heat flux. The physical problem is a hydrodynamically developed, thermally developing, three-dimensional steady state laminar flow of a Newtonian fluid inside a circular sector duct, insulated in the flat walls and subject to unknown wall heat flux at the curved wall. Results are presented for polynomial and sinusoidal trial functions, and the unknown parameters as well as surface heat fluxes are determined. Depending on the nature of the flow, on the position of experimental points the inverse problem sometimes could not be solved. Therefore, an identification condition is defined to specify a condition under which the inverse problem can be solved. Once the parameters have been computed it is possible to obtain the statistical significance of the inverse problem solution. Therefore, approximate confidence bounds based on standard statistical linear procedure, for the estimated parameters, are analyzed and presented.

NONLINEAR SURFACE-WAVE EXCITATIONS IN THE BENARD-MARANGONI SYSTEM

We examine the appearance of surface waves governed by Burgers and Korteweg-de Vries equations in a shallow viscous heated fluid. We consider waves triggered by a surface-tension variation induced by both temperature and concentration gradients. We also establish the range of parameters for which the above-mentioned equations appear.

Transient non-Darcy forced convective heat transfer from a flat plate embedded in a fluid-saturated porous medium

Transient non-Darcy forced convection on a flat plate embedded in a porous medium is investigated using the Forchheimer-extended Darcy law. A sudden uniform pressure gradient is applied along the flat plate, and at the same time, its wall temperature is suddenly raised to a high temperature. Both the momentum and energy equations are solved by retaining the unsteady terms. An exact velocity solution is obtained and substituted into the energy equation, which then is solved by means of a quasi-similarity transformation. The temperature field can be divided into the one-dimensional transient (downstream) region and the quasi-steady-state (upstream) region. Thus the transient local heat transfer coefficient can be described by connecting the quasi-steady-state solution and the one-dimensional transient solution. The non-Darcy porous inertia works to decrease the velocity level and the time required for reaching the steady-state velocity level. The porous-medium inertia delays covering of the plate by the steady-state thermal boundary layer. © 1990.

Numerical and experimental analysis of natural convection in a cavity with flush mounted heat sources on a side wall

This work considers a problem of interest in several technological applications such as the thermal control of electronic equipment. It is also important to study the heat transfer performance of these components under off-normal conditions, such as during failure of cooling fans. The effect of natural convection on the flow and heat transfer in a cavity with two flush mounted heat sources on the left vertical wall, simulating electronic components, is studied numerically and experimentally. The influence of the power distribution, spacing between the heat sources and cavity aspect ratio have been investigated. An analysis of the average Nusselt number of the two heat sources was performed to investigate the behavior of the heat transfer coefficients. The results obtained numerically and experimentally, after an error analysis, showed a good agreement.

Thermally developing forced convection of non-Newtonian fluids inside elliptical ducts

Laminar-forced convection inside tubes of various cross-section shapes is of interest in the design of a low Reynolds number heat exchanger apparatus. Heat transfer to thermally developing, hydrodynamically developed forced convection inside tubes of simple geometries such as a circular tube, parallel plate, or annular duct has been well studied in the literature and documented in various books, but for elliptical duct there are not much work done. The main assumptions used in this work are a non-Newtonian fluid, laminar flow, constant physical properties, and negligible axial heat diffusion (high Peclet number). Most of the previous research in elliptical ducts deal mainly with aspects of fully developed laminar flow forced convection, such as velocity profile, maximum velocity, pressure drop, and heat transfer quantities. In this work, we examine heat transfer in a hydrodynamically developed, thermally developing laminar forced convection flow of fluid inside an elliptical tube under a second kind of a boundary condition. To solve the thermally developing problem, we use the generalized integral transform technique (GITT), also known as Sturm-Liouville transform. Actually, such an integral transform is a generalization of the finite Fourier transform, where the sine and cosine functions are replaced by more general sets of orthogonal functions. The axes are algebraically transformed from the Cartesian coordinate system to the elliptical coordinate system in order to avoid the irregular shape of the elliptical duct wall. The GITT is then applied to transform and solve the problem and to obtain the once unknown temperature field. Afterward, it is possible to compute and present the quantities of practical interest, such as the bulk fluid temperature, the local Nusselt number, and the average Nusselt number for various cross-section aspect ratios.

Modelling study of the impact of deep convection on the utls air composition - Part I: Analysis of ozone precursors

The aim of this work is to study the local impact on the upper troposphere/lower stratosphere air composition of an extreme deep convective system. For this purpose, we performed a simulation of a convective cluster composed of many individual deep convective cells that occurred near Bauru (Brazil). The simulation is performed using the 3-D mesoscale model RAMS coupled on-line with a chemistry model. The comparisons with meteorological measurements show that the model produces meteorological fields generally consistent with the observations. The present paper (part I) is devoted to the analysis of the ozone precursors (CO, NO x and non-methane volatile organic compounds) and HO x in the UTLS. The simulation results show that the distribution of CO with altitude is closely related to the upward convective motions and consecutive outflow at the top of the convective cells leading to a bulge of CO between 7 km altitude and the tropopause (around 17km altitude). The model results for CO are consistent with satellite-borne measurements at 700 hPa. The simulation also indicates enhanced amounts of NO x up to 2 ppbv in the 7-17 km altitude layer mainly produced by the lightning associated with the intense convective activity. For insoluble non-methane volatile organic compounds, the convective activity tends to significantly increase their amount in the 7-17km layer by dynamical effects. During daytime in the presence of lightning NO x, this bulge is largely reduced in the upper part of the layer for reactive species (e.g. isoprene, ethene) because of their reactions with OH that is increased on average during daytime. Lightning NO x also impacts on the oxydizing capacity of the upper troposphere by reducing on average HO x, HO 2, H 2O 2 and organic hydroperoxides. During the simulation time, the impact of convection on the air composition of the lower stratosphere is negligible for all ozone precursors although several of the simulated convective cells nearly reach the tropopause. There is no significant transport from the upper troposphere to the lower stratosphere, the isentropic barrier not being crossed by convection. The impact of the increase of ozone precursors and HO x in the upper troposphere on the ozone budget in the LS is discussed in part II of this series of papers.