High order fast sweeping methods for efficiently solving steady state solutions of hyperbolic PDEs were not available yet on unstructured meshes. In this paper, we extend high order fast sweeping methods to unstructured triangular meshes by applying fixed-point iterative sweeping techniques to a fifth-order finite volume unstructured WENO scheme, for solving steady state solutions of hyperbolic conservation laws. An advantage of fixed-point fast sweeping methods which distinguishes them from other fast sweeping methods is that they are explicit and do not require inverse operation of nonlinear local systems. As in the first order fast sweeping methods on triangular meshes, we introduce multiple reference points to determine alternating sweeping directions on unstructured meshes. All the cells on the mesh are ordered according to their centroids' distances to those reference points, and the resulted orderings provide sweeping directions for iterations. To make the residue of the fast sweeping iterations converge to machine zero / round off errors, we follow the approach in our early work of developing the absolutely convergent fixed-point fast sweeping WENO methods on rectangular me
In this paper we consider the Moreau's sweeping processes driven by a time dependent prox-regular set $C(t)$ which is continuous in time with respect to the asymmetric distance $e$ called the excess, defined by $e(A,B) := \sup_{x \in A} d(x,B)$ for every pair of sets $A$, $B$ in a Hilbert space. As observed by J.J. Moreau in his pioneering works, the excess provides the natural topological framework for sweeping process. Assuming a uniform interior cone condition for $C(t)$, we prove that the associated sweeping process has a unique solution, thereby improving the existing result on continuous prox-regular sweeping processes in two directions: indeed, in the previous literature $C(t)$ was supposed to be continuous in time with respect to the symmetric Hausdorff distance instead of the excess and also its boundary $\partial C(t)$ was required to be continuous in time, an assumption which we completely drop. Therefore our result allows to consider a much wider class of continuously moving constraints.
In this paper, the study of nonsmooth optimal control problems (P) involving a controlled sweeping process with three main characteristics is launched. First, the sweeping sets are nonsmooth, time-dependent, and uniformly prox-regular. Second, the sweeping process is coupled with a controlled differential equation. Third, a joint-state endpoints constraint set S is present. This general model incorporates different important controlled submodels, such as a class of second order sweeping processes, and coupled evolution variational inequalities. A full form of the nonsmooth Pontryagin maximum principle for strong local minimizers in (P) is derived for bounded or unbounded moving sweeping sets satisfying local constraint qualifications (CQ) without any additional restriction. The existence and uniqueness of a Lipschitz solution for the Cauchy problem of our dynamic is established and the existence of an optimal solution for (P) is obtained. Two of the novelties in achieving the first goal are (i) the construction of a problem over truncated sweeping sets and truncated joint endpoints constraint set that has the same strong local minimizer as (P) and its (CQ) automatically holds, and
Standard sweep algorithms require an order of discrete points in Euclidean space, and rely on the property that, at a given point, all points in the halfspace below come earlier in this order. We are motivated by the problem of reconstructing a graph in $\mathbb{R}^d$ from vertex locations and degree information, which was addressed using standard sweep algorithms by Fasy et al. We generalize this to the reconstruction of general simplicial complexes. As our main ingredient, we introduce a generalized \emph{sweeping order} on $i$-simplices, maintaining the property that, at a given $i$-simplex $σ$, all $(i+1)$-dimensional cofaces of $σ$ in the halfspace below $σ$ have an $i$-dimensional face that appeared earlier in the order ("below" with respect to some direction perpendicular to $σ$). We then go on to incorporate computing such sweeping orders to reconstruct an unknown simplicial complex $K$, starting with only its vertex locations, i.e., its $0$-skeleton. Specifically, once we have found the $i$-skeleton of $K$, we compute a sweeping order for the $i$-simplices, and use it to reconstruct the $(i+1)$-skeleton of $K$ by querying the \emph{indegree}, or the number of $(i+1)$-simpl
Consider a given planar circular region, in which there is an unknown number of smart mobile evaders. We wish to detect evaders using a line formation of sweeping agents whose total sensing length is predetermined. We propose procedures for designing spiral sweeping protocols that ensure the successful completion of the task, thus deriving conditions on the sweeping speed of the linear formation and its path. Successful completion of the task implies that evaders with a given limit on their speed cannot escape the sweeping agents. A simpler task for the sweeping formation is the confinement of evaders to a desired region, such as their original domain. The feasibility of completing these tasks depends on geometric and dynamic constraints that impose a lower bound on the speed that the sweeping agents must have. This critical speed is derived to ensure the satisfaction of the confinement task. Increasing the speed above the lower bound enables the sweepers to complete the search task as well. We develop two spiral line formation search processes for smart evaders, that address current limitations in search against smart evaders. Additionally, we present a quantitative and qualitativ
This paper addresses, for the first time in the literature, optimal control problems for dynamic systems governed by a novel class of sweeping processes with time delay. We establish well-posedness of such processes, in the sense of the existence and uniqueness of feasible trajectories corresponding to feasible controls under fairly unrestrictive assumptions. Then we construct a well-posed family of discrete approximations and find efficient conditions under the discretized time-delayed sweeping process exhibits stability with respect to strong convergence of feasible and optimal solutions. This creates a bridge between optimization of continuous-time and discrete-time sweeping control systems and justifies the effective use of discrete approximations in deriving optimality conditions and numerical techniques to solve the original time-delayed sweeping control problems via discrete approximations.
Fixed-point fast sweeping methods are a class of explicit iterative methods developed in the literature to efficiently solve steady state solutions of hyperbolic partial differential equations (PDEs). As other types of fast sweeping schemes, fixed-point fast sweeping methods use the Gauss-Seidel iterations and alternating sweeping strategy to cover characteristics of hyperbolic PDEs in a certain direction simultaneously in each sweeping order. The resulting iterative schemes have fast convergence rate to steady state solutions. Moreover, an advantage of fixed-point fast sweeping methods over other types of fast sweeping methods is that they are explicit and do not involve inverse operation of any nonlinear local system. Hence they are robust and flexible, and have been combined with high order accurate weighted essentially non-oscillatory (WENO) schemes to solve various hyperbolic PDEs in the literature. For multidimensional nonlinear problems, high order fixed-point fast sweeping WENO methods still require quite large amount of computational costs. In this technical note, we apply sparse-grid techniques, an effective approximation tool for multidimensional problems, to fixed-point
In this paper, we study the well-posedness of state-dependent and state-independent sweeping processes driven by prox-regular sets and perturbed by a history-dependent operator. Our approach, based on an enhanced version of Gronwall's lemma and fixed-point arguments, provides an efficient framework for analyzing sweeping processes. In particular, our findings recover all existing results for the class of Volterra sweeping processes and provide new insights into history-dependent sweeping processes. Finally, we apply our theoretical results to establish the well-posedness of a viscoelastic model with long memory.
In recent research, the parallel performances of sweeping-type algorithms for high-frequency time-harmonic wave problems have been improved by departing from standard layer-type domain decomposition and introducing a new sweeping strategy on a checkerboard-type domain decomposition, where sweeps can be performed more flexibly. These sweeps can be done by a certain number of steps, each of which provides the necessary information from subdomains solved at the current iteration to their next neighboring subdomains. Although, subproblems in these subdomains can be solved concurrently at each step, the sequential nature of the process of the sweeping approaches still exists, which limits their potential for parallelization. We propose a block Jacobi sweeping preconditioner, which is an improved variant of sweeping-type preconditioners. The new feature of these improved variants can be interpreted as several partial sweeps, which can be thought of as sweeps that operate on a subset of the subdomains in parallel. We present several two-dimensional finite element results to study and compare the sweeping preconditioner and the block Jacobi sweeping preconditioner.
We study the problem of sweeping a pseudoline arrangement with $n$ $x$-monotone curves with a rope (an $x$-monotone curve that connects the points at infinity). The rope can move by flipping over a face of the arrangement, replacing parts of it from the lower to the upper chain of the face. Counting as length of the rope the number of edges, what rope-length can be needed in such a sweep? We show that all such arrangements can be swept with rope-length at most $2n-2$, and for some arrangements rope-length at least $7(n-2)/4+1$ is required. We also discuss some complexity issues around the problem of computing a sweep with the shortest rope-length.
In this paper, we study the well-posedness of integro-differential sweeping processes of Volterra type. Using new enhanced versions of Gronwall's inequality, a reparametrization technique, and a fixed point argument for history-dependent operators, we obtain the existence of solutions and provide a fully continuous dependence result for the integro-differential sweeping process. The paper ends with an application to projected dynamical systems.
The numerical method developed in [30] for optimal control problems involving sweeping processes with smooth sweeping set C is generalized to the case where C is nonsmooth, namely, C is the intersection of a finite number of sublevel sets of smooth functions. The novelty of this extension resides in producing for the general setting a different approach, since the one used for the smooth sweeping sets is not applicable here.
In this study, we investigate optimal control problems that involve sweeping processes with a drift term and mixed inequality constraints. Our goal is to establish necessary optimality conditions for these problems. We address the challenges that arise due to the combination of sweeping processes and inequality mixed constraints in two contexts: regular and non-regular. This requires working with different types of multipliers, such as finite positive Radon measures for the sweeping term and integrable functions for regular mixed constraints. For non-regular mixed constraints, the multipliers correspond to purely finitely additive set functions.
We investigate the spin tunnelling of spin-1 Bose-Einstein condensates in a linearly sweeping magnetic field with a mean-field treatment. We focus on the two typical alkali Bose atoms $^{87}Rb$ and $^{23}Na$ condensates and study their tunnelling dynamics according to different sweeping rates of external magnetic fields. In the adiabatic (i.e., slowly sweeping) and sudden (i.e., fast sweeping) limits, no tunnelling is observed. For the case of moderate sweeping rates, the tunnelling dynamics is found to be very sensitive on the sweeping rates with showing a chaotic-like tunnelling regime. With magnifying the regime, however, we find interestedly that the plottings become resolvable under a resolution of $10^{-4}$ G/s where the tunnelling probability with respect to the sweeping rate shows a regular periodic-like pattern. Moreover, a conserved quantity standing for the magnetization in experiments is found can dramatically affect the above picture of the spin tunnelling. Theoretically we have given a reasonable interpretation to the above findings and hope our studies would bring more attention to spin tunnelling experimentally.
Suppose that in a given planar circular region, there are some smart mobile evaders and we would like to find them using sweeping agents. We assume that each agent has a line sensor of length 2r. We propose procedures for designing cooperative sweeping processes that ensure the successful completion of the task, thereby deriving conditions on the sweeping velocity of the agents and their paths. Successful completion of the task means that evaders with a given limit on their velocity cannot escape the sweeping agents. A simpler task for the sweeping swarm is the confinement of the evaders to their initial domain. The feasibility of completing these tasks depends on geometric and dynamic constraints that impose a lower bound on the velocity that the sweeper swarm must have. This critical velocity is derived to ensure the satisfaction of the confinement task. Increasing the velocity above the lower bound enables the agents to complete the search task as well. We present results on the total search time as a function of the sweeping velocity of the swarm's agents given the initial conditions on the size of the search region and the maximal velocity of the evaders.
The aim of this paper is twofold. On one hand we prove that the Moreau's sweeping process driven by a uniformly prox-regular moving set with local bounded retraction has a unique solution provided that the coefficient of prox-regularity is larger than the size of any jump of the driving set. On the other hand we show how the case of local bounded retraction can be easily reduced to the $1$-Lipschitz continuous case: indeed we first solve the Lipschitz continuous case by means of the so called ``catching-up algorithm", and we reduce the local bounded retraction case to the Lipschitz one by using a reparametrization technique for functions with values in the family of prox-regular sets.
Fixed-point iterative sweeping methods were developed in the literature to efficiently solve steady state solutions of Hamilton-Jacobi equations and hyperbolic conservation laws. Similar as other fast sweeping schemes, the key components of this class of methods are the Gauss-Seidel iterations and alternating sweeping strategy to achieve fast convergence rate. Furthermore, good properties of fixed-point iterative sweeping methods include that they have explicit forms and do not involve inverse operation of nonlinear local systems, and they can be applied to general hyperbolic equations using any monotone numerical fluxes and high order approximations easily. In [L. Wu, Y.-T. Zhang, S. Zhang and C.-W. Shu, Commun. Comput. Phys., 20 (2016)], a fifth order fixed-point sweeping WENO scheme was designed and it was shown that the scheme converges much faster than the total variation diminishing (TVD) Runge-Kutta approach by stability improvement of high order schemes with a forward Euler time-marching. An open problem is that for some benchmark numerical examples, the iteration residue of the fixed-point sweeping WENO scheme hangs at a truncation error level instead of settling down to m
This paper introduces the recursive sweeping preconditioner for the numerical solution of the Helmholtz equation in 3D. This is based on the earlier work of the sweeping preconditioner with the moving perfectly matched layers (PMLs). The key idea is to apply the sweeping preconditioner recursively to the quasi-2D auxiliary problems introduced in the 3D sweeping preconditioner. Compared to the non-recursive 3D sweeping preconditioner, the setup cost of this new approach drops from $O(N^{4/3})$ to $O(N)$, the application cost per iteration drops from $O(N\log N)$ to $O(N)$, and the iteration count only increases mildly when combined with the standard GMRES solver. Several numerical examples are tested and the results are compared with the non-recursive sweeping preconditioner to demonstrate the efficiency of the new approach.
We consider a Caratheodory differential equation with a state-dependent convex constraint that changes BV-continuously in time (a perturbed BV-continuous state-dependent sweeping processes). By setting up an appropriate catching-up algorithm we prove solvability of the initial value problem. Then, for sweeping processes with $T$-periodic right-hand-sides, we prove the existence of at least one $T$-periodic solution. Finally, we further consider a $T$-periodic sweeping process which is close to an autonomous sweeping process with a constant constraint and prove the existence of a $T$-periodic solution specifically located near the boundary switched equilibrium of the autonomous sweeping process.
This paper introduces sweeping permutation automata, which move over an input string in alternating left-to-right and right-to-left sweeps and have a bijective transition function. It is proved that these automata recognize the same family of languages as the classical one-way permutation automata (Thierrin, "Permutation automata", Mathematical Systems Theory, 1968). An n-state two-way permutation automaton is transformed to a one-way permutation automaton with F(n)=\max_(k+l=n, m <= l) k (l \choose m) (k - 1 \choose l - m) (l - m)! states. This number of states is proved to be necessary in the worst case, and its growth rate is estimated as F(n) = n^(n/2 - (1 + \ln 2)/2 \cdot n/(\ln n) \cdot (1 + o(1))).