A moving point approach to model shallow ice sheets: a study case with radially-symmetrical ice sheets

Predicting the evolution of ice sheets requires numerical models able to accurately track the migration of ice sheet continental margins or grounding lines. We introduce a physically based moving point approach for the flow of ice sheets based on the conservation of local masses. This allows the ice sheet margins to be tracked explicitly and the waiting time behaviours to be modelled efficiently. A finite difference moving point scheme is derived and applied in a simplified context (continental radially-symmetrical shallow ice approximation). The scheme, which is inexpensive, is validated by comparing the results with moving-margin exact solutions and steady states. In both cases the scheme is able to track the position of the ice sheet margin with high precision.

A moving-point approach to model shallow ice sheets: a study case with radially symmetrical ice sheets

Predicting the evolution of ice sheets requires numerical models able to accurately track the migration of ice sheet continental margins or grounding lines. We introduce a physically based moving-point approach for the flow of ice sheets based on the conservation of local masses. This allows the ice sheet margins to be tracked explicitly. Our approach is also well suited to capture waiting-time behaviour efficiently. A finite-difference moving-point scheme is derived and applied in a simplified context (continental radially symmetrical shallow ice approximation). The scheme, which is inexpensive, is verified by comparing the results with steady states obtained from an analytic solution and with exact moving-margin transient solutions. In both cases the scheme is able to track the position of the ice sheet margin with high accuracy.

Thermal wave scattering by two overlapping and parallel cylinders

In this paper an analytical solution of the temperature of an opaque material containing two overlapping and parallel subsurface cylinders, illuminated by a modulated light beam, is presented. The method is based on the expansion of plane and cylindrical thermal waves in series of Bessel and Hankel functions. This model is addressed to the study of heat propagation in composite materials with interconnection between inclusions, as is the case of inverse opals and fiber reinforced composites. Measurements on calibrated samples using lock-in infrared thermography confirm the validity of the model.

Fock space for fermion-like lattices and the linear Glauber model

The concept of Fock space representation is developed to deal with stochastic spin lattices written in terms of fermion operators. A density operator is introduced in order to follow in parallel the developments of the case of bosons in the literature. Some general conceptual quantities for spin lattices are then derived, including the notion of generating function and path integral via Grassmann variables. The formalism is used to derive the Liouvillian of the d-dimensional Linear Glauber dynamics in the Fock-space representation. Then the time evolution equations for the magnetization and the two-point correlation function are derived in terms of the number operator. (C) 2008 Elsevier B.V. All rights reserved.

PARALLEL ALGORITHMS FOR MAXIMAL CLIQUES IN CIRCLE GRAPHS AND UNRESTRICTED DEPTH SEARCH

We present parallel algorithms on the BSP/CGM model, with p processors, to count and generate all the maximal cliques of a circle graph with n vertices and m edges. To count the number of all the maximal cliques, without actually generating them, our algorithm requires O(log p) communication rounds with O(nm/p) local computation time. We also present an algorithm to generate the first maximal clique in O(log p) communication rounds with O(nm/p) local computation, and to generate each one of the subsequent maximal cliques this algorithm requires O(log p) communication rounds with O(m/p) local computation. The maximal cliques generation algorithm is based on generating all maximal paths in a directed acyclic graph, and we present an algorithm for this problem that uses O(log p) communication rounds with O(m/p) local computation for each maximal path. We also show that the presented algorithms can be extended to the CREW PRAM model.

Cellular automaton model for molecular traffic jams

We consider the time evolution of an exactly solvable cellular automaton with random initial conditions both in the large-scale hydrodynamic limit and on the microscopic level. This model is a version of the totally asymmetric simple exclusion process with sublattice parallel update and thus may serve as a model for studying traffic jams in systems of self-driven particles. We study the emergence of shocks from the microscopic dynamics of the model. In particular, we introduce shock measures whose time evolution we can compute explicitly, both in the thermodynamic limit and for open boundaries where a boundary-induced phase transition driven by the motion of a shock occurs. The motion of the shock, which results from the collective dynamics of the exclusion particles, is a random walk with an internal degree of freedom that determines the jump direction. This type of hopping dynamics is reminiscent of some transport phenomena in biological systems.

