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

Most numerical schemes proposed for solving BGK models for rarefied gas dynamics are based on the discrete velocity approximation. Since such approach uses fixed velocity grids, one must secure a sufficiently large domain with fine velocity grids to resolve the structure of distribution functions. When one treats high Mach number problems, the computational cost becomes prohibitively expensive. In this paper, we propose a velocity adaptation technique in the semi-Lagrangian framework for BGK model. The velocity grid will be set locally in time and space, according to mean velocity and temperature. We apply a weighted minimization approach to impose conservation. We presented several numerical tests that illustrate the effectiveness of our proposed scheme.

相關內容

ACM/IEEE第23屆模型驅動工程語言和系統國際會議,是模型驅動軟件和系統工程的首要會議系列,由ACM-SIGSOFT和IEEE-TCSE支持組織。自1998年以來,模型涵蓋了建模的各個方面,從語言和方法到工具和應用程序。模特的參加者來自不同的背景,包括研究人員、學者、工程師和工業專業人士。MODELS 2019是一個論壇,參與者可以圍繞建模和模型驅動的軟件和系統交流前沿研究成果和創新實踐經驗。今年的版本將為建模社區提供進一步推進建模基礎的機會,并在網絡物理系統、嵌入式系統、社會技術系統、云計算、大數據、機器學習、安全、開源等新興領域提出建模的創新應用以及可持續性。 官網鏈接: · 變換 · 近似 · Performer · 分段 ·
2021 年 9 月 20 日

This paper is concerned with the development of suitable numerical method for the approximation of discontinuous solutions of parameter-dependent linear hyperbolic conservation laws. The objective is to reconstruct such approximation, for new instances of the parameter values for any time, from a transformation of pre-computed snapshots of the solution trajectories for new parameter values. In a finite volume setting, a Reconstruct-Transform-Average (RTA) algorithm inspired from the Reconstruct-Evolve-Average one of Godunov's method is proposed. It allows to perform, in three steps, a transformation of the snapshots with piecewise constant reconstruction. The method is fully detailed and analyzed for solving a parameter-dependent transport equation for which the spatial transformation is related to the characteristic intrinsic to the problem. Numerical results for transport equation and linear elastodynamics equations illustrate the good behavior of the proposed approach.

In this paper, we analyze the convergence of the operator-compressed multiscale finite element method (OC MsFEM) for Schr\"{o}dinger equations with general multiscale potentials in the semiclassical regime. In the OC MsFEM the multiscale basis functions are constructed by solving a constrained energy minimization. Under a mild assumption on the mesh size $H$, we prove the exponential decay of the multiscale basis functions so that localized multiscale basis functions can be constructed, which achieve the same accuracy as the global ones if the oversampling size $m = O(\log(1/H))$. We prove the first-order convergence in the energy norm and second-order convergence in the $L^2$ norm for the OC MsFEM and super convergence rates can be obtained if the solution possesses sufficiently high regularity. By analysing the regularity of the solution, we also derive the dependence of the error estimates on the small parameters of the Schr\"{o}dinger equation. We find that the OC MsFEM outperforms the finite element method (FEM) due to the super convergence behavior for high-regularity solutions and weaker dependence on the small parameters for low-regularity solutions in the presence of the multiscale potential. Finally, we present numerical results to demonstrate the accuracy and robustness of the OC MsFEM.

We show how to combine in a natural way (i.e. without any test nor switch) the conservative and non conservative formulations of an hyperbolic system that has a conservative form. This is inspired from two different class of schemes: the Residual Distribution one \cite{MR4090481}, and the Active Flux formulations \cite{AF1, AF3, AF4,AF5,RoeAF}. The solution is globally continuous, and as in the active flux method, described by a combination of point values and average values. Unlike the "classical" active flux methods, the meaning of the pointwise and cella averaged degrees of freedom is different, and hence follow different form of PDEs: it is a conservative version of the cell average, and a possibly non conservative one for the points. This new class of scheme is proved to satisfy a Lax-Wendroff like theorem. We also develop a method to perform non linear stability. We illustrate the behaviour on several benchmarks, some quite challenging.

In this paper, we propose a unified convergence analysis for a class of generic shuffling-type gradient methods for solving finite-sum optimization problems. Our analysis works with any sampling without replacement strategy and covers many known variants such as randomized reshuffling, deterministic or randomized single permutation, and cyclic and incremental gradient schemes. We focus on two different settings: strongly convex and nonconvex problems, but also discuss the non-strongly convex case. Our main contribution consists of new non-asymptotic and asymptotic convergence rates for a wide class of shuffling-type gradient methods in both nonconvex and convex settings. We also study uniformly randomized shuffling variants with different learning rates and model assumptions. While our rate in the nonconvex case is new and significantly improved over existing works under standard assumptions, the rate on the strongly convex one matches the existing best-known rates prior to this paper up to a constant factor without imposing a bounded gradient condition. Finally, we empirically illustrate our theoretical results via two numerical examples: nonconvex logistic regression and neural network training examples. As byproducts, our results suggest some appropriate choices for diminishing learning rates in certain shuffling variants.

