Avec l’évolution rapide des techniques de nanofabrication et les besoins croissants d’intégration et d’utilisation à moindre cout énergétique, l’étude et la manipulation de résonances électromagnétiques d’objets de faibles dimensions représentent des enjeux cruciaux. Un des objectifs de ce travail de thèse a donc été d’approfondir nos connaissances de l’interaction entre le champ électromagnétique et la matière. Dans ce but la réalisation conjointe d’expériences en champ proche optique et le développement numérique de modèles associés nous ont permis d’étudier différentes résonances électromagnétiques basées sur l’interaction lumière-matière au sein d’objets de dimensions sub-micrométriques. Dans une première partie, ce manuscrit présente les phénomènes mis en jeu en microscopie champ proche optique et décrit le fonctionnement du microscope SNOM utilise. La seconde partie est dédiée à l’étude de l’interaction d’une sonde champ proche avec une nano-cavité Fabry-Perot en régime non-linéaire. Dans un premier temps, on présente l’étude des non-linéarités de nano-cavités en Silicium à grand facteur de qualité et faible volume modal, démontrant ainsi l’obtention d’un régime de fonctionnement bistable pour de très faibles puissances. Dans un deuxième temps, on démontre la modulation possible via la sonde champ proche de ce régime de bistabilité. Enfin, on étudie dans la troisième partie de ce manuscrit les résonances de films minces métalliques. L’étude complète des modes propres et de l’excitation optique de cette structure a permis de connaitre avec précision les modes résonants et les angles de Brewster du film mince. Cette étude a été prolongée expérimentalement par l’étude en champ proche optique des modes résonants de demi-films minces métalliques. Lors de ces deux études, les mesures expérimentales ont été systématiquement accompagnées par l’analyse théorique et le développement numérique de modèles que les expériences réalisées ont permis de valider et de discuter. / Within the quick nano-fabrication techniques evolution and the increasing needs in low space and energy opto-electronic functions on-chip integration, studies and control of electromagnetic resonances of sub-micrometric objects are crucial issues. One aim of this work was then to improve our knowledge of the light-matter interaction. To this end, optical near-field experiments and complementary numerical development of models have been performed to study electromagnetic resonances due to light-matter interaction in sub-micrometric resonators. As an introduction, the first part of the manuscript presents the main phenomena involved in near-field optical microscopy and the used SNOM operation. The second part is dedicated to nearfield tip/Fabry-Perot micro-resonator interaction in the non-linear regime. In this part, we first present the study of high-Q/small volume Silicon nano-cavity non-linearities, showing possible bistable operation at very low incident power. We then demonstrate the cavity stability modulation obtain with the near-field tip. In the third part of the manuscript, we consider thin metallic film resonances. As a beginning, we realize an exhaustive study of the film eigenmodes and then discuss the resonant optical excitation of the film as well as precise determination of Brewster angle. These theoretical studies have been extended experimentally with the optical near-field determination of the resonant modes of semi-infinite thin films. For the different studies, experimental measures have been systematically associated with theoretical analyses and numerical developments of models, which have then been confirmed or discussed in light of the experimental results.

