Les aubes des turbines à haute pression des réacteurs d'avion subissent des chargements complexes dans un environnement réactif. Prédire leur durée de vie peut nécessiter une approche en tolérance aux dommages, basée sur la prédiction de la propagation d'une fissure supposée. Mais cette approche est confrontée au comportement non linéaire sous des chargements à amplitudes variables et au coût énorme des calculs elasto-plastiques des structures 3D complexes sur des millions des cycles. Dans ce cadre, un modèle incrémental de fissuration a été proposé. Ce modèle est basé sur la plasticité comme mécanisme principal de propagation de fissure par fatigue pure. Cette modélisation passe par une réduction de modèle de type POD. La plasticité en pointe de la fissure est alors modélisée par un nombre réduit de variables non locales et des variables internes. Un ensemble d'hypothèses doit être respecté pour garantir la validité de cette modélisation. Pour décliner ce modèle dans le cas d'un matériau anisotrope représentatif du comportement des monocristaux, une première étude a été faite sur le cas d'une élasticité cubique avec de la plasticité de Von-Mises. Une stratégie a été proposée pour identifier un modèle matériau basé sur les facteurs d'intensité non locaux. Cette stratégie comporte une détermination de la fonction critère basée sur les solutions élastiques en anisotrope. L'étude des directions d'écoulement plastique avec les variables non locales montre une forte dépendance à l'anisotropie élastique du modèle même avec une plasticité associée de Von-Mises. La stratégie comporte également une identification des variables internes.Dans la deuxième partie, le problème d'une fissure avec un modèle de plasticité cristalline a été traité. L'activation de différents systèmes de glissement a été alors prise en compte dans la modélisation. Finalement, différentes méthodologies ont été explorées en vue de transposer le modèle local de plasticité cristalline à l'échelle non locale de la région en pointe de la fissure. / The fatigue life prediction of high pressure turbine blades may require a damage tolerance approach based on the study of possible crack propagation. The nonlinear behavior of the material under complex nonproportional loadings and the high cost of running long and expensive elastic-plastic FE computations on complex 3D structures over millions of cycles are some major issues that may encounter this type of approach.Within this context, an incremental model was proposed based on plasticity as a main mechanism for fatigue crack growth.A model reduction strategy using the Proper Orthogonal Decomposition (POD) was used to reduce the cost of FEA. Based on a set of hypotheses, the number of the degrees of freedom of the problem is reduced drastically. The plasticity at the crack tip is finally described by a set of empirical equations of few nonlocal variables and some internal variables.In order to apply this modeling strategy to the case of anisotropic materials that represent the behavior of single crystals, a first study was done with cubic elasticity and a Von-Mises plasticity. Elastic and plastic reference fields, required to reduce the model, were determined. Then, a material model of the near crack tip region was proposed based on nonlocal intensity factors. A yield criterion function was proposed based on Hoenig's asymptotic solutions for anisotropic materials. The study of plastic flow directions with the nonlocal variables of the model shows a strong dependency on the cubic elasticity. A strategy to identify internal variables is proposed as well. In the second part, a crystal plasticity model was implemented. The activation of different slip systems was taken into account in the model reduction strategy. A kinematic basis was constructed for each slip system. Finally, a strategy was proposed to transpose the local crystal plasticity model to the nonlocal scale of the crack.

In order to understand biochemical networks we need to know not only how their parts work but also how they interact with each other. The goal of systems biology is to look at biological systems as a whole to understand how interactions of the parts can give rise to complex dynamics. In order to do this efficiently, new techniques have to be developed. This work shows how tools from mathematics are suitable to study problems in systems biology such as modeling, dynamics prediction, reverse engineering and many others. The advantage of using mathematical tools is that there is a large number of theory, algorithms and software available. This work focuses on how algebra can contribute to answer questions arising from systems biology. / Ph. D.

