[pt] A natureza física do amortecimento nos solos e sua formulação matemática têm sido objeto de intensa pesquisa ao longo do último meio século. Estes vários ensaios para explicar os mecanismos de perda de energia durante os episódios de vibração surgem da importância do amortecimento no problema de propagação de ondas para fins de engenharia. O objetivo deste estudo é propor uma maneira alternativa de analisar o fenômeno do amortecimento em solos, considerando o problema de vibração como um problema matemático poromecânico acoplado onde deslocamentos relativos e velocidades entre as duas fases que compõem o material geológico, fluido e sólido, gera forças de interação nas interfaces e pode desempenhar um papel importante na perda de energia vibracional dos solos. Presume-se que o problema seja resolvido na escala de poros, concentrando-se nas condições cinéticas e dinâmicas na interface entre o esqueleto sólido e o fluido de saturação. Variáveis adimensionais que misturam propriedades fluidas e sólidas permitem o acoplamento do problema, resultando no surgimento de uma massa aparente, coeficientes de amortecimento e rigidez que serão introduzidos posteriormente nas equações de movimento. A equação proposta de amortecimento será então justaposta a testes de laboratório e uma comparação com os modelos mais importantes da literatura será feita usando o software DEEPSOIL para análise de propagação de ondas para verificar sua validade. / [en] The physical nature of damping in soils and its mathematical formulation attempts have been the subject of intensive researches along the last half century. These several trials to explain the mechanisms of energy loss during vibration episodes arise from the importance of damping in the wave propagation problem for engineering purposes. The aim of this study is to propose an alternative way to analyze the damping phenomenon in soils by looking at the vibration problem as a coupled poromechanical mathematical problem where relative displacements and velocities in between the two phases that compose the geological material, fluid and solid, generates interaction forces at the interfaces and can play an important role in the vibrational energy loss of soils. The problem is assumed to be solved at the pore scale, focusing on the kinetic and dynamic conditions at the interface in between the solid skeleton and the saturating fluid. Dimensionless variables that mixes fluid and solid properties allows the coupling of the problem, resulting on the rising of an apparent mass, damping and stiffness coefficients that will be introduced later in the equations of motion. The proposed equation of damping will thus be juxtaposed to laboratory tests data and a comparison with the most important models in the literature will be done further using the software DEEPSOIL for wave propagation analysis to check its validity.

[pt] Os fluxos de detritos são perigosos riscos naturais, que afetam países com intensas precipitações e terrenos montanhosos. Tais eventos configuram alto perigo para a vida humana e danificação de infraestrutura, resultando em importantes perdas econômicas. O estudo de fluxos de detritos envolve um mecânismo complexo e suas técnicas de previsão são baseadas na calibração de modelos, que devem ser delimitados por tentativa e erro de eventos anteriores. Tais previsões são ferramentas valiosas para delimitar as potenciais áreas de risco e, dessa forma, projetar medidas de mitigação e convivência. O principal objetivo deste trabalho foi analisar o comportamento de quatro fluxos de detritos deflagrados por precipitações de alta intensidade em regiões tropicais utilizando modelagem numérica em 2D e 3D. Foram analisados os casos de Lajas e Llano de la Piedra na Costa Rica e os casos de Córrego D Antas e Hospital São Lucas no Rio de Janeiro. Os principais parâmetros utilizados, na avaliação do risco deste tipo de movimentos de massa, são: a distância percorrida, a área de impacto, a velocidade e profundidade do fluxo. Os casos foram calibrados utilizando a reologia de Voellmy. A definição dos parametros na calibração é vital, pois oferece a possibilidade de previsões de primeira ordem, feita sobre escorregamentos acontecendo em condições semelhantes. Os resultados da análise dinâmica mostram valores consistentes entre os valores observados e as modelagens numéricas em 2D e 3D para os principais parâmetros avaliados, corroborando o uso destas ferramentas para análises de risco e projeção de medidas de mitigação e convivência. / [en] Debris flows are dangerous natural hazards affecting countries with steep terrains and heavy rainfall. They are associated with serious risks to human lives and infrastructure leading to important economic losses and fatalities. Debris flows involve complex mechanics and prediction techniques that are based on the calibration of models that must be constrained by trial-and-error back-analysis of previous landslides. Such predictions are a valuable tool for outlining potential hazard areas and the development of mitigation strategies and design of protective structures. The main goal of this work was to analyze the behavior of four debris flows triggered by heavy rainfall in tropical regions with numerical modelling. The Voellmy rheology was used to calibrate the cases occurred in Lajas and Llano de la Piedra in Costa Rica, and Córrego D Antas and Hospital São Lucas in Rio de Janeiro. The main parameters used for landslide risk assessment are runout distance, potential impact area, flow velocity and flow depth. The definition of appropriate calibrating parameters is important because it provides the possibility of first order predictions to be made about the motion of future landslides happening under similar conditions. The results of the dynamic analysis showed that consistent values were obtained for the main parameters evaluated in the 2D and 3D runout models, verifying the usefulness of these tools for landslide risk assessment and the project of protection structures.