Comme la puissance de calcul des systèmes de calcul haute performance continue de croître, en utilisant un grand nombre de cœurs CPU ou d’unités de calcul spécialisées, les applications hautes performances destinées à la résolution des problèmes de très grande échelle sont de plus en plus sujettes à des pannes. En conséquence, la communauté de calcul haute performance a proposé de nombreuses contributions pour concevoir des applications tolérantes aux pannes. Cette étude porte sur une nouvelle classe d’algorithmes numériques de tolérance aux pannes au niveau de l’application qui ne nécessite pas de ressources supplémentaires, à savoir, des unités de calcul ou du temps de calcul additionnel, en l’absence de pannes. En supposant qu’un mécanisme distinct assure la détection des pannes, nous proposons des algorithmes numériques pour extraire des informations pertinentes à partir des données disponibles après une pannes. Après l’extraction de données, les données critiques manquantes sont régénérées grâce à des stratégies d’interpolation pour constituer des informations pertinentes pour redémarrer numériquement l’algorithme. Nous avons conçu ces méthodes appelées techniques d’Interpolation-restart pour des problèmes d’algèbre linéaire numérique tels que la résolution de systèmes linéaires ou des problèmes aux valeurs propres qui sont indispensables dans de nombreux noyaux scientifiques et applications d’ingénierie. La résolution de ces problèmes est souvent la partie dominante; en termes de temps de calcul, des applications scientifiques. Dans le cadre solveurs linéaires du sous-espace de Krylov, les entrées perdues de l’itération sont interpolées en utilisant les entrées disponibles sur les nœuds encore disponibles pour définir une nouvelle estimation de la solution initiale avant de redémarrer la méthode de Krylov. En particulier, nous considérons deux politiques d’interpolation qui préservent les propriétés numériques clés de solveurs linéaires bien connus, à savoir la décroissance monotone de la norme-A de l’erreur du gradient conjugué ou la décroissance monotone de la norme résiduelle de GMRES. Nous avons évalué l’impact du taux de pannes et l’impact de la quantité de données perdues sur la robustesse des stratégies de résilience conçues. Les expériences ont montré que nos stratégies numériques sont robustes même en présence de grandes fréquences de pannes, et de perte de grand volume de données. Dans le but de concevoir des solveurs résilients de résolution de problèmes aux valeurs propres, nous avons modifié les stratégies d’interpolation conçues pour les systèmes linéaires. Nous avons revisité les méthodes itératives de l’état de l’art pour la résolution des problèmes de valeurs propres creux à la lumière des stratégies d’Interpolation-restart. Pour chaque méthode considérée, nous avons adapté les stratégies d’Interpolation-restart pour régénérer autant d’informations spectrale que possible. Afin d’évaluer la performance de nos stratégies numériques, nous avons considéré un solveur parallèle hybride (direct/itérative) pleinement fonctionnel nommé MaPHyS pour la résolution des systèmes linéaires creux, et nous proposons des solutions numériques pour concevoir une version tolérante aux pannes du solveur. Le solveur étant hybride, nous nous concentrons dans cette étude sur l’étape de résolution itérative, qui est souvent l’étape dominante dans la pratique. Les solutions numériques proposées comportent deux volets. A chaque fois que cela est possible, nous exploitons la redondance de données entre les processus du solveur pour effectuer une régénération exacte des données en faisant des copies astucieuses dans les processus. D’autre part, les données perdues qui ne sont plus disponibles sur aucun processus sont régénérées grâce à un mécanisme d’interpolation. / As the computational power of high performance computing (HPC) systems continues to increase by using huge number of cores or specialized processing units, HPC applications are increasingly prone to faults. This study covers a new class of numerical fault tolerance algorithms at application level that does not require extra resources, i.e., computational unit or computing time, when no fault occurs. Assuming that a separate mechanism ensures fault detection, we propose numerical algorithms to extract relevant information from available data after a fault. After data extraction, well chosen part of missing data is regenerated through interpolation strategies to constitute meaningful inputs to numerically restart the algorithm. We have designed these methods called Interpolation-restart techniques for numerical linear algebra problems such as the solution of linear systems or eigen-problems that are the inner most numerical kernels in many scientific and engineering applications and also often ones of the most time consuming parts. In the framework of Krylov subspace linear solvers the lost entries of the iterate are interpolated using the available entries on the still alive nodes to define a new initial guess before restarting the Krylov method. In particular, we consider two interpolation policies that preserve key numerical properties of well-known linear solvers, namely the monotony decrease of the A-norm of the error of the conjugate gradient or the residual norm decrease of GMRES. We assess the impact of the fault rate and the amount of lost data on the robustness of the resulting linear solvers.For eigensolvers, we revisited state-of-the-art methods for solving large sparse eigenvalue problems namely the Arnoldi methods, subspace iteration methods and the Jacobi-Davidson method, in the light of Interpolation-restart strategies. For each considered eigensolver, we adapted the Interpolation-restart strategies to regenerate as much spectral information as possible. Through intensive experiments, we illustrate the qualitative numerical behavior of the resulting schemes when the number of faults and the amount of lost data are varied; and we demonstrate that they exhibit a numerical robustness close to that of fault-free calculations. In order to assess the efficiency of our numerical strategies, we have consideredan actual fully-featured parallel sparse hybrid (direct/iterative) linear solver, MaPHyS, and we proposed numerical remedies to design a resilient version of the solver. The solver being hybrid, we focus in this study on the iterative solution step, which is often the dominant step in practice. The numerical remedies we propose are twofold. Whenever possible, we exploit the natural data redundancy between processes from the solver toperform an exact recovery through clever copies over processes. Otherwise, data that has been lost and is not available anymore on any process is recovered through Interpolationrestart strategies. These numerical remedies have been implemented in the MaPHyS parallel solver so that we can assess their efficiency on a large number of processing units (up to 12; 288 CPU cores) for solving large-scale real-life problems.