Science and engineering problems frequently require solving a sequence of single linear systems or a sequence of dual linear systems. We develop algorithms that recycle Krylov subspaces and preconditioners from one system (or pair of systems) in the sequence to the next, leading to efficient solutions. Besides the benefit of only having to store few Lanczos vectors, using BiConjugate Gradients (BiCG) to solve dual linear systems may have application-specific advantages. For example, using BiCG to solve the dual linear systems arising in interpolatory model reduction provides a backward error formulation in the model reduction framework. Using BiCG to evaluate bilinear forms -- for example, in the variational Monte Carlo (VMC) algorithm for electronic structure calculations -- leads to a quadratic error bound. Since one of our focus areas is sequences of dual linear systems, we introduce recycling BiCG, a BiCG method that recycles two Krylov subspaces from one pair of dual linear systems to the next pair. The derivation of recycling BiCG also builds the foundation for developing recycling variants of other bi-Lanczos based methods like CGS, BiCGSTAB, BiCGSTAB2, BiCGSTAB(l), QMR, and TFQMR. We develop a generalized bi-Lanczos algorithm, where the two matrices of the bi-Lanczos procedure are not each other's conjugate transpose but satisfy this relation over the generated Krylov subspaces. This is sufficient for a short term recurrence. Next, we derive an augmented bi-Lanczos algorithm with recycling and show that this algorithm is a special case of generalized bi-Lanczos. The Petrov-Galerkin approximation that includes recycling in the iteration leads to modified two-term recurrences for the solution and residual updates. We generalize and extend the framework of our recycling BiCG to CGS, BiCGSTAB and BiCGSTAB2. We perform extensive numerical experiments and analyze the generated recycle space. We test all of our recycling algorithms on a discretized partial differential equation (PDE) of convection-diffusion type. This PDE problem provides well-known test cases that are easy to analyze further. We use recycling BiCG in the Iterative Rational Krylov Algorithm (IRKA) for interpolatory model reduction and in the VMC algorithm. For a model reduction problem, we show up to 70% savings in iterations, and we also demonstrate that solving the problem without recycling leads to (about) a 50% increase in runtime. Experiments with recycling BiCG for VMC gives promising results. We also present an algorithm that recycles preconditioners, leading to a dramatic reduction in the cost of VMC for large(r) systems. The main cost of the VMC method is in constructing a sequence of Slater matrices and computing the ratios of determinants for successive Slater matrices. Recent work has improved the scaling of constructing Slater matrices for insulators, so that the cost of constructing Slater matrices in these systems is now linear in the number of particles. However, the cost of computing determinant ratios remains cubic in the number of particles. With the long term aim of simulating much larger systems, we improve the scaling of computing determinant ratios in the VMC method for simulating insulators by using preconditioned iterative solvers. The main contribution here is the development of a method to efficiently compute for the Slater matrices a sequence of preconditioners that make the iterative solver converge rapidly. This involves cheap preconditioner updates, an effective reordering strategy, and a cheap method to monitor instability of ILUTP preconditioners. Using the resulting preconditioned iterative solvers to compute determinant ratios of consecutive Slater matrices reduces the scaling of the VMC algorithm from O(n^3) per sweep to roughly O(n^2), where n is the number of particles, and a sweep is a sequence of n steps, each attempting to move a distinct particle. We demonstrate experimentally that we can achieve the improved scaling without increasing statistical errors. / Ph. D.

For repeated transient solutions of geometrically nonlinear structures the numerical effort often poses a major obstacle. Thus, the introduction of a reduced order model, which takes the nonlinear effects into account and accelerates the calculations considerably, is often necessary.This work yields a method that allows for rapid, accurate and parameterisable solutions by means of a reduced model of the original structure. The structure is discretised and its dynamic equilibrium described by a matrix equation. The projection on a reduced basis is introduced to obtain the reduced model. A comprehensive numerical study on several common reduced bases shows that the simple introduction of a constant basis is not sufficient to account for the nonlinear behaviour. Three requirements for an rapid, accurate and parameterisable solution are derived. The solution algorithm has to take into account the nonlinear evolution of the solution, the solution has to be independent of the nonlinear finite element terms and the basis has to be adapted to external parameters.Three approaches are provided, each responding to one requirement. These approaches are assembled to the integrated method. The approaches are the update and augmentation of the basis, the polynomial formulation of the nonlinear terms and the interpolation of the basis. A Newmark-type time-marching algorithm provides the frame of the integrated method. The application of the integrated method on test-cases with geometrically nonlinear finite elements confirms that this method leads to the initial aim of a rapid, accurate and parameterisable transient solution.

