RADIAL BASIS FUNCTION COLLOCATION FOR THE CHEMICAL MASTER EQUATION

2010 ◽  
Vol 07 (03) ◽  
pp. 477-498 ◽  
Author(s):  
JINGWEI ZHANG ◽  
LAYNE T. WATSON ◽  
CHRISTOPHER A. BEATTIE ◽  
YANG CAO

The chemical master equation (CME), formulated from the Markov assumption of stochastic processes, offers an accurate description of general chemical reaction systems. This paper proposes a collocation method using radial basis functions to numerically approximate the solution to the CME. Numerical results for some systems biology problems show that the collocation approximation method has good potential for solving large-scale CMEs.

Acta Numerica ◽  
2015 ◽  
Vol 24 ◽  
pp. 215-258 ◽  
Author(s):  
Bengt Fornberg ◽  
Natasha Flyer

Finite differences provided the first numerical approach that permitted large-scale simulations in many applications areas, such as geophysical fluid dynamics. As accuracy and integration time requirements gradually increased, the focus shifted from finite differences to a variety of different spectral methods. During the last few years, radial basis functions, in particular in their ‘local’ RBF-FD form, have taken the major step from being mostly a curiosity approach for small-scale PDE ‘toy problems’ to becoming a major contender also for very large simulations on advanced distributed memory computer systems. Being entirely mesh-free, RBF-FD discretizations are also particularly easy to implement, even when local refinements are needed. This article gives some background to this development, and highlights some recent results.


Author(s):  
Kaveh Amouzgar ◽  
Niclas Stromberg

In this paper, an approach to generate surrogate models constructed by radial basis function networks (RBFN) with a priori bias is presented. RBFN as a weighted combination of radial basis functions only, might become singular and no interpolation is found. The standard approach to avoid this is to add a polynomial bias, where the bias is defined by imposing orthogonality conditions between the weights of the radial basis functions and the polynomial basis functions. Here, in the proposed a priori approach, the regression coefficients of the polynomial bias are simply calculated by using the normal equation without any need of the extra orthogonality prerequisite. In addition to the simplicity of this approach, the method has also proven to predict the actual functions more accurately compared to the RBFN with a posteriori bias. Several test functions, including Rosenbrock, Branin-Hoo, Goldstein-Price functions and two mathematical functions (one large scale), are used to evaluate the performance of the proposed method by conducting a comparison study and error analysis between the RBFN with a priori and a posteriori known biases. Furthermore, the aforementioned approaches are applied to an engineering design problem, that is modeling of the material properties of a three phase spherical graphite iron (SGI). The corresponding surrogate models are presented and compared.


Author(s):  
Jan Kamenik ◽  
Michele Stramacchia ◽  
David J. J. Toal ◽  
Andy J. Keane ◽  
Ron Bates

A novel infill criterion for so-called ensemble of surrogates-based optimization is proposed and applied in practice for an aerodynamic compressor rotor design optimization. The ensemble uses a combined approach based on different radial basis functions and aims to reduce prediction errors through weighted linear combinations of radial basis functions. The update strategy uses a new hybrid custom metric termed α, which incorporates information about each surrogate’s local agreement through correlation coefficients and also information about the global accuracy of each ensemble combination through the root-mean-square error. Surrogate models are searched using a hybrid optimizer, i.e., with a genetic algorithm and sequential quadratic programming, and proposed update points are evaluated using the high-fidelity black box function. The results are compared with established optimization approaches and the best design is analyzed further in terms of the flow physics. Results show that α-based ensemble of surrogates approaches are particularly efficient for large-scale cases, where other types of surrogates such as Kriging models are onerous to construct.


Sign in / Sign up

Export Citation Format

Share Document