On local Fourier analysis of multigrid methods for PDEs with jumping and random coefficients

In this paper, we propose a novel nonstandard local Fourier analysis (LFA) variant for accurately predicting the multigrid convergence of problems with random and jumping coefficients. This LFA method is based on a specific basis of the Fourier space rather than the commonly used Fourier modes. To show the utility of this analysis, we consider, as an example, a simple cell-centered multigrid method for solving a steady-state single phase flow problem in a random porous medium. We successfully demonstrate the predictive capability of the proposed LFA using a number of challenging benchmark problems. The information provided by this analysis could be used to estimate a priori the time needed for solving certain uncertainty quantification problems by means of a multigrid multilevel Monte Carlo method

A parametric acceleration of multilevel Monte Carlo convergence for nonlinear variably saturated flow

We present a multilevel Monte Carlo (MLMC) method for the uncertainty quantification of variably saturated porous media flow that is modeled using the Richards equation. We propose a stochastic extension for the empirical models that are typically employed to close the Richards equations. This is achieved by treating the soil parameters in these models as spatially correlated random fields with appropriately defined marginal distributions. As some of these parameters can only take values in a specific range, non-Gaussian models are utilized. The randomness in these parameters may result in path-wise highly nonlinear systems, so that a robust solver with respect to the random input is required. For this purpose, a solution method based on a combination of the modified Picard iteration and a cell-centered multigrid method for heterogeneous diffusion coefficients is utilized. Moreover, we propose a non-standard MLMC estimator to solve the resulting high-dimensional stochastic Richards equation. The improved efficiency of this multilevel estimator is achieved by parametric continuation that allows us to incorporate simpler nonlinear problems on coarser levels for variance reduction while the target strongly nonlinear problem is solved only on the finest level. Several numerical experiments are presented showing computational savings obtained by the new estimator compared with the original MC estimator

A multigrid multilevel Monte Carlo method for transport in the Darcy–Stokes system

A multilevel Monte Carlo (MLMC) method for Uncertainty Quantification (UQ) of advection-dominated contaminant transport in a coupled Darcy–Stokes flow system is described. In particular, we focus on high-dimensional epistemic uncertainty due to an unknown permeability field in the Darcy domain that is modelled as a lognormal random field. This paper explores different numerical strategies for the subproblems and suggests an optimal combination for the MLMC estimator. We propose a specific monolithic multigrid algorithm to efficiently solve the steady-state Darcy–Stokes flow with a highly heterogeneous diffusion coefficient. Furthermore, we describe an Alternating Direction Implicit (ADI) based time-stepping for the flux-limited quadratic upwinding discretization for the transport problem. Numerical experiments illustrating the multigrid convergence and cost of the MLMC estimator with respect to the smoothness of permeability field are presented

A parametric acceleration of multilevel Monte Carlo convergence for nonlinear variably saturated flow

We present a multilevel Monte Carlo (MLMC) method for the uncertainty quantification of variably saturated porous media flow that are modeled using the Richards’ equation. We propose a stochastic extension for the empirical models that are typically employed to close the Richards’ equations. This is achieved by treating the soil parameters in these models as spatially correlated random fields with appropriately defined marginal distributions. As some of these parameters can only take values in a specific range, non-Gaussian models are utilized. The randomness in these parameters may result in path-wise highly nonlinear systems, so that a robust solver with respect to the random input is required. For this purpose, a solution method based on a combination of the modified Picard iteration and a cell-centered multigrid method for heterogeneous diffusion coefficients is utilized. Moreover, we propose a non-standard MLMC estimator to solve the resulting high-dimensional stochastic Richards’ equation. The improved efficiency of this multilevel estimator is achieved by parametric continuation that allows us to incorporate simpler nonlinear problems on coarser levels for variance reduction while the target strongly nonlinear problem is solved only on the finest level. Several numerical experiments are presented showing computational savings obtained by the new estimator compared to the original MC estimator

Strategies and options for increasing and sustaining fisheries and aquaculture production to benefit poorer households in Asia [PDF in letter standard]