Motivé par l'étude de tissus biologiques, ce travail contribue aux développements d'outils numériques permettant de prédire la réponse mécanique de matériaux hétérogènes non linéaires dans lesquels les énergies d'interfaces deviennent prépondérantes. Ainsi, une méthode d'homogénéisation multi échelle combinée à une technique de réduction de modèle basée sur la décomposition orthogonale aux valeurs propres est proposée dans un cadre thermique et hyperélastique. Les énergies d'interfaces entre les différentes phases des composites sont décrites par un modèle d'interface cohérent et prises en compte numériquement par une approche liant la méthode des éléments finis étendus et la méthode level-set. Une étude de l'étalement d'une cellule vivante entre deux lamelles fixes est ensuite réalisée. Les deux modèles utilisés pour les simulations montrent que l'assemblage cortex d'actine-membrane plasmique ne joue qu'un rôle minime dans la réponse mécanique cellulaire / Motivated by the study of biological tissues, this work contributes to developing numerical tools to predict the mechanical response of nonlinear heterogeneous materials in which the energies of interfaces can no longer be ignored. First, a computational homogenization strategy combined with a model reduction technique based on the proper orthogonal decomposition is implemented in the cases of large elastic deformations and highly nonlinear conduction. The interfaces between the different phases of a composite are described by means of a coherent interface model and taken into account numerically by an extended finite element method in tandem with a level-set technique. Finally, experimental results of single cell spreading between two fixed parallel microplates are exploited through finite element modelling. Our two models show that the bilayer membrane and the actin cortex do not play a significant role in the cell mechanical response

<p>Disertacijom su opisane postojeće tehnike redukcije dinamičkih modela koje se koriste u teoriji upravljanja i postojeće tehnike za redukciju dinamičkih modela i ekvivalentiranje elektroenergetskih sistema. Predložen je nov pristup na fizici problema zasnovanoj redukciji dinamičkog modela elektroenergetskog sistema korišćenjem teorije balansnih realizacija. Takođe se predlaže korišćenje aproksimativnih bisimulacionih relacija za redukciju dinamičkih modela elektroenergetskog sistema. Postojeće tehnike i predloženi pristupi i algoritmi su primenjeni za redukciju dinamičkih modela dva razmatrana test sistema.</p> / <p>Dissertation describes the existing dynamic model reduction techniques used in control theory and existing techniques that are used for the reduction (equivalencing) of power system dynamic models. A new approach to physics-based reduction of power system dynamic model based on the balanced realization theory is proposed. Use of approximate bisimulation relations for reduction of power system dynamic models is also proposed. Existing techniques and proposed approaches and algorithms are applied to reduce the dynamic models of two considered test systems.<br />&nbsp;</p>

Cette thèse s’intéresse à la faisabilité de la mise en place de modèles d’ordre réduit pour l’étude des vibrations sous écoulement au sein de faisceaux tubulaires d’échangeurs de chaleur. Ces problématiques sont cruciales car les systèmes étudiés sont des éléments majeurs des centrales nucléaires civiles et des chaufferies embarquées dans les sous-marins.Après avoir rappelé le fonctionnement et les risques vibratoires existants au sein des échangeurs de chaleur, des calculs complets d’écoulement et de vibrations sous écoulement ont été effectués, d’abord pour un tube seul en milieu infini, puis pour un faisceau de tubes. Ces calculs ont été menés avec l’outil CFD Code_Saturne. La méthode de réduction de modèle POD (Proper Orthogonal De-composition) a été appliquée au cas des écoulements avec la structure laissée fixe.Les résultats obtenus montrent l’efficacité de la méthode pour ces configurations, moyennant l’introduction de méthodes de stabilisation pour l’écoulement au sein du faisceau. La méthode POD-multiphasique, permettant d’adapter la méthode POD à l’interaction fluide-structure, a ensuite été appliquée. Les grands déplacements d’un cylindre seul dans la zone d’accrochage (lock-in) ont été correctement reproduits par cette méthode de réduction de modèle. De même, on montre que les grands déplacements d’un cylindre en milieu confiné dans un faisceau de tubes sont fidèlement reconstruits.Enfin, l’extension de l’utilisation de la réduction de modèle aux études d’évolution paramétrique a été testée. Nous avons d’abord utilisé la technique considérant une base POD unique pour reproduire des écoulements à divers nombres de Reynolds autour d’un cylindre seul. Les résultats confirment la prédictivité bornée à une gamme de paramètres de cette méthode. Enfin, l’interpolation de bases POD pré-calculées pour une famille de paramètres donnés, utilisant les variétés de Grassmann et permettant de générer de nouvelles bases POD, a été testée sur des cas modèles. / The objective of this thesis is to study the ability of model reduction for investigations of flow-induced vibrations in heat exchangers tube bundle systems.These mechanisms are a cause of major concern because heat exchangers are key elements of nuclear power plants and on-board stoke-holds.In a first part, we give a recall on heat exchangers functioning and on vi-bratory problems to which they are prone. Then, complete calculations leaded with the CFD numerical code Code_Saturne are carried out, first for the flow around a single circular cylinder (fixed then elastically mounted) and then for the case of a tube bundle system submitted to cross-flow. Reduced-order method POD is ap-plied to the flow resolution with fixed structures. The obtained results show the efficiency of this technique for such configurations, using stabilization methods for the dynamical system resolution in the tube-bundle case.Multiphase-POD, which is a method enabling the adaptation of POD to fluid-structure interactions, is applied. Large displacements of a single cylinder elastically mounted under cross-flow, corresponding to the lock-in phenomenon,are well reproduced with this reduction technique. In the same way, large displace-ments of a confined moving tube in a bundle are shown to be faithfully recon-structed.Finally, the use of model reduction is extended to parametric studies. First,we propose to use the method which consists in projecting Navier-Stokes equations for several values of the Reynolds number on to a unique POD basis. The resultsobtained confirm the fact that POD predictability is limited to a range of parameter values. Then, a basis interpolation method, constructed using Grassmann mani-folds and allowing the construction of a POD basis from other pre-calculated basis,is applied to basic cases.

