Search

A second generation wavelet-based adaptive finite-difference Lattice Boltzmann method (FD-LBM) is developed in this paper. In this approach, the adaptive wavelet collocation method (AWCM) is firstly, to the best of our knowledge, incorporated into the FD-LBM. According to the grid refinement criterion based on the wavelet amplitudes of density distribution functions, an adaptive sparse grid is generated by the omission and addition of collocation points. On the sparse grid, the finite differences are used to approximate the derivatives. To eliminate the special treatments in using the FD-based derivative approximation near boundaries, the immersed boundary method (IBM) is also introduced into FD-LBM. By using the adaptive technique, the adaptive code requires much less grid points as compared to the uniform-mesh code. As a consequence, the computational efficiency can be improved. To justify the proposed method, a series of test cases, including fixed boundary cases and moving boundary cases, are invested. A good agreement between the present results and the data in previous literatures is obtained, which demonstrates the accuracy and effectiveness of the present AWCM-IB-LBM. [ABSTRACT FROM AUTHOR]

A parallel adaptive wavelet collocation method for solving a large class of Partial Differential Equations is presented. The parallelization is achieved by developing an asynchronous parallel wavelet transform, which allows one to perform parallel wavelet transform and derivative calculations with only one data synchronization at the highest level of resolution. The data are stored using tree-like structure with tree roots starting at a priori defined level of resolution. Both static and dynamic domain partitioning approaches are developed. For the dynamic domain partitioning, trees are considered to be the minimum quanta of data to be migrated between the processes. This allows fully automated and efficient handling of non-simply connected partitioning of a computational domain. Dynamic load balancing is achieved via domain repartitioning during the grid adaptation step and reassigning trees to the appropriate processes to ensure approximately the same number of grid points on each process. The parallel efficiency of the approach is discussed based on parallel adaptive wavelet-based Coherent Vortex Simulations of homogeneous turbulence with linear forcing at effective non-adaptive resolutions up to 2048 3 using as many as 2048 CPU cores. [ABSTRACT FROM AUTHOR]

Abstract: We address the problem of spherical deconvolution in a non-parametric statistical framework, where both the signal and the operator kernel are subject to measurement errors. After a preliminary treatment of the kernel, we apply a thresholding procedure to the signal in a second generation wavelet basis. Under standard assumptions on the kernel, we study the minimax performances of the resulting algorithm in terms of losses ( ) on Besov spaces on the sphere. We hereby extend the application of second generation spherical wavelets to the blind spherical deconvolution framework. It is important to stress that the procedure is adaptive with regard to both the target function sparsity and the kernel blurring effect. We end with the study of a concrete example, putting into evidence the improvement of our procedure on the recent blockwise SVD algorithm of Delattre et al. (2012). [Copyright &y& Elsevier]

Abstract: We consider the problem of estimating a density of probability from indirect data in the spherical convolution model. We aim at building an estimate of the unknown density as a linear combination of functions of an overcomplete dictionary. The procedure is devised through a well-calibrated -penalized criterion. The spherical deconvolution setting has been barely studied so far, and the two main approaches to this problem, namely the SVD and the hard thresholding ones considered only one basis at a time. The dictionary approach allows to combine various bases and thus enhances estimates sparsity. We provide an oracle inequality under global coherence assumptions. Moreover, the calibrated procedure that we put forward gives quite satisfying results in the numerical study when compared with other procedures. [Copyright &y& Elsevier]

Abstract: In this paper, a new method is presented to identify crack location and size, which is based on stress intensity factor suitable for pipe structure and finite element method of second generation wavelets (SGW-FEM). Pipe structure is dispersed into a series of nested thin-walled pipes. By making use of stress intensity factor of the thin-walled pipe, a new calculation method of crack equivalent stiffness is proposed to solve the stress intensity factor of the pipe structure. On this basis, finite element method of second generation wavelets is used to establish the dynamic model of cracked pipe. Then we combine forward problem with inverse problem in order to establish quantitative identification method of the crack based on frequency change, which provides a non-destructive testing technology with vibration for the pipe structure. The efficiency of the proposed method is verified by experiments. [Copyright &y& Elsevier]

Abstract: A dynamic adaptive numerical method for solving partial differential equations on the sphere is developed. The method is based on second generation spherical wavelets on almost uniform nested spherical triangular grids, and is an extension of the adaptive wavelet collocation method to curved manifolds. Wavelet decomposition is used for grid adaption and interpolation. An hierarchical finite difference scheme based on the wavelet multilevel decomposition is used to approximate Laplace–Beltrami, Jacobian and flux-divergence operators. The accuracy and efficiency of the method is demonstrated using linear and nonlinear examples relevant to geophysical flows. Although the present paper considers only the sphere, the strength of this new method is that it can be extended easily to other curved manifolds by considering appropriate coarse approximations to the desired manifold (here we used the icosahedral approximation to the sphere at the coarsest level). [Copyright &y& Elsevier]

Abstract: Dynamically adaptive numerical methods have been developed to efficiently solve differential equations whose solutions are intermittent in both space and time. These methods combine an adjustable time step with a spatial grid that adapts to spatial intermittency at a fixed time. The same time step is used for all spatial locations and all scales: this approach clearly does not fully exploit space–time intermittency. We propose an adaptive wavelet collocation method for solving highly intermittent problems (e.g. turbulence) on a simultaneous space–time computational domain which naturally adapts both the space and time resolution to match the solution. Besides generating a near optimal grid for the full space–time solution, this approach also allows the global time integration error to be controlled. The efficiency and accuracy of the method is demonstrated by applying it to several highly intermittent (1D + t)-dimensional and (2D + t)-dimensional test problems. In particular, we found that the space–time method uses roughly 18 times fewer space–time grid points and is roughly 4 times faster than a dynamically adaptive explicit time marching method, while achieving similar global accuracy. [Copyright &y& Elsevier]

Abstract: An adaptive multilevel wavelet collocation method for solving multi-dimensional elliptic problems with localized structures is described. The method is based on multi-dimensional second generation wavelets, and is an extension of the dynamically adaptive second generation wavelet collocation method for evolution problems [Int. J. Comp. Fluid Dyn. 17 (2003) 151]. Wavelet decomposition is used for grid adaptation and interpolation, while a hierarchical finite difference scheme, which takes advantage of wavelet multilevel decomposition, is used for derivative calculations. The multilevel structure of the wavelet approximation provides a natural way to obtain the solution on a near optimal grid. In order to accelerate the convergence of the solver, an iterative procedure analogous to the multigrid algorithm is developed. The overall computational complexity of the solver is , where is the number of adapted grid points. The accuracy and computational efficiency of the method are demonstrated for the solution of two- and three-dimensional elliptic test problems. [Copyright &y& Elsevier]