This research studies the green geopolymer incorporated with phase change material (PCM) for energy saving in buildings. First class F fly ash (FA) based-geopolymer binder was studied. In order to improve the mechanical properties, low calcium slag (SG) was added to the FA to produce geopolymer. The effect of different factors including SG content (at different relative amounts FA/SG = 0/100, 25/75, 50/50, 75/25 and 100/0), NaOH solution at different concentrations (7.5, 10 and 15 M), various curing times (1, 2, 4, 7, 14 and 28 days) and curing temperatures (25 (ambient), 45, 60, 75 and 90°C) was investigated. The unit weight and uniaxial compressive strength (UCS) of the geopolymer specimens were measured. Scanning electron microscopy/energy-dispersive X-ray spectroscopy (SEM/EDX) and X-ray diffraction (XRD) were also performed to characterize the microstructure and phase composition of the geopolymer specimens. The results show that the incorporation of SG not only improves the strength of the geopolymer specimens but also decreases the initial water content and thus the NaOH consumption at the same NaOH concentration required for geopolymer production. In addition, the inclusion of SG increases the unit weight of the geopolymer specimens, simply because SG has a much greater specific gravity than FA. The results also show that the strength of the FA/SG-based geopolymer develops rapidly within only 2 days and no obvious gain of the strength after 7 days. The optimum curing temperature (the curing temperature at which the maximum UCS is obtained) at a FA/SG ratio of 50/50 is around 75°C. Second, FA-based geopolymer concrete was synthesized and the effect of different factors including sodium silicate to sodium hydroxide (SS/SH) ratio, aggregate shape, water to fly ash (W/FA) ratio, curing time, water exposure and PCM inclusion on the compressive strength of the geopolymer concrete specimens cured at different ambient temperatures was studied. The results show that the UCS of the specimens increases with higher SS/SH and W/FA ratios up to a certain level and then starts to decrease at higher ratios. The results also indicate that a major portion of the strength of the specimens cured at ambient temperatures develops within the first four weeks. In addition the strength of the FA-based geopolymer concrete is slightly decreased with water exposure and PCM incorporation. Third, the mechanical and thermal properties of geopolymer mortar synthesized with FA and different amount of PCM were studied and the effect of incorporated PCM on the unit weight and UCS of geopolymer mortar was evaluated. SEM imaging was performed to identify the change of micro structure of the geopolymer mortar after incorporation of PCM. The thermal properties of the geopolymer mortar containing different amount of PCM were also characterized using differential scanning calorimetry (DSC) analysis. In addition model tests were performed using small cubicles built with geopolymer mortar slabs containing different amount of PCM to evaluate the effectiveness of geopolymer mortar wall with incorporated PCM in controlling the heat flow and internal temperature. The results indicate that both the unit weight and UCS of the geopolymer mortar decrease slightly after PCM is incorporated, mainly due to the small unit weight and low strength and stiffness of the PCM, respectively. However, the compressive strength of geopolymer mortar containing up to 20% PCM is still sufficiently high for applications in buildings. The results also show that the incorporation of PCM leads to substantial increase of heat capacity and decrease of thermal conductivity of the geopolymer mortar and is very effective in decreasing the temperature inside the cubicles. Finally, a numerical study on the thermal performance of geopolymer with incorporated PCM was carried out. In order to simulate the heat transfer through geopolymer containing PCM, a simplified method was first presented. The influence of phase transition was linked to the energy balance equation through variable specific heat capacity of the PCM-geopolymer. The thermal properties of the geopolymer containing PCM for the numerical analysis were determined using DSC and guarded heat flow (GHF) tests. The simplified method was validated based on the good agreement between the numerical and experimental results. With the validated model, the effect of various factors including the specific heat capacity, thermal conductivity and wall thickness on the thermal performance of PCM-geopolymer walls was studied. Then a modified numerical method was proposed for simulating the whole thermal transfer processes and the simulation results were used to conduct the economic evaluation of PCM-geopolymer walls for energy savings in buildings.