Les freins aéronautiques sont soumis à des instabilités vibratoires induites par le frottement. Il en résulte des vibrations qui présentent un risque pour la structure du frein et de l’atterrisseur et posent des problèmes d’intégration. Safran Landing Systems doit donc répondre à des spécifications avionneurs strictes sur les niveaux des vibrations générées par son équipement. Le respect de ces spécifications est actuellement contrôlé par la réalisation d’essais de freinage longs et coûteux. L’objectif de ces travaux de recherche est de reproduire numériquement ces phénomènes vibratoires via des outils intégrables au processus de conception d’un frein. Le crissement de frein, bien qu’il soit l’objet de recherches depuis le début du XXe siècle, demeure un phénomène assez mal compris, notamment dans l’aéronautique. Des vibrations instables apparaissent régulièrement sur l’ensemble de la plage fréquentielle 0-2 kHz. Au cours de la dernière décennie, une instabilité vibratoire vers 200 Hz dénommée whirl 2 s’est manifestée de manière récurrente et souvent critique sur la plupart des nouveaux freins développés. On cherche donc à mettre en place une méthode permettant de simuler l’apparition et les amplitudes des instabilités vibratoires, notamment du mode de whirl 2. Dans une première partie, on présente des analyses d’essais vibratoires réalisés en conditions opérationnelles et expérimentales. On décrit ensuite la modélisation par la méthode des éléments finis du frein instable au sens de Lyapunov. La stabilité du système linéarisé est étudiée et on montre une corrélation en fréquence et déformée entre le modèle et les essais. Ce modèle éléments finis est trop volumineux en l’état pour permettre la simulation d’amplitudes de vibrations non linéaires. On propose donc dans une seconde partie deux méthodes de réduction adaptées à l’architecture complexe d’un système de freinage aéronautique et permettant la prise en compte du frottement. La première est une méthode semi-analytique qui se révèle très performante jusqu’à 500 Hz. La seconde méthode de réduction mise en oeuvre est la double synthèse modale. Elle est implémentée dans sa version classique, puis une amélioration est proposée avec succès : la double synthèse modale complexe. La troisième partie est consacrée à l’étude de la dynamique non linéaire du whirl 2 par la réalisation d’intégrations temporelles. La simulation des amplitudes de vibration nécessite la prise en compte réaliste du comportement non linéaire du frein. Or, on fait d’abord le constat que, contrairement à une hypothèse communément admise, les non-linéarités de contact situées aux interfaces frottantes ne suffisent pas à expliquer à elles seules la saturation des amplitudes vibratoires constatée expérimentalement. La recherche des phénomènes physiques non linéaires influents nous amène a considérer l’interaction de la structure vibrante avec le circuit hydraulique de commande du frein. La modélisation du couplage hydrodynamique fournit alors des éléments de compréhension inédits et permet de formuler des règles de conception. Enfin on étudie l’impact du frottement sec dans les contacts périphériques des disques de freinage avec la structure. Ce phénomène, jusque là négligé, apparaît largement prépondérant. Des études d’influences, présentant une bonne corrélation avec les essais, permettent de mettre en évidence de manière robuste l’influence du design et des scénarios de freinage sur les amplitudes vibratoires. / These vibrations are a threat for the brake and landing-gear structural integrity and represent an issue in terms of integration. Thus Safran Landing Systems has to comply with aircraft manufacturers’ strict requirements on the vibration amplitude its product is likely to generate. Compliance to these requirements is assessed by long and costly braking test campaigns. The objective of the research presented here is to reproduce by simulation the brake dynamic instabilities with numerical tools that could be integrated in the design process. Brake squeal has been a research topic since the early XXth century. However it remains a rather ill-understood phenomenon, especially in aeronautics. Unstable vibrations regularly appear on the whole 0-2kHz frequency spectrum. In the last decade, an instability located around 200 Hz called whirl 2 persistently appeared on the newly developed wheel and brake assemblies, sometimes exhibiting critical vibration amplitudes. Consequently, Safran Landing Systems wishes to develop numerical tools able to simulate both the occurrence and the amplitudes associated with friction-induced instabilities, especially with the whirl 2 mode. In the first part of this report, an experimental analysis of the brake is conducted, on both laboratory and in operational set-ups. The modelling of the wheel and brake assembly using the finite element method is then described. The system stability in a Lyapunov’s sense is studied and shows good correlation in both frequencies and mode shapes with the experiments. This finite element model is too big to be used to perform the transient simulation of the nonlinear amplitudes. In the second part, two reduction methods, tailored to the complex aircraft brakes architectures, are thus presented. The first method is a semi-analytical. It shows excellent performances up to 500 Hz. The second reduction method is the double modal synthesis, implemented under its classical version. It is then successfully improved and called "complex double modal synthesis". The third part is dedicated to the study of the nonlinear dynamics of the whirl 2 through transient analyses. The nonlinear amplitudes simulation requires taking into account the relevant nonlinear brake behavior. However, it is first observed that, contrary to a commonly accepted hypothesis, the contact nonlinearities located at the friction interfaces cannot single-handedly account for the vibration amplitudes saturation observed in the tests. The need to identify the relevant physical phenomena leads then to consider the interaction between the squealing brake structure and its hydraulic command circuit. The modelling of the hydro-mechanical coupling provides an unprecedented insight and allows to prescribe design rules. Finally, we study the impact of dry friction in the peripheral contacts between the braking discs and the structure. This phenomenon, neglected until now, appears to have a major influence. Sensitivity studies exhibit a good correlation with tests, allowing to highlight, in a robust manner, the impact of brake design and braking scenarii on the nonlinear vibration amplitudes.

