We present 2-D numerical simulations of two phase flow in seafloor hydrothermal systems using the finite control volume numerical scheme FISHES. The FISHES code solves the coupled non-linear equations for mass, momentum, energy, and salt conservation in a NaCl-H2O fluid to model the seafloor hydrothermal processes. These simulations use homogeneous box geometries at a fixed seafloor pressure of 25 MPa with constant bottom temperature boundary conditions that represent a sub-axial magma chamber to explore the effects of permeability, maximum bottom temperature and system depth on the evolution of vent fluid temperature and salinity, and heat output. We also study the temporal and spatial variability in hydrothermal circulation. The two-phase simulation results show that permeability plays an important role in plume structure and heat output of hydrothermal systems, but it has little effect on vent fluid temperature and salinity, given the same bottom temperature. For some permeability values, multiple plumes can vent at the seafloor above the simulated magma chamber. Temporal variability of vent fluid temperature and salinity and the complexity of phase separation suggest that pressure and temperature conditions at the top of the axial magma chamber cannot be easily inferred from vent fluid temperature and salinity alone. Vapor and brine derived fluids can vent at the seafloor simultaneously, even from neighboring locations that are fed by the same plume. / Master of Science

A novel micromechanical model for predicting the failure of polymeric and ceramic matrices under an applied global tensile loading is presented. The model is based on a Weibull statistical type of approach and incorporates eleven different variables to simulate different micromechanical failure phenomenon. Through the variables, the model incorporates such phenomenon as fiber-fiber interaction; matrix-fiber interaction; damage due to processing; and local stress distribution at the interphase. The effects of the load sharing constants and shape parameter on failure are studied. A software package, used in the Windows environment, was also developed to perform a numerical analysis of the model. / Master of Science

[pt] A caracterização das zonas de dano no entorno das falhas geológicas é importante para a indústria de óleo e gás pois as estruturas geológicas presentes podem afetar tanto os processos de fluxo nos reservatórios quanto o comportamento geomecânico dos campos. Como aproximação inicial, a espessura das zonas de dano tem sido correlacionada diretamente com os rejeitos de falha. No entanto, a dispersão dos dados levantados em campo indica que outras variáveis, além do rejeito, também podem afetar à espessura da zona de dano, particularmente em rochas carbonáticas. Neste sentido, o objetivo desta dissertação é a avaliação do impacto das propriedades geomecânicas dos carbonatos na previsão da espessura da zona de dano em falhas geológicas dos campos do pré-sal Brasileiro. Em função das incertezas na definição das propriedades geomecânicas, estas foram inferidas a partir de correlações empíricas que foram estabelecidas com a porosidade. Em seguida, fazendo uso de modelagem numérica, se realizou um planejamento de experimentos, obtendo-se a espessura da zona de dano a partir de diferentes combinações de parâmetros geomecânicos. Logo, empregando o método da superfície de resposta, foram definidos modelos de regressão que preveem a espessura da zona de dano a partir do rejeito e da porosidade das rochas carbonáticas avaliadas. Os resultados indicam que os modelos de regressão propostos permitem estimar espessuras de zonas de dano consistentes com as observações em campo, e, portanto, podem ser utilizadas para uma caracterização preliminar das zonas de dano nos processos exploratórios da indústria de óleo e gás. / [en] The characterization of damage zones surrounding geological faults is important for the oil and gas industry, as the geological structures present can affect both the fluid flow in reservoirs and the geomechanical behavior of fields. As an initial approximation, the width of damage zones has been directly correlated with fault throw. However, the dispersion of data collected in the field indicates that other variables, besides fault throw, can also affect the width of the damage zone, particularly in carbonate rocks. In this sense, the objective of this dissertation is to evaluate the impact of geomechanical properties of carbonates on the prediction of the thickness of the damage zone in geological faults in the Brazilian pre-salt fields. Due to uncertainties in defining geomechanical properties, they were inferred from empirical correlations that were established with porosity. Then, using numerical modeling, a design of experiments was carried out, obtaining the width of the damage zone from different combinations of geomechanical parameters. Next, employing the response surface method, regression models were defined that predict the width of the damage zone from fault throw and porosity of the evaluated carbonate rocks. The results indicate that the proposed regression models allow for the estimation of damage zone widths consistent with field observations and therefore can be used for preliminary characterization of damage zones in exploration processes in the oil and gas industry.

