In this work we calculate the binding energies and root-mean-square radii for A−body nuclear bound state systems, where A ≥ 3. To study three−body systems, we employ the three−dimensional differential Faddeev equations with nucleon-nucleon semi-realistic potentials. The equations are solved numerically. For this purpose, the equations are transformed into an eigenvalue equation via the orthogonal collocation procedure using triquintic Hermite splines. The resulting eigenvalue equation is solved using the Restarted Arnoldi Algorithm. Ground state binding energies of the 3H nucleus are determined. For A > 3, the Potential Harmonic Expansion Method is employed. Using this method, the Schr¨odinger equation is transformed into coupled Faddeev-like equations. The Faddeevlike amplitudes are expanded on the potential harmonic basis. To transform the resulting coupled differential equations into an eigenvalue equation, we employ again the orthogonal collocation procedure followed by the Gauss-Jacobi quadrature. The corresponding eigenvalue equation is solved using the Renormalized Numerov Method to obtain ground state binding energies and root-mean-square radii of closed shell nuclei 4He, 8Be, 12C, 16O and 40Ca. / Physics / M. Sc. (Physics)

The developments of advanced non-linear control strategies have attracted a considerable research interests over the past decades especially in process control. Rather than an absolute reliance on mathematical models of process plants which often brings discrepancies especially owing to design errors and equipment degradation, non-linear models are however required because they provide improved prediction capabilities but they are very difficult to derive. In addition, the derivation of the global optimal solution gets more difficult especially when multivariable and non-linear systems are involved. Hence, this research investigates soft computing techniques for the implementation of a novel real time constrained non-linear model predictive controller (NMPC). The time-frequency localisation characteristics of wavelet neural network (WNN) were utilised for the non-linear models design using system identification approach from experimental data and improve upon the conventional artificial neural network (ANN) which is prone to low convergence rate and the difficulties in locating the global minimum point during training process. Salient features of particle swarm optimisation and a genetic algorithm (GA) were combined to optimise the network weights. Real time optimisation occurring at every sampling instant is achieved using a GA to deliver results both in simulations and real time implementation on coupled tank systems with further extension to a complex quadruple tank process in simulations. The results show the superiority of the novel WNN-NMPC approach in terms of the average controller energy and mean squared error over the conventional ANN-NMPC strategies and PID control strategy for both SISO and MIMO systems.

Large amplitude internal solitary waves in the coastal ocean are commonly modelled with the Korteweg-de Vries (KdV) equation or a closely related evolution equation. The characteristic feature of these models is the solitary wave solution, and it is well documented that these provide the basic paradigm for the interpretation of oceanic observations. However, often internal waves in the ocean survive for several inertial periods, and in that case, the KdV equation is supplemented with a linear non-local term representing the effects of background rotation, commonly called the Ostrovsky equation. This equation does not support solitary wave solutions, and instead a solitary-like initial condition collapses due to radiation of inertia-gravity waves, with instead the long-time outcome typically being an unsteady nonlinear wave packet. The KdV equation and the Ostrovsky equation are formulated on the assumption that only a single vertical mode is used. In this thesis we consider the situation when two vertical modes are used, due to a near-resonance between their respective linear long wave phase speeds. This phenomenon can be described by a pair of coupled Ostrovsky equations, which is derived asymptotically from the full set of Euler equations and solved numerically using a pseudo-spectral method. The derivation of a system of coupled Ostrovsky equations is an important extension of coupled KdV equations on the one hand, and a single Ostrovsky equation on the other hand. The analytic structure and dynamical behaviour of the system have been elucidated in two main cases. The first case is when there is no background shear flow, while the second case is when the background state contains current shear, and both cases lead to new solution types with rich dynamical behaviour. We demonstrate that solitary-like initial conditions typically collapse into two unsteady nonlinear wave packets, propagating with distinct speeds corresponding to the extremum value in the group velocities. However, a background shear flow allows for several types of dynamical behaviour, supporting both unsteady and steady nonlinear wave packets, propagating with the speeds which can be predicted from the linear dispersion relation. In addition, in some cases secondary wave packets are formed associated with certain resonances which also can be identified from the linear dispersion relation. Finally, as a by-product of this study it was shown that a background shear flow can lead to the anomalous version of the single Ostrovsky equation, which supports a steady wave packet.