Cette thèse porte sur la compréhension et le contrôle des interactions dynamiques entre les mécanismes physiques et biologiques en considérant un procédé alternatif de séparation membranaire pour les bioprocédés industriels. L’objectif premier est un apport de connaissances scientifiques liées à la maîtrise de la bioréaction en considérant l'hydrodynamique complexe et les verrous rétention-perméation. Une technologie de filtration dynamique, appelée Rotating and Vibrating Filtration (RVF), a été spécifiquement étudiée. Elle se compose de cellules de filtration en série comprenant deux membranes circulaires planes fixées sur des supports poreux au voisinage d'un agitateur à trois pales planes attachées à un arbre central. Ce dispositif mécanique simple fonctionne en continu et génère une contrainte de cisaillement élevée ainsi qu'une perturbation hydrodynamique dans un entrefer étroit (pale-membrane). Les verrous scientifiques et techniques qui motivent ce travail, sont la caractérisation et la quantification (i) des champs de vitesse locaux et instantanés, (2) des contraintes pariétales de cisaillement à la surface de la membrane et (3) l'impact mécanique sur les cellules microbiennes.Dans ce but, des expériences et des simulations numériques ont été réalisées pour étudier l'hydrodynamique à des échelles globales et locales, en régimes laminaire et turbulent avec des fluides newtoniens dans des environnements biotique et abiotique. Pour l'approche globale, la distribution des temps de séjour (RTD) et le bilan thermique ont été réalisés et comparés aux précédentes études globales (courbes de consommation de puissance et de frottement). Une étude analytique des fonctions de distribution a été effectuée et les moments statistiques ont été calculés et discutés. Une analyse systémique a été utilisée pour décrire les comportements hydrodynamiques du module RVF. En combinant la simulation des écoulements (CFD) et les observations (RTD), les conditions et les zones de dysfonctionnement des cellules de filtration sont éclairées. Pour l'approche locale, la vélocimétrie laser (PIV) a été réalisée dans les plans horizontaux et verticaux et comparée à la simulation numérique (CFD). Une étude préliminaire basée sur une synchronisation entre la prise d’image et la position de l’agitateur (résolution angulaire) a permis d’accéder aux champs de vitesse moyens. Une campagne de mesure PIV a été réalisée sans synchronisation afin d’appliquer une décomposition orthogonale aux valeurs propres (POD) pour 'identifier les composantes moyennes, organisées et turbulentes des champs de vitesse (énergie cinétique). Pour l'application aux bioprocédés, un travail exploratoire a caractérisé l'effet de la filtration dynamique sur des cellules procaryotes (E. coli) en quantifiant l'intégrité cellulaire ou leur dégradation en fonction du temps et de la vitesse de rotation. / This thesis focuses on the understanding and the control of dynamic interactions between physical and biological mechanisms considering an alternative membrane separation into industrial bioprocess. It aims to carry scientific knowledge related to the control of bioreaction considering complex hydrodynamics and retention-permeation locks specific to membrane separation. A dynamic filtration technology, called Rotating and Vibrating Filtration (RVF), was investigated. It consists of filtration cells in series including two flat disc membranes fixed onto porous substrates in the vicinity of a three-blade impeller attached to a central shaft. This simple mechanical device runs continuously and generates a high shear stress as well as a hydrodynamic perturbation in the narrow membrane-blade gap. Several scientific and technical locks motivating this work are to characterize and to quantify (i) the velocity fields locally and instantaneously, (2) the shear stresses at membrane surface and (3) the mechanical impact on microbial cells.To this end, experiments and numerical simulations have been performed to investigate the hydrodynamics at global and local scales under laminar and turbulent regimes with Newtonian fluids under biotic and abiotic environment. For global approach, investigation of Residence Time Distribution (RTD) and thermal balance was carried out and compared to the previous global study (power consumption and friction curves). Analytical study of distribution functions was conducted and statistical moments were calculated and discussed. A systemic analysis was used to describe the hydrodynamic behaviors of the RVF module. Combining Computational Fluid Dynamics (CFD) and RTD observations, it leads to demonstrate dysfunctioning conditions and area. For the local approach, Particle Image Velocimetry (PIV) was be carried out in both horizontal and vertical planes and compared to CFD simulation. PIV preliminary study was conducted with a trigger strategy to access through angle-resolved measurements to an averaged velocity field. PIV further study were performed with a non-trigger strategy and applied to Proper Orthogonal Decomposition (POD) analysis in order to identify the coherent structure of the flow by decomposing the organized and turbulent fluctuations. For the bioprocess application, an exploratory work characterized the effect of Dynamic Filtration on prokaryote cell population (Escherichia coli) by quantifying cell integrity or damage as a function of time and rotation speed during filtration process in turbulent regime.