[pt] Compósitos cimentícios estão ganhando cada vez mais relevância na indústria da construção civil. No entanto, as diretrizes para o projeto do material compósito e dos seus elementos estruturais são ainda incipientes, pois mecanismos de ponte de transferência de forças providos pelas fibras ainda estão sob investigação. Este trabalho apresenta uma estratégia de modelagem de elementos finitos que leva em consideração a estrutura de nível mesoestrutural do material cimentício reforçado com fibras. Desta forma, quatro fases do material são consideradas no modelo numérico: agregados graúdos, argamassa, zona de transição interfacial (ZTI) e fibras. A argamassa e os agregados são modelados usando elementos contínuos triangulares com comportamento linear-elástico. As fibras são incluídas usando elementos de treliça unidimensionais acopladas a elementos bidimensionais contínuos. Uma técnica de fragmentação de malha é usada para introduzir elementos de interface nas arestas dos elementos de argamassa e na interface entre agregados e argamassa para representar a ZTI. O método Take-and-Place, proposto por Wriggers e Moftah (2006), foi adotado neste estudo para incluir agregados no modelo. Primeiro, os agregados são gerados seguindo uma curva de Fuller, que define um empacotamento entre os agregados perfeitos. Na segunda fase, os agregados são introduzidos no modelo garantindo a não sobreposição entre eles. Finalmente, as fibras são adicionadas. Para validar a metodologia proposta, testes experimentais foram simulados com sucesso em um framework de simulação numérica – GeMA. Por fim, o trabalho explora a influência do empacotamento fibra-agregado na resposta mecânica e nos padrões de fraturamento de compósitos cimentícios fibrosos. / [en] Fiber Reinforced Concrete (FRC) materials are gaining more relevance in the construction industry. However, the guidelines for the design of the composite material and of structural elements thereof are incipient and the stress bridging mechanisms are still under investigation. This work presents a finite element modelling strategy that takes into account the material meso-level structure. Four phases of the FRC material are considered in the model: coarse aggregates, mortar, interfacial transition zone (ITZ), and fibers. The mortar and aggregates are modelled using triangular linear elements with linear–elastic behavior. Fibers are included using one-dimensional truss elements which are coupled to the matrix through the technique proposed by Congro (2021). Zero-thickness interface elements are introduced at the interface between mortar elements, and at the interface between aggregates and mortar to represent the ITZ. The Take-and-Place method, obtained from Wriggers and Moftah (2006), was adopted in this study to include aggregates in the model. First, the aggregates are generated following a Fuller s curve that means a perfect aggregate package. In the second phase, the aggregates are placed in the model without overlapping. Finally, fibers were added. A mesh fragmentation technique is used to introduced zero-thickness interface elements at the interface between mortar elements, and at the interface between aggregates and mortar to represent the ITZ. To validate the proposed methodology, direct tensile test models were successfully reproduced in finite element analyses performed in an in-house framework – GeMA. Based on the obtained results, the authors could explore the influence of the fibers-aggregate packing in the mechanical response of the composite material.