There exists substantial applications motivating the study of nonlinear longitudinal wave propagation in layered (or laminated) elastic waveguides, in particular within areas related to non-destructive testing, where there is a demand to understand, reinforce, and improve deformation properties of such structures. It has been shown [76] that long longitudinal waves in such structures can be accurately modelled by coupled regularised Boussinesq (cRB) equations, provided the bonding between layers is sufficiently soft. The work in this thesis firstly examines the initial-value problem (IVP) for the system of cRB equations in [76] on the infinite line, for localised or sufficiently rapidly decaying initial conditions. Using asymptotic multiple-scales expansions, a nonsecular weakly nonlinear solution of the IVP is constructed, up to the accuracy of the problem formulation. The asymptotic theory is supported with numerical simulations of the cRB equations. The weakly nonlinear solution for the equivalent IVP for a single regularised Boussinesq equation is then constructed; constituting an extension of the classical d'Alembert's formula for the leading order wave equation. The initial conditions are also extended to allow one to separately specify an O(1) and O(ε) part. Large classes of solutions are derived and several particular examples are explicitly analysed with numerical simulations. The weakly nonlinear solution is then improved by considering the IVP for a single regularised Boussinesq-type equation, in order to further develop the higher order terms in the solution. More specifically, it enables one to now correctly specify the higher order term's time dependence. Numerical simulations of the IVP are compared with several examples to justify the improvement of the solution. Finally an asymptotic procedure is developed to describe the class of radiating solitary wave solutions which exist as solutions to cRB equations under particular regimes of the parameters. The validity of the analytical solution is examined with numerical simulations of the cRB equations. Numerical simulations throughout this work are derived and implemented via developments of several finite difference schemes and pseudo-spectral methods, explained in detail in the appendices.

The performance of commercial Radio Frequency Identification (RFID) tags is primarily limited by present techniques used for tag antenna design. Currently, industry techniques rely on identifying the RFID tag application (books, clothing, etc.) and then building antenna prototypes of different configurations in order to satisfy minimum read range requirements. However, these techniques inherently lack an electromagnetic basis and are unable to provide a low cost solution to the tag antenna design process. RFID tag performance characteristics (read-range, chip-antenna impedance matching, surrounding environment) can be very complex, and a thorough understanding of the RFID tag antenna design may be gained through an electromagnetic approach in order to reduce the tag antenna size and the overall cost of the RFID system. The research presented in this thesis addresses RFID tag antenna design process for passive RFID tags. With the growing number of applications (inventory, supply-chain, pharmaceuticals, etc), the proposed RFID antenna design process demonstrates procedures to design tag antennas for such applications. Electrical/geometrical properties of the antennas designed were investigated with the help of computer electromagnetic simulations in order to achieve optimal tag performance criteria such as read range, chip-impedance matching, antenna efficiency, etc. Experimental results were performed on the proposed antenna designs to compliment computer simulations and analytical modelling.