Dans le dimensionnement des lanceurs spatiaux, la maîtrise de l'amortissement est une problématique majeure. Faute d'essais sur structure réelle très couteux avant la phase finale de qualification, la modélisation de l'amortissement peut conduire à un sur-dimensionnement de la structure alors que le but recherché est de diminuer le coût du lancement d'une fusée tout en garantissant le confort vibratoire de la charge utile.Nos contributions sont les suivantes. Premièrement, une méthode de prédiction par le calcul des niveaux vibratoires dans les structures de lanceurs en utilisant une stratégie d'essais virtuels qui permet de prédire les amortissements en basses fréquences, est proposée. Cette méthode est basée sur l'utilisation de méta-modèles construits à partir de plans d'expériences numériques à l'aide de modèles détaillés des liaisons. Ces méta-modèles peuvent être obtenus grâce à des calculs spécifiques utilisant une résolution 3D par éléments finis avec prise en compte du contact. En utilisant ces méta-modèles, l'amortissement modal dans un cycle de vibration peut être calculé comme étant le ratio entre l'énergie dissipée et l'énergie de déformation. L'approche utilisée donne une approximation précise et peu coûteuse de la solution. Le calcul non-linéaire global qui est inaccessible pour les structures complexes est rendu accessible en utilisant l'approche virtuelle basées sur les abaques.Deuxièmement, une validation des essais virtuels sur la structure du lanceur Ariane 5 a été élaborée en tenant compte des liaisons boulonnées entre les étages afin d'illustrer l'approche proposée. Lorsque la matrice d'amortissement généralisé n'est pas diagonale (car des dissipations localisées), ces méthodes modales ne permettent pas de calculer ou d'estimer les termes d'amortissement généralisé extra-diagonaux. La problématique posée est alors la quantification de l'erreur commise lorsque l'on néglige les termes extra-diagonaux dans le calcul des niveaux vibratoires ; avec un bon ratio précision / coût de calcul.Troisièmement, la validité de l'hypothèse de diagonalité de la matrice d'amortissement généralisée a été examinée et une méthode très peu coûteuse de quantification a posteriori de l'erreur d'estimation de l'amortissement modal par la méthodes des perturbations a été proposée.Finalement, la dernière contribution de cette thèse est la proposition d'un outil d'aide à la décision qui permet de quantifier l'impact des méconnaissances sur l'amortissement dans les liaisons sur le comportement global des lanceurs via l'utilisation de la méthode info-gap. / In the dimensioning of space launchers, controlling depreciation is a major problem. In the absence of very expensive real structural tests before the final qualification phase, damping modeling can lead to over-sizing of the structure while the aim is to reduce the cost of launching a rocket while guaranteeing the vibratory comfort of the payload.[...]

