GIS technologies and spatial organization patterns among hunter-gatherers: an archaeological case study from Southern Patagonia (southern South America)

This paper discusses some aspects of hunter-gatherer spatial organization in southern South Patagonia, in later times to 10,000 cal yr BP. Various methods of spatial analysis, elaborated with a Geographic Information System (GIS) were applied to the distributional pattern of archaeological sites with radiocarbon dates. The shift in the distributional pattern of chronological information was assessed in conjunction with other lines of evidence within a biogeographic framework. Accordingly, the varying degrees of occupation and integration of coastal and interior spaces in human spatial organization are explained in association with the adaptive strategies hunter-gatherers have used over time. Both are part of the same human response to changes in risk and uncertainty variability in the region in terms of resource availability and environmental dynamics.

Formation of helical ion chains

We study the nonequilibrium dynamics of the linear to zigzag structural phase transition exhibited by an ion chain confined in a trap with periodic boundary conditions. The transition is driven by reducing the transverse confinement at a finite quench rate, which can be accurately controlled. This results in the formation of zigzag domains oriented along different transverse planes. The twists between different domains can be stabilized by the topology of the trap and under laser cooling the system has a chance to relax to a helical chain with nonzero winding number. Molecular dynamics simulations are used to obtain a large sample of possible trajectories for different quench rates. The scaling of the average winding number with different quench rates is compared to the prediction of the Kibble-Zurek theory, and a good quantitative agreement is found.

Temperature fields in a channel partially filled with a porous material under local thermal non-equilibrium condition-an exact solution

This work examines analytically the forced convection in a channel partially filled with a porous material and subjected to constant wall heat flux. The Darcy–Brinkman–Forchheimer model is used to represent the fluid transport through the porous material. The local thermal non-equilibrium, two-equation model is further employed as the solid and fluid heat transport equations. Two fundamental models (models A and B) represent the thermal boundary conditions at the interface between the porous medium and the clear region. The governing equations of the problem are manipulated, and for each interface model, exact solutions, for the solid and fluid temperature fields, are developed. These solutions incorporate the porous material thickness, Biot number, fluid to solid thermal conductivity ratio and Darcy number as parameters. The results can be readily used to validate numerical simulations. They are, further, applicable to the analysis of enhanced heat transfer, using porous materials, in heat exchangers.

Monitoring and model updating of an FRP pedestrian truss bridge

This paper presents the results of a real bridge field experiment, carried out on a fiber reinforced polymer (FRP) pedestrian truss bridge of which nodes are reinforced with stainless steel plates. The aim of this paper is to identify the dynamic parameters of this bridge by using both conventional techniques and a model updating algorithm. In the field experiment, the bridge was instrumented with accelerometers at a number of locations on the bridge deck, recording both vertical and transverse vibrations. It was excited via jump tests at particular locations along its span and the resulting acceleration signals are used to identify dynamic parameters, such as the bridge mode shape, natural frequency and damping constant. Pedestrianinduced vibrations are also measured and utilized to identify dynamic parameters of the bridge. For a complete analysis of the bridge, a numerical model of the FRP bridge is created whose properties are calibrated utilizing a model updating algorithm. Comparable frequencies and mode shapes to those from the experiment were obtained by the FE models considering the reinforcement by increasing elastic modulus at every node of the bridge by stainless steel plate. Moreover, considering boundary conditions at both ends as fixed in the model resulted in modal properties comparable/similar to those from the experiment. This study also demonstrated that the effect of reinforcement and boundary conditions must be properly considered in an FE model to analyze real behavior of the FRP bridge.

ClimateWell TDC with District Heat

The PolySMART demonstration system SP1b has been modeled in TRNSYS and calibrated against monitored data. The system is an example of distributed cooling with centralized CHP, where the driving heat is delivered via the district heating network. The system pre-cools the cooling water for the head office of Borlänge municipality, for which the main cooling is supplied by a 200 kW compression chiller. The SP1b system thus provides pre-cooling. It consists of ClimateWell TDC with nominal capacity of 10 kW together with a dry cooler for recooling and heat exchangers in the cooling and driving circuits. The cooling system is only operated from 06:00 to 17:00 during working days, and the cooling season is generally from mid May to mid September. The nominal operating conditions of the main chiller are 12/15°C. The main aims of this simulation study were to: reduce the electricity consumption, and if possible to improve the thermal COP and capacity at the same time; and to study how the system would perform with different boundary conditions such as climate and load. The calibration of the system model was made in three stages: estimation of parameters based on manufacturer data and dimensions of the system; calibration of each circuit (pipes and heat exchangers) separately using steady state point; and finally calibration of the complete model in terms of thermal and electrical energy as well as running times, for a five day time series of data with one minute average data values. All the performance figures were with 3% of the measured values apart from the running time for the driving circuit that was 4% different. However, the performance figures for this base case system for the complete cooling season of mid-May to midSeptember were significantly better than those for the monitoring data. This was attributed to long periods when the monitored system was not in operation and due to a control parameter that hindered cold delivery at certain times. 