Récemment plusieurs récepteurs couplés aux protéines G (RCPGs) ont été caractérisés au niveau des membranes intracellulaires, dont la membrane nucléaire. Notre objectif était de déterminer si les sous-types de récepteurs β-adrénergiques (βAR) et leurs machineries de signalisation étaient fonctionnels et localisés à la membrane nucléaire des cardiomyocytes. Nous avons démontré la présence des β1AR et β3AR, mais pas du β2AR à la membrane nucléaire de myocytes ventriculaires adultes par immunobuvardage, par microscopie confocale, et par des essais fonctionnels. De plus, certains partenaires de signalisation comme les protéines GαS, Gαi, l’adénylate cyclase II, et V/VI y étaient également localisés. Les sous-types de βAR nucléaires étaient fonctionnels puisqu'ils pouvaient lier leurs ligands et activer leurs effecteurs. En utilisant des noyaux isolés, nous avons observé que l'agoniste non-sélectif isoprotérénol (ISO), et que le BRL37344, un ligand sélectif du β3AR, stimulaient l'initiation de la synthèse de l’ARN, contrairement à l'agoniste sélectif du β1AR, le xamotérol. Cette synthèse était abolie par la toxine pertussique (PTX). Cependant, la stimulation des récepteurs nucléaires de type B de l’endothéline (ETB) causaient une réduction de l'initiation de la synthèse d’ARN. Les voies de signalisations impliquées dans la régulation de la synthèse d’ARN par les RCPGs ont ensuite été étudiées en utilisant des noyaux isolés stimulés par des agonistes en présence ou absence de différents inhibiteurs des voies MAP Kinases (proteines kinases activées par mitogènes) et de la voie PI3K/PKB. Les protéines impliquées dans les voies de signalisation de p38, JNK, ERK MAP Kinase et PKB étaient présents dans les noyaux isolés. L'inhibition de PKB par la triciribine, inhibait la synthèse d’ARN. Nous avons ensuite pu mettre en évidence par qPCR que la stimulation par l’ISO entrainait une augmentation du niveau d'ARNr 18S ainsi qu’une diminution de l'expression d’ARNm de NFκB. En contraste, l’ET-1 n’avait aucun effet sur le niveau d’expression de l’ARNr 18S. Nous avons ensuite montré que la stimulation par l’ISO réduisait l’expression de plusieurs gènes impliqués dans l'activation de NFκB, tandis que l’inhibition de ERK1/2 et PKB renversait cet effet. Un microarray global nous a ensuite permis de démontrer que les βARs et les ETRs nucléaires régulaient un grand nombre de gènes distincts. Finalement, les βARs et ETRs nucléaires augmentaient aussi une production de NO de noyaux isolés, ce qui pouvait être inhibée par le LNAME. Ces résultats ont été confirmés dans des cardiomyocytes intacts en utilisant des analogues cagés et perméables d’ISO et de l'ET-1: l'augmentation de NO nucléaire détectée par DAF2-DA, causée par l'ET-1 et l'ISO, pouvait être prévenue par le LNAME. Finalement, l’augmentation de l’initiation de la transcription induite par l'ISO était aussi bloquée par le L-NAME ou par un inbitheur de PKG, le KT5823, suggérant que la voie NO-GC-PKG est impliquée dans la régulation de la transcription par les βAR. En conclusion, les βARs et les ETRs nucléaires utilisent des voies de signalisation différentes et exercent ainsi des effets distincts sur l’expression des gènes cardiaques. Ils représentent donc une avenue intéressante pour le développement de drogues pharmacologiques. / Recently several G protein-coupled receptors (GPCRs) have been shown to localize to intracellular membranes, in particular the nuclear membrane. As such, we sought to determine if the β-adrenergic receptor (βAR) subtypes and their associated signalling machinery are functionally localized to nuclear membranes. We demonstrated the presence of β1AR and the β3AR, but not the β2AR, in adult ventricular myocyte nuclei by western blotting, confocal microscopy and functional assays. Downstream signalling partners such as GαS, Gαi and adenylyl cyclase II and V/VI were also present. Nuclear-localized βARs were functional with respect to ligand binding and effector activation. In isolated nuclei, the non-selective βAR agonist isoproterenol (ISO) and the β3AR-selective ligand BRL37344, but not the β1AR-selective xamoterol, stimulated transcription initiation in a pertussis toxin (PTX)-sensitive manner. In contrast, stimulation of type B endothelin receptors (ETB), another GPCR family shown to be present on the nuclear membrane, decreased de novo RNA synthesis. To investigate the signalling pathway(s) involved in GPCR-mediated regulation of RNA synthesis, nuclei were isolated from intact adult rat hearts and treated with receptor agonists in the presence or absence of inhibitors of the PI3K/PKB and mitogen-activated protein kinase (MAPK) pathways. Components of p38, JNK, and ERK1/2 MAPK cascades as well as PKB were detected in nuclear preparations. Inhibition of PKB with triciribine converted the activation of the βAR from stimulatory to inhibitory with regards to transcription initiation. Analysis by qPCR indicated isoproterenol treatment increased 18S rRNA but decreased NFκB mRNA. In contrast, ET-1 had no effect on 18S rRNA expression. Further investigation using pathway-specific PCR arrays revealed that isoproterenol treatment also reduced the expression of several other genes involved in the activation of NFκB and that ERK1/2 and PKB inhibitors attenuated this effect. Subsequent genome-wide microarray analysis has revealed that nuclear βAR and ETB regulated a host of genes in an overlapping but distinct manner. Moreover, both ET-1 and ISO produced an L-NAME-sensitive increase in NO production in isolated cardiac nuclei. These observations were confirmed in intact cardiomyocytes using novel caged analogues of ISO and ET-1 and the cell-permable NO-sensitive fluorescent dye, DAF-2 DA. Briefly, both ET-1 and isoproterenol increased NO production, and this increase was prevented upon preincubation with L-NAME. Moreover, the ability of isoproterenol to increase transcription initiation in isolated nuclei was blocked by L-NAME or the PKG inhibitor KT5823, indicating the NO-GC-PKG pathway is involved in the regulation of gene expression by nuclear βARs. Hence, we have shown that βARs and ETRs in the nuclear membrane activate distinct signalling pathways, resulting in different effects on gene transcription and thus represent potentially important targets for drug development.