Cette thèse porte sur le contrôle acoustique actif (ANC) dans une cavité. L’objectif est d’atténuer l’effet d’une onde sonore perturbatrice en des points ou dans un volume. Ceci est réalisé à l’aide d’un contre-bruit généré, par exemple, par un haut-parleur. Cette étude requiert l’utilisation de modèles dynamiques rendant compte de l’évolution des pressions aux points d’intérêt en fonction des bruits exogènes. Ce modèle peut être obtenu par une identification fréquentielle des réponses point-à-point ou en utilisant le modèle physique sous jacent (équation des ondes). Dans ce dernier cas, la recherche d'un modèle de dimension finie est souvent un préalable à l’étude conceptuelle d'un système d’ANC. Les contributions de cette thèse portent donc sur l’élaboration de différents modèles simplifiés paramétrés par la position pour les systèmes acoustiques et sur la conception de lois de commande pour l’ANC. Le premier volet de la thèse est dédié à l’élaboration de différents modèles simplifiés de système de propagation acoustique au sein d’une cavité. Pour cela, les simplifications envisagées peuvent être de nature spatiale autant que fréquentielle. Nous montrons notamment qu'il est possible, sous certaines conditions, d’approximer le système 3D par un système 1D. Ceci a été mis en évidence expérimentalement sur le banc d’essai LS2NBox. Le second volet porte sur la conception de lois de commande. En premier lieu, les stratégies de commandes couramment utilisées pour l’ANC sont comparées. L'effet dela commande multi-objectif H en différents points voisins des points d'atténuation est analysé. La possibilité d’une annulation parfaitedu bruit en un point est aussi discutée. / This thesis deals with active noise control (ANC) in a cavity. The aim is to mitigate the effect of a disturbing sound wave at some points or in a volume. This is achieved using ananti-noise generated, for example, by a loudspeaker. This study requires the use of dynamic models that report changes in pressure at points of interest in response to exogenous noises. Such models can be obtained by frequency identification of point-to-point responses or by using the underlying physical model (wave equation). In the latter case, the search for a low-complexity model (finite dimensional model) is often a prerequisite for the conceptual study of an active control system. The contributions of this thesis concern the development of different simplified models parameterized by the spatial position for acoustic systems, and the design of control laws for noise attenuation. The first part of the thesis is dedicated to the development of various simplified models of acoustic propagation system within a cavity. For that, the simplifications envisaged can be of spatial nature as much as frequential. We show in particular that it is possible, under certain conditions, to approximate the 3D system by a 1D system. This has been demonstrated experimentally on the prototype system, LS2NBox. The second part of the thesis deals with the design of control laws. First, the control strategies commonly used for ANC are compared. The effect of multi-objective H control at different spatial positions close to the attenuation points is analyzed. The possibility of perfect noise cancellation at one point is also discussed.