Modeling and simulation of severe slugging in air-water pipeline-riser systems

A mathematical model, numerical simulations and stability and flow regime maps corresponding to severe slugging in pipeline riser systems, are presented. In the simulations air and water were used as flowing fluids. The mathematical model considers continuity equations for liquid and gas phases, with a simplified momentum equation for the mixture, neglecting inertia. A drift-flux model, evaluated for the local conditions in the riser, is used as a closure law. The developed model predicts the location of the liquid accumulation front in the pipeline and the liquid level in the riser, so it is possible to determine which type of severe slugging occurs in the system. The numerical procedure is convergent for different nodalizations. A comparison is made with experimental results corresponding to a catenary riser, showing very good results for slugging cycle and stability and flow regime maps. (c) 2010 Elsevier Ltd. All rights reserved.

3D CFD Simulation of Vortex-induced Vibration of Cylinder

The 3D flow around a circular cylinder free to oscillate transversely to the free stream was simulated using Computational Fluid Dynamics (CFD) and the Spalart-Allmaras Detached Eddy Simulation (DES) turbulence model for a Reynolds number Re = 10(4). Simulations were carried out for a small mass-damping parameter m*zeta = 0.00858, where m* = 3.3 and zeta = 0.0026. We found good agreement between the numerical results and experimental data. The simulations predicted the high observed amplitudes of the upper branch of vortex-induced vibrations for low mass-damping parameters.

Safety of a railway wheelset - derailment simulation with increasing lateral force