Cette thèse aborde le processus ponctuel déterminantal, un modèle probabiliste qui capture la répulsion entre les points d’un certain espace. Celle-ci est déterminée par une matrice de similarité, la matrice noyau du processus, qui spécifie quels points sont les plus similaires et donc moins susceptibles de figurer dans un même sous-ensemble. Contrairement à la sélection aléatoire uniforme, ce processus ponctuel privilégie les sous-ensembles qui contiennent des points diversifiés et hétérogènes. La notion de diversité acquiert une importante grandissante au sein de sciences comme la médecine, la sociologie, les sciences forensiques et les sciences comportementales. Le processus ponctuel déterminantal offre donc une alternative aux traditionnelles méthodes d’échantillonnage en tenant compte de la diversité des éléments choisis. Actuellement, il est déjà très utilisé en apprentissage automatique comme modèle de sélection de sous-ensembles. Son application en statistique est illustrée par trois articles. Le premier article aborde le partitionnement de données effectué par un algorithme répété un grand nombre de fois sur les mêmes données, le partitionnement par consensus. On montre qu’en utilisant le processus ponctuel déterminantal pour sélectionner les points initiaux de l’algorithme, la partition de données finale a une qualité supérieure à celle que l’on obtient en sélectionnant les points de façon uniforme. Le deuxième article étend la méthodologie du premier article aux données ayant un grand nombre d’observations. Ce cas impose un effort computationnel additionnel, étant donné que la sélection de points par le processus ponctuel déterminantal passe par la décomposition spectrale de la matrice de similarité qui, dans ce cas-ci, est de grande taille. On présente deux approches différentes pour résoudre ce problème. On montre que les résultats obtenus par ces deux approches sont meilleurs que ceux obtenus avec un partitionnement de données basé sur une sélection uniforme de points. Le troisième article présente le problème de sélection de variables en régression linéaire et logistique face à un nombre élevé de covariables par une approche bayésienne. La sélection de variables est faite en recourant aux méthodes de Monte Carlo par chaînes de Markov, en utilisant l’algorithme de Metropolis-Hastings. On montre qu’en choisissant le processus ponctuel déterminantal comme loi a priori de l’espace des modèles, le sous-ensemble final de variables est meilleur que celui que l’on obtient avec une loi a priori uniforme. / This thesis presents the determinantal point process, a probabilistic model that captures repulsion between points of a certain space. This repulsion is encompassed by a similarity matrix, the kernel matrix, which selects which points are more similar and then less likely to appear in the same subset. This point process gives more weight to subsets characterized by a larger diversity of its elements, which is not the case with the traditional uniform random sampling. Diversity has become a key concept in domains such as medicine, sociology, forensic sciences and behavioral sciences. The determinantal point process is considered a promising alternative to traditional sampling methods, since it takes into account the diversity of selected elements. It is already actively used in machine learning as a subset selection method. Its application in statistics is illustrated with three papers. The first paper presents the consensus clustering, which consists in running a clustering algorithm on the same data, a large number of times. To sample the initials points of the algorithm, we propose the determinantal point process as a sampling method instead of a uniform random sampling and show that the former option produces better clustering results. The second paper extends the methodology developed in the first paper to large-data. Such datasets impose a computational burden since sampling with the determinantal point process is based on the spectral decomposition of the large kernel matrix. We introduce two methods to deal with this issue. These methods also produce better clustering results than consensus clustering based on a uniform sampling of initial points. The third paper addresses the problem of variable selection for the linear model and the logistic regression, when the number of predictors is large. A Bayesian approach is adopted, using Markov Chain Monte Carlo methods with Metropolis-Hasting algorithm. We show that setting the determinantal point process as the prior distribution for the model space selects a better final model than the model selected by a uniform prior on the model space.

