In this paper, we combine the multiscale flnite element method to propose an algorithm for solving the non-stationary Stokes-Darcy model, where the permeability coefflcient in the Darcy region exhibits multiscale characteristics. Our algorithm involves two steps: first, conducting the parallel computation of multiscale basis functions in the Darcy region. Second, based on these multiscale basis functions, we employ an implicitexplicit scheme to solve the Stokes-Darcy equations. One signiflcant feature of the algorithm is that it solves problems on relatively coarse grids, thus signiflcantly reducing computational costs. Moreover, under the same coarse grid size, it exhibits higher accuracy compared to standard flnite element method. Under the assumption that the permeability coefflcient is periodic and independent of time, this paper demonstrates the stability and convergence of the algorithm. Finally, the rationality and effectiveness of the algorithm are verifled through three numerical experiments, with experimental results consistent with theoretical analysis.
相關內容
ACM/IEEE第23屆模型驅動工程語言和系統國際會議，是模型驅動軟件和系統工程的首要會議系列，由ACM-SIGSOFT和IEEE-TCSE支持組織。自1998年以來，模型涵蓋了建模的各個方面，從語言和方法到工具和應用程序。模特的參加者來自不同的背景，包括研究人員、學者、工程師和工業專業人士。MODELS 2019是一個論壇，參與者可以圍繞建模和模型驅動的軟件和系統交流前沿研究成果和創新實踐經驗。今年的版本將為建模社區提供進一步推進建模基礎的機會，并在網絡物理系統、嵌入式系統、社會技術系統、云計算、大數據、機器學習、安全、開源等新興領域提出建模的創新應用以及可持續性。
官網鏈接： ·
混合
·
估計/估計量
·
離散化
·
CC
·
In this paper, a novel multigrid method based on Newton iteration is proposed to solve nonlinear eigenvalue problems. Instead of handling the eigenvalue $\lambda$ and eigenfunction $u$ separately, we treat the eigenpair $(\lambda, u)$ as one element in a product space $\mathbb R \times H_0^1(\Omega)$. Then in the presented multigrid method, only one discrete linear boundary value problem needs to be solved for each level of the multigrid sequence. Because we avoid solving large-scale nonlinear eigenvalue problems directly, the overall efficiency is significantly improved. The optimal error estimate and linear computational complexity can be derived simultaneously. In addition, we also provide an improved multigrid method coupled with a mixing scheme to further guarantee the convergence and stability of the iteration scheme. More importantly, we prove convergence for the residuals after each iteration step. For nonlinear eigenvalue problems, such theoretical analysis is missing from the existing literatures on the mixing iteration scheme.
We present a novel family of particle discretisation methods for the nonlinear Landau collision operator. We exploit the metriplectic structure underlying the Vlasov-Maxwell-Landau system in order to obtain disretisation schemes that automatically preserve mass, momentum, and energy, warrant monotonic dissipation of entropy, and are thus guaranteed to respect the laws of thermodynamics. In contrast to recent works that used radial basis functions and similar methods for regularisation, here we use an auxiliary spline or finite element representation of the distribution function to this end. Discrete gradient methods are employed to guarantee the aforementioned properties in the time discrete domain as well.
We propose a local discontinuous Galerkin (LDG) method for fractional Korteweg-de Vries equation involving the fractional Laplacian with exponent $\alpha\in (1,2)$ in one and two space dimensions. By decomposing the fractional Laplacian into a first order derivative and a fractional integral, we prove $L^2$-stability of the semi-discrete LDG scheme incorporating suitable interface and boundary fluxes. We analyze the error estimate by considering linear convection term and utilizing the estimate, we derive the error estimate for general nonlinear flux and demonstrate an order of convergence $\mathcal{O}(h^{k+1/2})$. Moreover, the stability and error analysis have been extended to multiple space dimensional case. Additionally, we devise a fully discrete LDG scheme using the four-stage fourth-order Runge-Kutta method. We prove that the scheme is strongly stable under an appropriate time step constraint by establishing a \emph{three-step strong stability} estimate. Numerical illustrations are shown to demonstrate the efficiency of the scheme by obtaining an optimal order of convergence.
We detail the mathematical formulation of the line of "functional quantizer" modules developed by the Mathematics and Music Lab (MML) at Michigan Technological University, for the VCV Rack software modular synthesizer platform, which allow synthesizer players to tune oscillators to new musical scales based on mathematical functions. For example, we describe the recently-released MML Logarithmic Quantizer (LOG QNT) module that tunes synthesizer oscillators to the non-Pythagorean musical scale introduced by indie band The Apples in Stereo.
We introduce a novel method for non-convex optimization which is at the interface between the swarm-based gradient-descent (SBGD) [J. Lu et. al., ArXiv:2211.17157; E.Tadmor and A. Zenginoglu, Acta Applicandae Math., 190, 2024] and Simulated Annealing (SA) [V. Cerny, J. optimization theory and appl., 45:41-51, 1985; S.Kirkpatrick et. al., Science, 220(4598):671-680, 1983; S. Geman and C.-R. Hwang, SIAM J. Control and Optimization, 24(5):1031-1043, 1986]. We follow the methodology of SBGD in which a swarm of agents, each identified with a position, ${\mathbf x}$ and mass $m$, explores the ambient space. The agents proceed in gradient descent direction, and are subject to Brownian motion with annealing-rate dictated by a decreasing function of their mass. Thus, instead of the SA protocol for time-decreasing temperature, we let the swarm decide how to `cool down' agents, depending on their accumulated mass over time. The dynamics of masses is coupled with the dynamics of positions: agents at higher ground transfer (part of) their mass to those at lower ground. Consequently, the swarm is dynamically divided between heavier, cooler agents viewed as `leaders' and lighter, warmer agents viewed as `explorers'. Mean-field convergence analysis and benchmark optimizations demonstrate the effectiveness of the swarm-based method as a multi-dimensional global optimizer.
We introduce a nonconforming hybrid finite element method for the two-dimensional vector Laplacian, based on a primal variational principle for which conforming methods are known to be inconsistent. Consistency is ensured using penalty terms similar to those used to stabilize hybridizable discontinuous Galerkin (HDG) methods, with a carefully chosen penalty parameter due to Brenner, Li, and Sung [Math. Comp., 76 (2007), pp. 573-595]. Our method accommodates elements of arbitrarily high order and, like HDG methods, it may be implemented efficiently using static condensation. The lowest-order case recovers the $P_1$-nonconforming method of Brenner, Cui, Li, and Sung [Numer. Math., 109 (2008), pp. 509-533], and we show that higher-order convergence is achieved under appropriate regularity assumptions. The analysis makes novel use of a family of weighted Sobolev spaces, due to Kondrat'ev, for domains admitting corner singularities.
In this paper, we present a new computational framework using coupled and decoupled approximations for a Cahn-Hilliard-Navier-Stokes model with variable densities and degenerate mobility. In this sense, the coupled approximation is shown to conserve the mass of the fluid, preserve the point-wise bounds of the density and decrease an energy functional. In contrast, the decoupled scheme is presented as a more computationally efficient alternative but the discrete energy-decreasing property can not be assured. Both schemes are based on a finite element approximation for the Navier-Stokes fluid flow with discontinuous pressure and an upwind discontinuous Galerkin scheme for the Cahn-Hilliard part. Finally, several numerical experiments contrasting both approaches are conducted. In particular, results for a convergence test, a simple qualitative comparison and some well-known benchmark problems are shown.
In this work, a family of symmetric interpolation points are generated on the four-dimensional simplex (i.e. the pentatope). These points are optimized in order to minimize the Lebesgue constant. The process of generating these points closely follows that outlined by Warburton in "An explicit construction of interpolation nodes on the simplex," Journal of Engineering Mathematics, 2006. Here, Warburton generated optimal interpolation points on the triangle and tetrahedron by formulating explicit geometric warping and blending functions, and applying these functions to equidistant nodal distributions. The locations of the resulting points were Lebesgue-optimized. In our work, we extend this procedure to four dimensions, and construct interpolation points on the pentatope up to order ten. The Lebesgue constants of our nodal sets are calculated, and are shown to outperform those of equidistant nodal distributions.
In the present work, we examine and analyze an hp-version interior penalty discontinuous Galerkin finite element method for the numerical approximation of a steady fluid system on computational meshes consisting of polytopic elements on the boundary. This approach is based on the discontinuous Galerkin method, enriched by arbitrarily shaped elements techniques as has been introduced in [13]. In this framework, and employing extensions of trace, Markov-type, and H1/L2-type inverse estimates to arbitrary element shapes, we examine a stationary Stokes fluid system enabling the proof of the inf/sup condition and the hp- a priori error estimates, while we investigate the optimal convergence rates numerically. This approach recovers and integrates the flexibility and superiority of the discontinuous Galerkin methods for fluids whenever geometrical deformations are taking place by degenerating the edges, facets, of the polytopic elements only on the boundary, combined with the efficiency of the hp-version techniques based on arbitrarily shaped elements without requiring any mapping from a given reference frame.
In this paper, we introduce a conservative Crank-Nicolson-type finite difference schemes for the regularized logarithmic Schr\"{o}dinger equation (RLSE) with Dirac delta potential in 1D. The regularized logarithmic Schr\"{o}dinger equation with a small regularized parameter $0<\eps \ll 1$ is adopted to approximate the logarithmic Schr\"{o}dinger equation (LSE) with linear convergence rate $O(\eps)$. The numerical method can be used to avoid numerical blow-up and/or to suppress round-off error due to the logarithmic nonlinearity in LSE. Then, by using domain-decomposition technique, we can transform the original problem into an interface problem. Different treatments on the interface conditions lead to different discrete schemes and it turns out that a simple discrete approximation of the Dirac potential coincides with one of the conservative finite difference schemes. The optimal $H^1$ error estimates and the conservative properties of the finite difference schemes are investigated. The Crank-Nicolson finite difference methods enjoy the second-order convergence rate in time and space. Numerical examples are provided to support our analysis and show the accuracy and efficiency of the numerical method.
注冊地址： 北京市海淀區羊坊店路18號2幢3層301-191