Geothermal energy has become a focal point of the renewable energy revolution. Both shallow and deep types of geothermal energy have the potential to offset carbon emissions, reduce energy costs, and stimulate the economy. Before widespread geothermal exploration and exploitation can occur, both shallow and deep technologies require improvement by theoretical and experimental investigations. This thesis investigated one aspect of both shallow and deep geothermal energy technologies. First, a group of shallow geothermal energy piles was modeled numerically. The model was constructed, calibrated, and validated using available data collected from full-scale in-situ experimental energy piles. Following calibration, the model was parameterized to demonstrate the impact of construction specifications on energy pile performance and cross-sectional thermal stress distribution. The model confirmed the role of evenly spaced heat exchangers in optimal pile performance. Second, experimental methods were used to demonstrate the evolution of a fractured granite permeability as a function of mineral dissolution. Steady-state flow-through experiments were performed on artificially fractured granite cores constrained by 5 MPa pore pressure, 30 MPa confining pressure, and a 120°C temperature. Upstream pore pressures, effluent mineral concentrations, and X-Ray tomography confirmed the hypothesis that fracture asperities dissolve during the flow through experiment, resulting in fracture closure.

L’étude présentée ici entre dans la problématique générale des arcs électriques intervenant dans des applications industrielles telles que le soudage, le découpage, le traitement des déchets. Ce travail constitue une première phase de modélisation de cette étude générale. Le plasma est décrit par un ensemble d'équations de conservation de fluide et de l’électromagnétisme, complétés par des propriétés thermodynamiques et des coefficients de transport appropriés, en formant un système d’équations non linéaires fortement couplées. Ces équations sont écrites en supposant l’équilibre thermodynamique local, une symétrie cylindrique et un écoulement laminaire stationnaire. Ce système d’équations est résolu à l’aide du logiciel commercial FLUENT de type CFD fondé sur l’approche des volumes finis. Pour pouvoir utiliser la partie solveur nous avons résolu notre modèle en utilisant les routines UDF Users-Defined-Function. Dans une première partie, nous présentons la validation du modèle à deux dimensions et à 100 A dans l’argon par des résultats de la littérature. Cette comparaison laisse apparaître un accord satisfaisant sur les profils de température dans la colonne de plasma et des différences dans les zones proches des électrodes dues aux conditions aux limites. Dans une deuxième partie, nous présentons une étude expérimentale, à l’issue de laquelle on constate que les profils de température expérimentaux sont en accord avec ceux du modèle dans la zone de colonne positive. / The study presented here enters the general problem of electric arc involved in industrial applications such as welding, cutting, waste treatment. This work constitutes the first numerical phase of modeling of this general study. The plasma is described by a set of fluid conservation equations, electromagnetic equations complemented by suitable thermodynamic and transport properties, forming a strongly coupled non-linear system. These equations are written assuming local thermodynamic equilibrium, a cylindrical symmetry and steady laminar flow. This system of equations is solved using commercial software FLUENT CFD-type based on finite volume algorithm. To use the solver part we solved our model using UDF macro Users-Defined-Function. In the first part, we present the validation of the two dimensional model in and 100 A in argon with the literature results. A detailed analysis of various characteristic quantities is presented. This comparison reveals a good agreement of the temperature profiles in the column plasma and differences in electrode areas due to boundary conditions. In a second part, we present an experimental study, the experimental temperature profiles are consistent with the model in the column area which means that the model is validated in the column of plasma.

