Ce travail est consacré au développement de schémas numériques pour approcher les solutions de modèles d'écoulement type “depth averaged”. Dans un premier temps, nous détaillons la construction d'approches Volumes Finis pour le système Shallow Water avec termes sources sur maillages non structurés. En se basant sur une reformulation appropriée des équations, nous mettons en place un schéma équilibré et préservant la positivité de la hauteur d'eau, et suggérons des extensions MUSCL adaptées. La méthode est capable de gérer des topographies irrégulières et exhibe de fortes propriétés de stabilité. L'inclusion des termes de friction fait l'objet d'une analyse poussée, aboutissant à l'établissement d'une propriété type “Asymptotic Preserving” à travers l'amélioration d'un autre récent schéma Volumes Finis. La seconde composante de cette étude concerne les méthodes Elements Finis type Galerkin discontinu. Certaines des idées avancées dans le contexte Volumes Finis sont employées pour aborder le système Shallow Water surmaillages triangulaires. Des résultats numériques sont exposés et la méthode se révèle bien adaptée à la description d'une large variété d'écoulements. Partant de ces observations nous proposons finalement d'exploiter ces caractéristiques pour étendre l'approche à une nouvelle famille d'équations type Green-Nadghi. Des validations numériques sont également proposées pour valider le modèle numérique. / This work is devoted to the development of numerical schemes to approximatesolutions of depth averaged flow models. We first detail the construction of Finite Volume approaches for the Shallow Water system with source terms on unstructured meshes. Based on a suitable reformulation of the equations, we implement a well-balanced and positive preserving approach, and suggest adapted MUSCL extensions. The method is shown to handle irregular topography variations and demonstrates strong stabilities properties. The inclusion of friction terms is subject to a thorough analysis, leading to the establishment of some Asymptotic Preserving property through the enhancement of another recent Finite Volume scheme.The second aspect of this study concerns discontinuous Galerkin Finite Elementmethods. Some of the ideas advanced in the Finite Volume context areemployed to broach the Shallow Water system on triangular meshes. Numericalresults are exposed and the method turns out to be well suited to describe a large variety of flows. On these observations we finally propose to exploit its features to extend the approach to a new family of Green-Nadghi equations. Numerical experiments are also proposed to validate this numerical model.

Made available in DSpace on 2015-04-22T21:58:27Z (GMT). No. of bitstreams: 1 Elizabeth Tavares Pimentel.pdf: 1766165 bytes, checksum: 014e86903e1eaf04dbad226c2450a908 (MD5) Previous issue date: 2009-04-06 / From October 2007 to September 2008 a geothermal monitoring experiment was conducted at depths of 0.02 m, 0.5 m and 1.0 m to quantify the variations of temperature, thermal conductivity and the shallow geothermal heat flow at places with and without vegetation cover in the Humaitá city, Amazonas. The influence of the vegetation cover on the shallow geothermal system was observed in the sites studied. There were variations of monthly average values of temperature between the places with and without vegetation cover. During the "dry" period, this variation was up to 6.01ºC at the depth of 0.02 m, and 2.84ºC at the depth of 1.0 m. During the "rainy" period, however, the variation was up to 2.94ºC, at the depth of 0.02 m, and 2.51ºC at the depth of 1.0 m. The difference of the daily extreme values of temperature between sites with and without vegetation cover were 3.97ºC during the "rainy" period and 9.63ºC during the "dry" period, at the depth of 0.02 m. It was noticed that at 06:00 PM the magnitude of the temperature remained high compared to other times on the day studied. The values of the thermal conductivity were 0.54 W/mºC during the "dry" period, and 1.23 W/mºC during the "rainy" period. The values of the shallow geothermal flows at depths of 0.5 m and 1.0 m, were 2.51 W/m² and 0.64 W/m², respectively. These values are 10³ larger than the terrestrial heat flow in the region. The thermal variations at0.5 m to 1.0 m depth are influenced by external sources that reach the surface and cannot be neglected. The thermal variations recorded in this work are important and fundamental to better understanding the shallow geothermal structure in the southern, part of Amazonas state, and they also contribute as input to models that allow the mitigation or elimination of the effects caused by anthropogenic actions / No período de outubro de 2007 a setembro de 2008 foi realizado monitoramento geotermal, às profundidades de 0,02 m, 0,5 m e 1,0 m, em locais com e sem cobertura vegetal, na cidade de Humaitá (AM), a fim de quantificar as variações de temperatura, condutividade térmica e fluxo geotermal raso local. Constatou-se a influência da cobertura vegetal sobre o regime geotermal raso na região estudada. Houve variação dos valores médios mensais da temperatura nos locais cc e sc. No período "seco", esta variação foi de até 6,01ºC à profundidade de 0,02 m, e de 2,84ºC à profundidade de 1,0 m; já no período "chuvoso", a variação foi de até 2,94ºC a 0,02 m de profundidade e de 2,51ºC à profundidade de 1,0 m. Na profundidade de 0,02 m, a diferença entre os valores diários máximos, nos locais cc e sc, foi de 3,97ºC no período "chuvoso" e de 9,63ºC no período "seco". Às 18 h, as magnitudes da temperatura permaneceram elevadas em relação aos outros horários estudados. Os valores de condutividade térmica foram de 0,54 W/mºC no período "seco" e de 1,23 W/mºC no período "chuvoso". Os valores do fluxo geotermal raso, às profundidades de 0,5 m e 1,0 m, variaram até 2,51 W/m² e 0,64 W/m², respectivamente. Tais valores são da ordem de 10³ acima do valor do fluxo térmico terrestre profundo na região. As variações termais a 0,5 m e a 1,0 m de profundidade são influenciadas por fontes externas que atingem a superfície e não podem ser negligenciadas. As variações térmicas registradas neste trabalho são importantes e fundamentais para o melhor conhecimento da estrutura geotermal rasa na cidade de Humaitá (AM), como também, contribuem para a elaboração de modelos que possibilitem mitigar ou eliminar os efeitos causados por ações antrópicas