Ce mémoire porte sur l’étude du laplacien sur des surfaces de Riemann. En particulier, nous nous intéressons à ses valeurs propres qui représentent les notes que jouerait la surface si elle était un tambour. Les valeurs les plus étudiées sont la première valeur propre non nulle λ1 ainsi que sa multiplicité m1 (la dimension de l’espace propre). Notamment, Colin de Verdière conjecturait que m1 est toujours borné supérieurement par le nombre chromatique moins 1. Des travaux de Fortier Bourque et Petri ont montré que parmi toutes les surfaces hyperboliques de genre 3, c’est la quartique de Klein qui maximise la multiplicité et atteint la borne supérieure conjecturée par Colin de Verdière. Cette surface est la première d’une suite de surfaces hautement symétriques, les surfaces de Hurwitz. Nous montrons à l’aide de la formule des traces de Selberg que pour la prochaine surface dans la suite, la surface de Fricke–Macbeath F, nous avons m1(F) = 7. Une recherche indépendante menée par Chul-hee Lee arrive au même résultat à propos de la multiplicité. Le chapitre 1 introduit des notions géométriques comme la géométrie hyperbolique, les surfaces hyperboliques et triangulaires ainsi que le théorème de Hurwitz. Le chapitre 2 présente des concepts de base de théorie spectrale ainsi que des outils comme la formule des traces de Selberg et la théorie de la représentation. Le chapitre 3 est dédié à l’étude de la surface de Fricke–Macbeath et à la preuve de notre résultat principal à l’aide des outils des chapitres précédents. Dans le chapitre 4, nous discutons de nouvelles techniques de calcul de m1 qui ont été utilisées pour montrer l’existence de contre-exemples à la conjecture de Colin de Verdière dans des travaux conjoints avec Fortier Bourque, Gruda-Mediavilla et Petri. / This master’s thesis studies the Laplace operator on Riemann surfaces. We are especially interested in its eigenvalues, which correspond to the notes that the surface would play if it were a drum. In particular, the first non-zero eigenvalue λ1 and its multiplicity m1 (the dimension of the corresponding eigenspace) have been well studied. For instance, Colin de Verdière conjectured that m1 is bounded above by the chromatic number minus 1 based on a few examples. Later work by Fortier Bourque and Petri showed that among hyperbolic surfaces of genus 3, the Klein quartic maximizes the multiplicity, and attains the upper bound conjectured by Colin de Verdière. This surface is the first of a sequence of highly symmetrical surfaces named Hurwitz surfaces. We will show using the Selberg trace formula that for the next surface in the sequence, the Fricke–Macbeath surface F, we have m1(F) = 7. This result was also obtained independently by Chul-hee Lee. Chapter 1 introduces some geometric notions including hyperbolic geometry, hyperbolic surfaces, and triangular surfaces, followed by Hurwitz’s automorphism theorem. Chapter 2 covers some basic concepts in spectral theory as well as some useful tools like the Selberg trace formula and a bit of representation theory. Chapter 3 focuses on the study of the Fricke–Macbeath surface and the proof of our main result using the techniques introduced in previous chapters. Finally, Chapter 4 discusses new methods for computing m1 which were used to show the existence of counterexamples to Colin de Verdière’s conjecture in joint work with Fortier Bourque, Gruda-Mediavilla, and Petri.

Le dioxyde de carbone est certes le gaz à effet de serre le plus connu du grand public. Mais il existe dans l'atmosphère terrestre des gaz peu médiatisés, présents en infimes quantités, dont la dangerosité potentielle n'est pas moindre. Ainsi, la molécule de SF6 a un pouvoir de forçage radiatif environ 24 000 fois supérieur à celui du CO2. D'origine anthropogénique, utilisé majoritairement comme isolant électrique, ce gaz voit depuis quelques années sa concentration dans l'atmosphère augmenter au taux alarmant de 8% par an. Ces données expliquent pourquoi l'hexafluorure de soufre est l'une des cibles de la lutte contre le réchauffement climatique dans le protocole de Kyoto et justifient le regain d'intérêt qui lui est porté dans le domaine académique. Dans cette thèse, nous avons réalisé une étude à la fois expérimentale et théorique de certaines transitions faiblement permises en spectroscopie Raman relatives à l'hexafluorure de soufre. Les transitions étudiées sont d'une part les harmoniques 2ν3 et 2ν5 de la molécule, toutes deux actives en Raman en raison d'anharmonicités électriques ou mécaniques dans la molécule isolée. Le fait que les vibrations ν3 et ν5 sont respectivement considérées comme les plus représentatives des modes d'étirement et de fléchissement de SF6 fût une motivation importante dans ce choix. Une étude expérimentale de la bande ν3 " induite par les collisions ", correspondant aux paires SF6−SF6, a d'autre part été menée et s'ensuit dans ce rapport. D'autres études relatives aux molécules isolées sont également présentées. Les expériences ont été réalisées avec un montage Raman artisanal de haute sensibilité et au moyen d'un protocole particulièrement strict, perpétuant un savoir-faire développé au sein d'une équipe reconnue pour sa capacité à détecter et analyser des flux lumineux de très faible intensité.

