Optimal Estimates Research Articles

A novel discontinuous Galerkin projection scheme for the hydrodynamics of nematic liquid crystals

This paper is focused on the numerical approximations for the hydrodynamic model of nematic liquid crystals. Under the framework of a splitting projection method, we propose a novel interior penalty discontinuous Galerkin (DG) method for solving the coupled system, which is employed by combining the scalar auxiliary variables (SAV) approach, implicit-explicit (IMEX) treatments and a rotational pressure-correction method. One prominent feature of the developed scheme here is by introducing an additional stabilization term artificially in liquid crystal equation to balance the explicit treatment for the coupling term, so that the computations of vector field from velocity field are decoupled. Hence a linear, fully decoupled and unconditionally energy-stable DG scheme can be achieved in a fully discrete manner. We show that the resulting scheme is uniquely solvable and unconditional energy stable, and the optimal error estimates of the proposed scheme are further proved theoretically. Finally, numerical results are carried out to demonstrate the accuracy and energy stability of our scheme.

Soap bubbles and convex cones: optimal quantitative rigidity

We consider a class of recent rigidity results in a convex cone Σ ⊆ R N \Sigma \subseteq \mathbb {R}^N . These include overdetermined Serrin-type problems for a mixed boundary value problem relative to Σ \Sigma , Alexandrov’s soap bubble-type results relative to Σ \Sigma , and Heintze-Karcher’s inequality relative to Σ \Sigma . Each rigidity result is obtained here by means of a single integral identity and holds true under weak integral overdeterminations in possibly non-smooth cones. Optimal quantitative stability estimates are obtained in terms of an L 2 L^2 -pseudodistance. In particular, the optimal stability estimate for Heintze-Karcher’s inequality is new even in the classical case Σ = R N \Sigma = \mathbb {R}^N . Stability bounds in terms of the Hausdorff distance are also provided. Several new results are established and exploited, including a new Poincaré-type inequality for vector fields whose normal components vanish on a portion of the boundary and an explicit (possibly weighted) trace theory – relative to the cone Σ \Sigma – for harmonic functions satisfying a homogeneous Neumann condition on the portion of the boundary contained in ∂ Σ \partial \Sigma . We also introduce new notions of uniform interior and exterior sphere conditions relative to the cone Σ ⊆ R N \Sigma \subseteq \mathbb {R}^N , which allow to obtain (via barrier arguments) uniform lower and upper bounds for the gradient in the mixed boundary value-setting. In the particular case Σ = R N \Sigma =\mathbb {R}^N , these conditions return the classical uniform interior and exterior sphere conditions (together with the associated classical gradient bounds of the Dirichlet setting).

Error estimates and numerical simulations of a thermoviscoelastic contact problem with damage and long memory

This paper aims to investigate a thermal frictional contact model with damage and long memory effects. We consider a deformable body made of viscoelastic material and assume the process to be dynamic. The material is expected to adhere to the Kelvin–Voigt constitutive law, with damage and thermal effects incorporated. The variational formulation of the model results in a coupled system comprising a history-dependent hemivariational inequality governing the displacement field, a parabolic variational inequality describing the damage field and an evolution equation for the temperature field. In the analysis of this system, we initially introduce a fully discrete scheme, and then concentrate on deriving error estimates of numerical solutions. An optimal order error estimate is attained under some appropriate solution regularity assumptions. At the tail of this manuscript, numerical simulations are provided for the contact problem to validate the theoretical results.

Designing optimal protocols in Bayesian quantum parameter estimation with higher-order operations

Using quantum systems as sensors or probes has been shown to greatly improve the precision of parameter estimation by exploiting unique quantum features such as entanglement. A major task in quantum sensing is to design the optimal protocol, i.e., the most precise one. It has been solved for some specific instances of the problem, but in general even numerical methods are not known. Here, we focus on the single-shot Bayesian setting, where the goal is to find the optimal initial state of the probe (which can be entangled with an auxiliary system), the optimal measurement, and the optimal estimator function. We leverage the formalism of higher-order operations to develop a method based on semidefinite programming that finds a protocol that is close to the optimal one with arbitrary precision. Crucially, our method is not restricted to any specific quantum evolution, cost function, or prior distribution, and thus can be applied to any estimation problem. Moreover, it can be applied to both single-parameter or multiparameter estimation tasks. We demonstrate our method with three examples, consisting of unitary phase estimation, thermometry in a bosonic bath, and multiparameter estimation of an SU(2) transformation. Exploiting our methods, we extend several results from the literature. For example, in the thermometry case, we find the optimal protocol at any finite time and quantify the usefulness of entanglement. Published by the American Physical Society 2024