We develop and analyse explicit in time schemes for the computation of compressible flows, based on staggered in space. Upwinding is performed equation by equation only with respect to the velocity. The pressure gradient is built as the transpose of the natural divergence. For the barotropic Euler equations, the velocity convection is built to obtain a discrete kinetic energy balance, with residual terms which are non-negative under a CFL condition. We then show that, in 1D, if a sequence of discrete solutions converges to some limit, then this limit is the weak entropy solution. For the full Euler equations, we choose to solve the internal energy balance since a discretization of the total energy is rather unnatural on staggered meshes. Under CFL-like conditions, the density and internal energy are kept positive, and the total energy cannot grow. To obtain correct weak solutions with shocks satisfying the Rankine-Hugoniot conditions, we establish a kinetic energy identity at the discrete level, then choose the source term of the internal energy equation to recover the total energy balance at the limit. More precisely speaking, we prove that in 1D, if we assume the L∞ and BV-stability and the convergence of the scheme, passing to the limit in the discrete kinetic and internal energy equations, we show that the limit of the sequence of solutions is a weak solution. Finally, we consider the computation of radial flows, governed by Euler equations in axisymetrical (2D) or spherical (3D) coordinates, and obtain similar results to the previous sections. In all chapters, we show numerical tests to illustrate for theoretical results. / We develop and analyse explicit in time schemes for the computation of compressible flows, based on staggered in space. Upwinding is performed equation by equation only with respect to the velocity. The pressure gradient is built as the transpose of the natural divergence. For the barotropic Euler equations, the velocity convection is built to obtain a discrete kinetic energy balance, with residual terms which are non-negative under a CFL condition. We then show that, in 1D, if a sequence of discrete solutions converges to some limit, then this limit is the weak entropy solution. For the full Euler equations, we choose to solve the internal energy balance since a discretization of the total energy is rather unnatural on staggered meshes. Under CFL-like conditions, the density and internal energy are kept positive, and the total energy cannot grow. To obtain correct weak solutions with shocks satisfying the Rankine-Hugoniot conditions, we establish a kinetic energy identity at the discrete level, then choose the source term of the internal energy equation to recover the total energy balance at the limit. More precisely speaking, we prove that in 1D, if we assume the L∞ and BV-stability and the convergence of the scheme, passing to the limit in the discrete kinetic and internal energy equations, we show that the limit of the sequence of solutions is a weak solution. Finally, we consider the computation of radial flows, governed by Euler equations in axisymetrical (2D) or spherical (3D) coordinates, and obtain similar results to the previous sections. In all chapters, we show numerical tests to illustrate for theoretical results.