L’objectif de ce travail est de développer des outils numériques utilisables dans la détermination de manière exacte des propriétés modales des structures sandwichs viscoélastiques composites au vue de la conception des structures sandwichs viscoélastiques légères mais à haut pouvoir amortissant. Pour cela nous avons tout d’abord développé un outil générique implémenté en Matlab pour la détermination des propriétés modales en vibration libre des plaques sandwichs viscoélastiques dont les faces sont en stratifié de plusieurs couches orientées dans diverses directions. L’intérêt de cet outil, basé sur une formulation éléments finis, réside dans sa capacité à prendre en compte l’anisotropie des couches composites, la non linéarité matérielle de la couche viscoélastique traduit par diverses lois viscoélastiques dépendant de la fréquence ainsi que diverses conditions aux limites. La résolution du problème aux valeurs propres non linéaires complexes se fait par le couplage entre la technique d’homotopie, la méthode asymptotique numérique et la différentiation automatique. Ensuite pour permettre une étude continue des effets d’un paramètre de modélisation sur les propriétés modales des sandwichs viscoélastiques, nous avons proposé une méthode générique de résolution de problème résiduel non linéaire aux valeurs propres complexes possédant en plus de la dépendance en fréquence introduite par la couche viscoélastique du coeur, la dépendance du paramètre de modélisation qui décrit un intervalle d’étude bien spécifique. Cette résolution est basée sur la méthode asymptotique numérique, la différentiation automatique, la technique d’homotopie et la continuation et prend en compte diverses lois viscoélastiques. Nous proposons après cela, deux formulations distinctes pour étudier les effets, sur les propriétés amortissantes, de deux paramètres de modélisation qui sont importants dans la conception de sandwichs viscoélastiques à haut pouvoir amortissement. Le premier est l’orientation des fibres des composites dans la référence du sandwich et le second est l’épaisseur des couches qui lorsqu’elles sont bien définies permettent d’obtenir non seulement des structures sandwichs à haut pouvoir amortissant mais très légères. Les équations fortement non linéaires aux valeurs propres complexes obtenues dans ces formulations sont résolues par la nouvelle méthode de résolution d’équation résiduelle développée. Des comparaisons avec des résultats discrets sont faites ainsi que les temps de calcul pour montrer non seulement l’utilité de ces deux formulations mais également celle de la méthode de résolution d’équations résiduelles non linéaires aux valeurs propres complexes à double dépendance / Modeling and design of composite viscoelastic structures with high damping powerThe aim of this thesis is to develop numerical tools to determine accurately damping properties of composite sandwich structures for the design of lightweight viscoelastic sandwichs structures with high damping power. In a first step, we developed a generic tool implemented in Matlab for determining damping properties in free vibration of viscoelastic sandwich plates with laminate faces composed of multilayers. The advantage of this tool, which is based on a finite element formulation, is its ability to take into account the anisotropy of composite layers, the material non-linearity of the viscoelastic core induiced by the frequency-dependent viscoelastic laws and various boundary conditions . The nonlinear complex eigenvalues problem is solved by coupling homotopy technic, asymptotic numerical method and automatic differentiation. Then for the continuous study of a modeling parameter on damping properties of viscoelastic sandwichs, we proposed a generic method to solve the nonlinear residual complex eigenvalues problem which has in addition to the frequency dependence introduced by the viscoelastic core, a modeling parameter dependence that describes a very specific study interval. This resolution is based on asymptotic numerical method, automatic differentiation, homotopy technique and continuation technic and takes into account various viscoelastic laws. We propose after that, two separate formulations to study effects on the damping properties according to two modeling parameters which are important in the design of high viscoelastic sandwichs with high damping power. The first is laminate fibers orientation in the sandwich reference and the second is layers thickness which when they are well defined allow to obtain not only sandwich structures with high damping power but also very light. The highly nonlinear complex eigenvalues problems obtained in these formulations are solved by the new method of resolution of eigenvalue residual problem with two nonlinearity developed before. Comparisons with discrete results and computation time are made to show the usefulness of these two formulations and of the new method of solving nonlinear complex eigenvalues residual problem of two dependances