On the homogenization of some linear problems in domains weakly connected by a system of traps

The aim of the paper is to study the asymptotic behaviour of solutions of second‐order elliptic and parabolic equations, arising in modelling of flow in cavernous porous media, in a domain Ω^ε weakly connected by a system of traps ??^ε, where ε is the parameter that characterizes the scale of the microstructure. Namely, we consider a strongly perforated domain Ω^ε ?Ω a bounded open set of ?³ such that Ω^ε =Ω₁^ε ∪Ω₂^ε ∪??^ε ∪ W^ε, where Ω₁^ε, Ω₂^ε are non‐intersecting subdomains strongly connected with respect to Ω, ??^ε is a system of traps and meas W^ε → 0 as ε → 0. Without any periodicity assumption, for a large range of perforated media and by means of variational homogenization, we find the homogenized models. The effective coefficients are described in terms of local energy characteristics of the domain Ω^ε associated with the problem under consideration. The resulting homogenized problem in the parabolic case is a vector model with memory terms. An example is presented to illustrate the methodology. Copyright © 2007 John Wiley & Sons, Ltd.     

ΓD‐convergence for a class of quasilinear elliptic equations in thin structures

We study the asymptotic behaviour of the solution of a stationary quasilinear elliptic problem posed in a domain Ω^(ε) of asymptotically degenerating measure, i.e. meas Ω^(ε) → 0 as ε → 0, where ε is the parameter that characterizes the scale of the microstructure. We obtain the convergence of the solution and the homogenized model of the problem is constructed using the notion of convergence in domains of degenerating measure. Proofs are given using the method of local characteristics of the medium Ω^(ε) associated with our problem in a variational form. Copyright © 2005 John Wiley & Sons, Ltd.     

Asymptotic limit of initial boundary value problems for conservation laws with relaxational extensions

We study the boundary layer effect in the small relaxation limit to the equilibrium scalar conservation laws in one space dimension for the relaxation system proposed in [6]. First, it is shown that for initial and boundary data satisfying a strict version of the subcharacteristic condition, there exists a unique global (in time) solution, (u^ε, v^ε), to the relaxation system (1.4) for each ε > 0. The spatial total variation of (u^ε, v^ε) is shown to be bounded independently of ε, and consequently, a subsequence of (u^ε, v^ε) converges to a limit (u, v) as ε → 0⁺. Furthermore, u(x, t) is a weak solution to the scalar conservation law (1.5) and v = f(u). Next, we prove that for data that are suitably small perturbations of a nontransonic state, the relaxation limit function satisfies the boundary-entropy condition (2.11). Finally, the weak solutions to (1.5) with the boundary-entropy condition (2.11) is shown to be unique. Consequently, the relaxation limit of solutions to (1.4) is unique, and the whole sequence converges to the unique limit. One consequence of our analysis shows that the boundary layer occurs only in the u-component in the sense that v^ε(0, ·) converges strongly to γ ○ v = f(γ ○ u), the trace of f(u) on the t-axis. © 1998 John Wiley & Sons, Inc.     

Homogenization of eigenvalue problem for Laplace–Beltrami operator on Riemannian manifold with complicated ‘bubble‐like’ microstructure

We study the asymptotic behavior of the eigenvalues and the eigenfunctions of the Laplace–Beltrami operator on a Riemannian manifold M^ε depending on a small parameter ε>0 and whose structure becomes complicated as ε→0. Under a few assumptions on scales of M^ε we obtain the homogenized eigenvalue problem. In addition we study the behavior of the heat equation on M^ε and investigate the large time behavior of the homogenized equation. Copyright © 2009 John Wiley & Sons, Ltd.     

RADICAL THEORY FOR SEMIRINGS

Abstract

In this article we continue investigations on a Kurosh-Amitsur radical theory for a universal class U of hemirings as introduced by O.M. Olson et al. We give some necessary and sufficient conditions that such a universal class U consists of all hemirings. Further we consider special and weakly special subclasses M of U which yield hereditary radical classes P = um of U. In this context we correct some statements in the papers of Olson et al. Moreover, a problem posed there concerning the equality of two radicals ?_∞(S) and ?^ε(S) and two similar ideals β_∞ (S) and β^ε(S) is widely solved. We prove ?_∞(S) ? ?^ε(S) = β^∞(S) = β^ε(S) and give necessary and sufficient conditions for equality in the first inclusion. This yields in particular that the weakly special class M_ε(U) is always semisimple, a result which is not true for the special class M_∞(U).     

Averaging of a 3D primitive equations with oscillating external forces