Synthesis of NiO-MgO-ZrO(2) catalysts and their performance in reforming of model biogas

Catalysts containing NiO/MgO/ZrO(2) mixtures were synthesized by the polymerization method in a single step. They were characterized by X-ray diffraction (XRD), temperature programmed reduction (TPR) and physisorption of N(2) (BET) and then tested in the reforming of a model biogas (1.5CH4:1CO(2)) in the presence of air (1.5CH(4) + 1CO(2) + 0.25O(2)) at 750 degrees C for 6h. It was observed that the catalyst Ni20MZ performed better in catalytic processes than the well known catalysts, Ni/ZrO(2) and Ni/MgO, synthesized under the same conditions. The formation of solid solutions, MgO-ZrO(2) and NiO-MgO, increased the rate of conversion of reactants (CH(4) and CO(2)) into synthesis gas (H(2) + CO). The formation of oxygen vacancies (in samples containing ZrO(2) and MgO) seems to promote removal of the coke deposited on the nickel surface. The values of the H(2)/CO ratio were generally found to be slightly lower than stoichiometric, owing to the reverse water gas shift reaction occurring in parallel. (C) 2011 Elsevier B.V. All rights reserved.

Hydropower Energy Simulation Using Mike 11 Model; A Case Study In South Germany's Small Run-Of-River Hydropower Plants

Renewable energy production is a basic supplement to stabilize rapidly increasing global energy demand and skyrocketing energy price as well as to balance the fluctuation of supply from non-renewable energy sources at electrical grid hubs. The European energy traders, government and private company energy providers and other stakeholders have been, since recently, a major beneficiary, customer and clients of Hydropower simulation solutions. The relationship between rainfall-runoff model outputs and energy productions of hydropower plants has not been clearly studied. In this research, association of rainfall, catchment characteristics, river network and runoff with energy production of a particular hydropower station is examined. The essence of this study is to justify the correspondence between runoff extracted from calibrated catchment and energy production of hydropower plant located at a catchment outlet; to employ a unique technique to convert runoff to energy based on statistical and graphical trend analysis of the two, and to provide environment for energy forecast. For rainfall-runoff model setup and calibration, MIKE 11 NAM model is applied, meanwhile MIKE 11 SO model is used to track, adopt and set a control strategy at hydropower location for runoff-energy correlation. The model is tested at two selected micro run-of-river hydropower plants located in South Germany. Two consecutive calibration is compromised to test the model; one for rainfall-runoff model and other for energy simulation. Calibration results and supporting verification plots of two case studies indicated that simulated discharge and energy production is comparable with the measured discharge and energy production respectively.

Temperature Dynamics Investigation At Small And Shallow Lakes Using Hydrodynamic Model

A three-dimensional time-dependent hydrodynamic and heat transport model of Lake Binaba, a shallow and small dam reservoir in Ghana, emphasizing the simulation of dynamics and thermal structure has been developed. Most numerical studies of temperature dynamics in reservoirs are based on one- or two-dimensional models. These models are not applicable for reservoirs characterized with complex flow pattern and unsteady heat exchange between the atmosphere and water surface. Continuity, momentum and temperature transport equations have been solved. Proper assignment of boundary conditions, especially surface heat fluxes, has been found crucial in simulating the lake’s hydrothermal dynamics. This model is based on the Reynolds Average Navier-Stokes equations, using a Boussinesq approach, with a standard k − ε turbulence closure to solve the flow field. The thermal model includes a heat source term, which takes into account the short wave radiation and also heat convection at the free surface, which is function of air temperatures, wind velocity and stability conditions of atmospheric boundary layer over the water surface. The governing equations of the model have been solved by OpenFOAM; an open source, freely available CFD toolbox. As its core, OpenFOAM has a set of efficient C++ modules that are used to build solvers. It uses collocated, polyhedral numerics that can be applied on unstructured meshes and can be easily extended to run in parallel. A new solver has been developed to solve the hydrothermal model of lake. The simulated temperature was compared against a 15 days field data set. Simulated and measured temperature profiles in the probe locations show reasonable agreement. The model might be able to compute total heat storage of water bodies to estimate evaporation from water surface.

