Modelling ionospheric total electron content (TEC) is an important area of interest for radio wave propagation, geodesy, surveying, the understanding of space weather dynamics and error correction in relation to Global Navigation Satellite Systems (GNNS) applications. With the utilisation of improved ionosonde technology coupled with the use of GNSS, the response of technological systems due to changes in the ionosphere during both quiet and disturbed conditions can be historically inferred. TEC values are usually derived from GNSS measurements using mathematically intensive algorithms. However, the techniques used to estimate these TEC values depend heavily on the availability of near-real time GNSS data, and therefore, are sometimes unable to generate complete datasets. This thesis investigated possibilities for the modelling of TEC values derived from the South African Global Positioning System (GPS)receiver network using linear regression methods and artificial neural networks (NNs). GPS TEC values were derived using the Adjusted Spherical Harmonic Analysis (ASHA) algorithm. Considering TEC and the factors that influence its variability as “dependent and independent variables” respectively, the capabilities of linear regression methods and NNs for TEC modelling were first investigated using a small dataset from two GPS receiver stations. NN and regression models were separately developed and used to reproduce TEC fluctuations at different stations not included in the models’ development. For this purpose, TEC was modelled as a function of diurnal variation, seasonal variation, solar and magnetic activities. Comparative analysis showed that NN models provide predictions of GPS TEC that were an improvement on those predicted by the regression models developed. A separate study to empirically investigate the effects of solar wind on GPS TEC was carried out. Quantitative results indicated that solar wind does not have a significant influence on TEC variability. The final TEC simulation model developed makes use of the NN technique to find the relationship between historical TEC data variations and factors that are known to influence TEC variability (such as solar and magnetic activities, diurnal and seasonal variations and the geographical locations of the respective GPS stations) for the purposes of regional TEC modelling and mapping. The NN technique in conjunction with interpolation and extrapolation methods makes it possible to construct ionospheric TEC maps and to analyse the spatial and temporal TEC behaviour over Southern Africa. For independent validation, modelled TEC values were compared to ionosonde TEC and the International Reference Ionosphere (IRI) generated TEC values during both quiet and disturbed conditions. This thesis provides a comprehensive guide on the development of TEC models for predicting ionospheric variability over the South African region, and forms a significant contribution to ionospheric modelling efforts in Africa.

Le sujet de thèse concerne l’analyse de la propagation des ondes dans les structures complexes et leurs exploitations pour le contrôle semiactif et le contrôle de santé de structures intelligentes. Les structures composites munies de patches piézoélectriques sont la cible principale des investigations. Les patches piézoélectriques sont disposés avec une périodicité. Des travaux précédents ont montré l’intérêt de ce type de configuration pour l’amortissement actif de modes de structures en basses fréquences. L’objectif principal de cette thèse est l’extension de ces constatations dans une bande de fréquences plus large : basses et moyennes fréquences. La maîtrise des paramètres de propagation et de diffusion des ondes est la finalité recherchée. Dans ce cadre, les travaux proposés se baseront sur une technique particulière développée au sein de l’équipe Dynamique des Systèmes et des Structures: la technique WFE (Wave Finite Element), Ondes par éléments finis. Cette approche, construite à l’aide d’un modèle éléments finis d’une cellule représentative de l’essentiel des paramètres de propagation et de diffusion des ondes dans les structures. Elle a été validée sur des cas simples de structures, principalement isotrope monodimensionnel. La modélisation dans ce cas des sandwichs plaques composites munies de couches piézoélectriques sera opérée. Des simulations numériques poussées seront effectuées afin de cerner le cadre d’application de la WFE pour ce type de structures. Des optimisations pourront être réalisées avec ces outils numériques afin d’obtenir des paramètres géométriques et électriques optimaux dans la conception des structures intelligentes. Les travaux de cette thèse sont intégrés dans le projet CALIOP en collaborant avec le laboratoire de Mécanique Appliquée R.Chaléat de l’Institut FEMTOSTet G.W. Woodruff School of Mechanical Engineering de Georgia Institute of Technology. / The analysis of wave propagation in complex structures and its application for the semi-active control of smart structures and health monitoring of these structures are dealt with in this thesis. The design of composite structures with shunted piezoelectric patches is one of the main objectives of all the investigations. This kind of smart composite structures is equipped with periodically distributed shunted piezoelectric patches. Former studies have shown the great interest of such a configuration for the active damping of structural modes at low frequencies. This thesis is focused on the extension of all these interesting characteristics of the smart structures to a larger frequency band: low and medium frequencies. The mastering of the propagation parameters and energy diffusion characteristics is targeted. In this context, the proposed work is based on techniques specifically developed in the research team "Dynamics of Systems and Structures"(D2S): the Wave Finite Element (WFE) method and Diffusion Matrix Model(DMM). The WFE approach is constructed via the finite element model of a unit cell, representative of the waveguide structure. It enables the calculation of essential wave propagation parameters like wavenumbers. The DMM, associated with the WFE approach, enables the calculation of energy diffusion characteristics like reflection and transmission coefficients of specific wave modes. These approaches are extended to consider shunted piezoelectric elements and then to evaluate the performance of shunted piezoelectric patches on the control of wave propagation in the aforementioned smart composite structures. Intensive optimizations can be carried out, with these tools, so as to obtain optimal geometric and electric parameters in the design of these smart structures. The present work is integrated in the CALIOP project in cooperation with the Laboratory of Applied Mechanics R.Chaléat at FEMTO-ST Institute and the G.W. Woodruff School of Mechanical Engineering of Georgia Institute of Technology.