Frisch–Waugh–Lovell theorem-type results for the k-Class and 2SGMM estimators

The Frisch–Waugh–Lovell (FWL) theorem shows that for the least squares estimator, parameter estimates from full and partial models are identically same. I show that in linear regression models with a mix of exogenous and endogenous regressors, FWL theorem-type results hold for the k-class estimators (including LIML) and the two-step optimal GMM estimator.

Translation-Invariant Equations With at Least Four Variables

Abstract We prove that every subset of $\{1,\dots , N\}$ that does not contain any solutions to a translation invariant equation $a_{1}x_{1}+\dots +a_{k}x_{k}=0$ with $k\ge 4$ has at most $$ \begin{align*} & \exp\big(-c(\log N)^{1/5+o(1)}\big)N\end{align*} $$ elements, for some $c>0$. This theorem improves upon previous estimates. Additionally, our method has the potential to yield an optimal estimate for this problem that matches Behrend’s classical lower bound. Our approach relies on a new result on almost-periodicity of convolutions.

An energy-stable variable-step L1 scheme for time-fractional Navier–Stokes equations

We propose a structure-preserving scheme and its error analysis for time-fractional Navier–Stokes equations (TFNSEs) with periodic boundary conditions. The equations are first rewritten as an equivalent system by eliminating the pressure explicitly. Then, the spatial and temporal discretization are done by the Fourier spectral method and variable-step L1 scheme, respectively. It is proved that the fully-discrete scheme is energy-stable and divergence-free. The energy is an asymptotically compatible one since it recovers the classical energy when α→1. Moreover, optimal error estimates are presented very technically by the obtained boundedness of the numerical solutions and some Sobolev inequalities. To our knowledge, they are the first results of the construction and analysis of structure-preserving schemes for TFNSEs. Several interesting numerical examples are given to confirm the theoretical results at last.

MAXIMUM PROFIT ANALYSIS USING LINEAR PROGRAMMING SIMPLEX METHOD AND POM-QM SOFTWARE AT UKM PIE BU SRI

In current conditions, many people are competing to set up businesses to meet their daily living needs. One of the various ways to improve people's welfare is through small and medium enterprises (SMEs). As many businesses develop and are accompanied by intense competition, many problems will arise and will affect the profits of SMEs. The problem that occurs at UKM Pie Bu Sri is the problem of calculating the optimum profit obtained every day from the production of its activities. The aim of this research is to help calculate the optimum profit for UKM Pie Bu Sri so that it can be done accurately and quickly. In order to achieve this goal, the Simplex Linear Programming Method and QM for Windows Software were used to find the optimum profit estimate obtained in each production activity carried out by UKM Pie Bu Sri within a period of one day. The results of the maximization calculations show that in order to achieve maximum profits, Mrs. Sri must make 60 Brownie Pies and produce 25 Fruit Pies with a profit of IDR 99,500 per day.

Two-grid methods for nonlinear pseudo-parabolic integro-differential equations by finite element method

In this paper, two efficient two-grid finite element algorithms are proposed for solving two-dimensional nonlinear pseudo-parabolic integro-differential equations. Firstly, we obtain optimal error estimates of the fully discrete finite element method using a temporal-spatial error splitting technique in H1 and Lp norms. Then the two-grid technique to improve computation efficiency of the proposed finite element method. Error estimates in H1 and Lp norms of two-grid solutions are presented. Theoretical analysis shows that the two-grid algorithms maintain asymptotically optimal accuracy. Finally, numerical examples are provided to support our theoretical results and demonstrate the effectiveness of these methods.

Precise error bounds for numerical approximations of fractional HJB equations