In this article, we consider a non-autonomous three-dimensional primitive equations of the ocean with a singularly oscillating external force g ^ε?=?g ₀(t)?+?ε^?ρ g ₁(t/ε) depending on a small parameter ε?>?0 and ρ?∈?[0,?1) together with the averaged system with the external force g ₀(t), formally corresponding to the case ε?=?0. Under suitable assumptions on the external force, we prove as in [V.V. Chepyzhov, V. Pata, and M.M.I. Vishik, Averaging of 2D Navier–Stokes equations with singularly oscillating forces, Nonlinearity, 22 (2009), pp. 351–370] the boundness of the uniform global attractor 𝒜^ε as well as the convergence of the attractors 𝒜^ε of the singular systems to the attractor 𝒜⁰ of the averaged system as ε?→?0⁺. When the external force is small enough and the viscosity is large enough, the convergence rate is controlled by Kε^(1?ρ). Let us note that the main difference between this work and that of Chepyzhov et al. (2009) is that the non-linearity involved in the three-dimensional primitive equation is stronger than the one in the two-dimensional Navier–Stokes equations considered in Chepyzhov et al. (2009), which makes the analysis of the problem studied in this article more involved.     

Stochastic Differential Games with Multiple Modes and a Small Parameter

Abstract

Two persons stochastic differential games with multiple modes where the system is driven by a wideband noise is considered. The state of the system at time t is given by a pair (x ^ε(t), θ^ε(t)), where θ^ε(t) takes values in S = {1, 2,…, N} and ε is a small parameter. The discrete component θ^ε(t) describes various modes of the system. The continuous component x ^ε(t) is governed by a “controlled process” with drift vector which depends on the discrete component θ^ε (t). Thus, the state of the system, x ^ε(t), switches from one random path to another at random times as the mode θ^ε(t) changes. The discrete component θ^ε (t) is a “controlled Markov chain” with transition rate matrix depending on the continuous component. Both zero-sum and nonzero-sum games will be considered. In a zero-sum game, player I is trying to maximize certain expected payoff over his/her admissible strategies, where as player II is trying to minimize the same over his/her admissible strategies. This kind of game typically occurs in a pursuit-evation problem where an interceptor tries to destroy a specific target. Due to swift manuaring of the evador and the corresponding reaction by the interceptor, the trajectory keep switching rapidly at random times. We will show that the process (x ^ε(t), θ^ε(t)) converges to a process whose evolution is given by a “controlled diffusion process” with switching random paths. We will also establish the existence of randomized δ -optimal strategies for both players.     

Long time approximations for solutions of wave equations associated with the Steklov spectral homogenization problems

The interest in the use of quasimodes, or almost frequencies and almost eigenfunctions, to describe asymptotics for low‐frequency and high‐frequency vibrations in certain singularly perturbed spectral problems, which depend on a small parameter ε, has been recently highlighted in many papers. In this paper we deal with the low frequencies for a Steklov‐type eigenvalue homogenization problem: we consider harmonic functions in a bounded domain of ?², and strongly alternating boundary conditions of the Dirichlet and Steklov type on a part of the boundary. The spectral parameter appears in the boundary condition on small segments T^ε of size O(ε) periodically distributed along the boundary; ε also measures the periodicity of the structure. We consider associated second‐order evolution problems on spaces of traces that depend on ε, and we provide estimates for the time t in which standing waves, constructed from quasimodes, approach their solutions u^ε(t) as ε→0. Copyright © 2009 John Wiley & Sons, Ltd.     

Removable singularities of weak solutions to the navier-stokes equations

Consider the Navier-Stokes equations in Ω×(0,T), where Ω is a domain in R³. We show that there is an absolute constant ε₀ such that ever, y weak solution u with the property that Sup_tε(a,b)|u(t)|_L(D)≤ε0 is necessarily of class C^∞ in the space-time variables on any compact suhset of D × (a,b) , where D?? and 0 a<b<T. As an application. we prove that if the weak solution u behaves around (x_o, t_o) εΩ×(o,T) 1ike u(x, t) = o(|x - xo|^-1) as x→x ₀ uniforlnly in t in some neighbourliood of t_o, then (x_o,t_o) is actually a removable singularity of u.     

Rigorous results in selection of steady needle crystals

This paper concerns the existence of solutions to a steady needle crystal growth problem in a one-sided model. We rigorously prove that for small nonzero anisotropy γ, analytic symmetric needle crystal solutions exist in the limit of surface tension ε² if only if the stokes constant S for a relatively simple nonlinear differential equation is zero. This Stokes constant S depends on the parameter β=2^9/7γε^−8/7 and earlier numerical calculations by a number of investigators have shown this to be zero for a discrete set of values of β. It is also proven that for γ=0, there can be no symmetric needle crystal solution in the considered space.The methodology consists of two steps. First, the original problem is reduced to a weak half-strip problem for any γ in a compact set of [0,1) by relaxation of the symmetry condition. The weak problem is shown to have a unique solution in the function space considered for any γ∈[0,γ_m] for some γ_m>0. When a symmetry is invoked, the weak problem is shown equivalent to the original needle crystal problem. Next, by considering the behavior of the solution in neighborhood of an appropriate complex turning point for γ∈(0,γ_m], we extract an exponentially small term in ε as ε→0⁺ that generally violates the symmetric condition. We prove that the symmetry condition is satisfied for small ε when the parameter β is constrained appropriately.     