L'exploitation, le contrôle et la signalisation des systèmes de métros modernes et en particulier les métros sans conducteur, reposent sur deux familles de systèmes de communication radio sans fil : les transmissions vitales pour la signalisation et le contrôle des trains à faible débit et les transmissions non vitales très haut débit pour la vidéo surveillance embarquée, le télé-diagnostic,l'information des clients, etc. Lors de la mise en œuvre de tels systèmes de transmission, l'industriel doit donc garantir à l'exploitant de métro des performances concernant les débits et la qualité de la transmission en termes de perte de paquets d'information. Dès lors, appréhender le comportement du canal de propagation est un élément clé pour prédire et améliorer les performances des systèmes de transmission.L'objectif de ces travaux de thèse est double et concerne :• D'une part, la réduction des temps de calcul d'un simulateur de canal à tracé de rayons 3D qui prédit avec précision le comportement du canal de propagation, mais reste coûteux en temps de calculs pour les simulations dynamiques. Nous proposons dans cette thèse une méthode qui s'appuie sur trois critères de visibilité pour simplifier la description de lagéométrie de l'environnement sans pénaliser la prédiction des paramètres caractéristiques du canal de propagation ;• D'autre part, la mise en œuvre d'un nouveau modèle de canal de propagation MIMO dynamique en tunnel rectiligne de section rectangulaire qui permettra d'optimiser des systèmes de transmission multi antennaires (MIMO) pour des applications de transmissions sans fil pour les métros en tunnel. Ce modèle s'inspire globalement du modèle de canal utilisé / The exploitation, control and signalling systems for metro and especially modern driverless metros, are based on two families of radio communication systems Wireless : vital transmissions for signalling and train control low-flow and nonvital transmissions with very high bandwidth for video-tracking systems, telediagnosis, customer information, etc.. During the implementation ofsuch transmission systems, the industrial must therefore ensure the metro operator performance on throughput and transmission quality in terms of packet loss information. Therefore, understanding the behavior of the propagation channel is a key to predict and improve performance of transmission systems.The aim of this thesis is twin and concerns :• On the one hand, the computation time reduction of a 3D ray tracing simulator that accurately predicts the behavior of the propagation channel, but it is expensive in times of calculation for dynamic simulations. We propose here a new method based on three visibility criteria to simplify the geometry description of the environment without degrading the characteristic parameters prediction of the propagation channel ;• On the other hand, the construction of a new model for MIMO propagation channel dynamics in straight tunnel of rectangular section which will optimize transmission systems multi-antennary (MIMO) applications for wireless transmission in metros tunnel. This model draws broadly from the channel model used in the standard WINNER and is fed by the results extracted from the 3D ray tracing channel simulator. The results obtained in this thesis are encouraging and offer many opportunities.

[pt] O esforço de caracterização da influência da vegetação sobre enlaces de radiocomunicação é de grande importância para o dimensionamento de modernos sistemas de comunicação. Baseado em uma extensa campanha de medições, este trabalho pretende caracterizar e propor um modelo de predição de perda de propagação sobre espaços vegetados em banda estreita na faixa de 700 MHz adequado às características morfológicas brasileiras. / [en] The characterization effort of vegetation influence on radiocommunication link fading has a great importance on the design of modern communication systems. Based on an intensive narrow band measurement campaign, this work intends to characterize and propose a model for predicting propagation loss over vegetated areas at the UHF range of 700 MHz.

