Numerical solution of the PTT constitutive equation for unsteady three-dimensional free surface flows

This work deals with the development of a numerical technique for simulating three-dimensional viscoelastic free surface flows using the PTT (Phan-Thien-Tanner) nonlinear constitutive equation. In particular, we are interested in flows possessing moving free surfaces. The equations describing the numerical technique are solved by the finite difference method on a staggered grid. The fluid is modelled by a Marker-and-Cell type method and an accurate representation of the fluid surface is employed. The full free surface stress conditions are considered. The PTT equation is solved by a high order method, which requires the calculation of the extra-stress tensor on the mesh contours. To validate the numerical technique developed in this work flow predictions for fully developed pipe flow are compared with an analytic solution from the literature. Then, results of complex free surface flows using the FIT equation such as the transient extrudate swell problem and a jet flowing onto a rigid plate are presented. An investigation of the effects of the parameters epsilon and xi on the extrudate swell and jet buckling problems is reported. (C) 2010 Elsevier B.V. All rights reserved.

An unstructured CVFEM and moving interface algorithm for non-Newtonian Hele-Shaw flows in injection molding

Purpose - The purpose of this paper is to develop a novel unstructured simulation approach for injection molding processes described by the Hele-Shaw model. Design/methodology/approach - The scheme involves dual dynamic meshes with active and inactive cells determined from an initial background pointset. The quasi-static pressure solution in each timestep for this evolving unstructured mesh system is approximated using a control volume finite element method formulation coupled to a corresponding modified volume of fluid method. The flow is considered to be isothermal and non-Newtonian. Findings - Supporting numerical tests and performance studies for polystyrene described by Carreau, Cross, Ellis and Power-law fluid models are conducted. Results for the present method are shown to be comparable to those from other methods for both Newtonian fluid and polystyrene fluid injected in different mold geometries. Research limitations/implications - With respect to the methodology, the background pointset infers a mesh that is dynamically reconstructed here, and there are a number of efficiency issues and improvements that would be relevant to industrial applications. For instance, one can use the pointset to construct special bases and invoke a so-called ""meshless"" scheme using the basis. This would require some interesting strategies to deal with the dynamic point enrichment of the moving front that could benefit from the present front treatment strategy. There are also issues related to mass conservation and fill-time errors that might be addressed by introducing suitable projections. The general question of ""rate of convergence"" of these schemes requires analysis. Numerical results here suggest first-order accuracy and are consistent with the approximations made, but theoretical results are not available yet for these methods. Originality/value - This novel unstructured simulation approach involves dual meshes with active and inactive cells determined from an initial background pointset: local active dual patches are constructed ""on-the-fly"" for each ""active point"" to form a dynamic virtual mesh of active elements that evolves with the moving interface.

Numerical Solution of the Upper-Convected Maxwell Model for Three-Dimensional Free Surface Flows

This work presents a finite difference technique for simulating three-dimensional free surface flows governed by the Upper-Convected Maxwell (UCM) constitutive equation. A Marker-and-Cell approach is employed to represent the fluid free surface and formulations for calculating the non-Newtonian stress tensor on solid boundaries are developed. The complete free surface stress conditions are employed. The momentum equation is solved by an implicit technique while the UCM constitutive equation is integrated by the explicit Euler method. The resulting equations are solved by the finite difference method on a 3D-staggered grid. By using an exact solution for fully developed flow inside a pipe, validation and convergence results are provided. Numerical results include the simulation of the transient extrudate swell and the comparison between jet buckling of UCM and Newtonian fluids.

A Robust, Fully Adaptive Hybrid Level-Set/Front-Tracking Method for Two-Phase Flows with an Accurate Surface Tension Computation