[pt] Colunas de aço estaiadas e protendidas são conhecidas como excelente solução em escoramento de grandes estruturas, como colunas de coberturas de lonas tensionadas, etc. Este trabalho apresenta uma série de ensaios experimentais tridimensionais em escala real desenvolvidos para determinação do comportamento estrutural de colunas de aço estaiadas e protendidas. Foram estudadas as variações no nível de protensão e na rigidez dos estais, além de descrever como os ensaios experimentais das três colunas de aço estaiadas são realizados: compreendendo os materiais utilizados; um novo sistema de medição de força nos estais; os passos e dificuldades na montagem das três colunas e o sistema de protensão aplicado. Existem alguns parâmetros que influenciam diretamente na resistência dessas colunas estaiadas, como por exemplo, a altura da coluna, o diâmetro externo, entre outros. Devido ao comportamento complexo deste tipo de colunas e ao grande esforço computacional para simulação do comportamento estrutural, através de uma análise paramétrica, optou-se por utilizar um projeto de experimentos junto com redes neurais a fim de extrapolar e obter novos resultados para carga crítica do sistema estrutural sem a necessidade de análise por programas de elementos finitos. De forma a complementar a tese, realizou-se um estudo do comportamento do sistema estrutural sujeito a ações dinâmicas através do programa de elementos finitos ANSYS com o objetivo de determinar as frequências naturais associadas aos seus modos de vibração. Também foi estudada a aplicação de um carregamento súbito para determinação do fator de amplificação dinâmico da coluna de aço estaiada e protendida. / [en] Prestressed steel columns are known as an efficient structural solution for great variety of temporary or permanent supporting systems for large span spatial frames and tensile surface structures. This work presents of full-scale three-dimensional tests carriedout for the assessment of structural behaviour of prestressed stayed steel columns. It was studied the effect prestress force level, stiffness of column braces and stays. Test setup and a new force measuring system for the column stays is fully described. Prestressed stayed steel columns have their strength dependant of parameters like: length, hollow section diameter, brace length and stiffness and axial stiffness of stays. Due to the complex behaviour of such columns that demands great computational effort for numerical simulations required for a parametric analysis it was used an experiment design tool coupled with neural network techniques employed to generate new data for the prestressed column buckling load. A study of the dynamic behaviour of prestressed columns using the finite element package ANSYS was carried-out in order to determine the column natural frequencies and their associated vibration modes. It was also studied the application of sudden loads to determine the dynamic amplification factor of this type of prestressed stayed steel column.

[pt] Esta pesquisa de tese teve como objetivo explorar a análise de mídias sociais, especialmente as disponíveis em comunidades online de sites de redes sociais, a fim de realizar estudos sociais sobre questões de saúde. Com base em uma abordagem prática foi definido um processo para realizar esses estudos. Este processo contou com ferramentas computacionais adaptados para fornecer apoio em tarefas específicas, tais como recuperação de conteúdo, seleção e análise. Duas ferramentas que se destacam são apresentadas por causa de sua utilidade e a complexidade do processo em que a sua construção se baseou. Para o benefício da análise de comunidades online, o Mapa de Associação de Comunidades é um processo desenvolvido para apoiar especialistas em compreender os interesses dos usuários com base em suas associações dentro de suas comunidades. A outra ferramenta visa auxiliar analistas a selecionar discussões de fóruns online a serem analisados manualmente com técnicas de pesquisa qualitativa, por exemplo, análise de conteúdo e do discurso. Esta ferramenta, TorchSR, foi criada baseada em aprendizado de máquina não supervisionado, usando agrupamento hierárquico, para dar suporte na resolução do problema de seleção de conteúdo. Um estudo de caso exploratório mostra que esta ferramenta ajuda na resolução do problema. O processo proposto foi utilizado em dois estudos sobre questões relevantes de saúde (hepatite C e o abuso de drogas), que resultou em descobertas relevantes sobre saúde pública. Em conclusão, este trabalho apresenta a aplicação prática de ciência social computacional no campo da saúde, através do desenvolvimento de um processo e ferramentas utilizadas para apoiar os analistas e melhorar a sua aplicação. / [en] This research thesis is aiming to exploit valuable social media, especially those available in online communities of social network sites, in order to perform social studies about healthcare issues. Based on a practical approach, a process was defined to conduct such studies. This process relied on tailored computational tools to provide support for specific tasks such as contente retrieval, selection, and analysis. Two tools that stand out are presented because of their utility and the complexity of the process in which their development was based on. The first tool, for the benefit of online community analysis, is the Community Association Map, a process developed to support experts in understanding users’ interests based on their associations within their communities. Our second tool (TorchSR) aims to aid analysts in the selection of discussions from online forums to be manually analyzed by (qualitative) research techniques (e.g. content and discourse analysis). This task, which was defined as solving the content selection problem, was tackled with a tool based on unsupervised machine learning techniques, such as hierarchical clustering. An exploratory study case shows that TorchSR helps analysts in dealing with the problem. The proposed process was employed in two studies about relevant healthcare issues (i.e. hepatitis C and drug abuse) which resulted in interesting findings in the field of public health. In conclusion, this thesis presents a practical application of computational social science to the field of health, through development of a process and tools used to support analysts and improve its application.