Ninety-eight percent of water taken from the Mississippi River Shallow Alluvial Aquifer, hereafter referred to as “the aquifer” or “MRVA,” is used by the agricultural industry for irrigation. Mississippi Delta agriculture is increasingly using more water from the MRVA and the aquifer has been losing about 300,000 acreeet per year. This research expands on previous work in which a model was developed that simulates the effects of climatic variability, crop acreage changes, and specific irrigation methods on consequent variations in the water volume of the MRVA. This study corrects an identified problem by replacing total growing season precipitation with an irrigation demand driver based on evaporation and crop coefficients and changing the time scale from the entire growing season to a daily resolution. The calculated irrigation demand, as a climatological driver for the model, captures effective precipitation more precisely than the initial growing season precipitation driver. Predictive equations resulting from regression analyses of measured versus calculated irrigation water use showed R2 and correlations of 0.33 and 0.57, 0.77 and 0.88, 0.71 and 0.84, and 0.68 and 0.82 for cotton, corn, soybeans and rice, respectively. Ninetyive percent of the predicted values fall within a range of + or - about 23,000 acreeet, an error of about 10-percent. The study also adds an additional conservation strategy through the use of surface water from onarm reservoirs in lieu of groundwater. Analyses show that climate could provide the entire water need of the plants in 70-percent of the years for corn, 65-percent of the years for soybeans and cotton, and even 5-percent of the years for rice. Storing precipitation in onarm structures is an effective way to reduce reliance of Delta producers on groundwater. If producers adopted, at a minimum, the 97.5:2.5 ratio suggested management practice, this minimal management strategy could potentially conserve 48-percent, 35-percent and 42-percent of groundwater for cotton, corn and soybeans, respectively. Even in extreme drought years such as 2007, cotton, corn and soybeans produced under the 97.5:2.5 management strategy could conserve 32-percent, 46-percent and 38-percent of groundwater, respectively.

Convective snow bands develop commonly over the open water surface of lakes or seas when cold airgets advected from a continent. Enhanced heat and moisture fluxes from the comparatively warm waterbody trigger shallow convection and an unstable boundary layer builds up. Relatively strong wind canorganize this convection into wind-parallel quasi-stationary cloud bands with moving individual cells.Depending on various factors like the horizontal wind, the vertical shear or the shape of the coast, thosecloud bands can form of different strength and structure. When the air mass meets the coast orographicforcing causes horizontal convergence and vertical lifting intensifies the precipitation at the coast. If thewind direction stays constant for several days a single snow band would accumulate its precipitation ina very restricted region and cause locally a significant increase in snow depth. This process leads in thecold season repeatedly to severe precipitation events at the Swedish east coast. Large amounts of snowalong with strong wind speeds can cause serious problems for traffic and infrastructure.Two different cases of convective snow bands in the Baltic Sea area were selected to simulate theassociated atmospheric conditions with a total of five different model systems. The atmosphere climatemodel RCA has been used independently at default settings as well as with increased resolution on avertical and a horizontal scale and furthermore coupled either to the ice-ocean model NEMO or the wavemodel component WAM.Comparing all models the crucial parameters like wind, temperature, heat fluxes, and precipitationvary generally in a reasonable range. However, the model systems show systematical differences amongthemselves. The strongest 10 meter wind speeds can be observed for both RCA models with increasedresolution. The RCA-WAM simulation shows its wind enhancement during the snow band event witha time shift to the other models by several hours. The mean directional wind shear above the Gulf ofBothnia, the snow band’s region of origin, is for all models small. The warmest sea surface temperaturesare reached by the RCA-NEMO simulation, which as a result also stands out for its most intense heatfluxes in both sensible and latent heat. Both high resolution RCA models as well as RCA-NEMO givethe most remarkable local precipitation rates. The original RCA and RCA-WAM simulate significantlyless snowfall. Local comparison with SMHI station measurements show that the models represent thetrend of wind, temperature and precipitation evolution well. However, all models decelerate the air masstoo rapidly when meeting the coast. Moreover, it remains a challenge to simulate the exact time andlocation of the extreme precipitation.The coupling of the atmosphere model with the ice-ocean model as well as the increased resolution ofthe atmospheric component have been observed to show great improvements in the model performanceand are suggested for future research work to be used in combination with each other for the regionalmodelling of convective snow bands in the Baltic Sea area.