La turbulence optique et son impact sur les images obtenues à partir d'instruments de mesure est un phénomène connu dans le domaine de l'astronomie. Des modèles issus de la théorie de Kolmogorov, développée pour une turbulence dynamique, ainsi que des méthodes de correction, telles que l'optique adaptative, existent pour l'analyse et la compensation des effets de cette turbulence optique. L'analyse de cette dernière dans les milieux confinés est cependant plus limitée. Les sources susceptibles de générer une turbulence optique dans ces espaces sont pourtant multiples et peuvent avoir un impact non négligeable sur les mesures des instruments installés à proximité. Ce mémoire constitue une première approche de la caractérisation de la turbulence optique dans un espace confiné. Après l'introduction des phénomènes étudiés et des outils mathématiques employés, ce mémoire présente les résultats issus d'analyses de coupoles de télescopes et de caractérisations de salles blanches employées par l'entreprise Thales Alenia Space pour l'intégration et le test d'instruments optiques. Ces résultats sont obtenus avec l'instrument INTENSE, développé durant la thèse pour la caractérisation de la turbulence optique locale via la mesure des fluctuations des angles d'arrivée de multiples faisceaux lasers. En prévision de futurs améliorations de l'analyse de la turbulence en espace confiné, un chapitre du mémoire est dédié aux travaux réalisés sur une méthode d'extraction du profil de l'énergie de la turbulence et à son application sur l'instrument INTENSE. Les conclusions et les perspectives des travaux réalisés pendant la thèse sont présentées à la fin du mémoire. / Optical turbulence and its impact on measured images is a well-known phenomenon in astronomy. Models based on the Kolmogorov theory, elaborated for a dynamical turbulence description, and methods, such as Adaptive Optics, were both developed so as to understand and correct the degradations caused by this turbulence. Analysis of the same phenomenon in indoor situation was however less investigated. The local air volume is nonetheless prone to optical perturbations sources which could have non negligible impacts on the measurements of instruments installed at proximity. This document introduces a first approach of indoor optical turbulence characterization. After the introduction of the studied phenomenon and the mathematical tools employed, this thesis present optical turbulence characterizations inside Thales Alenia Space clean rooms used for optical instrument integration and testing. Analyses inside telescope domes are also shown in this document. All the results were obtained using the INTENSE instrument which was developed during the thesis for optical turbulence characterizations using angle of arrival fluctuations of laser beams. In anticipation for future ameliorations of optical turbulence analysis methods, a chapter of this thesis is dedicated to the work made on a turbulence energy profile extraction and its application on the INTENSE instrument. Conclusions and perspectives of the work made during this thesis are presented at the end of the document.

Les ondes de surface se propagent le long de la surface terrestre par-delà l'horizon et ces propriétés de propagation sont des atouts pour diverses applications industrielles (civiles ou militaires) opérant dans la bande HF. C'est le cas notamment des radars transhorizon, des systèmes de communications, de la détection de tsunami, ou encore l'étude de la courantologie. Actuellement, les aériens utilisés dans ce contexte ne sont pas optimisés pour privilégier la propagation par ondes de surface. C'est dans ce contexte que TDF, la Direction Générale de l'Armement (DGA) et l'Institut d'Electronique et de Télécommunications de Rennes (IETR) se sont associés pour proposer cette thèse. L'objectif principal de la thèse est de concevoir des antennes en bande HF de dimensions réduites (<3 mètres) réduisant le rayonnement vers le ciel afin favoriser l'onde de surface. Une première approche a permis de concevoir deux antennes élémentaires, l'une omnidirectionnelle et l'autre directive, permettant de réduire le rayonnement vers le ciel pour les sites compris entre -45° et +45°. Les solutions proposées sont de tailles très réduites grâce à l'utilisation d'une antenne à polarisation verticale électriquement courte dans le plan vertical et de parasites horizontaux alimentés par couplage EM. Les études expérimentales réalisées à échelle réelle sont en bon accord avec les performances théoriques obtenus à l'aide de logiciels d'analyse électromagnétique commerciaux. Une seconde approche a permis de concevoir des réseaux d'antennes électriquement courtes dans le plan vertical selon une disposition linéaire des aériens. La réduction du rayonnement vers le ciel est réalisée grâce aux réseaux de type Endfire/Broadside. Enfin, une troisième étude a permis de concevoir deux systèmes antennaires agiles en diagrammes tout en conservant les propriétés de réduction de rayonnement vers le ciel. Le premier système à l'aide d'antennes électriquement courtes dans le plan vertical chargées par des réactances, et le second en appliquant une pondération en amplitude et en phase d'un réseau à deux dimensions. / The surface-waves propagate beyond the horizon. This phenomenon is an attractive feature for many applications operating in the HF band both in civil or military field. By example, we can quote the High Frequency Surface Waves Radar (HFSWR), communicating systems, or the study of the currentology. Currently, the antennas used in the field of surface-waves propagation are not fully optimized and radiate also in sky wave. In order to optimize the radiation pattern of the antennas used in these systems, the company TDF, the Direction Générale de l’Armement (DGA) and the Institute of Electronic and Telecommunications of Rennes (IETR) have proposed this projetct. The main objective of this work is to design electrically small antennas in the vertical plane (< 3 meters) in the HF band while reducing the radiation toward the sky. The first concept is based on vertically polarized antennas electrically small in the vertical plane coupled to horizontally radiators placed above the antenna. As a result, the radiation pattern in the horizontal plane is omnidirectional or directive with a radiation toward the sky reduced for the elevations range from -45° to +45°. Prototypes have been realized and validated by measurements. Then, in the second study we have developed antennas array in order to reduce the radiation toward the sky with another way that the horizontal radiators. This method allows reducing significantly the radiation toward the sky according to the number of elements in the array. Finally, in the third study we have studied electronically scanning antennas while maintaining the reduction of the radiation toward the sky. The first system is based on an array of radiator loaded by reactance, and the second is based on a planar array with an appropriate feeding.