Master of Science / Department of Civil Engineering / Dunja Peric / Piles have been used for many years in civil infrastructure as foundations for buildings, bridges, and retaining walls. Energy piles are thermo-active foundation systems that use geothermal energy for heating and cooling of buildings. Ground source heat is a very attractive, economical, efficient and sustainable alternative to current heating practices. Unlike the air temperature, the temperature below the Earth’s surface remains relatively constant throughout the year, somewhere between 10oC to 15oC below a depth of 6 m to 9 m (Kelly, 2011). This provides an opportunity for construction of thermo-active foundation systems with embedded geothermal loops. The main purpose of such thermo-active system is to transfer deep ground heat to a building through the fluid circulating within the geothermal loop. It is because these thermo-active foundation systems enable heat exchange between the deep ground and the building that is called the heat exchanger pile (HEP). The thermal energy supplied by a HEP can then supplement air-pump-based heating/cooling system. Although heat exchanger piles have been successfully implemented in Europe and Asia, their usage in U.S. remains uncommon. One reason for this might be currently limited understanding of the associated soil-structure interaction, thus unfavorably affecting the design procedures. To this end, a study was undertaken to investigate the predictive capabilities of computational models and to gain a better understanding of the load-transfer mechanisms of energy piles. Thus, coupled thermo-hydro-mechanical computational modeling of a single actual end bearing HEP was carried out for different loading scenarios including thermal and mechanical loads by using the finite element code ABAQUS/Standard 6.13-2. The results of the analyses of the heat exchanger pile with two different types of layered soil profile are presented: isotropic and anisotropic. The computational model was validated and verified successfully against field test results for all considered loading scenarios. Additional analyses were performed to gain a deeper insight into the effects of soil layering and on the behavior of energy piles. It was found that changes in the soil stiffness affected primarily the head displacement and vertical stresses and strains in the pile.

<div>Two-phase flow and heat transfer has a wide variety of applications ranging from nuclear power plants to computer chip cooling. The efficient designs of these systems require a clear understanding of the mechanisms by which two-phase flows enhance heat transfer. With the rapid growth in computing power, Computational Fluid Dynamics is becoming an increasingly reliable predictive tool to understand the physics underlying two-phase flow and heat transfer. We identify the two chief phenomena</div><div>affecting heat transfer in two-phase flows as being the improved convective effect in bubbly flows, and the phase change phenomenon. We examine three key aspects of</div><div>bubbly flows in the present work namely: a) The flow of bubbles near vertical walls, b) the heat transfer associated with a non-condensable bubble rising near a vertical wall, and c) the heat transfer associated with boiling and condensation involving bubbles.</div><div><br></div><div>The first part involves studying the rise velocity of a layer of bubbles rising near a vertical wall. We derive a scaling between the rise velocity based Reynold’s number and the Archimedes number. The second part involves examining the flow pattern around a single bubble rising under the buoyancy effect in a shear flow near a heated wall, and how it affects the heat transfer from the wall. We study the dependence of the fractional improvement in Nusselt number at the wall on various non-dimensional parameters such as the Archimedes number, the Laplace number and the shear rate. Our study shows the existence of an optimum dimensionless shear rate for heat transfer enhancement and a strong dependence between the flow pattern around the bubble and its associated heat transfer enhancement. The third part involves building a numerical model to study flow boiling in micro-channels. We validate the proposed model with two benchmark problems and two experimental studies. The validated numerical tool is then used to understand the effect of varying the micro-channel inlet flow rate on its heat transfer characteristics. This numerical tool is further developed to include a stagnant micro-layer model that can simulate nucleate boiling. We then use it to study the flow boiling characteristics of a line of bubbles undergoing boiling and lift-off in a shear flow. In the end, based on existing literature in the field, we propose future tasks to be undertaken in the area of numerical two-phase flow.<br></div><div><br></div>