The last three decades have wi tnessed dramatic changes in the structure of supply and demand for fish, especially in Asia. This WorldFish research study sponsored by the Asian Development Bank focussed on nine developing countries û Bangladesh, China, India, Indonesia, Malaysia, the Philippines, Sri Lanka, Thailand, and Vietnam, all active players in the transformation of global fish supply and demand. The study, broken into five components and reported here, considered: 1) the profile of key aquaculture technologies and fishing practices; 2) analysis of policies, institutions and support services; 3) socioeconomic profile of major stakeholders in the fisheries sector; 4) projections of fish demand and supply in the nine Asian countries; and 5) formulation of national action plans based on the findings and recommendations of the study.Research, Fisheries, Economic analysis, Aquaculture, Fish consumption, Trade, Policies, Socioeconomic aspects, Technology, Fishery products, Asia, China, People's Rep., Bangladesh, India, Indonesia, Malaysia, Philippines, Malaysia, Sri Lanka, Thailand, Vietnam,

Effects of a nanoscopic filler on the structure and dynamics of a simulated polymer melt and the relationship to ultra-thin films

We perform molecular dynamics simulations of an idealized polymer melt surrounding a nanoscopic filler particle to probe the effects of a filler on the local melt structure and dynamics. We show that the glass transition temperature $T_g$ of the melt can be shifted to either higher or lower temperatures by appropriately tuning the interactions between polymer and filler. A gradual change of the polymer dynamics approaching the filler surface causes the change in the glass transition. We also find that while the bulk structure of the polymers changes little, the polymers close to the surface tend to be elongated and flattened, independent of the type of interaction we study. Consequently, the dynamics appear strongly influenced by the interactions, while the melt structure is only altered by the geometric constraints imposed by the presence of the filler. Our findings show a strong similarity to those obtained for ultra-thin polymer films (thickness $\lesssim 100$ nm) suggesting that both ultra-thin films and filled-polymer systems might be understood in the same context

Scale invariance in coarsening of binary and ternary fluids

Phase separation in binary and ternary fluids is studied using a two dimensional Lattice Gas Automata. The lengths, given by the the first zero crossing point of the correlation function and the total interface length is shown to exhibit power law dependence on time. In binary mixtures, our data clearly indicate the existence of a regime having more than one length scale where the coarsening process proceeds through the rupture and reassociation of domains. In ternary fluids; in the case of symmetric mixtures there exists a regime with a single length scale having dynamic exponent 1/2, while in asymmetric mixtures our data establish the break down of scale invariance.Comment: 20 pages, 13 figure

Low-rank representation of extended image volumes: Applications to imaging and velocity continuation

In this paper, we present a cost-effective low-rank factorized form of computing the full subsurface extended image volumes without explicitly calculating the receiver wavefields. The propose approach is computationally feasible, which exploits the low-rank structure of full subsurface extended image volumes organized as a matrix, thus avoiding the customary loop over sources. Using carefully selected stylized examples, we show how conventional migration images can be extracted from the low-rank factorized form. We also present a velocity continuation procedure that uses the same low-rank form and allows us to map the extended image for one velocity model to another velocity model without recomputing all the source wavefields

Stellar evolution and modelling stars

In this chapter I give an overall description of the structure and evolution of stars of different masses, and review the main ingredients included in state-of-the-art calculations aiming at reproducing observational features. I give particular emphasis to processes where large uncertainties still exist as they have strong impact on stellar properties derived from large compilations of tracks and isochrones, and are therefore of fundamental importance in many fields of astrophysics.Comment: Lecture presented at the IVth Azores International Advanced School in Space Sciences on "Asteroseismology and Exoplanets: Listening to the Stars and Searching for New Worlds" (arXiv:1709.00645), which took place in Horta, Azores Islands, Portugal in July 201

Strategies and options for increasing and sustaining fisheries and aquaculture production to benefit poorer households in Asia

The last three decades have witnessed dramatic changes in the structure of supply and demand for fish, especially in Asia. This WorldFish research study sponsored by the Asian Development Bank focussed on nine developing countries – Bangladesh, China, India, Indonesia, Malaysia, the Philippines, Sri Lanka, Thailand, and Vietnam, all active players in the transformation of global fish supply and demand. The study, broken into five components and reported here, considered: 1) the profile of key aquaculture technologies and fishing practices; 2) analysis of policies, institutions and support services; 3) socioeconomic profile of major stakeholders in the fisheries sector; 4) projections of fish demand and supply in the nine Asian countries; and 5) formulation of national action plans based on the findings and recommendations of the study