Abstract We prove precise rates of convergence for monotone approximation schemes of fractional and nonlocal Hamilton–Jacobi–Bellman equations. We consider diffusion-corrected difference-quadrature schemes from the literature and new approximations based on powers of discrete Laplacians, approximations that are (formally) fractional order and second-order methods. It is well known in numerical analysis that convergence rates depend on the regularity of solutions, and here we consider cases with varying solution regularity: (i) strongly degenerate problems with Lipschitz solutions and (ii) weakly nondegenerate problems where we show that solutions have bounded fractional derivatives of order $\sigma \in (1,2)$. Our main results are optimal error estimates with convergence rates that capture precisely both the fractional order of the schemes and the fractional regularity of the solutions. For strongly degenerate equations, these rates improve earlier results. For weakly nondegenerate problems of order greater than one, the results are new. Here we show improved rates compared to the strongly degenerate case, rates that are always better than $\mathcal{O}\big (h^{\frac{1}{2}}\big )$.

An efficient weak Galerkin FEM for third-order singularly perturbed convection-diffusion differential equations on layer-adapted meshes

In this article, we study the weak Galerkin finite element method to solve a class of a third order singularly perturbed convection-diffusion differential equations. Using some knowledge on the exact solution, we prove a robust uniform convergence of order O(N−(k−1/2)) on the layer-adapted meshes including Bakhvalov-Shishkin type, and Bakhvalov-type and almost optimal uniform error estimates of order O((N−1ln⁡N)(k−1/2)) on Shishkin-type mesh with respect to the perturbation parameter in the energy norm using high-order piecewise discontinuous polynomials of degree k. Here N is the number mesh intervals. We conduct numerical examples to support our theoretical results.

Prediction problem for continuous time stochastic processes with periodically correlated increments observed with noise

We propose solution of the problem of the mean square optimal estimation of linear functionals which depend on the unobserved values of a continuous time stochastic process with periodically correlated increments based on observations of this process with periodically stationary noise. To solve the problem, we transform the processes to the sequences of stochastic functions which form an infinite dimensional vector stationary sequences. In the case of known spectral densities of these sequences, we obtain formulas for calculating values of the mean square errors and the spectral characteristics of the optimal estimates of the functionals. Formulas determining the least favorable spectral densities and the minimax (robust) spectral characteristics of the optimal linear estimates of functionals are derived in the case where the sets of admissible spectral densities are given.

A second-order Strang splitting scheme for the generalized Allen–Cahn type phase-field crystal model with FCC ordering structure

In this paper, we consider the generalized Allen–Cahn-type phase-field crystal model with face-centered-cubic ordering structure (PFC-FCC). Due to the combined complexity of the eighth-order spatial derivative and inherent nonlinearity, it poses a significant challenge to design a numerical scheme of high accuracy, stability, and efficiency to solve the PFC-FCC model. Endeavoring towards this objective, we adopt operator splitting method to address the PFC-FCC model. Based on Fourier spectral method for spatial discretization and SSP-RK method for temporal discretization, we propose an effective and easy-to-implement second-order scheme. We are also able to proffer an analytical proof for discrete mass conservation, and engage a discussion on an optimal error estimate for the second-order scheme. In addition, the energy stability of the numerical scheme is derived. Ultimately, a series of numerical experiments specifically designed to substantiate its accuracy, efficiency, and capability for phase transition are performed.

The forward self-similar solution of fractional incompressible Navier-Stokes system: The critical case

In this paper, we study the regularity and pointwise estimates of forward self-similar solutions of fractional Navier-Stokes system under the critical case. By employing a Caffarelli,Kohn and Nirenberg-type iteration, L∞ estimates of the self-similar solution's profile are established, which is a key ingredient to ensure that the global weighted energy estimate procedure used in [28] is performed under the critical case. As a product, its natural pointwise bounds are recovered. Moreover, to obtain the optimal spatial decay estimate of self-similar solution's profile, a new technique is required due to lack of the related regularity theory.

A structure-preserving local discontinuous Galerkin method for the stochastic KdV equation

This paper proposes a local discontinuous Galerkin (LDG) method for the stochastic Korteweg-de Vries (KdV) equation with multi-dimensional multiplicative noise. In the mean square sense, we show that the numerical method is L2 stable and it preserves energy conservation and energy dissipation. If the degree of the polynomial is n, the optimal error estimate in the mean square sense can reach as n+1. Finally, structure-preserving and convergence are verified by numerical experiments.

Distributed Cooperative Quantum Learning for Discrete-Time Multiagent Source Exploration With Information Prompts.