In this work we analyse a PDE-ODE problem modelling the evolution of a Glioblastoma, which includes an anisotropic nonlinear diffusion term with a diffusion velocity increasing with respect to vasculature. First, we prove the existence of global in time weak-strong solutions using a regularization technique via an artificial diffusion in the ODE-system and a fixed point argument. In addition, stability results of the critical points are given under some constraints on parameters. Finally, we design a fully discrete finite element scheme for the model which preserves the pointwise and energy estimates of the continuous problem.

In the first part of this work, we develop a novel scheme for solving nonparametric regression problems. That is the approximation of possibly low regular and noised functions from the knowledge of their approximate values given at some random points. Our proposed scheme is based on the use of the pseudo-inverse of a random projection matrix, combined with some specific properties of the Jacobi polynomials system, as well as some properties of positive definite random matrices. This scheme has the advantages to be stable, robust, accurate and fairly fast in terms of execution time. Moreover and unlike most of the existing nonparametric regression estimators, no extra regularization step is required by our proposed estimator. Although, this estimator is initially designed to work with random sampling set of uni-variate i.i.d. random variables following a Beta distribution, we show that it is still work for a wide range of sampling distribution laws. Moreover, we briefly describe how our estimator can be adapted in order to handle the multivariate case of random sampling sets. In the second part of this work, we extend the random pseudo-inverse scheme technique to build a stable and accurate estimator for solving linear functional regression (LFR) problems. A dyadic decomposition approach is used to construct this last stable estimator for the LFR problem. The performance of the two proposed estimators are illustrated by various numerical simulations. In particular, a real dataset is used to illustrate the performance of our nonparametric regression estimator.

In this paper, we study the connection between entropic optimal transport and entropy power inequality (EPI). First, we prove an HWI-type inequality making use of the infinitesimal displacement convexity of optimal transport map. Second, we derive two Talagrand-type inequalities using the saturation of EPI that corresponds to a numerical term in our expression. We evaluate for a wide variety of distributions this term whereas for Gaussian and i.i.d. Cauchy distributions this term is found in explicit form. We show that our results extend previous results of Gaussian Talagrand inequality for Sinkhorn distance to the strongly log-concave case.

We study the optimization version of the equal cardinality set partition problem (where the absolute difference between the equal sized partitions' sums are minimized). While this problem is NP-hard and requires exponential complexity to solve in general, we have formulated a weaker version of this NP-hard problem, where the goal is to find a locally optimal solution. The local optimality considered in our work is under any swap between the opposing partitions' element pairs. To this end, we designed an algorithm which can produce such a locally optimal solution in $O(N^2)$ time and $O(N)$ space. Our approach does not require positive or integer inputs and works equally well under arbitrary input precisions. Thus, it is widely applicable in different problem scenarios.

This work contributes to the limited literature on estimating the diffusivity or drift coefficient of nonlinear SPDEs driven by additive noise. Assuming that the solution is measured locally in space and over a finite time interval, we show that the augmented maximum likelihood estimator introduced in Altmeyer, Reiss (2020) retains its asymptotic properties when used for semilinear SPDEs that satisfy some abstract, and verifiable, conditions. The proofs of asymptotic results are based on splitting the solution in linear and nonlinear parts and fine regularity properties in $L^p$-spaces. The obtained general results are applied to particular classes of equations, including stochastic reaction-diffusion equations. The stochastic Burgers equation, as an example with first order nonlinearity, is an interesting borderline case of the general results, and is treated by a Wiener chaos expansion. We conclude with numerical examples that validate the theoretical results.

A nonlinear multigrid solver for two-phase flow and transport in a mixed fractional-flow velocity-pressure-saturation formulation is proposed. The solver, which is under the framework of the full approximation scheme (FAS), extends our previous work on nonlinear multigrid for heterogeneous diffusion problems. The coarse spaces in the multigrid hierarchy are constructed by first aggregating degrees of freedom, and then solving some local flow problems. The mixed formulation and the choice of coarse spaces allow us to assemble the coarse problems without visiting finer levels during the solving phase, which is crucial for the scalability of multigrid methods. Specifically, a natural generalization of the upwind flux can be evaluated directly on coarse levels using the precomputed coarse flux basis vectors. The resulting solver is applicable to problems discretized on general unstructured grids. The performance of the proposed nonlinear multigrid solver in comparison with the standard single level Newton's method is demonstrated through challenging numerical examples. It is observed that the proposed solver is robust for highly nonlinear problems and clearly outperforms Newton's method in the case of high Courant-Friedrichs-Lewy (CFL) numbers.

北京阿比特科技有限公司