亚洲男人的天堂2018av,欧美草比,久久久久久免费视频精选,国色天香在线看免费,久久久久亚洲av成人片仓井空

A quasi-second order scheme is developed to obtain approximate solutions of the shallow water equationswith bathymetry. The scheme is based on a staggered finite volume scheme for the space discretization:the scalar unknowns are located in the discretisation cells while the vector unknowns are located on theedges (in 2D) or faces (in 3D) of the mesh. A MUSCL-like interpolation for the discrete convectionoperators in the water height and momentum equations is performed in order to improve the precisionof the scheme. The time discretization is performed either by a first order segregated forward Eulerscheme in time or by the second order Heun scheme. Both schemes are shown to preserve the waterheight positivity under a CFL condition and an important state equilibrium known as the lake at rest.Using some recent Lax-Wendroff type results for staggered grids, these schemes are shown to be Lax-consistent with the weak formulation of the continuous equations; besides, the forward Euler schemeis shown to be consistent with a weak entropy inequality. Numerical results confirm the efficiency andaccuracy of the schemes.

相關內容

In this paper we study the numerical method for approximating the random periodic solution of semiliear stochastic evolution equations. The main challenge lies in proving a convergence over an infinite time horizon while simulating infinite-dimensional objects. We first show the existence and uniqueness of the random periodic solution to the equation as the limit of the pull-back flows of the equation, and observe that its mild form is well-defined in the intersection of a family of decreasing Hilbert spaces. Then we propose a Galerkin-type exponential integrator scheme and establish its convergence rate of the strong error to the mild solution, where the order of convergence directly depends on the space (among the family of Hilbert spaces) for the initial point to live. We finally conclude with the best order of convergence that is arbitrarily close to 0.5.

Continuous determinantal point processes (DPPs) are a class of repulsive point processes on $\mathbb{R}^d$ with many statistical applications. Although an explicit expression of their density is known, it is too complicated to be used directly for maximum likelihood estimation. In the stationary case, an approximation using Fourier series has been suggested, but it is limited to rectangular observation windows and no theoretical results support it. In this contribution, we investigate a different way to approximate the likelihood by looking at its asymptotic behaviour when the observation window grows towards $\mathbb{R}^d$. This new approximation is not limited to rectangular windows, is faster to compute than the previous one, does not require any tuning parameter, and some theoretical justifications are provided. It moreover provides an explicit formula for estimating the asymptotic variance of the associated estimator. The performances are assessed in a simulation study on standard parametric models on $\mathbb{R}^d$ and compare favourably to common alternative estimation methods for continuous DPPs.

We consider point sources in hyperbolic equations discretized by finite differences. If the source is stationary, appropriate source discretization has been shown to preserve the accuracy of the finite difference method. Moving point sources, however, pose two challenges that do not appear in the stationary case. First, the discrete source must not excite modes that propagate with the source velocity. Second, the discrete source spectrum amplitude must be independent of the source position. We derive a source discretization that meets these requirements and prove design-order convergence of the numerical solution for the one-dimensional advection equation. Numerical experiments indicate design-order convergence also for the acoustic wave equation in two dimensions. The source discretization covers on the order of $\sqrt{N}$ grid points on an $N$-point grid and is applicable for source trajectories that do not touch domain boundaries.

In this paper, a numerical scheme for a nonlinear McKendrick-von Foerster equation with diffusion in age (MV-D) with the Dirichlet boundary condition is proposed. The main idea to derive the scheme is to use the discretization based on the method of characteristics to the convection part, and the finite difference method to the rest of the terms. The nonlocal terms are dealt with the quadrature methods. As a result, an implicit scheme is obtained for the boundary value problem under consideration. The consistency and the convergence of the proposed numerical scheme is established. Moreover, numerical simulations are presented to validate the theoretical results.

Stochastic Hamiltonian partial differential equations, which possess the multi-symplectic conservation law, are an important and fairly large class of systems. The multi-symplectic methods inheriting the geometric features of stochastic Hamiltonian partial differential equations provide numerical approximations with better numerical stability, and are of vital significance for obtaining correct numerical results. In this paper, we propose three novel multi-symplectic methods for stochastic Hamiltonian partial differential equations based on the local radial basis function collocation method, the splitting technique, and the partitioned Runge-Kutta method. Concrete numerical methods are presented for nonlinear stochastic wave equations, stochastic nonlinear Schr\"odinger equations, stochastic Korteweg-de Vries equations and stochastic Maxwell equations. We take stochastic wave equations as examples to perform numerical experiments, which indicate the validity of the proposed methods.