Numerical Simulation Of Ultrasonic Wave Propagation In Elastically Anisotropic Media

The ultrasonic non-destructive testing of components may encounter considerable difficulties to interpret some inspections results mainly in anisotropic crystalline structures. A numerical method for the simulation of elastic wave propagation in homogeneous elastically anisotropic media, based on the general finite element approach, is used to help this interpretation. The successful modeling of elastic field associated with NDE is based on the generation of a realistic pulsed ultrasonic wave, which is launched from a piezoelectric transducer into the material under inspection. The values of elastic constants are great interest information that provide the application of equations analytical models, until small and medium complexity problems through programs of numerical analysis as finite elements and/or boundary elements. The aim of this work is the comparison between the results of numerical solution of an ultrasonic wave, which is obtained from transient excitation pulse that can be specified by either force or displacement variation across the aperture of the transducer, and the results obtained from a experiment that was realized in an aluminum block in the IEN Ultrasonic Laboratory. The wave propagation can be simulated using all the characteristics of the material used in the experiment evaluation associated to boundary conditions and from these results, the comparison can be made.

Driving Flows in Laboratory Astrophysical Plasma Jets: The Mochi.LabJet Experiment

Thesis (Master's)--University of Washington, 2016-08

Stochastic Modeling of Reversible Biochemical Reaction-Diffusion Systems and High-Resolution Shock-Capturing Methods for Fluid Interfaces

Thesis (Ph.D.)--University of Washington, 2016-08

The Roles of Land and Orography on Precipitation and Ocean Circulation in Global Climate Models

Thesis (Ph.D.)--University of Washington, 2016-08

Study of Tsunami-Induced Fluid and Debris Load on Bridges using the Material Point Method

Thesis (Ph.D.)--University of Washington, 2016-08

Evaluation of a versatile energy auditing tool

Energy auditing can be an important contribution for identification and assessment of energy conservation measures (ECMs) in buildings. Numerous tools and software have been developed, with varying degree of precision and complexity and different areas of use.   This paper evaluates PHPP as a versatile, easy-to-use energy auditing tool and gives examples of how it has been compared to a dynamic simulation tool, within the EU-project iNSPiRe. PHPP is a monthly balance energy calculation tool based on EN13790. It is intended for assisting the design of Passive Houses and energy renovation projects and as guidance in the choice of appropriate ECMs.   PHPP was compared against the transient simulation software TRNSYS for a single family house and a multi-family house. It should be mentioned that dynamic building simulations might strongly depend on the model assumptions and simplifications compared to reality, such as ideal heating or real heat emission system. Setting common boundary conditions for both PHPP and TRNSYS, the ideal heating and cooling loads and demands were compared on monthly and annual basis for seven European locations and buildings with different floor area, S/V ratio, U-values and glazed area of the external walls.   The results show that PHPP can be used to assess the heating demand of single-zone buildings and the reduction of heating demand with ECMs with good precision. The estimation of cooling demand is also acceptable if an appropriate shading factor is applied in PHPP. In general, PHPP intentionally overestimates heating and cooling loads, to be on the safe side for system sizing. Overall, the agreement with TRNSYS is better in cases with higher quality of the envelope as in cold climates and for good energy standards. As an energy auditing tool intended for pre-design it is a good, versatile and easy-to-use alternative to more complex simulation tools.

Computer simulations of soft particles in flow

Understanding the dynamics of blood cells is a crucial element to discover biological mechanisms, to develop new efficient drugs, design sophisticated microfluidic devices, for diagnostics. In this work, we focus on the dynamics of red blood cells in microvascular flow. Microvascular blood flow resistance has a strong impact on cardiovascular function and tissue perfusion. The flow resistance in microcirculation is governed by flow behavior of blood through a complex network of vessels, where the distribution of red blood cells across vessel cross-sections may be significantly distorted at vessel bifurcations and junctions. We investigate the development of blood flow and its resistance starting from a dispersed configuration of red blood cells in simulations for different hematocrits, flow rates, vessel diameters, and aggregation interactions between red blood cells. Initially dispersed red blood cells migrate toward the vessel center leading to the formation of a cell-free layer near the wall and to a decrease of the flow resistance. The development of cell-free layer appears to be nearly universal when scaled with a characteristic shear rate of the flow, which allows an estimation of the length of a vessel required for full flow development, $l_c \approx 25D$, with vessel diameter $D$. Thus, the potential effect of red blood cell dispersion at vessel bifurcations and junctions on the flow resistance may be significant in vessels which are shorter or comparable to the length $l_c$. The presence of aggregation interactions between red blood cells lead in general to a reduction of blood flow resistance. The development of the cell-free layer thickness looks similar for both cases with and without aggregation interactions. Although, attractive interactions result in a larger cell-free layer plateau values. However, because the aggregation forces are short-ranged at high enough shear rates ($\bar{\dot{\gamma}} \gtrsim 50~\text{s}^{-1}$) aggregation of red blood cells does not bring a significant change to the blood flow properties. Also, we develop a simple theoretical model which is able to describe the converged cell-free-layer thickness with respect to flow rate assuming steady-state flow. The model is based on the balance between a lift force on red blood cells due to cell-wall hydrodynamic interactions and shear-induced effective pressure due to cell-cell interactions in flow. We expect that these results can also be used to better understand the flow behavior of other suspensions of deformable particles such as vesicles, capsules, and cells. Finally, we investigate segregation phenomena in blood as a two-component suspension under Poiseuille flow, consisting of red blood cells and target cells. The spatial distribution of particles in blood flow is very important. For example, in case of nanoparticle drug delivery, the particles need to come closer to microvessel walls, in order to adhere and bring the drug to a target position within the microvasculature. Here we consider that segregation can be described as a competition between shear-induced diffusion and the lift force that pushes every soft particle in a flow away from the wall. In order to investigate the segregation, on one hand, we have 2D DPD simulations of red blood cells and target cell of different sizes, on the other hand the Fokker-Planck equation for steady state. For the equation we measure force profile, particle distribution and diffusion constant across the channel. We compare simulation results with those from the Fokker-Planck equation and find a very good correspondence between the two approaches. Moreover, we investigate the diffusion behavior of target particles for different hematocrit values and shear rates. Our simulation results indicate that diffusion constant increases with increasing hematocrit and depends linearly on shear rate. The third part of the study describes development of a simulation model of complex vascular geometries. The development of the model is important to reproduce vascular systems of small pieces of tissues which might be gotten from MRI or microscope images. The simulation model of the complex vascular systems might be divided into three parts: modeling the geometry, developing in- and outflow boundary conditions, and simulation domain decomposition for an efficient computation. We have found that for the in- and outflow boundary conditions it is better to use the SDPD fluid than DPD one because of the density fluctuations along the channel of the latter. During the flow in a straight channel, it is difficult to control the density of the DPD fluid. However, the SDPD fluid has not that shortcoming even in more complex channels with many branches and in- and outflows because the force acting on particles is calculated also depending on the local density of the fluid.

Estudo numérico do escoamento em descarregadores por orifício

Atualmente, no estudo do comportamento de descarregadores de cheias por orifício, é necessário conhecer as pressões ao longo do descarregador, para várias alturas de água na albufeira. Este tipo de estudos tem sido usualmente realizado recorrendo a ensaios experimentais. No entanto, a utilização de modelos numéricos para a simulação do escoamento em estruturas hidráulicas encontra-se numa fase emergente. Neste sentido, a presente dissertação pretende apresentar um estudo numérico relativo ao escoamento em descarregadores de cheias por orifício. Por forma a efetuar o estudo numérico, foi utlizado o programa comercial de CFD FLOW-3D®, reproduzindo um modelo reduzido construído no Laboratório Nacional de Engenharia Civil (LNEC) correspondente a um descarregador por orifício. Para a realização das simulações numéricas foi necessário definir o objeto de estudo, a malha de cálculo, as condições de fronteira e as propriedades do fluido e objeto de modo a reproduzir as situações ensaiadas experimentalmente no LNEC. A proximidade dos resultados experimentais e numéricos para vários níveis de água na albufeira permitiu validar o modelo numérico para este tipo de escoamentos no interior do orifício.

Estudo teórico da tomografia por impedância elétrica

O presente trabalho descreve um estudo sobre a metodologia matemática para a solução do problema direto e inverso na Tomografia por Impedância Elétrica. Este estudo foi motivado pela necessidade de compreender o problema inverso e sua utilidade na formação de imagens por Tomografia por Impedância Elétrica. O entendimento deste estudo possibilitou constatar, através de equações e programas, a identificação das estruturas internas que constituem um corpo. Para isto, primeiramente, é preciso conhecer os potencias elétricos adquiridos nas fronteiras do corpo. Estes potenciais são adquiridos pela aplicação de uma corrente elétrica e resolvidos matematicamente pelo problema direto através da equação de Laplace. O Método dos Elementos Finitos em conjunção com as equações oriundas do eletromagnetismo é utilizado para resolver o problema direto. O software EIDORS, contudo, através dos conceitos de problema direto e inverso, reconstrói imagens de Tomografia por Impedância Elétrica que possibilitam visualizar e comparar diferentes métodos de resolução do problema inverso para reconstrução de estruturas internas. Os métodos de Tikhonov, Noser, Laplace, Hiperparamétrico e Variação Total foram utilizados para obter uma solução aproximada (regularizada) para o problema de identificação. Na Tomografia por Impedância Elétrica, com as condições de contorno preestabelecidas de corrente elétricas e regiões definidas, o método hiperparamétrico apresentou uma solução aproximada mais adequada para reconstrução da imagem.