Les maladies inflammatoires de l’intestin (MIIs) comprennent la colite ulcéreuse (CU) et la maladie de Crohn (MC). Elles résultent en l’inflammation chronique du tractus gastro-intestinal. Des études d’association pan-génomiques ont associé le locus 14q31, incluant les gènes galactosylceramidase (GALC) et G protein-coupled receptor 65 (GPR65), à ces pathologies. Nous avons déterminé que le variant le plus associé aux MIIs de cette région (rs8005161; p=2,35x10-14) est parfaitement corrélé (r2=1) à un variant codant dans le gène GPR65 (rs3742704: Ile231Leu). GPR65 code pour un récepteur couplé à des protéines G (RCPG) senseur de pH. Nous avons observé que GPR65 est exprimé dans les tissus lymphoïdes et mucoïdes en plus des lignées cellulaires immunitaires et des cellules immunes primaires humaines. Son expression augmente significativement dans les biopsies inflammées de patients atteints de la CU comparativement à des biopsies non-inflammées ou à des biopsies d’individus sains. GPR65, lorsqu’activée par un pH acide, stimule l’accumulation d’AMPc, la formation de fibres de stress et l’activation de la voie RhoA. GPR65 semble donc être un bon candidat pour l’étude de l’étiologie des MIIs. Nos objectifs étaient donc de définir les voies découlant de l’activation de GPR65 et d’évaluer l’impact de GPR65*231Leu sur ces voies. Nous avons utilisé des HEK293 exprimant de façon stable l’une ou l’autre des deux allèles de GPR65 et déficientes pour Gαs/olf, Gαq/11 ou Gα12/13 puisqu’il est connu que Gαs/olf activent l’adénylate cyclase, Gαq/11 mènent au relâchement de Ca2+ intracellulaire et peuvent activer certaines isoformes de l’adénylate cyclase et Gα12/13 régulent le remodelage du cytosquelette d’actine. Nous avons démontré que l’accumulation d’AMPc dépendante de l’activation de GPR65 par un pH acide est causée, au moins en partie, aux voies Gαs/olf et peu ou pas aux voies Gαq/11. Il ne semble pas y avoir un effet du variant GPR65*231Leu sur ces voies. Cependant, nous avons observé que le variant GPR65*231Leu réduit la formation de fibres de stress et inhibe l’accumulation d’actine filamenteuse (F-actine) dépendantes de l’activation de GPR65 par un pH acide. Des données préliminaires nous montrent que la formation de fibres de stress est causée, au moins en partie, aux voies G12/13. En conclusion, nous avons démontré que GPR65 active les voies Gαs/olf et Gα12/13 et que le variant codant associé aux MIIs altère le remodelage de l’actine. GPR65 pourrait donc potentiellement avoir un rôle dans la pathologie des MIIs. / Inflammatory bowel diseases, mainly comprising of ulcerative colitis (UC) and Crohn’s disease (CD), result in the chronic inflammation of the gastrointestinal tract. Genome-wide association studies have associated the 14q31 locus, including the genes galactosylceramidase (GALC) and G protein-coupled receptor 65 (GPR65), to those phenotypes. The most associated variant in this region for IBD (rs8005161; p=2,35x10-14) is correlated (r2=1) to a missense coding variant of GPR65 (rs3742704: Ile231Leu). GPR65 encodes a pH-sensing G protein-coupled receptor (GPCR). We observed that GPR65 is expressed in lymphoid and mucosal tissues as well as in immune cell lines and human primary immune cells. We also found that its expression is significantly increased in inflamed biopsies from UC patients compared to non-inflamed biopsies and biopsies from healthy controls. Upon activation by low pH, GPR65 stimulates accumulation of cAMP, formation of stress fibers and activation of the RhoA pathway. Thus GPR65 is a good candidate causal gene for the IBD pathology. Our objective, therefore, was to define pathways downstream of activation of GPR65 and to evaluate the impact GPR65*231Leu on these pathways. We used HEK 293 cells stably expressing one or the other allele of GPR65 and deficient for either Gαs/olf, Gαq/11 or Gα12/13 as it is known that Gαs/olf activates adenylyl cyclase, Gαq/11 leads to the release of intracellular Ca2+, which can also activate certain isoforms of adenylyl cyclase, and that Gα12/13 regulates actin cytoskeletal remodeling. We demonstrated that cAMP accumulation upon activation of GPR65 is, at least partly, due to the Gαs/olf pathway and only slightly or not at all to the Gαq/11 pathway. The coding variant, however, does not appear to have an effect on that pathway. In contrast, we observed that GPR65*231Leu variant reduces the GPR65-dependent stress fiber formation and inhibits the increase of filamentous actin (F-actin) content versus free globular-actin (G-actin) upon activation by low pH. Also, preliminary data showed that the stress fiber formation in HEK293 cells is due to the G12/13 pathways. In conclusion, we demonstrate that GPR65 activates both Gαs and Gα12/13, the coding variant alters the actin remodeling pathway and that GPR65 potentially has a causal role in IBD.