O principal objetivo deste trabalho é prever a evolução de margens divergentes desde o início da extensão litosférica, levando-se em consideração a interação entre processos superficiais e tectônicos. Para isto, foi desenvolvido um modelo numérico que acopla isostasia flexural, efeitos térmicos, estiramento litosférico e processos superficiais. A isostasia flexural é simulada através de uma placa elástica fina sobre um fluido invíscido, representando o comportamento flexural da litosfera flutuando sobre a astenosfera. Durante a simulação, a estrutura térmica da litosfera evolui como um resultado da advecção e difusão do calor no interior da Terra. Considera-se que o estiramento da litosfera é acomodado por falhas planas na crosta superior e deformação dúctil na crosta inferior e manto. O modelo de processos superficiais descreve como a paisagem é erodida e como os sedimentos são transportados e depositados nas bacias sedimentares. Através desse modelo numérico, é mostrado que o estiramento litosférico tem uma profunda influência na evolução da migração de escarpas em margens divergentes. Os resultados sugerem que escarpas limitadas por falhas criadas em flancos de rifts por descarregamento mecânico e resposta flexural têm pouca chance de \"sobreviver\" através de recuo erosivo se a crosta inferior sob o flanco do rift foi substancialmente estirada. Nessa configuração, o divisor de drenagem que persiste através do tempo é criado em direção ao continente em uma posição que depende da rigidez flexural da crosta superior. Esse cenário ocorre quando a topografia pré-rift mergulha para o continente, caso contrário a evolução da escarpa é guiada pelo divisor de águas interior pré-existente. Esses conceitos são aplicados no estudo das margens do sudeste da Australia e do sudeste do Brasil, onde o cenário de retração de escarpas através de recuo erosivo mostrou-se improvável. O mesmo modelo numérico foi utilizado para estudar como a passagem de uma anomalia térmica sob a litosfera pode afetar a evolução pós-rift de bacias sedimentares em margens divergentes. Os resultados numéricos mostram que a velocidade da litosfera em relação à anomalia térmica e a rigidez flexural da litosfera oceânica e continental afetam a evolução de bacias sedimentares devido ao soerguimento da superfície relacionado com a expansão térmica da litosfera. Como exemplo, é estudada a possível influência de uma anomalia térmica (Pluma de Trindade?) na evolução das bacias de Campos e Espírito Santo, na margem sudeste brasileira. / The purpose of this work is to predict the evolution of divergent margins since the onset of lithospheric extension, taking into account the interaction between surface and tectonic processes. For this, a numerical model was developed to study the coupling of flexural isostasy, thermal effects, stretching of the lithosphere and surface processes. The flexural isostasy is simulated through a thin elastic plate overlying an inviscid fluid, representing the flexural behavior of the lithosphere floating on the asthenosphere. During the simulation, the thermal structure of the lithosphere evolves as a result of advection and diffusion of heat in the Earths interior. The stretching of the lithosphere is assumed to be accommodated by planar faults in the upper crust and ductile flow in the lower crust and mantle. The surface processes model describes how the landscape is eroded and how the sediments are transported and deposited in the sedimentary basins. The results from this numerical model show that the amount of lithospheric stretching has a profound influence on the evolution of escarpment migration in divergent margins. These results suggest that fault-bounded escarpments created at rift flanks by mechanical unloading and flexural rebound have little potential to survive as retreating escarpments if the lower crust under the rift flank is substantially stretched. In this configuration, a drainage divide that persists through time is created landward in a position that depends on the flexural rigidity of the upper crust. This scenario occurs when the pre-rift topography dips landward, otherwise the evolution of the escarpment is guided by the pre-existing inland drainage divide. These concepts are applied to study the margins of Southeastern Australia and Southeastern Brazil, where the retreating escarpment scenario showed to be unlikely. The same numerical model is used to study how the passage of a thermal anomaly under the lithosphere can affect the post-rift evolution of sedimentary basins in divergent margins. The numerical results show that the velocity of the lithosphere relative to the thermal anomaly and the flexural rigidity of the continental and oceanic lithospheres affect the evolution of sedimentary basins due to surface uplift related to thermal expansion of the lithosphere. As an example, the model is applied to assess the possible influence of a thermal anomaly (Trindade Plume?) on the evolution of the Campos and Esp rito Santo Basins, in Southeastern Brazilian margin.