Survey propagation (SP) has recently been discovered as an efficient algorithm in solving classes of hard constraint-satisfaction problems (CSP). Powerful as it is, SP is still a heuristic algorithm, and further understanding its algorithmic nature, improving its effectiveness and extending its applicability are highly desirable. Prior to the work in this thesis, Maneva et al. introduced a Markov Random Field (MRF) formalism for k-SAT problems, on which SP may be viewed as a special case of the well-known belief propagation (BP) algorithm. This result had sometimes been interpreted to an understanding that “SP is BP” and allows a rigorous extension of SP to a “weighted” version, or a family of algorithms, for k-SAT problems. SP has also been generalized, in a non-weighted fashion, for solving non-binary CSPs. Such generalization is however presented using statistical physics language and somewhat difficult to access by more general audience. This thesis generalizes SP both in terms of its applicability to non-binary problems and in terms of introducing “weights” and extending SP to a family of algorithms. Under a generic formulation of CSPs, we first present an understanding of non-weighted SP for arbitrary CSPs in terms of “probabilistic token passing” (PTP). We then show that this probabilistic interpretation of non-weighted SP makes it naturally generalizable to a weighted version, which we call weighted PTP. Another main contribution of this thesis is a disproof of the folk belief that “SP is BP”. We show that the fact that SP is a special case of BP for k-SAT problems is rather incidental. For more general CSPs, SP and generalized SP do not reduce from BP. We also established the conditions under which generalized SP may reduce as special cases of BP. To explore the benefit of generalizing SP to a wide family and for arbitrary, particularly non-binary, problems, we devised a simple weighted PTP based algorithm for solving 3-COL problems. Experimental results, compared against an existing non-weighted SP based algorithm, reveal the potential performance gain that generalized SP may bring.