No description available.

Ce travail de recherche propose un modèle numérique pour simuler les procédés par infusion de résine en utilisant la méthode des éléments finis. Ce modèle permet de représenter l'écoulement d'une résine liquide dans des préformes poreuses subissant de grandes déformations. Dans cette étude, une modélisation macroscopique est utilisée. Au niveau du procédé, une zone de résine liquide est déposée sur les préformes. Ces dernières étant considérées comme un milieu poreux. Les équations de Stokes et de Darcy sont utilisées pour modéliser l'écoulement de la résine respectivement dans le drainant et dans les préformes. L'originalité du modèle réside dans le fait qu'un seul maillage est utilisé pour les deux milieux. La discrétisation est réalisée avec des éléments mixtes : dans Stokes, des éléments P1+/P1 sont utilisés et dans Darcy, des éléments P1/P1 stabilisés avec une formulation multi-échelle sont employés. Des fonctions distances signées sont utilisées pour représenter l'interface entre Stokes-Darcy et pour représenter le front de résine. Concernant la déformation des préformes, une formulation Lagrangienne réactualisée est utilisée. Dans cette formulation Lagrangienne, le comportement des préformes humides est représenté à l'aide du modèle de Terzaghi dans lequel les préformes sèches ont un comportement élastique non-linéaire. La perméabilité est reliée à la porosité via la relation de Carman-Kozeny. Celle-ci est déterminée à partir de l'équation de conservation de la masse. Ce modèle a été implémenté dans ZéBuLoN. Plusieurs simulations numériques d'infusion de résine sont présentées à la fin de ce manuscrit. / This work proposes a numerical model to simulate the manufacturing processes by resin infusion using the finite element method. This model allows to represent the resin flow into porous preforms, which are themselves subject to large deformations. In this study, a macroscopic description is used. The preforms are considered as a porous medium. The Stokes and the Darcy equations are used respectively to describe the resin flow into the liquid zone and into the preforms.The originality of the model consists in using one single unstructured mesh. The discretization is ensured by using a mixed velocity-pressure formulation. Indeed, a P1/P1 formulation is employed throughout the entire discretized domain, stabilized in the Darcy region with a multi-scale formulation and in the Stokes subdomain with a hierarchical-based bubble, i.e. a P1+/P1 finite element. Signed distance functions are used both to represent the Stokes-Darcy interface and to capture the moving flow front. Concerning the deformations of the preforms, an updated Lagrangian scheme is used. In the Lagrangian formulation, the behavior of the wet preforms is represented by using the Terzaghi model in which the dry preforms have a non-linear elastic behavior. The permeability depends on the porosity through the Carman-Kozeny relationship. This model has been implemented in Zset. Several numerical simulations of manufacturing processes by resin infusion are presented at the end of this manuscript.