We present a variable time step, fully adaptive in space, hybrid method for the accurate simulation of incompressible two-phase flows in the presence of surface tension in two dimensions. The method is based on the hybrid level set/front-tracking approach proposed in [H. D. Ceniceros and A. M. Roma, J. Comput. Phys., 205, 391400, 2005]. Geometric, interfacial quantities are computed from front-tracking via the immersed-boundary setting while the signed distance (level set) function, which is evaluated fast and to machine precision, is used as a fluid indicator. The surface tension force is obtained by employing the mixed Eulerian/Lagrangian representation introduced in [S. Shin, S. I. Abdel-Khalik, V. Daru and D. Juric, J. Comput. Phys., 203, 493-516, 2005] whose success for greatly reducing parasitic currents has been demonstrated. The use of our accurate fluid indicator together with effective Lagrangian marker control enhance this parasitic current reduction by several orders of magnitude. To resolve accurately and efficiently sharp gradients and salient flow features we employ dynamic, adaptive mesh refinements. This spatial adaption is used in concert with a dynamic control of the distribution of the Lagrangian nodes along the fluid interface and a variable time step, linearly implicit time integration scheme. We present numerical examples designed to test the capabilities and performance of the proposed approach as well as three applications: the long-time evolution of a fluid interface undergoing Rayleigh-Taylor instability, an example of bubble ascending dynamics, and a drop impacting on a free interface whose dynamics we compare with both existing numerical and experimental data.

Discontinuous Galerkin approximation of two-phase flows in heterogeneous porous media with discontinuous capillary pressures

We design and investigate a sequential discontinuous Galerkin method to approximate two-phase immiscible incompressible flows in heterogeneous porous media with discontinuous capillary pressures. The nonlinear interface conditions are enforced weakly through an adequate design of the penalties on interelement jumps of the pressure and the saturation. An accurate reconstruction of the total velocity is considered in the Raviart-Thomas(-Nedelec) finite element spaces, together with diffusivity-dependent weighted averages to cope with degeneracies in the saturation equation and with media heterogeneities. The proposed method is assessed on one-dimensional test cases exhibiting rough solutions, degeneracies, and capillary barriers. Stable and accurate solutions are obtained without limiters. (C) 2010 Elsevier B.V. All rights reserved.

Accurate velocity reconstruction for Discontinuous Galerkin approximations of two-phase porous media flows

We consider Discontinuous Galerkin approximations of two-phase, immiscible porous media flows in the global pressure/fractional flow formulation with capillary pressure. A sequential approach is used with a backward Euler step for the saturation equation, equal-order interpolation for the pressure and the saturation, and without any limiters. An accurate total velocity field is recovered from the global pressure equation to be used in the saturation equation. Numerical experiments show the advantages of the proposed reconstruction. To cite this article: A. Ern et al., C R. Acad. Sci. Paris, Ser. 1347 (2009). (C) 2009 Academie des sciences. Published by Elsevier Masson SAS. All rights reserved.

Russian Arctic Ports : Problems with the development and management of tourist flows

The potential changes to the territory of the Russian Arctic open up unique possibilities for the development of tourism. More favourable transport opportunities along the Northern Sea Route (NSR) create opportunities for tourism development based on the utilisation of the extensive areas of sea shores and river basins. A major challenge for the Russian Arctic sea and river ports is their strong cargo transport orientation originated by natural resource extraction industries. A careful assessment of the prospects of current and future tourism development is presented here based on the development of regions located along the shores of the Arctic ocean (including Murmansk and Arkhangelsk oblast, Nenets Autonomous okrug (AO), Yamal-Nenets AO, Taymyr AO, Republic of Sakha, Chykotsky AO). An evaluation of the present development of tourism in maritime cities suggests that a considerable qualitative and quantitative increase of tourism activities organised by domestic tourism firms is made virtually impossible. There are several factors contributing to this. The previously established Soviet system of state support for the investments into the port facilities as well as the sea fleet were not effectively replaced by creation of new structures. The necessary investments for reconstruction could be contributed by the federal government but the priorities are not set towards the increased passenger transportation. Having in mind, increased environmental pressures in this highly sensitive area it is especially vital to establish a well-functioning monitoring and rescue system in the situation of ever increasing risks which come not only from the increased transports along the NSR, but also from the exploitation of the offshore oil and gas reserves in the Arctic seas. The capacity and knowledge established in Nordic countries (Norway, Finland) concerning cruise tourism should not be underestimated and the already functioning cooperation in Barents Region should expand towards this particular segment of the tourism industry. The current stage of economic development in Russia makes it clear that tourism development is not able to compete with the well-needed increase in the cargo transportation, which means that Russia’s fleet is going to be utilised by other industries. However, opening up this area to both local and international visitors could contribute to the economic prosperity of these remote areas and if carefully managed could sustain already existing maritime cities along the shores of the Arctic Ocean.