Pile-supported T-walls, which are concrete floodwalls that are shaped like an inverted "T" and supported by batter piles, are commonly used by the United States Army Corps of Engineers (USACE) to protect low-lying portions of New Orleans and other areas. The design of a T-wall in southern Louisiana is complex, as the structure needs to resist both 1) large settlements caused by fill placed beneath and beside the T-wall before T-wall construction or by fill placed beside the T-wall after T-wall construction, and 2) large lateral flood loads that are imposed during a hurricane. As a result of these loading conditions, large bending moments can develop in the batter piles and these moments need to be accounted for as part of the T wall design. The goal of this research is to develop a more complete understanding of the pile bending moments in T wall systems, specifically for cross sections where large settlements may occur. As a first step towards this goal, Rensselaer Polytechnic Institute (RPI) performed a series of eight centrifuge tests to investigate and physically model the effects of settlement-induced bending moments on pile-supported T-walls. The centrifuge tests were evaluated and interpreted, and in order to better capture uncertainty, upper and lower bounds were estimated for the interpreted data. The centrifuge results offered some valuable insights on their own, but were ultimately used as the basis for validating and calibrating corresponding numerical models. The numerical models were developed following a rigorous modeling approach and using rational and reasonable assumptions based on widely-accepted and well-justified procedures. The numerical model results were in good agreement with the centrifuge results without the need for significant calibration or modifications. This good agreement indicates that similar numerical models can be developed to reliably analyze actual T-wall cross sections. Detailed recommendations were developed for using numerical models to analyze pile-supported T walls, and an example problem is presented herein that illustrates the application of this approach. These same techniques were then used to perform a parametric study to analyze the combined and separate effects of flood loading for a wide range of different T-wall cross sections. The range was selected in collaboration with the USACE in order to reasonably cover cross sections and conditions that 1) are typically encountered in practice, and 2) were expected to generate both upper and lower bound pile bending moments. In total, 3,648 cross sections were analyzed, and 29,184 sets of analysis results were generated since each cross section was analyzed for eight different loading conditions. Summary results are provided to show the influence of the loading conditions and parameters on T-wall response, including the influence of flood loading, new fill symmetry, pile fixity, number of piles, subsurface profile, pile batter, pile type, levee slope, T-wall elevation, and the presence of existing levee fill. In addition, the key results for all of the analyses are provided in the appendices and in an electronic database. Based on the parametric study results, a simplified analysis procedure was developed that can be used to calculated maximum pile bending moments for T walls installed directly on foundation soils due to settlements. In this procedure, the loads from new fill placed during or after T-wall construction are distributed onto the pile, and the pile response is analyzed using traditional p-y curves and a beam on elastic foundation formulation. This procedure shows good agreement with the numerical model results for a range of conditions. To demonstrate the application of the procedure, the same example problem that is analyzed numerically is reanalyzed using the simplified analysis procedure. Due to the complexity of the problem, it was not possible to modify this procedure or develop a similar procedure for T-walls installed on top of new or existing levees. Overall, this research demonstrates that numerical models can be used to calculate the bending moments that can develop in pile-supported T-walls due to settlements and flood loading, provides valuable insights into the behavior of T-walls and the influence of various parameters on T-wall response, presents a large database of T-wall analysis results, and recommends a simplified analysis procedure that can be used in some cases to calculate pile bending moments due to settlements. / Ph. D.

Coronary stents have become a very popular treatment for cardiovascular disease, historically the leading cause of death in the United States. Stents, while successful in the short term, are subject to high failure rates (up to 24% in the first six months) due to wall regrowth and clotting, probably due to a combination of abnormal mechanical stresses and disruption of the arterial blood flow. The goal of this research was to develop recommendations concerning ways in which stent design might be improved, focusing on the problem of pressure wave reflections. A one-dimensional finite-difference model was developed to predict these reflections, and effects of variations in stent and vessel properties were examined, including stent stiffness, length, and compliance transition region, as well as vessel radius and wall thickness. The model was solved using a combination of Weighted Essentially Non-Oscillatory (WENO) and Runge-Kutta methods. Over 100 cases were tested. Results showed that reasonable variations in these parameters could induce changes in reflection magnitude of up to ±50%. It was also discovered that the relationship between each of these properties and the resulting wave reflection could be described simply, and the effect of all of them together could in fact be encompassed by a single non-dimensional parameter. This parameter was titled"Stent Authority," and several variations were proposed. It is believed this parameter is a novel way of relating the energy imposed upon the arterial wall by the stent, to the fraction of the incident pressure energy which is reflected from the stented region. / Master of Science