On the vibrations of a plate with a concentrated mass and very small thickness

We consider the vibrations of an elastic plate that contains a small region whose size depends on a small parameter ε. The density is of order O(ε^–m) in the small region, the concentrated mass, and it is of order O(1) outside; m is a positive parameter. The thickness plate h being fixed, we describe the asymptotic behaviour, as ε→O, of the eigenvalues and eigenfunctions of the corresponding spectral problem, depending on the value of m: Low‐ and high‐frequency vibrations are studied for m>2. We also consider the case where the thickness plate h depends on ε; then, different values of m are singled out. Copyright © 2003 John Wiley & Sons, Ltd.     

Fast Thermoviscous Convection

This paper studies the asymptotic structure of convection in an infinite Prandtl number fluid with strongly temperature-dependent viscosity, in the limit where the dimensionless activation energy 1/ε is large, and the Rayleigh number R, defined (essentially) with the basal viscosity and the prescribed temperature drop, is also large. We find that the Nusselt number N is given by N~CεR^1/5, where C depends on the aspect ratio a. The relative error in this result is O(R^?1/10ε^?1/4, ε^1/2, R^?2/5ε^?2, R^?2/20ε^?1/24), so that we cannot hope to find accurate confirmation of this result at moderate Rayleigh numbers, though it should serve as a useful indicator of the relative importance of R and ε. For the above result to be valid, we require R ? 1/ε⁵ ?1. More important, however, is the asymptotic structure of the flow: there is a cold (hence rigid) lid with sloping base, beneath which a rapid, essentially isoviscous, convection takes place. This convection is driven by plumes at the sides, which generate vorticity due to thermal buoyancy, as in the constant viscosity case (Roberts, 1979). However, the slope of the lid base is sufficient to cause a large shear stress to be generated in the thermal boundary layer which joins the lid to the isoviscous region underneath (though a large velocity is not generated); consequently, the layer does not “see” the shear stress exerted by the interior flow (at leading order), and therefore the thermal boundary layer structure is totally self-determined: it even has a similarity structure (as a consequence). This fact makes it easy to analyse the problem, since the boundary layer uncouples from the rest of the flow. In addition, we find an alternative scaling (in which the lid base is “almost” flat), but it seems that the resulting boundary layer equations have no solution, though this is certainly open to debate: the results quoted above are not for this case. When a free slip boundary condition is applied at the top surface, one finds that there exists a thin “skin” at the top of the lid which is a stress boundary layer. The shear stress changes rapidly to zero, and there exists a huge longitudinal stress (compressive/tensile) in this skin. For earthlike parameters, this stress far exceeds the fracture strength of silicate rocks.     

A mixed problem for the Laplace operator in a domain with moderately close holes

We investigate the behavior of the solution of a mixed problem in a domain with two moderately close holes. We introduce a positive parameter ε and we define a perforated domain Ω_ε obtained by making two small perforations in an open set. Both the size and the distance of the cavities tend to 0 as ε → 0. For ε small, we denote by u_ε the solution of a mixed problem for the Laplace equation in Ω_ε. We describe what happens to u_ε as ε → 0 in terms of real analytic maps and we compute an asymptotic expansion.     

Multidimensional Oscillations for Non-linear Hyperbolic Mixed Problem: The Justified Non-linear Geometric Optics Method

The mixed-Neumann problem for the non-linear wave equation □u−a(u)(∣∂_tu)∣²−∣∇u∣² = f_ε(z) is studied. The function f_ε(z) = ∑_k∈Kf_k(z,ε⁻¹ϕ_k(z),ε), ε∈[0,1], K is finite, f_k(z,θ_k,ε) are 2π-periodic with respect to θ_k. The existence of solution u_ε on a domain z = (t,x,y)∈[0,T]×ℝ₊×ℝ^d, d = 1 or 2, is proved when ε is sufficiently small; T does not depend on ε. By the non-linear geometric optics method the asymptotic (with respect to ε→0) solution ũ _ε is constructed. The estimation for the rest ε²r_ε = u_ε−ũ_ε is derived and the limit r_ε, ε→0, is studied.     

Homogenization of a non‐linear degenerate parabolic problem in a highly heterogeneous periodic medium