In light of optimization theory and swarm evolutionary schemes, under multiple single-integrator mobile agents equipped with sensors and prompters, this article addresses a discrete-time multiagent source exploration problem with information prompts. Regarding information prompts as constraints on the unknown target, by virtue of penalty function skills (PFSs) and sequential unconstrained minimization techniques (SUMTs), the agents are driven toward the source under the guidance of the control strategy. In two cases of available and unavailable gradient information, a quantum potential well, an average optimal position estimator (AOPE), and a global optimal position estimator (GOPE) are introduced into swarm evolutionary schemes with a periodically oscillating weight, such that distributed cooperative quantum learning (DCQL) policy is proposed as a control strategy under communication restrictions, where AOPE and GOPE are developed relying on distributed consensus theory. In particular, when the gradient is unavailable, we put forth an adaptive generalized Bernstein neural network (AGBNN) to replace it based on excellent properties of Bernstein polynomials and adaptive approaches. Further, a performance analysis for the proposed policy is executed on the convergence and computational complexity, which ensures the accuracy and efficiency of the source exploration in theory. Ultimately, a simulation test is carried out, and the results validate the practicability and effectiveness of the offered method.

Probing the galactic and extragalactic gravitational wave backgrounds with space-based interferometers

We employ the formalism developed in [1] and [2] to study the prospect of detecting an anisotropic Stochastic Gravitational Wave Background (SGWB) with the Laser Interferometer Space Antenna (LISA) alone, and combined with the proposed space-based interferometer Taiji. Previous analyses have been performed in the frequency domain only. Here, we study the detectability of the individual coefficients of the expansion of the SGWB in spherical harmonics, by taking into account the specific motion of the satellites. This requires the use of time-dependent response functions, which we include in our analysis to obtain an optimal estimate of the anisotropic signal. We focus on two applications. Firstly, the reconstruction of the anisotropic galactic signal without assuming any prior knowledge of its spatial distribution. We find that both LISA and LISA with Taiji cannot put tight constraints on the harmonic coefficients for realistic models of the galactic SGWB. We then focus on the discrimination between a galactic signal of known morphology but unknown overall amplitude and an isotropic extragalactic SGWB component of astrophysical origin. In this case, we find that the two surveys can confirm, at a confidence level ≳ 3σ, the existence of both the galactic and extragalactic background if both have amplitudes as predicted in standard models. We also find that, in the LISA-only case, the analysis in the frequency domain (under the assumption of a time average of data taken homogeneously across the year) provides a nearly identical determination of the two amplitudes as compared to the optimal analysis.

Global well-posedness and large time behavior for the 3D magneto-micropolar equations

This paper focuses on the global well-posedness and optimal decay estimates to the 3D magneto-micropolar equations in Besov spaces. We obtain the global existence and uniqueness of solutions in Besov spaces. Moreover, the decay estimates of the global solutions are established provided that the initial data belongs to the negative Besov space.

A parareal exponential integrator finite element method for semilinear parabolic equations

AbstractIn this article, we present a parareal exponential finite element method, with the help of variational formulation and parareal framework, for solving semilinear parabolic equations in rectangular domains. The model equation is first discretized in space using the finite element method with continuous piecewise multilinear rectangular basis functions, producing the semi‐discrete system. We then discretize the temporal direction using the explicit exponential Runge–Kutta approach accompanied by the parareal framework, resulting in the fully‐discrete numerical scheme. To further improve computational speed, we design a fast solver for our method based on tensor product spectral decomposition and fast Fourier transform. Under certain regularity assumption, we successfully derive optimal error estimates for the proposed parallel‐based method with respect to ‐norm. Extensive numerical experiments in two and three dimensions are also carried out to validate the theoretical results and demonstrate the performance of our method.

On the Cauchy problem for acoustic waves in hereditary fluids: Decay properties and inviscid limits

This paper considers the viscous/inviscid Moore–Gibson–Thompson (MGT) equations with memory of type I in the whole space . For one thing, associating with a new condition on initial data, we derive the optimal estimates and the optimal leading term of the acoustic velocity potential for large time, where we analyze different contributions from viscous, thermally relaxing, as well as hereditary fluids on large time asymptotic behavior for the acoustic waves models. For another, by using the multiscale analysis and energy methods in the Fourier space, we demonstrate the inviscid limits (i.e., as the diffusivity of sound tends to zero), which match our Wentzel‐Kramers‐Brillouin (WKB) expansion of the solution. Finally, we give a further application of our results on large time behavior for the nonlinear Jordan‐MGT equation in viscous hereditary fluids.