A water-filled geomembrane tube is considered for the purpose of temporary flood protection. With proper design, this tube can be a cheap and efficient breakwater, temporary levee, or cofferdam. This thesis considers a single tube resting on clay and sand foundations. A finite difference program, FLAC, is used in the numerical analyses. The tube is assumed to be infinitely long, and it is modeled two-dimensionally. Beam elements are used to model the tube. The tube is inflated with water. The hydrostatic pressure in the tube is converted to point loads and applied at the beam nodes in the direction perpendicular to the chord connecting two adjacent nodes. Two of FLAC's built-in soil models are used: elastic and Mohr-Coulomb. The Mohr-Coulomb model is used in all the cases except the preliminary analyses, in which the elastic soil model is used. The Mohr-Coulomb soil model is able to model the soil's nonlinear stress-strain and path-dependent deformation behavior. A tube without external water is placed on clay with various shear strengths to study how the clay consistency affects the height and the stresses in the tube. A tube with external water on one side is placed on medium dense sand. A wooden block is placed on the side opposite the floodwater. Three types of block geometry and two sizes are studied. The floodwater level is increased until the system fails. Three failure modes, rolling, sliding, and piping, are studied. The effect of pore pressure on these failure modes is examined. The influence of a filter placed under part of the tube and block is also investigated. The tube's tensile forces, shear forces, moments, and settlements are included. Soil stresses and pore pressures at the soil-tube interfaces are computed. The cross-section of the tube at various external water levels, and the pore pressures in the soil, are calculated. These results are compared with experimental results that were obtained by graduate students in geotechnical engineering at Virginia Tech. / Master of Science

Full Waveform Inversion (FWI) is a powerful seismic imaging technique used to reconstruct high-resolution velocity models of the subsurface. It relies on the inversion of seismic data acquired from multiple sources and receivers to estimate the mechanical properties of geologic materials and can be used to detect anomalous subsurface conditions. The accuracy of FWI results is influenced by various factors related to the workflow used for its implementation. This includes the survey parameters, the mathematical framework of the inversion, and the complexity of the subsurface conditions modeled during the inversion process. Therefore, it is crucial to have a fundamental understanding of the interplay between these factors and their impact on the accuracy of the reconstructed model, particularly given the effects of these factors on computational costs. This is an area that has been understudied within the context of near-surface geotechnical applications for anomaly detection, which is an application that presents unique challenges relative to seismic exploration for hydrocarbons where FWI has been more fully developed. One key aspect that has not received sufficient attention is the impact of survey parameters on the accuracy of FWI results. The lack of formal research in this topic may lead to near-surface FWI studies that use more seismic sources than required for subsurface feature reconstruction, which results in data collection and computational inefficiencies. The selection of misfit function and starting model are also essential factors influencing the reliability of the reconstructed model. The physics employed for forward modeling can also affect the ability to simulate wave propagation in the domain of interest. These factors have significant implications for near-surface applications of FWI, and further research is required to explore their interplay and improve FWI workflow.Given the gaps in the current implementation of FWI for geotechnical applications, this research will explore the role of parameterization and workflow on FWI results when applied to anomaly detection in karst conditions. This will include selection of an FWI workflow that can improve the feasibility of fieldwork and reduce the processing time. The research will investigate four key factors of the FWI workflow (i.e., survey design, initial model, misfit function, and forward modeling physics) for detection of sinkholes using numerical and field testing in different subsurface conditions. Overall, the outcomes of this research will help practitioners with more appropriate choices in the FWI process and consequently promote its high potential in near-surface applications. / Civil Engineering