[pt] Freqüentemente as análises de estabilidade de taludes são feitas por métodos determinísticos, nos quais é obtido um Fator de Segurança (FS). Estes métodos não quantificam as incertezas existentes nas variáveis de entrada (parâmetros de resistência) na análise. Tampouco mostram detalhes sobre qual variável afeta mais o resultado. Os métodos probabilísticos permitem superar estes problemas. A presente pesquisa tem como objetivo comparar dois métodos probabilísticos geralmente utilizados (Estimativas Pontuais e Segundo Momento). Três projetos utilizados comumente na geotecnia são analisados: barragem de rejeitos, talude de solo e aterro sobre argila mole.Um aspecto importante na análise probabilística é a quantificação adequada do desvio padrão. É sabido que as propriedades dos solos mostram uma correlação no espaço, pelo que o desvio padrão calculado por métodos clássicos da estatística é superestimado em comparação ao valor real no campo. El-Ramly (2001) desenvolveu uma metodologia para o cálculo da probabilidade de ruptura considerando a correlação espacial das propriedades do solo. Esta metodologia é de difícil aplicação prática. Serão, portanto, avaliadas na presente pesquisa as técnicas geralmente usadas de probabilidade e estabilidade de taludes, juntamente com um fator de correção proposto por Vanmarcke (1977a). Estas técnicas de probabilidade de estabilidade de taludes são factíveis de serem utilizadas junto com este fator de correção. Verificou-se que o Método de Segundo Momento é de mais fácil utilização, e portanto adequado para emprego em projetos geotécnicos. / [en] The analysis of slope stability is often determined by deterministic methods, in which a Factor of Security (FS) is obtained. These methods do not quantify the uncertainty in the input variables (strength parameters). Neither these methods show details which variable affects mostly the outcome results. The probabilistic methods allow overcoming these restrictions. This study aims at comparing two probabilistic methods of general use (Point Estimates and First Order Second Moment). Three projects commonly used in geotechnical engineering are analyzed: tailings dam, slope soil and embankment on soft clay. An important aspect of the probabilistic analysis is the proper quantification of the standard deviation. It is known that the soil properties show a correlation in space, so the standard deviation, calculated by traditional methods of statistics, is overestimated when compared to the real field value. El-Ramly (2001) developed a methodology for calculating the probability of failure considering the spatial correlation of the soil properties. This methodology is difficult to apply in practice. The present study will therefore evaluate the techniques generally used in probability of failure of slopes. These probability techniques applied to slope stability can be used together with a correction factor proposed by Vanmarcke (1977b). The Second Moment Method was found to be easier to use, and therefore more suitable for geotechnical projects.