Rule versus interaction function : evaluating regional aggregations of commuting flows in Sweden

The aim of this paper is to evaluate the performance of two divergent methods for delineating commuting regions, also called labour market areas, in a situation that the base spatial units differ largely in size as a result of an irregular population distribution. Commuting patterns in Sweden have been analyzed with geographical information system technology by delineating commuting regions using two regionalization methods. One, a rule-based method, uses one-way commuting flows to delineate local labour market areas in a top-down procedure based on the selection of predefined employment centres. The other method, the interaction-based Intramax analysis, uses two-way flows in a bottom-up procedure based on numerical taxonomy principles. A comparison of these methods will expose a number of strengths and weaknesses. For both methods, the same data source has been used. The performance of both methods has been evaluated for the country as a whole using resident employed population, self-containment levels and job ratios for criteria. A more detailed evaluation has been done in the Goteborg metropolitan area by comparing regional patterns with the commuting fields of a number of urban centres in this area. It is concluded that both methods could benefit from the inclusion of additional control measures to identify improper allocations of municipalities.

The Nine Circles of Surveillance Hell: An Institutional View of Information Flows and Information Threats in Libraries

In the past few years, libraries have started to design public programs that educate patrons about different tools and techniques to protect personal privacy. But do end user solutions provide adequate safeguards against surveillance by corporate and government actors? What does a comprehensive plan for privacy entail in order that libraries live up to their privacy values? In this paper, the authors discuss the complexity of surveillance architecture that the library institution might confront when seeking to defend the privacy rights of patrons. This architecture consists of three main parts: physical or material aspects, logical characteristics, and social factors of information and communication flows in the library setting. For each category, the authors will present short case studies that are culled from practitioner experience, research, and public discourse. The case studies probe the challenges faced by the library—not only when making hardware and software choices, but also choices related to staffing and program design. The paper shows that privacy choices intersect not only with free speech and chilling effects, but also with questions that concern intellectual property, organizational development, civic engagement, technological innovation, public infrastructure, and more. The paper ends with discussion of what libraries will require in order to sustain and improve efforts to serve as stewards of privacy in the 21st century.

Stochastic Assessment Of Environmental Flows In Semiarid Environments

The regimen of environmental flows (EF) must be included as terms of environmental demand in the management of water resources. Even though there are numerous methods for the computation of EF, the criteria applied at different steps in the calculation process are quite subjective whereas the results are fixed values that must be meet by water planners. This study presents a friendly-user tool for the assessment of the probability of compliance of a certain EF scenario with the natural regimen in a semiarid area in southern Spain. 250 replications of a 25-yr period of different hydrological variables (rainfall, minimum and maximum flows, ...) were obtained at the study site from the combination of Monte Carlo technique and local hydrological relationships. Several assumptions are made such as the independence of annual rainfall from year to year and the variability of occurrence of the meteorological agents, mainly precipitation as the main source of uncertainty. Inputs to the tool are easily selected from a first menu and comprise measured rainfall data, EF values and the hydrological relationships for at least a 20-yr period. The outputs are the probabilities of compliance of the different components of the EF for the study period. From this, local optimization can be applied to establish EF components with a certain level of compliance in the study period. Different options for graphic output and analysis of results are included in terms of graphs and tables in several formats. This methodology turned out to be a useful tool for the implementation of an uncertainty analysis within the scope of environmental flows in water management and allowed the simulation of the impacts of several water resource development scenarios in the study site.