[pt] A água de um aquífero não é totalmente pura, ela contém substâncias constituintes dissolvidas como sulfatos e cloretos, líquidos e constituintes orgânicos, inorgânicos dissolvidos e patógenos. É importante a remediação da água contaminada com o controle ou remoção da fonte de contaminação, evitando a contínua propagação do poluente no aquífero. A captura da pluma de contaminação pode ser realizada com poços de bombeamento, formando zonas de captura, impedindo sua propagação para áreas limpas e removendo o poluente. Esta extração deve ser com o menor volume de água possível pois esta deverá ser descontaminada e utilizada ou descartada. Tendo esta direção de pesquisa, este trabalho baseia-se no desenvolvimento de um software, que possa simular o fluxo de água em um aquífero confinado e a formação de zonas de captura para a extração de poluentes por poços; na investigação de diferentes operações de bombeamento minimizando as zonas de estagnação. O software desenvolvido utiliza elementos finitos 2d; variação no tempo das condições de fluxo monofásico; obtém a solução no tempo pelo método da integração temporal. Os estudos numéricos consideraram o aquífero homogêneo ou heterogêneo e isotrópico, e o fluxo de água segue a lei de Darcy. Foram analisadas zonas de captura sob diferentes condições hidrogeológicas; com poços em diferentes posições, com diferentes vazões e com bombeamento constante e; a utilização de poços de injeção reinjetando água extraída. Foi constatado que a modelagem numérica é importante na avaliação de estratégias de bombeamento. / [en] The aquifer water is not totally pure, it contains dissolved constituents from natural solids sources such as sulfates and chlorides. Also, there are organic liquids, organic and inorganic constituents dissolved or pathogens from an anthropogenic source. It is important the contaminated water remediation with the control or removal of the source of contamination, avoiding the continuous propagation of the contamination plume in the aquifer and the reduction of contaminant concentration levels. Control and capture of the contamination plume can be done by pumping wells positioned at the edge of the plume, forming the capture zones to prevent the plume propagation to the not contaminated areas and to remove the pollutant from the aquifer. This extraction should be with a few water volumes, as this should be decontaminated and used or discarded. With this research direction, this work is based on the development of a software to simulate the water flow in a confined aquifer and the formation of the capture zone by pumping wells. Also, it is based on the investigation of different pumping operation features to maximize the contaminant extraction with the minimizing of the withdrawn water volume. The developed software uses bi and three-dimensional finite elements; time variation of monophasic flow conditions; obtains the solution in time by the method of direct temporal integration. The numerical studies considered homogeneous or heterogeneous and isotropic aquifer, and the groundwater flow follows Darcy s law. The formation of the capture zone under different hydrogeological conditions was analyzed; also, with pupping wells in different positions, with different flows conditions and with constant pumping and; the use of injection wells reinjecting extracted water. It was verified that the numerical modeling is important to the evaluation of pumping strategies that provide the greatest removal of solute mass with the lowest water volume extraction.