O presente trabalho tem como objetivo estudar a dinâmica do manto superior em margens continentais através do uso de modelos numéricos que simulam processos convectivos e condutivos no tempo geológico, avaliando-se como a incorporação da convecçãao no estudo da história de subsidência de margens continentais difere do resultado obtido através de modelos puramente condutivos. Como primeiros testes, foram realizadas comparações dos resultados numéricos com soluções analíticas para diferentes valores do número de Rayleigh, verificando-se a validade das soluções computacionais. Também foi feita uma análise da estabilidade da litosfera no tempo geológico para diferentes perfis de viscosidade, servindo como base para a escolha dos parâmetros reológicos do manto para os modelos no contexto de margens divergentes. A partir dos cenários numéricos que melhor reproduziram a estrutura da litosfera terrestre, novos cenários foram criados para simular a evolução térmica e isostática de margens continentais. Como exemplo, utilizou-se dados geofísicos e geológicos extraídos da literatura para a bacia sedimentar do Golfo do Leão, no sudeste da Françaa, com o objetivo de comparar a evolução geodinâmica do presente modelo numérico com outros modelos publicados na literatura. Observou-se que o efeito convectivo astenosférico preserva a estrutura térmica aquecida da margem estirada por mais tempo em comparação com o modelo puramente condutivo. Isso implica que, possivelmente, outros fatores também devem ser levados em consideração como o efeito da geometria tridimensional da margem do Golfo do Leão que pode contribuir para um aumento da subsidência da margem em relação ao modelo obtido no presente trabalho. Adicionalmente, constatou-se que a convecção mantélica pode induzir tensões na base da litosfera que a deslocam dinamicamente ao longo do tempo geológico, podendo influenciar a evolução estratigráfica das bacias sedimentares marginais. São apresentados cerca de 60 cenários geodinâmicos mostrando como a variação da estrutura reológica do manto influencia a evolução térmica da litosfera e consequentemente, a história de subsidência da margem. / This work aims to study the dynamics of the upper mantle in continental margins by using numerical models that simulate convective and conductive processes in geological time scale. It was evaluated the contribution of convection and conduction for subsidence history of sedimentary basins. As first tests, simple numerical scenarios with different Rayleigh number were compared with analytic solutions, verifying the validate of the computational solutions. These numerical experiments were followed by the analysis of the lithospheric stability in the geological time scale for different values of viscosity. These experiments were used as a base for the choice of the rheological parameters of the mantle for the models in the context of divergent margins. From the numerical scenarios that better reproduced the lithospheric structure of the Earth, new scenarios were created to simulate the thermal and isostatic evolution of continental margins. As an example, geophysical and geological data extracted from the literature for the sedimentary basin of the Gulf of Lion, Southeastern France, were compared with the results of different geodynamic models published in the literature and with the numerical scenarios obtained in the present work. We observed that the effect of the astenospheric convection preserves the thermal structure of the stretched margin for a long time in comparison with purely conductive models. This implies that, possibly, other processes must be taken into account, such as the effect of the three-dimensional geometry of the Gulf of Lion margin that may contribute to a higher subsidence of the margin than the one obtained in the present work. Additionally, it was observed that mantle convection may induce stress at the base of the lithosphere that dynamically moves it in the geological time, and may influence the stratigraphic evolution of sedimentary basins. It is presented about 60 scenarios showing how the variation of the rheological structure of the mantle is taken into account in the thermal evolution of the lithosphere and consequently in the subsidence history of the margin.