International Telemetering Conference Proceedings / October 22-25, 2001 / Riviera Hotel and Convention Center, Las Vegas, Nevada / Undersea warfare (USW) was perceived as a large-area, deep-water operation in the past therefore Fleet USW training ranges were designed to meet these requirements. Currently the bigger threat is the likelihood of regional conflict throughout the world by aggressive nations in littoral waters. The U.S. Navy must stand ready to respond to these regional conflicts when national interests are threatened. Consequently, naval forces must train to operate in the littoral environments where such regional conflicts are likely to occur. The West Cost Shallow Water Undersea Warfare Training Range (WC SWUWTR) is being developed to provide this training.

Many applications in engineering practice can be described by thehyperbolic partial differential equations (PDEs). Numerical modeling of this typeof equations often involves large gradients or shocks, which makes it achallenging task for conventional numerical methods to accurately simulate suchsystems. Thus developing accurate and efficient shock capturing numericalschemes becomes important for the study of hyperbolic equations.In this dissertation, a detailed study of the numerical methods for linearand nonlinear unsteady hyperbolic equations was carried out. A new finitedifference shock capturing scheme of finite volume style was developed. Thisscheme is based on the high order Pad?? type compact central finite differencemethod with the weighted essentially non-oscillatory (WENO) reconstruction toeliminate non-physical oscillations near the discontinuities while maintain stablesolution in the smooth areas. The unconditionally stable semi-implicit Crank-Nicolson (CN) scheme is used for time integration.The theoretical development was conducted based on one-dimensionalhomogeneous scalar equation and system equations. Discussions were alsoextended to include source terms and to deal with problems of higher dimension.For the treatment of source terms, Strang splitting was used. For multidimensionalequations, the ?? -form Douglas-Gunn alternating direction implicit(ADI) method was employed. To compare the performance of the scheme withENO type interpolation, the current numerical framework was also applied usingENO reconstruction.The numerical schemes were tested on 1-D and 2-D benchmark problems,as well as published experimental results. The simulated results show thecapability of the proposed scheme to resolve discontinuities while maintainingaccuracy in smooth regions. Comparisons with the experimental results validatethe method for dam break problems. It is concluded that the proposed scheme isa useful tool for solving hyperbolic equations in general, and from engineeringapplication perspective it provides a new way of modeling open channel flows.

本研究旨在探討台灣高中學生學習英語字彙時所使用之學習策略，除了解各種不同字彙學習策略的使用頻率情形，找出較常和較少使用之字彙學習策略，探求字彙學習策略的使用和字彙成就表現的關係外，並找出字彙成就高者和字彙成就低者對學習策略使用之不同處，藉此，俾能提供老師有效教學及學生有效學習字彙之依據。 本研究以271位高三學生為調查對象，研究工具為字彙學習策略問卷和字彙測驗。分析方法主要採量化分析，包括描述性統計、皮爾遜積差相關分析及獨立樣本t檢定。 本研究結果如下：(一)最常使用的是認知策略，最少使用的是社會策略。(二)學生偏愛機械記憶或與拼字及發音有關之策略。(三)學生的策略使用和其字彙成就表現有明顯相關。(四)字彙學習成就較高者與字彙成就學習較低者在策略使用上亦有明顯差異。高字彙學習成就者較喜重複念誦記憶單字的策略，並傾向於有上下文的情境中學習單字；低字彙學習成就者偏愛重複書寫的策略，並傾向學習孤立的單字。 / This study investigated the vocabulary learning strategies used by senior high students in Taiwan. The aim was (a) to investigate the frequencies of different vocabulary learning strategies use; (b) to find out the most and the least frequently used strategies; (c) to explore the relationship between strategy use and vocabulary size; and (d) to identify the differences in strategy use between good and poor learners. A total of 271 senior high students in northern Taiwan participated in the study. A vocabulary learning strategies questionnaire and a Vocabulary Levels Test were administered to the participants as instruments. The collected data were analyzed by using SPSS version 10.0, including descriptive statistics, Pearson product-moment correlations, and independent-samples t-tests. The results are summarized as follows. (1) Cognitive strategies were reported as the most frequently used strategies while metacognitive and social strategies the least frequently used strategies. (2) Students favored strategies related to “rote repetition” or “the form of a word.” (3) Students’ strategy use was significantly correlated with their vocabulary size. (4) There was a significant difference in strategy use between good and poor learners. Good learners tended to learn words in context while poor learners tended to learn words in isolation. The findings of the study suggest that teachers should introduce students to a variety of vocabulary learning strategies, especially deep strategies and strategies related to context.