We consider the homogenization of a time‐dependent heat transfer problem in a highly heterogeneous periodic medium made of two connected components having comparable heat capacities and conductivities, separated by a third material with thickness of the same order ε as the basic periodicity cell but having a much lower conductivity such that the resulting interstitial heat flow is scaled by a factor λ tending to zero with a rate λ=λ(ε). The heat flux vectors a_j, j=1,2,3 are non‐linear, monotone functions of the temperature gradient. The heat capacities c_j(x) are positive, but may vanish at some subsets such that the problem can be degenerate (parabolic–elliptic). We show that the critical value of the problem is δ=lim_ε→0ε^p/λ and identify the homogenized problem depending on whether δ is zero, strictly positive finite or infinite. Copyright © 2005 John Wiley & Sons, Ltd.     

Convergence of a fourth-order singular perturbation of the n-dimensional radially symmetric Monge–Ampère equation

This paper concerns with the convergence analysis of a fourth-order singular perturbation of the Dirichlet Monge–Ampère problem in the n-dimensional radial symmetric case. A detailed study of the fourth- order problem is presented. In particular, various a priori estimates with explicit dependence on the perturbation parameter ε are derived, and a crucial convexity property is also proved for the solution of the fourth-order problem. Using these estimates and the convexity property, we prove that the solution of the perturbed problem converges uniformly and compactly to the unique convex viscosity solution of the Dirichlet Monge–Ampère problem. Rates of convergence in the H^k-norm for k = 0, 1, 2 are also established.     

Infeasible Mehrotra-Type Predictor-Corrector Interior-Point Algorithm for the Cartesian P *(κ)-LCP Over Symmetric Cones

We propose an infeasible Mehrotra-type predictor-corrector algorithm with a new center parameter updating scheme for Cartesian P _*(κ)-linear complementarity problem over symmetric cones. Based on the Nesterov-Todd direction, we show that the iteration-complexity bound of the proposed algorithm is 𝒪((1 + κ)³ r ²log ε^?1), where r is the rank of the associated Euclidean Jordan algebras and κ is the handicap of the problem and ε > 0 is the required precision. Some numerical results are reported as well.     

Construction of Solutions and L^1-error Estimates of Viscous Methods for Scalar Conservation Laws with Boundary

This paper is concerned with an initial boundary value problem for strictly convex conservation laws whose weak entropy solution is in the piecewise smooth solution class consisting of finitely many discontinuities. By the structure of the weak entropy solution of the corresponding initial value problem and the boundary entropy condition developed by Bardos-Leroux Nedelec, we give a construction method to the weak entropy solution of the initial boundary value problem. Compared with the initial value problem, the weak entropy solution of the initial boundary value problem includes the following new interaction type： an expansion wave collides with the boundary and the boundary reflects a new shock wave which is tangent to the boundary. According to the structure and some global estimates of the weak entropy solution, we derive the global L^1-error estimate for viscous methods to this initial boundary value problem by using the matching travelling wave solutions method. If the inviscid solution includes the interaction that an expansion wave collides with the boundary and the boundary reflects a new shock wave which is tangent to the boundary, or the inviscid solution includes some shock wave which is tangent to the boundary, then the error of the viscosity solution to the inviscid solution is bounded by O（ε^1/2） in L^1-norm; otherwise, as in the initial value problem, the L^1-error bound is O（ε｜ In ε｜）.     

The scaling window of the 2‐SAT transition

We consider the random 2‐satisfiability (2‐SAT) problem, in which each instance is a formula that is the conjunction of m clauses of the form x∨y, chosen uniformly at random from among all 2‐clauses on n Boolean variables and their negations. As m and n tend to infinity in the ratio m/n→α, the problem is known to have a phase transition at α_c=1, below which the probability that the formula is satisfiable tends to one and above which it tends to zero. We determine the finite‐size scaling about this transition, namely the scaling of the maximal window W(n, δ)=(α_?(n,δ), α₊(n,δ)) such that the probability of satisfiability is greater than 1?δ for α<α_? and is less than δ for α>α₊. We show that W(n,δ)=(1?Θ(n^?1/3), 1+Θ(n^?1/3)), where the constants implicit in Θ depend on δ. We also determine the rates at which the probability of satisfiability approaches one and zero at the boundaries of the window. Namely, for m=(1+ε)n, where ε may depend on n as long as |ε| is sufficiently small and |ε|n^1/3 is sufficiently large, we show that the probability of satisfiability decays like exp(?Θ(nε³)) above the window, and goes to one like 1?Θ(n^?1|ε|^?3 below the window. We prove these results by defining an order parameter for the transition and establishing its scaling behavior in n both inside and outside the window. Using this order parameter, we prove that the 2‐SAT phase transition is continuous with an order parameter critical exponent of 1. We also determine the values of two other critical exponents, showing that the exponents of 2‐SAT are identical to those of the random graph. © 2001 John Wiley & Sons, Inc. Random Struct. Alg., 18: 201–256 2001