[EN] In the first part of the current thesis, two fundamental numerical models (Fin2D-W and Fin1D-MB) for analyzing the air-side performance of minichannel evaporators were developed and verified. The Fin2D-W model applies a comprehensive two-dimensional scheme to discretize the evaporator. On the other hand, the Fin1D-MB model is based on the one-dimensional fin theory in conjunction with the moving boundaries technique along the fin height. The first objective of the two presented models is to identify and quantify the most influential phenomena encountered in the process of cooling and dehumidification. The second objective is to study the impact of the classical modeling assumptions on the air-side performance of minichannel evaporators. Different comparative studies between the traditional Effectiveness-NTU approach and the proposed numerical models were implemented to achieve the mentioned goals. The results revealed that the modeling assumptions which have the most significant impacts on the heat and mass transfer rates are: the uniform air properties along the fin height, adiabatic-fin-tip at half the height, and negligence of partial dehumidification scenarios. These widely used assumptions resulted in substantial deviations in total heat transfer rate, up to 52%, between the Effectiveness-NTU approach and Fin2D-W model. In the second part of the thesis, the Fin1D-MB model was integrated into the IMST-ART® simulation tool to evaluate the global performance of minichannel evaporators (air- and refrigerant-side). The Fin1D-MB model was selected because of its simplicity, calculation speed, and reasonable solution accuracy relative to the Fin2D-W model. The validation of the complete Fin1D-MB model was conducted against many experimental data and numerical models available in the literature. The validation process was achieved for different heat exchanger geometries, refrigerants, and operating conditions. The results showed that for the R134a minichannel evaporators studied, the Fin1D-MB model successfully predicted the Inlet refrigerant and outlet air temperatures, cooling capacity, and refrigerant-side pressure drop within error bands of ±0.5 ºC, ±5%, and ±20%, respectively. For the CO2 (R744) minichannel evaporator studied, the presented model estimated the cooling capacity and outlet air temperature within error bands of ±10% and ±1.0 ºC, respectively. Regarding the CO2 pressure drop, the Fin1D-MB model generally underpredicted the pressure drop values compared to the experimental data, with a maximum deviation of 11 kPa. / [ES] En la primera parte de la tesis actual, dos modelos numéricos fundamentales (Fin2D-W y Fin1D-MB) para analizar el lado del aire de los evaporadores de minicanales se han desarrollado y verificado. El modelo Fin2D-W aplica un esquema detallado de dos dimensiones para discretizar el evaporador mientras que el modelo Fin1D-MB se basa en la teoría de la aleta unidimensional junto con la técnica de fronteras móviles para el lado del aire. El primer objetivo de los dos modelos presentados es identificar y cuantificar los fenómenos más influyentes encontrados en el proceso de enfriamiento y deshumidificación. El segundo objetivo es estudiar el impacto de las hipótesis comúnmente usadas en el modelado de la transmisión de calor del aire de los evaporadores de minicanales. Se implementaron diferentes estudios comparativos entre el enfoque tradicional Effectiveness-NTU y los modelos numéricos propuestos para alcanzar los objetivos mencionados. Los resultados muestran que las hipótesis que provocan una mayor desviación con respecto a la solución detallada en la transferencia de calor y masa son: propiedades de aire uniforme a lo largo de la altura de la aleta, extremo adiabático de aleta a mitad de su longitud, y no contemplar el supuesto de deshumidificación parcial en la aleta. Estas hipótesis ampliamente utilizadas han resultado en errores importantes en la transferencia de calor total, hasta un 52%, entre el enfoque Effectiveness-NTU y el modelo Fin2D-W. En la segunda parte de la tesis, el modelo Fin1D-MB se integró en la herramienta de simulación IMST-ART® para evaluar el rendimiento global de los evaporadores de minicanales (en el lado del aire y del refrigerante). El modelo Fin1D-MB se seleccionó gracias a su simplicidad, velocidad de cálculo, y solución de una precisión razonable relativa al modelo Fin2D-W. Se realizó una validación del modelo completo Fin1D-MB con la ayuda de datos experimentales y modelos numéricos ya disponibles en la literatura. El modelo se ha validado para diferentes geometrías de intercambiadores de calor, refrigerantes y condiciones de funcionamiento. Los resultados han mostrado que para los evaporadores de minicanales funcionando con el refrigerante R134a, el modelo Fin1D-MB predice de manera correcta las temperaturas de entrada del refrigerante y de salida del aire, la capacidad de enfriamiento, y la caída de presión del lado de refrigerante dentro de las bandas de error de ±0.5 ºC, ±5%, y ±20%, respectivamente. Para el evaporador de minicanales con CO2 (R744) estudiado, el modelo estima la capacidad de refrigeración y la temperatura de salida del aire dentro de las bandas de error de ±10% y ±1.0 ºC, respectivamente. En cuanto a la caída de presión de CO2, el modelo Fin1D-MB generalmente predice a la baja los valores de la caída de presión en comparación con los datos experimentales, con una desviación máxima de 11 kPa. / [CA] A la primera part de la tesi actual, dos models numèrics fonamentals (Fin2D-W i Fin1D-MB) per analitzar el costat de l'aire dels evaporadors de minicanals s'han desenvolupat i verificat. Al model Fin2D-W s'aplica un esquema detallat de dues dimensions per discretitzar l'evaporador mentre que al model Fin1D-MB es basa en la teoria d'aleta unidimensional juntament amb la tècnica de frontera mòbil per al costat de l'aire. El primer objectiu dels dos models presentats és identificar i quantificar els fenòmens més influents trobats en el procés de refredament i deshumidificació. El segon objectiu és estudiar l'impacte de les hipòtesis comunament utilitzades en el modelatge de la transmissió de calor de l'aire dels evaporadors de minicanals. Es van implementar diferents estudis comparatius entre l'enfocament tradicional Effectiveness-NTU i els models numèrics proposats per assolir els objectius esmentats. Els resultats mostren que les hipòtesis que provoquen una major desviació respecte a la solució detallada a la transferència de calor i massa són: propietats d'aire uniforme al llarg de l'altura de l'aleta, extrem adiabàtic d'aleta a la meitat de la seua longitud, i no contemplar el supòsit de deshumidificació parcial en l'aleta. Aquestes hipòtesis àmpliament utilitzades donen errors importants en la transferència de calor total, fins a un 52%, entre l'enfocament Effectiveness-NTU i el model Fin2D-W. A la segona part de la tesi, el model Fin1D-MB es va integrar en l'eina de simulació IMST-ART® per avaluar el rendiment global dels evaporadors de minicanals (al costat de l'aire i del refrigerant). El model Fin1D-MB es va seleccionar gràcies a la seva simplicitat, velocitat de càlcul, i solució d'una precisió raonable relativa al model Fin2D-W. Es va realitzar una validació del model complet Fin1D-MB amb l'ajuda de dades experimentals i models numèrics ja disponibles a la literatura. El model s'ha validat per a diferents geometries d'intercanviadors de calor, refrigerants i condicions de funcionament. Els resultats mostren que per als evaporadors de minicanals funcionant amb el refrigerant R134a, el model Fin1D-MB prediu de manera correcta les temperatures d'entrada del refrigerant i de sortida de l'aire, la capacitat de refreda-ment, i la caiguda de pressió del costat de refrigerant dins de les bandes d'error de ±0.5 ºC, ±5%, i ±20%, respectivament. Per a l'evaporador de minicanals amb CO2 (R744) estudiat, el model estima la capacitat de refrigeració i la temperatura de sortida de l'aire dins de les bandes d'error de ±10% i ±1.0 ºC, respectivament. Pel que fa a la caiguda de pressió de CO2, el model Fin1D-MB generalment prediu a la baixa els valors de la caiguda de pressió en comparació amb les dades experimentals, amb una desviació màxima d'11 kPa. / Hassan, AHA. (2016). Development and Validation of a Minichannel Evaporator Model under Dehumidification [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/71357

Los modelos numéricos son herramientas fundamentales para realizar predicciones de muchos problemas enfrentados por ingenieros geotécnicos y geoambientales. Sin embargo, para que estos modelos puedan realizar predicciones confiables, los parámetros de entrada del modelo deben ser estimados considerando el efecto escala. En este contexto, esta tesis se concentra en las reglas del cambio de escala de los parámetros de flujo y transporte de masa en un suelo tropical a través de estudios numéricos, de laboratorio y de campo. Esta está organizada en cuatro partes. Primero, la heterogeneidad, correlación y correlación cruzada entre los parámetros de transporte de solutos (dispersividad, ¿, y coeficiente de partición, Kd) y las propiedades del suelo fueron estudiadas en detalle. En esta parte fue verificado que la conductividad hidráulica (K) y los parámetros de transporte de solutos son altamente heterogéneos, mientras que las propiedades del suelo no lo son. La correlación espacial de ¿ y K con variables estadísticamente significativas fue estudiada. Este resultado probablemente podrá mejorar la estimación en casos de estudios de pequeña escala debido a que solo fue observada correlaciones de hasta 2,5 m. Este estudio fue un primer intento de evaluar la variación espacial en el coeficiente de correlación de los parámetros de transporte de un soluto reactivo y de un no reactivo, indicando las variables más relevantes y aquella que debería ser incluida en estudios futuros. En la segunda parte, el efecto escala en K, dispersividad y coeficiente de partición de potasio y clorito fue estudiado experimentalmente a través de experimentos de laboratorio y de campo. El objetivo de esta parte fue contribuir a la discusión sobre el efecto escala en K, ¿ y Kd, y entender como estos parámetros se comportan con el cambio de escala de medición. La dispersividad tiende a aumentar con la altura de la muestra de manera exponencial. El coeficiente de partición tiende a aumentar con la altura, el diámetro y el volumen de la muestra. Estas diferencias encontradas en los parámetros de acuerdo con la escala de medición deben ser considerados cuando estos valores sean usados posteriormente como datos de entrada de modelos numéricos; de otra manera, las respuestas pueden ser malinterpretadas. Tercero, análisis estocásticos tridimensionales de cambio de escala de la conductividad hidráulica fueron realizados usando los métodos de promedios simples y de Laplace con piel para una variedad de tamaños de bloques usando mediciones reales de K. En esta parte son demostrados los errores que pueden ser introducidos al usar métodos determinísticos de cambio de escala usando promedios simples de las mediciones de K sin llevar en consideración la correlación espacial. La aplicación muestra que la heterogeneidad de K puede ser incorporada en la práctica diaria del modelador geotécnico. Los aspectos que considerar durante un proceso de cambio de escala también son discutidos. Finalmente, la dependencia del exponente de la norma-p como función del tamaño del bloque fue analizada. En la última parte, una aplicación de cambio de escala estocástico del coeficiente de dispersión hidrodinámica D y del factor de retardo R fue realizada usando datos reales con el objetivo de reducir la falta de casos de investigación experimental de cambio de escala de parámetros de transporte de solutos reactivos. El cambio de escala de D fue realizado usando el método de macrodispersión. El método de promedio simple de norma-p fue usado para realizar el cambio de escala de R. Una buena propagación de incertidumbres fue alcanzada. Métodos simples de cambio de escala pueden ser introducidos en la práctica del modelaje usando programas comerciales de transporte y conseguir reproducir el transporte en escala gruesa, pero puede requerir correcciones con el objetivo de reducir el efecto de suavizado de la heterogeneidad causado por el / Numerical models are becoming fundamental tools to predict a range of complex problems faced by geotechnical and geo-environmental engineers. However, to render the model reliable for future predictions, the model input parameters must be determined with consideration of the scale effects. In this context, this thesis focuses on upscaling of water flow and mass transport in a tropical soil by means of numerical, laboratory and field studies. This thesis is organized in four parts. First, the heterogeneity, correlation and cross-correlation between solute transport parameters (dispersivity, ¿, and partition coefficient, Kd) and soil properties were studied in detail. In this part, it was verified that the hydraulic conductivity (K) and solute transport parameters are highly heterogeneous, while soil properties not. Spatial correlation of ¿, K, and statistically significant variables were studied, and it would probably improve the estimation only in a small-scale study, since the spatial correlation were only observed up to 2.5 m. This study was a first attempt to evaluate the spatial variation in the correlation coefficient of transport parameters of a reactive and a nonreactive solute, indicating the more relevant variables and the one that should be included in future studies. In the second part, scale effect on K, dispersivity and partition coefficient of potassium and chloride is studied experimentally by means of laboratory and field experiments. The purpose of was to contribute to the discussion about scale effects on K, ¿ and Kd and understanding how these parameters behave with the change in the scale of measurement. Results shows that K increases with scale, regardless of the method of measurement. Dispersivity trends to increases exponentially with the sample height. Partition coefficient, tend to increase with sample length, diameter and volume. These differences in the parameters according to the scale of measurement must be considered when these observations are later used as input to numerical models, otherwise the responses can be misrepresented. Third, stochastic analysis of three-dimensional hydraulic conductivity upscaling was performed using a simple average and the Laplacian-with-skin methods for a variety of block sizes using real K measurements. In this part it was demonstrated the errors that can be introduced by using a deterministic upscaling using simple averages of the measured K without accounting for the spatial correlation. The application shows that K heterogeneity can be incorporated in the daily practice of the geotechnical modeler. The aspects to consider when performing the upscaling were also discussed. Finally, the dependence of the exponent of the p-norm as a function of the block size was analyzed. In the last part, an application of stochastic upscaling of hydrodynamic dispersion coefficient (D) and retardation factor (R) was performed using real data aiming to reduce the lack in experimental upscaling of reactive solute transport research. Upscaling of D was done using macrodispersion method. Simple average method based on p-norm was used to perform R upscaling. A good propagation of the uncertainties was achieved. Simple upscaling methods can be incorporated to the modeling practice using commercial transport codes and properly reproduce de transport at coarse scale but may require corrections to reduce smoothing of the heterogeneity caused by the upscaling procedure. / Els models numèrics s'estan constituint en eines fonamentals per a realitzar prediccions d'una àmplia gamma de problemes enfrontats per enginyers geotècnics i geoambientales. No obstant açò, perquè aquests models puguen realitzar prediccions fiables, els paràmetres d'entrada del model han de considerar l'efecte escala. En aquest context, aquesta tesi es concentra en les regles del canvi d'escala dels paràmetres de flux i transport de massa en un sòl tropical a través d'estudis numèrics, de laboratori i de camp. Aquesta tesi està organitzada en quatre parts. Primer, l'heterogeneïtat, correlació i correlació creuada entre els paràmetres de transport de soluts (dispersivitat, ¿, i coeficient de partició, Kd) i les propietats del sòl van ser estudiades detalladament. En aquesta part va ser verificat que la conductivitat hidràulica (K) i els paràmetres de transport de soluts són altament heterogenis, mentre que les propietats del sòl no ho són. La correlació espacial de ¿ i K amb variables estadísticament significatives va ser estudiada. Aquest resultat probablement podrà millorar l'estimació en casos d'estudis de xicoteta escala a causa que solament va ser observada correlacions de fins a 2,5 m. Aquest estudi va ser un primer intent d'avaluar la variació espacial en el coeficient de correlació dels paràmetres de transport d'un solut reactiu i d'un no reactiu, indicant les variables més rellevants i aquelles que haurien de ser inclosas en estudis futurs. En la segona part, l'efecte escala en K, dispersivitat i coeficient de partició de potassi i clorito va ser estudiat experimentalment a través d'experiments de laboratori i de camp. L'objectiu d'aquesta part va ser contribuir a la discussió sobre l'efecte escala en K, ¿ i Kd, i entendre com aquests paràmetres es comporten amb el canvi d'escala de mesurament. La dispersivitat tendeix a augmentar amb l'altura de la mostra, és a dir, amb la longitud del transport, de manera exponencial. El coeficient de partició tendeix a augmentar amb l'altura, el diàmetre i el volum de la mostra. Aquestes diferències en els paràmetres d'acord amb l'escala de mesurament han de ser considerats quan aquests valors siguen usats posteriorment com a dades d'entrada de models numèrics; d'una altra manera, les respostes poden ser malament interpretades. Tercer, anàlisis estocàstiques tridimensionals de canvi d'escala de la conductivitat hidràulica van ser realitzats usant els mètodes de mitjanes simples i de Laplace amb pell per a una varietat de grandàries de blocs usant mesuraments reals de K. En aquesta part són demostrats els errors que poden ser introduïts en usar mètodes determinístics de canvi d'escala usant mitjanes simples dels mesuraments de K sense tindre en consideració la correlació espacial. L'aplicació mostra que l'heterogeneïtat de K pot ser incorporada en la pràctica diària del modelador geotècnic. Els aspectes a considerar durant un procés de canvi d'escala també són discutits. Finalment, la dependència de l'exponent de la norma-p com a funció de la grandària del bloc va ser analitzada. En l'última part, una aplicació de canvi d'escala estocàstic del coeficient de dispersió hidrodinámica D i del factor de retard R va ser realitzada usant dades reals amb l'objectiu de reduir la falta de casos de recerca experimental de canvi d'escala de paràmetres de transport de soluts reactius. El canvi d'escala de D va ser realitzat usant el mètode de macrodispersió. El mètode de mitjana simple de norma-p va ser usat per a realitzar el canvi d'escala de R. Una bona propagació d'incerteses va ser aconseguida. Mètodes simples de canvi d'escala poden ser introduïts en la pràctica de la modelació usant programes comercials de transport i aconseguir reproduir el transport en escala gruixuda, però pot requerir correccions amb l'objectiu de reduir l'efecte de suavitzat de l'heterogeneïtat causat pel procés de canvi d'escala. / Almeida De Godoy, V. (2018). Upscaling of water flow and mass transport in a tropical soil: numerical, laboratory and field studies [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/102405