The Green's function approach of Giles and Pierce is used to build the lift and drag based analytic adjoint solutions for the two-dimensional incompressible Euler equations around irrotational base flows. The drag-based adjoint solution turns out to have a very simple closed form in terms of the flow variables and is smooth throughout the flow domain, while the lift-based solution is singular at rear stagnation points and sharp trailing edges owing to the Kutta condition. This singularity is propagated to the whole dividing streamline (comprising the incoming stagnation streamline and the wall) upstream of the rear singularity (trailing edge or rear stagnation point) by the sensitivity of the Kutta condition to changes in the stagnation pressure.

In this article, we present and analyze a stabilizer-free $C^0$ weak Galerkin (SF-C0WG) method for solving the biharmonic problem. The SF-C0WG method is formulated in terms of cell unknowns which are $C^0$ continuous piecewise polynomials of degree $k+2$ with $k\geq 0$ and in terms of face unknowns which are discontinuous piecewise polynomials of degree $k+1$. The formulation of this SF-C0WG method is without the stabilized or penalty term and is as simple as the $C^1$ conforming finite element scheme of the biharmonic problem. Optimal order error estimates in a discrete $H^2$-like norm and the $H^1$ norm for $k\geq 0$ are established for the corresponding WG finite element solutions. Error estimates in the $L^2$ norm are also derived with an optimal order of convergence for $k>0$ and sub-optimal order of convergence for $k=0$. Numerical experiments are shown to confirm the theoretical results.

We couple the L1 discretization for Caputo derivative in time with spectral Galerkin method in space to devise a scheme that solves quasilinear subdiffusion equations. Both the diffusivity and the source are allowed to be nonlinear functions of the solution. We prove method's stability and convergence with spectral accuracy in space. The temporal order depends on solution's regularity in time. Further, we support our results with numerical simulations that utilize parallelism for spatial discretization. Moreover, as a side result we find asymptotic exact values of error constants along with their remainders for discretizations of Caputo derivative and fractional integrals. These constants are the smallest possible which improves the previously established results from the literature.

In this paper, we propose an efficient proper orthogonal decomposition based reduced-order model(POD-ROM) for nonstationary Stokes equations, which combines the classical projection method with POD technique. This new scheme mainly owns two advantages: the first one is low computational costs since the classical projection method decouples the reduced-order velocity variable and reduced-order pressure variable, and POD technique further improves the computational efficiency; the second advantage consists of circumventing the verification of classical LBB/inf-sup condition for mixed POD spaces with the help of pressure stabilized Petrov-Galerkin(PSPG)-type projection method, where the pressure stabilization term is inherent which allows the use of non inf-sup stable elements without adding extra stabilization terms. We first obtain the convergence of PSPG-type finite element projection scheme, and then analyze the proposed projection POD-ROM's stability and convergence. Numerical experiments validate out theoretical results.

In this paper, we propose a cell-free scheme for unmanned aerial vehicle (UAV) base stations (BSs) to manage the severe intercell interference between terrestrial users and UAV-BSs of neighboring cells. Since the cell-free scheme requires enormous bandwidth for backhauling, we propose to use the sub-terahertz (sub-THz) band for the backhaul links between UAV-BSs and central processing unit (CPU). Also, because the sub-THz band requires a reliable line-of-sight link, we propose to use a high altitude platform station (HAPS) as a CPU. At the first time-slot of the proposed scheme, users send their messages to UAVs at the sub-6 GHz band. The UAVs then apply match-filtering and power allocation. At the second time-slot, at each UAV, orthogonal resource blocks are allocated for each user at the sub-THz band, and the signals are sent to the HAPS after analog beamforming. In the HAPS receiver, after analog beamforming, the message of each user is decoded. We formulate an optimization problem that maximizes the minimum signal-to-interference-plus-noise ratio of users by finding the optimum allocated power as well as the optimum locations of UAVs. Simulation results demonstrate the superiority of the proposed scheme compared with aerial cellular and terrestrial cell-free baseline schemes.

北京阿比特科技有限公司