Brazilian international and inter-state trade flows: an exploratory analysis using the gravity model

Recent efforts toward a world with freer trade, like WTO/GATT or regional Preferential Trade Agreements(PTAs), were put in doubt after McCallum's(1995) finding of a large border effect between US and Canadian provinces. Since then, there has been a great amount of research on this topic employing the gravity equation. This dissertation has two goals. The first goal is to review comprehensively the recent literature about the gravity equation, including its usages, econometric specifications, and the efforts to provide it with microeconomic foundations. The second goal is the estimation of the Brazilian border effect (or 'home-bias trade puzzle') using inter-state and international trade flow data. It is used a pooled cross-section Tobit model. The lowest border effect estimated was 15, which implies that Brazilian states trade among themselves 15 times more than they trade with foreign countries. Further research using industry disaggregated data is needed to qualify the estimated border effect with respect to which part of that effect can be attributed to actual trade costs and which part is the outcome of the endogenous location problem of the firm.

The effect of guia exame's ratings on the brazilian fund industry: an analysis of net-worth flows

This paper infers the impact the publication Guia Exame (the guide) has on the Brazilian fund industry, more specifically on the ability the concerned funds develop on attracting new investment. The impact is measured using the event-study analysis based on the variation of net worth subsequently to the event of being rated, according to the methodology applied by the guide to rank the funds. We used five years of fund ratings according to Guia Exame (2000-2004) and analyzed the changes of these funds net worth. We also compared the event amongst different categories of funds. The results found confirm the expected effects according to star rankings and asset manager size in all years.

Capital flows and monetary control under a domestic currency substitution regime: the recent Brazilian experience

The 90s have witnessed a resumption in capital flows to Latin America. due to the conjugation of low interest rates in the US and economic reforms in most LA countries. In Brazil. however. substantial capital flows have becn induced by the extremely high domestic interest rates practiced by the Central Bank as a measure of last reson given the absence of successful stabilization policies. These very high interest rates were needed to prevent capital flight in a context of a surprisingly stable inflation rate above 20% a month. and keep interest bearing govemment securities preferable to foreign assets as money substitutes. We carefully describe how this domestic currency substitution regime (interest bearing govemment securities are substituted for MIas cash holdings) requires the Central Bank to renounce aoy control over monerary aggregates. In this domestic currency substitution regime. hyperinflation is the most likely outcome of an isolated (i.e.. without fiscal adjusanents) attempt by the Brazilian Central Bank to control money.

Balance of payment crises and capital flows : the role of liquidity

A model of externaI CrISIS is deveIoped focusing on the interaction between Iiquidity creation by financiaI intermediaries and foreign exchange collapses. The intermediaries' role of transforming maturities is shown to result in larger movements of capital and a higher probability of crisis. This resembles the observed cycle in capital fiows: large infiows, crisis and abrupt outfiows. The mo deI highlights how adverse productivity and international interest rate shocks can be magnified by the behavior of individual foreign investors linked together through their deposits in the intermediaries. An eventual collapse of the exchange rate can link investors' behavior even further. The basic model is then extended, quite naturally, to study the effects of capital fiow contagion between countries.

The effect of guia exame’s ratings on the brazilian fund industry: an analysis of net-worth flows

This paper infers the impact the publication “Guia Exame” (the guide) has on the Brazilian fund industry, more specifically on the ability the concerned funds develop on attracting new investment. The impact is measured using the event-study analysis based on the variation of net worth subsequently to the event of being rated, according to the methodology applied by the guide to rank the funds. We used five years of fund ratings according to Guia Exame (2000-2004) and analyzed the changes of these funds’ net worth. We also compared the event amongst different categories of funds. The results found confirm the expected effects according to star rankings and asset manager size in all years.