L'évolution importante des réseaux sans fil tend à fournir les supports nécessaires aux applications ubiquitaires émergentes dans les réseaux Mesh sans fil. Les réseaux mesh comprennent des nœuds stationnaires qui remplissent la fonction de routage et appelés routeurs Mesh sans fil (WMR) et qui consitutent le réseau backbone sans fil (WBMN) ainsi que des nœuds clients mesh sans fil (WMCs). Alors que les WMCs sont limités en termes de fonctions et de resources radio, les WMRS fournissent des fonctions de pont et de passerelle afin de connecter les réseaux WMNs aux autres réseaux comme les réseaux cellulaires, les réseaux IEEE 802.11, les réseaux IEEE 802.16, ou tout simplement à Internet. Par conséquent, les réseaux WMRs sont construits à partir sur la base d'équipement radio de communication rapide et/ou multi-radio et multi-canaux. Les routeurs WMRs sont supposes être auto-organisés, auto-configurés et constituant un réseau WNMN robuste ce qui nécessite de soutenir des volumes importants de trafic de données et sur de longues périodes. Cependant, répondre à attentes élevées en termes de services nécessite le développement d'approches décentralisés pour le control dynamique des puissances de transmission (DTPC). La présente thèse se focalise sur le problème DTPC pour les deux cas de réseaux; utilisant un canal de communication et multicanaux. Pour les réseaux exploitant un seul canal, le problème est formulé en termes de problème d'optimisation où l'objectif est de minimiser en même temps des critères convexes associés aux liens et aux réseaux. Afin de résoudre ce problème, des modèles et des algorithmes, appelés MATA (multiple access transmission aware), ont été proposés. Pour les réseaux WBMNs utilisant des systèmes multi-radio et multi-canaux (MRMC), le réseau est modélisé par un ensemble de graphes appelés UCGs (unified channel graphs), chacun consistant les utilisateurs connectés au point d'accès en utilisant le même canal fréquentiel. Pour chaque ensemble UCG, le problème à résoudre un problème quadratique et stochastique soumis aux contraintes des états des liens dynamiques LSI (Link State Information) de tous les UCGs. Le protocole PMMUP, mutli-radio et minimisant la consummation énegétique, est propose au niveau de la couche liaison. Algorithmes d'estimation prédictive base sur ce protocol sont proposes pour résoudre les problèmes d'optimisation associés aux UGCs. Les problèmes énergétiques, les instabilités des queues et les interférences, sont formulés en termes de problèmes de commande optimale couplée, appelés SPWC (singularly-perturbed weakly-coupled). Pour résoudre les problèmes SPWC caractérisant le problème de commande optimal des energies de transmission, l'algorithme HORA (aeneralized higher-order recursive algorithm) qui permet d'obtenir les solutions stables pour les équations de Riccati a été développé. Les performes des modèles et algorithmes proposés dans le cadre de la présente thèse ont été évalués tant sur le plan théorique qu'en simulation. Différentes simulations ont été effectuées sur un large ensemble de topologies réseaux générés aléatoirement. Les résultats de simulation et analytiques Simulation confirment l'efficacité des algorithmes proposés par rapport à la majorité des techniques existantes / The remarkable evolution of wireless networks into the next generation to provide ubiquitous and seamless broadband applications has recently triggered the emergence of Wireless Mesh Networks (WMNs). The WMNs comprise stationary Wireless Mesh Routers (WMRs) forming Wireless Backbone Mesh Networks (WBMNs) and mobile Wireless Mesh Clients (WMCs) forming the WMN access. While WMCs are limited in function and radio resources, the WMRs are expected to support heavy duty applications : that is, WMRs have gateway and bridge functions to integrate WMNs with other networks such as the Internet, cellular, IEEE 802.11, IEEE 802.15, IEEE 802.16, sensor networks, et cetera. Consequently, WMRs are constructed from fast switching radios or multiple radio devices operating on multiple frequency channels. WMRs are expected to be self-organized, self-configured and constitute a reliable and robust WBMN which needs to sustain high traffic volumes and long “online” time. However, meeting such stringent service expectations requires the development of decentralized dynamic transmission power control (DTPC) approaches. This thesis addresses the DTPC problem for both single and multiple channel WBMNs. For single channel networks, the problem is formulated as the minimization of both the link-centric and network-centric convex cost function. In order to solve this issue, multiple access transmission aware (MATA) models and algorithms are proposed. For multi-radio multi-channel (MRMC) WBMNs, the network is modelled as sets of unified channel graphs (UCGs), each consisting of interconnected active network users communicating on the same frequency channel. For each UCG set, the minimization of stochastic quadratic cost functions are developed subject to the dynamic Link State Information (LSI) equations from all UCGs. An energy-efficient multi-radio unification protocol (PMMUP) is then suggested at the Link-Layer (LL). Predictive estimation algorithms based on this protocol are proposed to solve such objective functions. To address transmission energy and packet instabilities, and interference across multiple channels, singularly-perturbed weakly-coupled (SPWC) control problems are formulated. In order to solve the SPWC transmission power control problem, a generalized higher-order recursive algorithm (HORA) that obtains the Riccati Stabilizing Solutions to the control problem is developed. The performance behaviours of the proposed models and algorithms are evaluated both analytically and through computer simulations. Several simulations are performed on a large number of randomly generated topologies. Simulation and analytical results confirm the efficacy of the proposed algorithms compared to the most recently studied techniques