Le travail réalisé pendant cette thèse a permis d'apporter de nouvelles contraintes sur le développement et la distribution de la déformation dans le manteau supérieur et plus particulièrement au niveau des grandes limites de plaques décrochantes. Grâce à l'apport de l'expérience USArray et d'une dizaine d'années d'enregistrements sismologiques supplémentaires, nous avons pu étudier, de manière précise, les variations d'anisotropie dans le voisinage de la Faille de San Andreas. Nous avons confirmé et étendu l'observation de deux couches anisotropes sous cette limite de plaque. On y observe une première couche localisée dans la lithosphère marquant la déformation induite à la limite de plaque, et une autre, asthénosphérique, cohérente avec l'anisotropie observée loin de la faille et d'origine plus discutée. Nous avons montré que la zone de déformation associée aux failles de San Andreas, Calaveras et d'Hayward a, vraisemblablement, une largeur d'au moins 40 kilomètres en base de lithosphère, sous chacune de ces failles. Nous avons ensuite procédé à la modélisation thermomécanique (ADELI) de la migration d'une limite de plaques décrochante couplée à une modélisation du développement de fabriques cristallographiques par une approche viscoplastique auto-cohérente (VPSC). Ceci nous a permis d'y observer le développement de la déformation et les conséquences des possibles interactions entre la déformation décrochante en surface et le cisaillement en base de lithosphère dû au déplacement horizontal des plaques. Les propriétés élastiques déduites des fabriques cristallographiques modélisées montrent que de telles interactions existent et provoquent, sous la limite de plaques, une rotation des orientations cristallographiques avec la profondeur. Le signal associé à ces rotations progressives n'est toutefois pas cohérent avec la présence de deux couches d'anisotropie comme proposée sous la faille de San Andreas. Nous pensons par conséquent qu'il existe, sous la Californie, une zone de découplage entre la lithosphère et l'asthénosphère, permettant d'individualiser une déformation lithosphérique d'une déformation asthénosphérique. Nous estimons, en outre, que l'anisotropie observée dans l'asthénosphère sous la Californie ne peut être expliquée seulement par le cisaillement induit par le déplacement de la lithosphère Nord Amérique. En effet, les propriétés anisotropes obtenues par modélisation à partir d'une plaque se déplaçant dans une direction et une vitesse proche de celle de la plaque Amérique du Nord montrent qu'on ne peut espérer guère plus que quelques dixièmes de seconde de délai au bout de 10 Ma de déplacement. Les déphasages mesurés en Californie étant de l'ordre de 1,5 s, il est donc nécessaire d'invoquer la présence d'écoulements mantelliques actifs sous cette région / This work provides new constraints on the development and on the distribution of the deformation in the upper mantle and particularly beneath transform plate boundaries. USArray experiment and the remarkable increase of the dataset in California for the past ten years allowed us to scrutinize the lateral variations of the anisotropy in the vicinity of the San Andreas Fault zone. We have confirmed and increased the detection of two layers of anisotropy beneath this plate boundary. The first layer, located in the lithosphere, is related to the deformation induced at the fault, and the other one, located in the asthenosphere, is coherent with the anisotropy observed far from it, its origin is however less clear. We show that the deformation zone associated both to the San Andreas, Calaveras and Hayward Faults, is likely 40 km wide at 70 km depth. We then performed numerical thermomechanical modeling (ADELI) of the displacement of a transform plate boundary associated with the computation of the development of crystallographic fabrics using a viscoplastic self-consistent approach (VPSC). We analyzed the distribution of the deformation in the model ant looked after the possible interactions at depth between deformation caused at surface by the strike-slip dynamic of the fault and the shearing at the base of the lithosphere caused by the horizontal displacement of the plates. Elastic properties derived from the crystallographic fabrics modeled, show that such interactions exist and induce, beneath the fault zone, a progressive rotation of the crystallographic fabrics with depth. Seismological signature of these smooth rotations is however not relevant with the presence of two anisotropic layers as proposed beneath California. We thus consider that a decoupling zone exists between the lithosphere and the asthenosphere beneath the California to account for the sharp separation between a lithospheric and an asthenospheric deformation. We furthermore estimate that anisotropy observed far form the San Andreas Fault in California cannot be explained only by the drag of the asthenosphere by the North America lithosphere as proposed in our article. Indeed, we can only expect few tenths of second of splitting delay from the anisotropic properties derived from the numerical modeling of a plate moving in the same direction and in the same velocity than the North American lithosphere only for 10 Ma of displacement. As delays observed in California rather reach 1.5 s, anisotropy in this region thus requires the existence of an active asthenospheric flow to be explained.