No contexto da necessidade de atracação de embarcações de maior porte em terminais portuários brasileiros, tanto privativos como públicos, considerando um aumento das exportações, intervenções de engenharia que possam significar ganhos nas dimensões das máximas embarcações permitidas podem representar substanciosos benefícios econômicos. Assim, este trabalho apresenta uma abordagem do conceito de Fundo Náutico, definido como profundidade até a qual as embarcações podem navegar sem que haja efeitos adversos na manobrabilidade das embarcações e danos no casco do navio, cujo objetivo é aumentar o calado máximo de embarcações em espaços náuticos cujo fundo apresente camadas de lama fluida. Devido às suas características reológicas, a lama fluida, de modo geral, permite que as embarcações naveguem com reduzida ou até negativa folga sob a quilha, respeitando-se o conceito estabelecido pelo Fundo Náutico. Exposta a importância sobre este conceito e suas peculiaridades relacionadas às características da lama fluida, como características reológicas e métodos para sua determinação, discutem-se as variações na espessura da camada de lama fluida e as variáveis ambientais que as condicionam, a partir de análise de levantamentos batimétricos de dupla frequência, dados densimétricos e aplicação de modelagem numérica, para a região da bacia de evolução do Complexo Portuário de Itajaí. Para esta região, observaram-se camadas de lama fluida entre 0,5 e 2,5 m, sendo que os resultados obtidos na modelagem indicam que esta variabilidade está associada ao regime de descarga sólida do rio Itajaí-açu. Além disto, a partir de uma análise comparativa entre os dados batimétricos e densimétricos, foi possível estabelecer que o fundo náutico na região pode estar associado a um valor de massa específica entre 1150 e 1200 kg/m³. / Considering the demand for mooring larger ships at Brazilian port terminals, both private and public, in a scenario of growing exports, engineering interventions that can provide improvements in the vessel\'s maximum allowed dimensions can represent significant profits. Hence, this work presents an approach of Nautical Bottom, defined as the maximum depth in which ships can navigate without significant adverse effects in ship control and maneuverability without physical damage, with the goal of raising the maximum ship draft allowed in nautical spaces with fluid mud beds. Due to its rheological properties, fluid mud, in general, allows for vessels navigate with low or negative under keel clearance, respecting the established Nautical Bottom concepts. Once highlighted the importance of these concepts, this study presents an analysis about fluid mud and Nautical Bottom concept analysis at Itajaí Port Complex (Santa Catarina, Brazil), one of the most important port areas in Brazil. This analysis was based on bathymetric surveys, density measurements and numerical modeling. In addition, fluid mud layer thickness at port areas can possibly vary according to hydrodynamics and sedimentologic variations. Fluid mud layer thickness varies from 0.5 to 2.5 meters at Itajai Port Complex turning basin. From numerical modeling results, it is possible to establish that the thickness variations of fluid mud layers are associated with fluvial sediment transport rates. Also, from the comparison between low frequencies bathymetric surveys and specific density measurements, it was possible to conclude that the reference specific density at the interest area is about 1150-1200 kg/m³, being very similar to reference densities adopted worldwide.

O Princeton Ocean Model foi adaptado ao Canal de São Sebastião (CSS) para estudar as variações sazonais de sua circulação e estrutura termohalina. Três grades numéricas foram aninhadas. A de menor resolução na Plataforma Continental Sudeste (PCSE) e a de maior resolução no Canal de São Sebastião (CSS) com uma grade de resolução intermediária na região adjacente ao canal (PCI). O modelo numérico partiu de condições termohalinas médias sazonais e teve como forçantes, ventos, fluxos de calor, de sal e de radiação de ondas curtas mensais. Nessas condições médias, o modelo representou razoavelmente bem as condições típicas de primavera, verão, outono e inverno, preenchendo o fundo do CSS com a Água Central do Atlântico Sul na primavera e no verão. No outono e no inverno esta massa de água não se encontra no canal, porém, seus sinais mais fracos são obtidos no outono sendo que nesta estação são encontrados os sinais mais fortes da Água Tropical. O modelo aponta o sul do canal, no fundo, ao lado da Ilha de São Sebastião como entrada preferencial de águas mais frias e a simulação numérica da passagem de uma frente fria pela Plataforma Continental Sudeste sugere a rápida resposta das águas do canal com o recuo para o largo da Água Central do Atlântico Sul e pronto retorno assim que a frente deixa a plataforma. Ventos de nordeste na grade da PCSE são imprescindíveis para que a Água Central do Atlântico Sul penetre o Canal de São Sebastião; todos os experimentos com ventos de outras direções nesta grade ou ventos de nordeste somente nas grades média e do CSS não colocaram esta massa de água no canal. A circulação de fundo obtida no CSS é basicamente para nordeste e associada à intrusão da ACAS forçada em primeira instância pelo vento de nordeste na PCSE e em um segundo momento pela força do gradiente de pressão (com destaque para a componente baroclínica), sempre maior na entrada sul do que na entrada norte e sempre maior no verão do que no inverno. A circulação superficial é para sudoeste com relaxamento no outono, intensificando-se em direção ao verão com máximo nesta estação. / The Princeton Ocean Model was adapted to São Sebastião Channel (SSC) so as to study the seasonal changes of its circulation and thermohaline structure. Three numerical grids were nested. The coarse grid on Southeast Continental Shelf (SCS) and the fine grid on São Sebastião Channel with an intermediary grid on the region adjacent to the channel. The numerical model started from average seasonal thermohaline conditions and average monthly data forcing such as wind, heat flux, salt flux and short wave radiation flux. Within these average conditions, the numerical model simulated reasonably the typical conditions of spring, summer, autumn and winter, filling the bottom channel with South Atlantic Central Water (SACW) in spring and in summer. In autumn and in winter this mass of water is not found, nevertheless its weaker signs appear in autumn, which season one finds the stronger signs of Tropical Water (never more than 50%). The model points out the bottom south entrance of the channel, next to the São Sebastião Island (SSI) as the natural gate of colder water and the numerical simulation of a cold front through SCS suggests the quick answer of the water channel with the falling back offshore of the SACW and immediate return when the cold front is gone. Northeasterly winds on the SCS grid are essential so that the SACW enters the SSC; all the experiments with another direction winds in this grid or northeasterly winds only in the intermediate and in the fine grid failed to get in the SACW in the channel. The bottom circulation obtained in SSC is essentially to northeast and associated to the intrusion of the SACW forced in first instance by the northeasterly winds in SCS and in second instance by the pressure force (with emphasis on the baroclinic component) always bigger at the south entrance than in the north entrance and always bigger in summer than in the winter. The superficial currents are southeastward with weakening in autumn and intensification towards summer with maximum in this season.