[pt] Uma técnica bem conhecida para resolver problemas dependentes do tempo é a formulação, desses problemas, no domínio da frequência por meio da transformada de Laplace ou Fourier, com as consequêntes expressões apropriadas desses resultados utilizando inversões numéricas. Embora de fácil implementação, tais inversões numéricas, são computacionalmente dispendiosas quando resultados mais exatos são desejados e suscetíveis a instabilidades num´ericas. Para problemas de tipo difusão, o algoritmo de Gaver-Stehfest parece ser satisfatório. Problemas gerais de dinâmica demandam algoritmos mais robustos, usualmente baseados em expansões em séries de Fourier tal como foi proposto inicialmente por Dubner e Abate. Algoritmos de outros tipos já são implementados em softwares matemáticos tais como Matlab e Mathematica. A livraria de Fortran possui um algoritmo proposto por Crump e aperfei¸coado por de Hoog e colegas. Mais recentemente, foi proposto resolver problemas transientes de potencial e elasticidade pelo uso de uma técnica avançada de superposição modal que é aplicado a modelos de elementos finitos e elementos de contorno baseados em equilíbrio. O método começa com a formulação no domínio da frequência a qual leva a uma matriz de rigidez efetiva, simétrica-complexa (quando amortecimento viscoso é considerado), expressa como uma série de potências de frequência com matrizes generalizadas de rigidez, amortecimento e massa. Após a solução do problema de autovalor não linear associado, obtém-se uma solução modal avançada do problema, a qual permite uma rápida solução no domínio do tempo obtendo as expressões imediatamente de qualquer resultado de interesse. O objetivo deste trabalho é comparar o desempenho computacional da técnica de superposição modal avançada com as técnicas baseadas em transformadas inversas numéricas de Laplace como aplicações a problemas generais de grande porte. A bibliografia relevante é revista e as principais diferenças conceituais desses métodos são brevemente tratados. Todos os algoritmos são implementados em Fortran com o intuito de garantir uma base comum de comparação. Alguns resultados iniciais são mostrados, conclusões mais definitivas so poderão ser obtidas após uma grande série de simulações numéricas. / [en] An established technique to solve time-dependent problems is the formulation of a complete frequency-domain analysis via Laplace or Fourier transforms, with subsequent ad hoc expression of results by numerical inversion. Although usually easy to implement, such a transform inversion is computationally intensive, if accurate results are desired, and liable to numerical instabilities. For diffusion-type problems, the Gaver-Stehfest algorithm seems well suited. General dynamics problems demand more robust algorithms usually based on Fourier series expansions, as firstly proposed by Dubner and Abate. Algorithms of either kind are already implemented in mathematical languages such as Matlab and Mathematica. The Fortran library has a Fourier-series algorithm proposed by Crump and improved by de Hoog et al. More recently, it has been proposed to solve transient problems of potential and elasticity by using an advanced mode superposition technique that applies to equilibrium-based finite element and boundary element models. One starts with a frequency-domain formulation that leads to a complex-symmetric (if viscous damping is included), effective stiffness matrix expressed as a frequency power series with generalized stiffness, dumping and mass matrices. After solution of the associated complex-symmetric, non-linear eigenvalue problem, one arrives at an advanced modal solution of the problem, which leads to the straightforward solution in the time domain and the immediate expression of any results of interest. Aim of the present research work is to compare the computational efficiency of the proposed advanced modal analysis with the techniques based on numerical inverse transforms, as applied to general, large scale problems. The relevant literature is reviewed and the main conceptual differences of the investigated methods are briefly outlined. All algorithms are implemented in Fortran so as to assure a common basis of comparison. Some initial results are displayed, as more definitive conclusions can only be expected after a large series of numerical simulations.

[pt] A análise de requisitos não funcionais (RNF) é um desafio e vem sendo explorado na literatura científica. Tal iniciativa deve-se ao fato da existência do problema de se verificar o uso das operacionalizações desse tipo de requisito no software construído. Nessa tese apresenta-se um método, com técnicas e ferramentas de apoio, que analisam se um software está em conformidade com padrões de RNFs estabelecidos em catálogo como alternativa para o problema de análise de RNF. A estratégia adotada nessa tese utiliza agentes autônomos para análise de conformidade de software em relação a operacionalizações de RNF. Para isso, utiliza uma base de conhecimentos de padrões persistidos em um catálogo. Os resultados parciais são indicativos de que a proposta de solução é aplicável. A avaliação da validade dá-se por demonstração de que um método parcialmente automatizado é eficaz na identificação de conformidades. Um diferencial do trabalho apresentado é a ligação dos RNFs a sua efetiva implementação. Para demonstração da tese aplicou-se e customizou-se uma técnica de padrões de RNFs, baseados em orientação a metas, em estudos de caso de exemplos do cotidiano prático de software. Apresentamos também a construção de um framework de agentes, que operam sob notações XML para identificar conformidades de software em relação a um catálogo de RNF. / [en] The analysis of non-functional requirements (NFR) is a challenge and has been explored in the literature. This initiative is due to the fact of the existence of the problem of analysis the use of the NFRs operationalization in software. In this thesis we present a method, with supporting tools and techniques, that checks, if a software complies with standards of non-functional requirements as described in a catalog, as an alternative to the NFR analysis problem. The strategy adopted in this thesis uses autonomous agents to check software compliance regarding the operationalization of an NFR, by using a knowledge base of patterns persisted in a catalog. Initial results show that the proposed solution is applicable. The evaluation of the validity is given by the demonstration that a partially automated method is effective in identifying compliance. This work differs form others by linking NFRs to their effective implementation. A method based on patterns NFRs was used in common software, as to show the application of the proposed strategy. An agent based framework, working with XML descriptions, for checking software compliance with respect to a NFR catalog was built.