Originally developed in the 1970s for the optimal control of systems governed by partial differential equations, the adjoint method has found several successful applications, e.g., in meteorology with large-scale 3D or 4D atmospheric data assimilation schemes, for carbon cycle data assimilation in biogeochemistry and climate research, or in oceanographic modelling with efficient adjoint codes of ocean general circulation models. Despite the variety of applications in these research fields, adjoint methods have only very recently drawn attention from the ocean acoustics community. In ocean acoustic tomography and geoacoustic inversion, where the inverse problem is to recover unknown acoustic properties of the water column and the seabed from acoustic transmission data, the solution approaches are typically based on travel time inversion or standard matched-field processing in combination with metaheuristics for global optimization. In order to complement the adjoint schemes already in use in meteorology and oceanography with an ocean acoustic component, this thesis is concerned with the development of the adjoint of a full-field acoustic propagation model for shallow water environments. In view of the increasing importance of global ocean observing systems such as the European Seas Observatory Network, the Arctic Ocean Observing System and Maritime Rapid Environmental Assessment (MREA) systems for defence and security applications, the adjoint of an ocean acoustic propagation model can become an integral part of a coupled oceanographic and acoustic data assimilation scheme in the future. Given the acoustic pressure field measured on a vertical hydrophone array and a modelled replica field that is calculated for a specific parametrization of the environment, the developed adjoint model backpropagates the mismatch (residual) between the measured and predicted field from the receiver array towards the source. The backpropagated error field is then converted into an estimate of the exact gradient of the objective function with respect to any of the relevant physical parameters of the environment including the sound speed structure in the water column and densities, compressional/shear sound speeds, and attenuations of the sediment layers and the sub-bottom halfspace. The resulting environmental gradients can be used in combination with gradient descent methods such as conjugate gradient, or Newton-type optimization methods tolocate the error surface minimum via a series of iterations. This is particularly attractive for monitoring slowly varying environments, where the gradient information can be used to track the environmental parameters continuously over time and space. In shallow water environments, where an accurate treatment of the acoustic interaction with the bottom is of outmost importance for a correct prediction of the sound field, and field data are often recorded on non-fully populated arrays, there is an inherent need for observation over a broad range of frequencies. For this purpose, the adjoint-based approach is generalized for a joint optimization across multiple frequencies and special attention is devoted to regularization methods that incorporate additional information about the desired solution in order to stabilize the optimization process. Starting with an analytical formulation of the multiple-frequency adjoint approach for parabolic-type approximations, the adjoint method is progressively tailored in the course of the thesis towards a realistic wide-angle parabolic equation propagation model and the treatment of fully nonlocal impedance boundary conditions. A semi-automatic adjoint generation via modular graph approach enables the direct inversion of both the geoacoustic parameters embedded in the discrete nonlocal boundary condition and the acoustic properties of the water column. Several case studies based on environmental data obtained in Mediterranean shallow waters are used in the thesis to assess the capabilities of adjoint-based acoustic inversion for different experimental configurations, particularly taking into account sparse array geometries and partial depth coverage of the water column. The numerical implementation of the approach is found to be robust, provided that the initial guesses are not too far from the desired solution, and accurate, and converges in a small number of iterations. During the multi-frequency optimization process, the evolution of the control parameters displays a parameter hierarchy which clearly relates to the relative sensitivity of the acoustic pressure field to the physical parameters. The actual validation of the adjoint-generated environmental gradients for acoustic monitoring of a shallow water environment is based on acoustic and oceanographic data from the Yellow Shark '94 and the MREA '07 sea trials, conducted in the Tyrrhenian Sea, south of the island of Elba. Starting from an initial guess of the environmental control parameters, either obtained through acoustic inversion with global search or supported by archival in-situ data, the adjoint method provides an efficient means to adjust local changes with a couple of iterations and monitor the environmental properties over a series of inversions. In this thesis the adjoint-based approach is used, e.g., to fine-tune up to eight bottom geoacoustic parameters of a shallow-water environment and to track the time-varying sound speed profile in the water column. In the same way the approach can be extended to track the spatial water column and bottom structure using a mobile network of sparse arrays. Work is currently being focused on the inclusion of the adjoint approach into hybrid optimization schemes or ensemble predictions, as an essential building block in a combined ocean acoustic data assimilation framework and the subsequent validation of the acoustic monitoring capabilities with long-term experimental data in shallow water environments.