O município de São José do Rio Preto (SJRP), localizado na porção noroeste do estado de São Paulo, é extremamente dependente da água subterrânea, tanto para abastecimento público quanto privado, justificando-se o grande número de poços existentes, que excede a mais de dois mil, sobretudo no Sistema Aquífero Bauru (SAB). Embora este recurso seja fundamental para o abastecimento da cidade, 79% dos poços não possuem outorga. O aumento do número de poços, sem o devido controle, explotando o aquífero, pode justificar e intensificar os problemas de forte abatimento da superfície potenciométrica na região central da área urbana do município. Nesse sentido, esta pesquisa objetivou desenvolver um aplicativo computacional, na forma de modelagem numérica de simulação de fluxos subterrâneos, para subsidiar a explotação sustentável do SAB. Para o sucesso do modelo numérico, foi realizada toda a caracterização hidrogeológica e definido o modelo conceitual de fluxo das águas subterrâneas. O modelo numérico foi calibrado em regime estacionário para duas situações: pré-explotação e bombeamento atual; e, posteriormente, foram gerados três diferentes cenários futuros de bombeamento: a) explotação duas vezes maior do que a atual (Q = 133.206 m³/dia); b) explotação três vezes maior do que a atual (Q = 199.809 m³/dia); c) explotação apenas dos poços públicos (Q = 39.662 m³/dia). Houve médias de rebaixamento na área urbana de 15 m, 42 m e 8 m, respectivamente. Rebaixamentos mais proeminentes do nível freático ocorrem na área central do município, onde se concentra 80% dos poços. Dessa forma, recomendam-se algumas ações direcionadas ao poço, à proteção do aquífero e ao usuário. / The city of São José do Rio Preto (SJRP), located in the northwestern portion of the state of São Paulo, is highly dependent of groundwater for both public and private supply, that justifies the large number of existing wells, which exceeds more than 2000, especially in the Bauru Aquifer System (BAS). Although this resource is essential for the city supply, 79% of the wells do not have grants. The increase of wells, without the correct control, exploring the aquifer can justify and intensify the problems of strong relief in the central region of the urban area by city. Then, this research aimed to develop a computational application how numerical simulation modeling of groundwater flow, to support the rational exploitation of the BAS. For the success of the numerical model, the entire hydrogeological characterization was performed and defined the conceptual model of groundwater flow. The numerical model was calibrated in steady state for two situations: pre-exploitation and current pumping, and then were generated three different future scenarios of pumping: a) exploitation twice larger than the current (Q = 133,206 m³/day); b) exploitation three times larger than the current (Q = 199 809 m³/day); and c) exploitation of only public wells (Q = 39,662 m³/day). There were averages of debasement in the urban area of 15 m, 42 m and 8 m, respectively. Most prominent debasements of groundwater level occurred in the central area of the city, which the concentrates is 80 % of the wells. Therefore, some actions are recommended directed to the well, to protect the aquifer and the user.