The motivation for this research is to make a comparison between dynamic results of a free railway wheelset derailment and safety limits. For this purpose, a numerical simulation of a wheelset derailment submitted to increasing lateral force is used to compare with the safety limit, using different criteria. A simplified wheelset model is used to simulate derailments with different adhesion conditions. The contact force components, including the longitudinal and spin effects, are identified in a steady-state condition on the verge of a derailment. The contact force ratios are used in a three-dimensional (3D) analytical formula to calculate the safety limits. Simulation results obtained with two contact methods were compared with the published results and the safety limit was identified with the two criteria. Results confirm Nadal`s conservative aspect and show that safety 3D analytical formula presents slightly higher safety limits for lower friction coefficients and smaller limits for high friction, in comparison with the simulation results with Fastsim.

Simulation of the Barkhausen Noise using random field Ising model with long-range interaction

We present a method to simulate the Magnetic Barkhausen Noise using the Random Field Ising Model with magnetic long-range interaction. The method allows calculating the magnetic flux density behavior in particular sections of the lattice reticule. The results show an internal demagnetizing effect that proceeds from the magnetic long-range interactions. This demagnetizing effect induces the appearing of a magnetic pattern in the region of magnetic avalanches. When compared with the traditional method, the proposed numerical procedure neatly reduces computational costs of simulation. (c) 2008 Published by Elsevier B.V.

Comparison between detailed model simulation and artificial neural network for forecasting building energy consumption

There are several ways to attempt to model a building and its heat gains from external sources as well as internal ones in order to evaluate a proper operation, audit retrofit actions, and forecast energy consumption. Different techniques, varying from simple regression to models that are based on physical principles, can be used for simulation. A frequent hypothesis for all these models is that the input variables should be based on realistic data when they are available, otherwise the evaluation of energy consumption might be highly under or over estimated. In this paper, a comparison is made between a simple model based on artificial neural network (ANN) and a model that is based on physical principles (EnergyPlus) as an auditing and predicting tool in order to forecast building energy consumption. The Administration Building of the University of Sao Paulo is used as a case study. The building energy consumption profiles are collected as well as the campus meteorological data. Results show that both models are suitable for energy consumption forecast. Additionally, a parametric analysis is carried out for the considered building on EnergyPlus in order to evaluate the influence of several parameters such as the building profile occupation and weather data on such forecasting. (C) 2008 Elsevier B.V. All rights reserved.

Matrix Method for Acoustic Levitation Simulation

A matrix method is presented for simulating acoustic levitators. A typical acoustic levitator consists of an ultrasonic transducer and a reflector. The matrix method is used to determine the potential for acoustic radiation force that acts on a small sphere in the standing wave field produced by the levitator. The method is based on the Rayleigh integral and it takes into account the multiple reflections that occur between the transducer and the reflector. The potential for acoustic radiation force obtained by the matrix method is validated by comparing the matrix method results with those obtained by the finite element method when using an axisymmetric model of a single-axis acoustic levitator. After validation, the method is applied in the simulation of a noncontact manipulation system consisting of two 37.9-kHz Langevin-type transducers and a plane reflector. The manipulation system allows control of the horizontal position of a small levitated sphere from -6 mm to 6 mm, which is done by changing the phase difference between the two transducers. The horizontal position of the sphere predicted by the matrix method agrees with the horizontal positions measured experimentally with a charge-coupled device camera. The main advantage of the matrix method is that it allows simulation of non-symmetric acoustic levitators without requiring much computational effort.

Cement stabilization of runoff residuals: A study of stabilization/solidification of urban rainfall-runoff residuals in type 1 Portland cement by XRD and Si-29 NMR analysis

Urban rainfall-runoff residuals contain metals such as Cr, Zn, Cu, As, Pb and Cd and are thus reasonable candidates for treatment using Portland cement-based solidification-stabilization (S/S). This research is a study of S/S of urban storm water runoff solid residuals in Portland cement with quicklime and sodium bentonite additives. The solidified residuals were analyzed after 28 days of hydration time using X-ray powder diffraction (XRD) and solid-state Si-29 nuclear magnetic resonance (NMR) spectroscopy. X-ray diffraction (XRD) results indicate that the main cement hydration products are ettringite, calcium hydroxide and hydrated calcium silicates. Zinc hydroxide and lead and zinc silicates are also present due to the reactions of the waste compounds with the cement and its hydration products. Si-29 NMR analysis shows that the coarse fraction of the waste apparently does not interfere with cement hydration, but the fine fraction retards silica polymerization.

The impact of El Nino Southern Oscillation phases on off-season maize yield for a subtropical region of Brazil

In recent years, maize has become one of the main alternative crops for the autumn winter growing season in the central-western and southeastern regions of Brazil. However, water deficits, sub-optimal temperatures and low solar radiation levels are common problems that are experienced during this growing season by local farmers. One methodology to assess the impact of variable weather conditions on crop production is the use of crop simulation models. The goal of this study was to evaluate the effect of climate variability on maize yield for a subtropical region of Brazil. Specific objectives for this study were (1) to analyse the effect of El Nino Southern Oscillation (ENSO) on precipitation and air temperature for four locations in the state of Sao Paulo and (2) to analyse the impact of ENSO on maize grown off-season for the same four locations using a crop simulation model. For each site, historical weather data were categorised as belonging to one of three phases of ENSO: El Nino (warm sea surface temperature anomalies in the Pacific), La Nina (cool sea surface temperature anomalies) or neutral, based on an index derived from observed sea surface temperature anomalies. During El Nino, there is a tendency for an increase in the rainfall amount during May for the four selected locations, and also during April, mainly in three of the locations, resulting in an increase in simulated maize yield planted between February 15 and March 15. In general, there was a decrease in the simulated yield for maize grown off-season during neutral years. This study showed how a crop model can be used to assess the impact of climate variability on the yield of maize grown off-season in a subtropical region of Brazil. The outcomes of this study can be very useful for both policy makers and local farmers for agricultural planning and decision making. Copyright (C) 2009 Royal Meteorological Society

Simulation of irrigated Thanzania grass growth based on photothermal units, nitrogen fertilization and water availability

Simulation of irrigated Thanzania grass growth based on photothermal units, nitrogen fertilization and water availability. The mathematical model to predict the forage yield using photothennal units was utilized with success in Elephant grass, Thanzania and Brachiaria niziziensis in the absence of water stress and nitrogen stress. The aim of this study was to propose models to estimate the forage yield of Thanzania grass under different irrigation (25, 50,75, 100 e 125% of ETc) and nitrogen level in various regions of Brazil. As such, models were developed to estimate the dry matter production of Panicum maximum Jacq. frass cv Thanzania in different irrigation and nitrogen levels, using photothermal units. The models were adjusted to doses of 0, 30, 60, 110 and 270 kg of N ha(-1), doses were divided in applications after each evaluation, with a rest cycle of 35 days. The adjusted model presented good performance in predicting dry matter production of Thanzania grass, with r(2) = 0.9999. The results made it possible to verify that the proposed model can be used to predict forage production in different regions of Brazil. It can be estimated, with good precision. The production of Thanzania grass dry matter can be accurately estimated in specific places (in function of latitude and time of year), with the maximum and minimum temperature values.

THE MONTHLY RAINFALL IN THE RIO DE JANEIRO STATE, BRAZIL: SEASONALITY AND TREND

The objective of this work was to carry a descriptive analysis in the monthly precipitation of rainfall stations from Rio de Janeiro State, Brazil, using data of position and dispersion and graphical analyses, and to verify the presence of seasonality and trend in these data, with a study about the application of models of time series. The descriptive statistics was to characterize the general behavior of the series in three stations selected which present consistent historical series. The methodology of analysis of variance in randomized blocks and the determination of models of multiple linear regression, considering years and months as predictors variables, disclosed the presence of seasonality, what allowed to infer on the occurrence of repetitive natural phenomena throughout the time and absence of trend in the data. It was applied the methodology of multiple linear regression to removal the seasonality of these time series. The original data had been deducted from the estimates made by the adjusted model and the analysis of variance in randomized blocks for the residues of regression was preceded again. With the results obtained it was possible to conclude that the monthly rainfall present seasonality and they don`t present trend, the analysis of multiple regression was efficient in the removal of the seasonality, and the rainfall can be studied by means of time series.

Distinct conformational properties determined by implicit and explicit representation of protein-solvent interactions. An analytical and computer simulation study