Water-induced soil erosion and shallow landslides interact with each other and need to be studied in an integrated approach to understand hillslope sediment yields. The principal aim of this thesis was to study and model soil erosion and shallow landslides in an integrated way. The thesis presents results from laboratory and catchment-scale studies and modelling. A laboratory flume under a rainfall simulator was used for shallow landslide and soil erosion experiments using sandy and silty loess soils. In the experiments, landslide initiation, retrogressions and slip surface depths were measured and monitored directly or by using video camera recordings. Sediment and runoff were collected from the flume outlet every minute during landslides and every 10 minutes before and after landslides. Changes in the soil slope, after landslides, were recorded. Initially, six experiments including two repetitions were conducted using sandy soils at a 30º and 10º compound slope configuration, but with different soil profile depths. The experimental results showed that total and landslide-driven sediment yields were affected by the original soil profile depth; the greater the depth, the higher the sediment yield. Later, twelve other experiments were conducted on different slopes using silty loess soils. The experimental observations were used to validate an integrated modelling approach which includes WEPP for runoff and soil erosion modelling, a slope stability model for simulating shallow landslides, and a simple soil redistribution model for runout distance prediction. The model predictions were in good alignment with the observations. In all (sandy and silty loess) experiments, peak sediment discharges were related to the landslide events, proximity to the outlet and landslide volume. The post-failure soil erosion rate decreased as a function of changes in the slope profile. The GeoWEPP-SLIP modelling approach was proposed for catchment-scale modelling. The approach simulates soil erosion using the Hillslope and Flowpath methods in WEPP, predicts shallow landslides using a slope stability model coupled with the WEPP’s hillslope hydrology and finally uses a simple rule-based soil redistribution model to predict runout distance and post-failure topography. A case study application of the model to the Bowenvale research catchment (300 ha) showed that the model predictions were in good agreement with the observed values. However, the Hillslope method over-predicted the outlet sediment yield due to the computational weighting involved in the method. The Hillslope method predicted consistent values of sediment yield and soil erosion regardless to the changes in topography and land-cover in the post-failure scenarios. The Flowpath method, on the other hand, predicted higher values of sediment yield in the post-failure vegetation removal scenario. The effects of DEM resolution on the approach were evaluated using four different resolutions. Statistical analyses for all methods and resolutions were performed by comparing the predicted versus measured runoff and sediment yield from the catchment outlet and the spatial distribution of shallow landslides. Results showed that changes in resolution did not significantly alter the sediment yield and runoff between the pre- and post-failure scenarios at the catchment outlet using the Hillslope method. However, the Flowpath method predicted higher hillslope sediment yields at a coarser resolution level. Similarly, larger landslide areas and volumes were predicted for coarser resolutions whereas deposition volume decreased with the increase in grid-cell size due to changes in slope and flowpath distributions. The research conducted in the laboratory and catchment presented in this thesis helped understand the interactions between shallow landslides and soil erosion in an integrated approach.