[pt] Esta dissertação tem como objetivo propor uma abordagem pouco utilizada na área financeira para prever cancelamentos em planos de previdência privada. Métodos e modelos provenientes da análise de sobrevivência foram utilizados para prever o cancelamento, além de estimar o risco associado de o cliente cancelar. Os modelos propostos, modelo de regressão paramétrico e modelo de riscos proporcionais de Cox, foram estimados utilizando-se uma base de dados de clientes acompanhados durante um período de 6 anos e meio de uma seguradora nacional. Os modelos mostraram-se equivalentes, tendo capacidade de generalização de 58 por cento. / [en] This paper aims to propose an approach rarely used in finance to predict cancellations in private pension plans. Methods and models from survival analysiswere used to predict the cancellation and estimate the risk associated with the cliente cancellation. The proposed models, regression model and parametric proportional hazards model of Cox, were estimated using a database of customers followed for a period of 6 and a half years of a national insurer. The models showed equivalence and generalizability of 58 per cent.

[pt] Reservatórios carbonáticos contêm entre 50 e 60 por cento de petróleo e gás do mundo. No entanto, estas rochas têm apresentado problemas operacionais consideráveis durante as etapas de perfuração e produção. Tais problemas são originados pelas dificuldades na caracterização adequada destas rochas em função da complexa distribuição espacial de sua micro e macro estrutura. Enquanto ambas as escalas possuem importância no entendimento do comportamento de carbonatos, o presente trabalho trata da análise da microestrutura. Nesse âmbito, a proposta deste trabalho integra o processamento e análises de imagens digitais adquiridas mediante microtomografia de raios-X, elaboração de malhas de elementos finitos e simulação numérica de forma a prever propriedades elásticas, com o objetivo de correlacionar a microestrutura e o módulo de Young. Foram utilizadas imagens digitais de amostras de afloramento (travertinos) considerados como possíveis análogas a rochas carbonáticas de reservatório. A metodologia implementada permitiu obter sub-amostras que incluem uma gama de microestruturas e porosidades numa única amostra para assim simular numericamente o módulo de Young. Os resultados dos ensaios virtuais foram comparados com os resultados de ensaios reais executados em amostras da mesma rocha e dados encontrados na literatura. A aplicação da metodologia e os resultados obtidos indicam o potencial e as limitações atuais desta técnica. Conclui-se neste estudo que o valor do módulo de Young das amostras simuladas numericamente é afetado pela distribuição espacial, conectividade dos poros e pela microporosidade da rocha. Os valores encontrados no ensaio virtual tendem a serem superiores aos obtidos em ensaios reais devido à dificuldade de inserir, no modelo de micro estrutura, os contatos entre os grãos assim como regiões mais compressíveis não determinadas na aquisição das imagens. / [en] Carbonate reservoirs contain between 50 and 60 per cent of oil and gas in the world. However, these rocks have presented considerable operational problems during drilling and production steps. Such problems are caused by the difficulties in proper characterization of these rocks due the complex spatial distribution of their micro and macrostructure. Although both scales have importance in understanding the behavior of carbonates, the present work deals with the analysis of the microstructure. In this context, the proposal of this work integrates the processing and analysis of digital images acquired by x-ray microtomography, finite element mesh generation and numerical simulations to predict elastic properties in order to correlate the microstructure and the Young s modulus. We used digital images of samples of outcrop (travertine) considered as a possible analogous to carbonate rocks of reservoir. The methodology implemented allowed get sub-samples that include a range of microstructures and porosities in a single sample to determine numerically the Young s modulus. The virtual test results were compared with laboratory test results performed on samples from the same rock and data found in the literature. The application of the methodology and the results obtained indicate the potential and current limitations of this technique. This study concluded that the value of Young s modulus of numerically simulated samples is affected by spatial distribution, pore connectivity and microporosity of these rocks. The values found in the virtual test tend to be higher than those obtained in laboratory tests due to the difficulty of inserting, in the model of micro structure, contacts between the grains as well as more compressible regions not certain on the acquisition of the images.