Search results for "Modeling and simulation"
showing 10 items of 1561 documents
Electrophoretic properties of charged colloidal suspensions: Application of a hybrid MD/LB method
2006
Abstract Electrophoretic properties of charged colloidal suspensions are investigated using a hybrid simulation method. In this method, the colloidal particles are propagated via Newton’s equations of motion using molecular dynamics (MD), while they are coupled to a structureless solvent that is modelled by the Lattice-Boltzmann (LB) method.
Numerical Study of Two Sparse AMG-methods
2003
A sparse algebraic multigrid method is studied as a cheap and accurate way to compute approximations of Schur complements of matrices arising from the discretization of some symmetric and positive definite partial differential operators. The construction of such a multigrid is discussed and numerical experiments are used to verify the properties of the method.
Well-balanced bicharacteristic-based scheme for multilayer shallow water flows including wet/dry fronts
2013
The aim of this paper is to present a new well-balanced finite volume scheme for two-dimensional multilayer shallow water flows including wet/dry fronts. The ideas, presented here for the two-layer model, can be generalized to a multilayer case in a straightforward way. The method developed here is constructed in the framework of the Finite Volume Evolution Galerkin (FVEG) schemes. The FVEG methods couple a finite volume formulation with evolution operators. The latter are constructed using the bicharacteristics of multidimensional hyperbolic systems. However, in the case of multilayer shallow water flows the required eigenstructure of the underlying equations is not readily available. Thus…
Discrete multiresolution based on hermite interpolation: computing derivatives
2004
Abstract Harten’s framework for multiresolution representation of data has been extended by Warming and Beam in [SIAM J. Sci. Comp. 22 (2000) 1269] to include Hermite interpolation. It needs the point-values of the derivative, which are usually unavailable, so they have to be approximated. In this work we show that the way in which the derivatives are approximated is crucial for the success of the method, and we present a new way to compute them that makes the scheme adequate for non-smooth data.
Comparison between adaptive and uniform discontinuous Galerkin simulations in dry 2D bubble experiments
2013
Accepted by the Journal of Computational Physics Adaptive mesh refinement generally aims to increase computational efficiency without compromising the accuracy of the numerical solution. However it is an open question in which regions the spatial resolution can actually be coarsened without affecting the accuracy of the result. This question is investigated for a specific example of dry atmospheric convection, namely the simulation of warm air bubbles. For this purpose a novel numerical model is developed that is tailored towards this specific application. The compressible Euler equations are solved with a Discontinuous Galerkin method. Time integration is done with an IMEXmethod and the dy…
Constructing adaptive generalized polynomial chaos method to measure the uncertainty in continuous models: A computational approach
2015
Due to errors in measurements and inherent variability in the quantities of interest, models based on random differential equations give more realistic results than their deterministic counterpart. The generalized polynomial chaos (gPC) is a powerful technique used to approximate the solution of these equations when the random inputs follow standard probability distributions. But in many cases these random inputs do not have a standard probability distribution. In this paper, we present a step-by-step constructive methodology to implement directly a useful version of adaptive gPC for arbitrary distributions, extending the applicability of the gPC. The paper mainly focuses on the computation…
A semi-Lagrangian AMR scheme for 2D transport problems in conservation form
2013
In this paper, we construct a semi-Lagrangian (SL) Adaptive-Mesh-Refinement (AMR) solver for 1D and 2D transport problems in conservation form. First, we describe the a-la-Harten AMR framework: the adaptation process selects a hierarchical set of grids with different resolutions depending on the features of the integrand function, using as criteria the point value prediction via interpolation from coarser meshes, and the appearance of large gradients. We integrate in time by reconstructing at the feet of the characteristics through the Point-Value Weighted Essentially Non-Oscillatory (PV-WENO) interpolator. We propose, then, an extension to the 2D setting by making the time integration dime…
Fractional differential equations solved by using Mellin transform
2014
In this paper, the solution of the multi-order differential equations, by using Mellin Transform, is proposed. It is shown that the problem related to the shift of the real part of the argument of the transformed function, arising when the Mellin integral operates on the fractional derivatives, may be overcame. Then, the solution may be found for any fractional differential equation involving multi-order fractional derivatives (or integrals). The solution is found in the Mellin domain, by solving a linear set of algebraic equations, whose inverse transform gives the solution of the fractional differential equation at hands.
Mellin transform approach for the solution of coupled systems of fractional differential equations
2015
In this paper, the solution of a multi-order, multi-degree-of-freedom fractional differential equation is addressed by using the Mellin integral transform. By taking advantage of a technique that relates the transformed function, in points of the complex plane differing in the value of their real part, the solution is found in the Mellin domain by solving a linear set of algebraic equations. The approximate solution of the differential (or integral) equation is restored, in the time domain, by using the inverse Mellin transform in its discretized form.
GPU accelerated Monte Carlo simulation of the 2D and 3D Ising model
2009
The compute unified device architecture (CUDA) is a programming approach for performing scientific calculations on a graphics processing unit (GPU) as a data-parallel computing device. The programming interface allows to implement algorithms using extensions to standard C language. With continuously increased number of cores in combination with a high memory bandwidth, a recent GPU offers incredible resources for general purpose computing. First, we apply this new technology to Monte Carlo simulations of the two dimensional ferromagnetic square lattice Ising model. By implementing a variant of the checkerboard algorithm, results are obtained up to 60 times faster on the GPU than on a curren…