Pigeons (Columba livia) are a suitable experimental model for Neospora caninum infection in birds

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

On the dynamic behaviour of a mass supported by a parallel combination of a spring and an elastically connected damper

This paper presents a consistent and concise analysis of the free and forced vibration of a mass supported by a parallel combination of a spring and an elastically supported damper (a Zener model). The results are presented in a compact form and the physical behaviour of the system is emphasised. This system is very similar to the conventional single-degree-of freedom system (sdof)-(Voigt model), but the dynamics can be quite different depending on the system parameters. The usefulness of the additional spring in series with the damper is investigated, and optimum damping values for the system subject to different types of excitation are determined and compared.There are three roots to the characteristic equation for the Zener model; two are complex conjugates and the third is purely real. It is shown that it is not possible to achieve critical damping of the complex roots unless the additional stiffness is at least eight times that of the main spring. For a harmonically excited system, there are some possible advantages in using the additional spring when the transmitted force to the base is of interest, but when the displacement response of the system is of interest then the benefits are marginal. It is shown that the additional spring affords no advantages when the system is excited by white noise. (c) 2007 Elsevier Ltd. All rights reserved.

Path dependence of the quark nonlocal condensate within the instanton model

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

A viscosity model for fluorochloroindate glasses

Viscosity of fluoride glasses is generally Arrhenian between glass transition temperature and crystallization temperature. This dependence on temperature is not observed in some special compositions which have two regions with different activation energies. The viscosity of glasses 40InF(3)-20ZnF(2)-25BaF(2)-10SrF(2)-5LiCl and 40InF(3)-20ZnF(2)-15SrF(2)-20BaF(2)-5NaCl was obtained by the parallel plate method. A theoretical model has been proposed to explain the viscosity data. (C) 1999 Elsevier B.V. B.V. All rights reserved.

An implicit BEM formulation to model strong discontinuities in solids

A boundary element method (BEM) formulation to predict the behavior of solids exhibiting displacement (strong) discontinuity is presented. In this formulation, the effects of the displacement jump of a discontinuity interface embedded in an internal cell are reproduced by an equivalent strain field over the cell. To compute the stresses, this equivalent strain field is assumed as the inelastic part of the total strain. As a consequence, the non-linear BEM integral equations that result from the proposed approach are similar to those of the implicit BEM based on initial strains. Since discontinuity interfaces can be introduced inside the cell independently on the cell boundaries, the proposed BEM formulation, combined with a tracking scheme to trace the discontinuity path during the analysis, allows for arbitrary discontinuity propagation using a fixed mesh. A simple technique to track the crack path is outlined. This technique is based on the construction of a polygonal line formed by segments inside the cells, in which the assumed failure criterion is reached. Two experimental concrete fracture tests were analyzed to assess the performance of the proposed formulation.

Fluorescent lamp model based on equivalent resistances, considering the effects of dimming operation

This paper presents a new methodology for the determination of fluorescent lamp models based on equivalent resistances. One important feature of the proposed methodology is concerned with the inclusion of the filaments into the model, considering the effects of dimming operation on the equivalent resistances. The classical Series-Resonant Parallel-Loaded Half-Bridge inverter is used as the power stage of the ballast. Moreover, the variation of the inverter's switching frequency is the dimming technique assumed for the analyses. Results obtained with a F32T8 lamp indicate that the accuracy of the model is very satisfactory. Thus, the lamp models obtained with the proposed methodology have the potential to serve as an important tool for ballast designers, considering the necessity for evaluating the lamp/ballast compatibility, according to issues concerned to the operating conditions of the electrodes' filaments.