In the protein folding problem, solvent-mediated forces are commonly represented by intra-chain pairwise contact energy. Although this approximation has proven to be useful in several circumstances, it is limited in some other aspects of the problem. Here we show that it is possible to achieve two models to represent the chain-solvent system. one of them with implicit and other with explicit solvent, such that both reproduce the same thermodynamic results. Firstly, lattice models treated by analytical methods, were used to show that the implicit and explicitly representation of solvent effects can be energetically equivalent only if local solvent properties are time and spatially invariant. Following, applying the same reasoning Used for the lattice models, two inter-consistent Monte Carlo off-lattice models for implicit and explicit solvent are constructed, being that now in the latter the solvent properties are allowed to fluctuate. Then, it is shown that the chain configurational evolution as well as the globule equilibrium conformation are significantly distinct for implicit and explicit solvent systems. Actually, strongly contrasting with the implicit solvent version, the explicit solvent model predicts: (i) a malleable globule, in agreement with the estimated large protein-volume fluctuations; (ii) thermal conformational stability, resembling the conformational hear resistance of globular proteins, in which radii of gyration are practically insensitive to thermal effects over a relatively wide range of temperatures; and (iii) smaller radii of gyration at higher temperatures, indicating that the chain conformational entropy in the unfolded state is significantly smaller than that estimated from random coil configurations. Finally, we comment on the meaning of these results with respect to the understanding of the folding process. (C) 2009 Elsevier B.V. All rights reserved.

XSophe-Sophe-XeprView (R). A computer simulation software suite (v. 1.1.3) for the analysis of continuous wave EPR spectra

The XSophe-Sophe-XeprView((R)) computer simulation software suite enables scientists to easily determine spin Hamiltonian parameters from isotropic, randomly oriented and single crystal continuous wave electron paramagnetic resonance (CW EPR) spectra from radicals and isolated paramagnetic metal ion centers or clusters found in metalloproteins, chemical systems and materials science. XSophe provides an X-windows graphical user interface to the Sophe programme and allows: creation of multiple input files, local and remote execution of Sophe, the display of sophelog (output from Sophe) and input parameters/files. Sophe is a sophisticated computer simulation software programme employing a number of innovative technologies including; the Sydney OPera HousE (SOPHE) partition and interpolation schemes, a field segmentation algorithm, the mosaic misorientation linewidth model, parallelization and spectral optimisation. In conjunction with the SOPHE partition scheme and the field segmentation algorithm, the SOPHE interpolation scheme and the mosaic misorientation linewidth model greatly increase the speed of simulations for most spin systems. Employing brute force matrix diagonalization in the simulation of an EPR spectrum from a high spin Cr(III) complex with the spin Hamiltonian parameters g(e) = 2.00, D = 0.10 cm(-1), E/D = 0.25, A(x) = 120.0, A(y) = 120.0, A(z) = 240.0 x 10(-4) cm(-1) requires a SOPHE grid size of N = 400 (to produce a good signal to noise ratio) and takes 229.47 s. In contrast the use of either the SOPHE interpolation scheme or the mosaic misorientation linewidth model requires a SOPHE grid size of only N = 18 and takes 44.08 and 0.79 s, respectively. Results from Sophe are transferred via the Common Object Request Broker Architecture (CORBA) to XSophe and subsequently to XeprView((R)) where the simulated CW EPR spectra (1D and 2D) can be compared to the experimental spectra. Energy level diagrams, transition roadmaps and transition surfaces aid the interpretation of complicated randomly oriented CW EPR spectra and can be viewed with a web browser and an OpenInventor scene graph viewer.

User guide for shock and blast simulation with the OctVCE code (version 3.5+), Department of Mechanical Engineering Report 2007/13

OctVCE is a cartesian cell CFD code produced especially for numerical simulations of shock and blast wave interactions with complex geometries. Virtual Cell Embedding (VCE) was chosen as its cartesian cell kernel as it is simple to code and sufficient for practical engineering design problems. This also makes the code much more ‘user-friendly’ than structured grid approaches as the gridding process is done automatically. The CFD methodology relies on a finite-volume formulation of the unsteady Euler equations and is solved using a standard explicit Godonov (MUSCL) scheme. Both octree-based adaptive mesh refinement and shared-memory parallel processing capability have also been incorporated. For further details on the theory behind the code, see the companion report 2007/12.

Flow and axial dispersion simulation for traveling axisymmetric Taylor vortices

Numerical experiments using a finite difference method were carried out to determine the motion of axisymmetric Taylor vortices for narrow-gap Taylor vortex flow. When a pressure gradient is imposed on the flow the vortices are observed to move with an axial speed of 1.16 +/- 0.005 times the mean axial flow velocity. The method of Brenner was used to calculate the long-time axial spread of material in the flow. For flows where there is no pressure gradient, the axial dispersion scales with the square root of the molecular diffusion, in agreement with the results of Rosen-bluth et al. for high Peclet number dispersion in spatially periodic flows with a roll structure. When a pressure gradient is imposed the dispersion increases by an amount approximately equal to 6.5 x 10(-4) (W) over bar(2)d(2)/D-m, where (W) over bar is the average axial velocity in the annulus, analogous to Taylor dispersion for laminar flow in an empty tube.