L’un des éléments clé pour atteindre de hauts rendements dans les centrales solaires à récepteur volumétrique est l’absorbeur. Sa structure poreuse permet au rayonnement solaire de pénétrer à l’intérieur, où l'énergie absorbée doit être transférée efficacement par convection au fluide qui la traverse. Dans ce travail, deux types d’absorbeurs innovants sont présentés et analysés : l’un constitué par une série d’éléments empilés avec une structure en forme de grille, l’autre par des éléments similaires à des stores vénitiens. A titre de référence, un absorbeur classique à nid d’abeilles est aussi évalué. La propagation du rayonnement solaire au sein des absorbeurs est modélisée au moyen de la technique de lancer de rayons, basée sur la méthode de Monte Carlo. Leur comportement thermo-hydraulique est simulé par la méthode des éléments finis. Les caractéristiques géométriques des deux absorbeurs proposés améliorent le transfert thermique par convection par rapport aux absorbeurs alvéolaires et la modification des leurs principaux paramètres géométriques nous a permis d’augmenter la longueur d’extinction du rayonnement solaire. Cependant, l’accroissement de leur surface frontale apparente augmente les pertes par réflexion. A l’issu des résultats théoriques, l’absorbeur à stores vénitiens a été retenu pour l’analyser expérimentalement. Ses performances thermiques sont comparées avec celles d’un absorbeur alvéolaire. Ces résultats montrent également que la structure de l’absorbeur proposé intensifie les échanges thermiques vers le fluide. De plus, ce type d’absorbeur atteint un meilleur comportement thermique à de hauts flux radiatifs et à des débits élevés. / One of the key elements for achieving a high efficiency in solar power plants with volumetric receiver is the absorber. Its porous structure allows the solar radiation to penetrate inside it, where the absorbed energy should be transferred efficiently by convection to the fluid which crosses through it. In this work, two types of innovative absorbers are presented and analyzed: one consisting of a set of stacked elements with a grid-like structure, the other with elements similar to venetian blinds. As a reference, a conventional honeycomb absorber is also evaluated. The solar radiation propagation within the absorbers is modeled through the ray tracing technique, based on the Monte Carlo method. Their thermohydraulic behavior is simulated by the finite element method. The geometrical characteristics of the two proposed absorbers improve convective heat transfer compared to honeycomb absorbers and the modification of their main geometric parameters allowed us to increase the extinction length of solar radiation. However, the increase of their apparent frontal surface rises up reflection losses. Based on the theoretical results, the venetian blind absorber was selected to analyze it experimentally. Its thermal performance is compared with that of a honeycomb absorber. These results also indicate that the structure of the proposed absorber intensifies the heat exchange to the fluid. Moreover, this kind of absorber reaches a better thermal behavior at high heat flux and at high flow rates.

Les trois principaux phénomènes acoustiques propres au milieu forestier nécessitant d'être pris en compte sont (1) l'absorption due à la présence d'un sol multi-couche, (2) la diffusion multiple due à la présence d'obstacles tels que les troncs, (3) les effets micro-météorologiques rattachés aux variations des gradients de vitesse de vent et de température. Parmi les méthodes numériques de référence, la méthode des lignes de transmission (TLM), semble particulièrement adaptée pour la modélisation de la propagation acoustique en présence de forêt, à condition de procéder à de nouveaux développements. La première nécessité pour l'adaptation de la méthode TLM aux simulations acoustiques sur de grandes distances est la définition de couches absorbantes, permettant de tronquer efficacement le domaine d'étude, sans introduire de réflexions parasites. La formulation ainsi développée dans le cadre de la thèse est rigoureusement équivalente à l'équation de propagation des ondes amorties, et se traduit dans la méthode TLM par l'introduction et l'optimisation d'un terme de dissipation. L'étape suivante a consisté à vérifier la capacité de la méthode TLM à modéliser les phénomènes de diffusion par des cylindres. L’une des originalités introduites dans cette thèse réside dans le placement des éléments diffuseurs, à partir de lois de distribution aléatoire et de Gibbs, permettant ainsi de définir des répartitions proches de celles rencontrées en forêt. À titre d'application de la méthode développée dans le cadre de la thèse, une étude paramétrique a été réalisée afin de définir les conditions pour lesquelles une forêt peut également être considérée comme un dispositif de protection. / The prediction of sound propagation in presence of forest remains a major challenge for the outdoor sound propagation community. Reference numerical models such as the Transmission Line Matrix (TLM) method can be developed in order to accurately predict each acoustical phenomenon that takes place inside forest. The first need for the TLM method is an efficient theory-based absorbing layer formulation that enables the truncation of the numerical domain. The two proposed absorbing layer formulations are based on the approximation of the perfectly matched layer theory. The most efficient proposed formulation is shown to be equivalent to wave propagation in a lossy media, which, in the TLM method formulation, is introduced using an additional dissipation term. Then, the ability of the TLM method for the simulation of scattering is studied comparing the numerical results to both analytical solutions and measurements on scale models. Lastly, the attenuation of acoustic levels by a simplified forest is numerically studied using several arrangements of cylinders placed normal to either reflecting or absorbing ground. It is observed that randomly spaced arrangements are more inclined to attenuate acoustic waves than periodic arrangements. Moreover, the sensitivity to the density, the length of the array and the ground absorption is tested. The main trend shows that the density and the distribution are two important parameters for the attenuation. In future work, it can be interesting to look at the sensitivity of each parameter. This study could then be used to relate the morphology (i.e. distribution, density, length) of a forest to the acoustical properties of the forest.