An intrinsic Proper Generalized Decomposition for parametric symmetric elliptic problems

We introduce in this paper a technique for the reduced order approximation of parametric symmetric elliptic partial differential equations. For any given dimension, we prove the existence of an optimal subspace of at most that dimension which realizes the best approximation in mean of the error with respect to the parameter in the quadratic norm associated to the elliptic operator, between the exact solution and the Galerkin solution calculated on the subspace. This is analogous to the best approximation property of the Proper Orthogonal Decomposition (POD) subspaces, excepting that in our case the norm is parameter-depending, and then the POD optimal sub-spaces cannot be characterized by means of a spectral problem. We apply a deflation technique to build a series of approximating solutions on finite-dimensional optimal subspaces, directly in the on-line step. We prove that the partial sums converge to the continuous solutions, in mean quadratic elliptic norm.Comment: 18 page

The Influence of Quadrature Errors on Isogeometric Mortar Methods

Mortar methods have recently been shown to be well suited for isogeometric analysis. We review the recent mathematical analysis and then investigate the variational crime introduced by quadrature formulas for the coupling integrals. Motivated by finite element observations, we consider a quadrature rule purely based on the slave mesh as well as a method using quadrature rules based on the slave mesh and on the master mesh, resulting in a non-symmetric saddle point problem. While in the first case reduced convergence rates can be observed, in the second case the influence of the variational crime is less significant

Multiscale computing with the multiscale modeling library and runtime environment

We introduce a software tool to simulate multiscale models: The Multiscale Coupling Library and Environment 2 (MUSCLE 2). MUSCLE 2 is a component-based modeling tool inspired by the multiscale modeling and simulation framework, with an easy-to-use API which supports Java, C++, C, and Fortran. We present MUSCLE 2's runtime features, such as its distributed computing capabilities, and its benefits to multiscale modelers. We also describe two multiscale models that use MUSCLE 2 to do distributed multiscale computing: An in-stent restenosis and a canal system model. We conclude that MUSCLE 2 is a notable improvement over the previous version of MUSCLE, and that it allows users to more flexibly deploy simulations of multiscale models, while improving their performance. © 2013 The Authors. Published by Elsevier B.V

ITERATED QUASI-REVERSIBILITY METHOD APPLIED TO ELLIPTIC AND PARABOLIC DATA COMPLETION PROBLEMS

International audienceWe study the iterated quasi-reversibility method to regularize ill-posed elliptic and parabolic problems: data completion problems for Poisson's and heat equations. We define an abstract setting to treat both equations at once. We demonstrate the convergence of the regularized solution to the exact one, and propose a strategy to deal with noise on the data. We present numerical experiments for both problems: a two-dimensional corrosion detection problem and the one-dimensional heat equation with lateral data. In both cases, the method prove to be efficient even with highly corrupted data

Light induced self-written waveguides interactions in photopolymer media

We present experimental and theoretical study of the interaction of Light Induced Self-Written (LISW) waveguides in photopolymers. We show that the diffusion of the monomer controls the refractive index distribution. Consequently it influences the interaction between the LISW channels allowing the observation of anti-crossing behavior or the propagation of an array of non interacting LISW waveguides

Reduced basis isogeometric mortar approximations for eigenvalue problems in vibroacoustics

We simulate the vibration of a violin bridge in a multi-query context using reduced basis techniques. The mathematical model is based on an eigenvalue problem for the orthotropic linear elasticity equation. In addition to the nine material parameters, a geometrical thickness parameter is considered. This parameter enters as a 10th material parameter into the system by a mapping onto a parameter independent reference domain. The detailed simulation is carried out by isogeometric mortar methods. Weakly coupled patch-wise tensorial structured isogeometric elements are of special interest for complex geometries with piecewise smooth but curvilinear boundaries. To obtain locality in the detailed system, we use the saddle point approach and do not apply static condensation techniques. However within the reduced basis context, it is natural to eliminate the Lagrange multiplier and formulate a reduced eigenvalue problem for a symmetric positive definite matrix. The selection of the snapshots is controlled by a multi-query greedy strategy taking into account an error indicator allowing for multiple eigenvalues

Performance of distributed multiscale simulations

Multiscale simulations model phenomena across natural scales using monolithic or component-based code, running on local or distributed resources. In this work, we investigate the performance of distributed multiscale computing of component-based models, guided by six multiscale applications with different characteristics and from several disciplines. Three modes of distributed multiscale computing are identified: supplementing local dependencies with large-scale resources, load distribution over multiple resources, and load balancing of small- and large-scale resources. We find that the first mode has the apparent benefit of increasing simulation speed, and the second mode can increase simulation speed if local resources are limited. Depending on resource reservation and model coupling topology, the third mode may result in a reduction of resource consumption

Weakly Consistent Regularisation Methods for Ill-Posed Problems

This Chapter takes its origin in the lecture notes for a 9 h course at the Institut Henri Poincaré in September 2016. The course was divided in three parts. In the first part, which is not included herein, the aim was to first recall some basic aspects of stabilised finite element methods for convection-diffusion problems. We focus entirely on the second and third parts which were dedicated to ill-posed problems and their approximation using stabilised finite element methods. First we introduce the concept of conditional stability. Then we consider the elliptic Cauchy-problem and a data assimilation problem in a unified setting and show how stabilised finite element methods may be used to derive error estimates that are consistent with the stability properties of the problem and the approximation properties of the finite element space. Finally, we extend the result to a data assimilation problem subject to the heat equation

Gustatory Perception and Fat Body Energy Metabolism Are Jointly Affected by Vitellogenin and Juvenile Hormone in Honey Bees

Honey bees (Apis mellifera) provide a system for studying social and food-related behavior. A caste of workers performs age-related tasks: young bees (nurses) usually feed the brood and other adult bees inside the nest, while older bees (foragers) forage outside for pollen, a protein/lipid source, or nectar, a carbohydrate source. The workers' transition from nursing to foraging and their foraging preferences correlate with differences in gustatory perception, metabolic gene expression, and endocrine physiology including the endocrine factors vitellogenin (Vg) and juvenile hormone (JH). However, the understanding of connections among social behavior, energy metabolism, and endocrine factors is incomplete. We used RNA interference (RNAi) to perturb the gene network of Vg and JH to learn more about these connections through effects on gustation, gene transcripts, and physiology. The RNAi perturbation was achieved by single and double knockdown of the genes ultraspiracle (usp) and vg, which encode a putative JH receptor and Vg, respectively. The double knockdown enhanced gustatory perception and elevated hemolymph glucose, trehalose, and JH. We also observed transcriptional responses in insulin like peptide 1 (ilp1), the adipokinetic hormone receptor (AKHR), and cGMP-dependent protein kinase (PKG, or “foraging gene” Amfor). Our study demonstrates that the Vg–JH regulatory module controls changes in carbohydrate metabolism, but not lipid metabolism, when worker bees shift from nursing to foraging. The module is also placed upstream of ilp1, AKHR, and PKG for the first time. As insulin, adipokinetic hormone (AKH), and PKG pathways influence metabolism and gustation in many animals, we propose that honey bees have conserved pathways in carbohydrate metabolism and conserved connections between energy metabolism and gustatory perception. Thus, perhaps the bee can make general contributions to the understanding of food-related behavior and metabolic disorders