Global Semigroup of Conservative Solutions of the Nonlinear Variational Wave Equation

We prove the existence of a global semigroup for conservative solutions of the nonlinear variational wave equation u _tt − c(u)(c(u)u _x ) _x  = 0. We allow for initial data u| _{t = 0} and u _t | _t=0 that contain measures. We assume that 0 < k^-1 \leqq c(u) \leqq k{0 < \kappa^{-1} \leqq c(u) \leqq \kappa}. Solutions of this equation may experience concentration of the energy density (u_t²+c(u)²u_x²)dx{(u_t^2+c(u)^2u_x^2){\rm d}x} into sets of measure zero. The solution is constructed by introducing new variables related to the characteristics, whereby singularities in the energy density become manageable. Furthermore, we prove that the energy may focus only on a set of times of zero measure or at points where c′(u) vanishes. A new numerical method for constructing conservative solutions is provided and illustrated with examples.     

Asymptotic Variational Wave Equations

We investigate the equation (u _t +(f(u)) _x ) _x =f ^{′ ′}(u) (u _x )²/2 where f(u) is a given smooth function. Typically f(u)=u ²/2 or u ³/3. This equation models unidirectional and weakly nonlinear waves for the variational wave equation u _tt − c(u) (c(u)u _x ) _x =0 which models some liquid crystals with a natural sinusoidal c. The equation itself is also the Euler–Lagrange equation of a variational problem. Two natural classes of solutions can be associated with this equation. A conservative solution will preserve its energy in time, while a dissipative weak solution loses energy at the time when singularities appear. Conservative solutions are globally defined, forward and backward in time, and preserve interesting geometric features, such as the Hamiltonian structure. On the other hand, dissipative solutions appear to be more natural from the physical point of view.We establish the well-posedness of the Cauchy problem within the class of conservative solutions, for initial data having finite energy and assuming that the flux function f has a Lipschitz continuous second-order derivative. In the case where f is convex, the Cauchy problem is well posed also within the class of dissipative solutions. However, when f is not convex, we show that the dissipative solutions do not depend continuously on the initial data.     

A New Approach to the Existence of Weak Solutions of the Steady Navier–Stokes System with Inhomogeneous Boundary Data in Domains with Noncompact Boundaries

We prove the existence of a weak solution to the steady Navier–Stokes problem in a three dimensional domain Ω, whose boundary ∂Ω consists of M unbounded components Γ¹, . . . , Γ ^M and N − M bounded components Γ ^M+1, . . . , Γ ^N . We use the inhomogeneous Dirichlet boundary condition on ∂Ω. The prescribed velocity profile α on ∂Ω is assumed to have an L ³-extension to Ω with the gradient in L ²(Ω)^3×3. We assume that the fluxes of α through the bounded components Γ ^M+1, . . . , Γ ^N of ∂Ω are “sufficiently small”, but we impose no restriction on the size of fluxes through the unbounded components Γ¹, . . . , Γ ^M .     

The Evolution of Stress Intensity Factors and the Propagation of Cracks in Elastic Media

When a crack Γ _s propagates in an elastic medium the stress intensity factors evolve with the tip x(s) of Γ _s . In this paper we derive formulae which describe the evolution of these stress intensity factors for a homogeneous isotropic elastic medium under plane strain conditions. Denoting by ψ=ψ(x,s) the stress potential (ψ is biharmonic and has zero traction along the crack Γ _s ) and by κ(s) the curvature of the crack at the tip x(s), we prove that the stress intensity factors A ₁(s), A ₂(s), as functions of s, satisfy:

Differential Independence of <Emphasis Type="Italic">Γ</Emphasis> and <Emphasis Type="Italic">ζ</Emphasis>

In a celebrated theorem H?lder proved that the Euler Γ-function is differential transcendental, i.e. Γ(z) is not a solution of any (non-trivial) algebraic ordinary differential equation with coefficients that are complex numbers; and we extend his methods to the Riemann ζ-function. Moreover, we conjecture that Γ and ζ are differential independent, i.e. Γ(z) is not a solution of any such algebraic differential equation—even allowing coefficients that are differential polynomials in ζ(z). However, we are able to demonstrate only the partial result that Γ(z) and ζ(sin 2πz) are differential independent.     

Static Solutions of the Vlasov-Einstein System

Static solutions of the SO(3)-symmetric Vlasov-Einstein system are studied via a variational approach. For the constitutive relation of the Emden-Fowler type φ(E,F)≡E ^{σ+ 1} F ^k we prove the existence of such solutions of sufficiently small mass-energy, provided 0<σ < k+3/2. These solutions are local minimizers of the energy-Casimir functional, subjected to a variational barrier. Accepted July 16, 2000?Published online January 22, 2001     

Mechanics of adhesive contact on a power-law graded elastic half-space

We consider adhesive contact between a rigid sphere of radius R and a graded elastic half-space with Young's modulus varying with depth according to a power law E=E₀(z/c₀)^k (0<k<1) while Poisson's ratio ν remaining a constant. Closed-form analytical solutions are established for the critical force, the critical radius of contact area and the critical interfacial stress at pull-off. We highlight that the pull-off force has a simple solution of P_cr=−(k+3)πRΔγ/2 where Δγ is the work of adhesion and make further discussions with respect to three interesting limits: the classical JKR solution when k=0, the Gibson solid when k→1 and ν=0.5, and the strength limit in which the interfacial stress reaches the theoretical strength of adhesion at pull-off.     

Stabilization using fractional-order PI and PID controllers

This paper presents a solution to the problem of stabilizing a given fractional dynamic system using fractional-order PIλ and PI^λD^μ controllers. It is based on plotting the global stability region in the (k _p, k _i)-plane for the PI^λ controller and in the (k _p , k _i , k _d)-space for the PIλDμ controller. Analytical expressions are derived for the purpose of describing the stability domain boundaries which are described by real root boundary, infinite root boundary and complex root boundary. Thus, the complete set of stabilizing parameters of the fractional-order controller is obtained. The algorithm has a simple and reliable result which is illustrated by several examples, and hence is practically useful in the analysis and design of fractional-order control systems.     

Boundary Layers and KPP Fronts in a Cellular Flow

We study an eigenvalue problem associated with a reaction-diffusion-advection equation of the KPP type in a cellular flow. We obtain upper and lower bounds on the eigenvalues in the regime of a large flow amplitude A ≪ 1. It follows that the minimal pulsating traveling front speed c _*(A) satisfies the upper and lower bounds C ₁ A ^1/4≦ c _*(A)≦ C ₂ A ^1/4. Physically, the speed enhancement is related to the boundary layer structure of the associated eigenfunction – accordingly, we establish an “averaging along the streamlines” principle for the unique positive eigenfunction.     

Concerning the W k,p -Inviscid Limit for 3-D Flows Under a Slip Boundary Condition

We consider the 3-D evolutionary Navier–Stokes equations with a Navier slip-type boundary condition, see (1.2), and study the problem of the strong convergence of the solutions, as the viscosity goes to zero, to the solution of the Euler equations under the zero-flux boundary condition. We prove here, in the flat boundary case, convergence in Sobolev spaces W ^k, p (Ω), for arbitrarily large k and p (for previous results see Xiao and Xin in Comm Pure Appl Math 60:1027–1055, 2007 and Beir?o da Veiga and Crispo in J Math Fluid Mech, 2009, doi:). However this problem is still open for non-flat, arbitrarily smooth, boundaries. The main obstacle consists in some boundary integrals, which vanish on flat portions of the boundary. However, if we drop the convective terms (Stokes problem), the inviscid, strong limit result holds, as shown below. The cause of this different behavior is quite subtle. As a by-product, we set up a very elementary approach to the regularity theory, in L ^p -spaces, for solutions to the Navier–Stokes equations under slip type boundary conditions.     

Another Thin-Film Limit of Micromagnetics

We consider the variational problem of micromagnetics for soft, relatively small thin films with no applied magnetic field. In terms of the film thickness t, the diameter l and the magnetic exchange length w, we study the asymptotic behavior in the small-aspect-ratio limit t/l→0, when either (a) w²/l²≫(t/l)| log (t/l)| or (b) w²/l²∼(t/l)| log (t/l)|. Our analysis builds on prior work by Gioia & James and Carbou. The limiting variational problem is much simpler than 3D micromagnetics; in particular it is two-dimensional and local, with no small parameters. The contribution of shape anisotropy reduces, in this limit, to a constant times the boundary integral of (m·n)².     

On resolution to Wu＇s conjecture on Cauchy function＇s exterior singularities

This is a series of studies on Wu’s conjecture and on its resolution to be presented herein. Both are devoted to expound all the comprehensive properties of Cauchy’s function f(z) (z = x + iy) and its integral J[f(z) ] ≡(2πi) -1 C f(t)(t z) -1dt taken along the unit circle as contour C,inside which(the open domain D+) f(z) is regular but has singularities distributed in open domain Doutside C. Resolution is given to the inverse problem that the singularities of f(z) can be determined in analytical form in terms of the values f(t) of f(z) numerically prescribed on C(|t| = 1) ,as so enunciated by Wu’s conjecture. The case of a single singularity is solved using complex algebra and analysis to acquire the solution structure for a standard reference. Multiple singularities are resolved by reducing them to a single one by elimination in principle,for which purpose a general asymptotic method is developed here for resolution to the conjecture by induction,and essential singularities are treated with employing the generalized Hilbert transforms. These new methods are applicable to relevant problems in mathematics,engineering and technology in analogy with resolving the inverse problem presented here.     

Variational Analysis of the Asymptotics of the <Emphasis Type="Italic">XY</Emphasis> Model

In this paper we consider the XY (N-dimensional possibly anisotropic) spin-type model and, by comparison with a Ginzburg–Landau-type functional, we perform a variational analysis in the limit when the number of particles diverges. In particular we show how the appearance of vortex-like singularities can be described by properly scaling the energy of the system through a Γ-convergence procedure. We also address the problem in the case of long-range interactions and solve it in two dimensions.     

A unified approach to closed-form solutions of moving heat-source problems

A closed-form model for the computation of temperature distribution in an infinitely extended isotropic body with a time-dependent moving-heat sources is discussed. The temperature solutions are presented for the sources of the forms: (i) ₀₁(t)=₀ exp(−λt), (ii) ₀₂(t) =₀(t/t ^*)exp(−λt), and ₀₃(t)=₀[1+a cos(ωt)], where λ and ω are real parameters and t ^* characterizes the limiting time. The reduced (or dimensionless) temperature solutions are presented in terms of the generalized representation of an incomplete gamma function Γ(α,x;b) and its decomposition C _Γ and S _Γ. The solutions are presented for moving, -point, -line, and -plane heat sources. It is also demonstrated that the present analysis covers the classical temperature solutions of a constant strength source under quasi-steady state situations. Received on 13 June 1997     

Crack Initiation in Brittle Materials

In this paper we study the crack initiation in a hyper-elastic body governed by a Griffith-type energy. We prove that, during a load process through a time-dependent boundary datum of the type t → t g(x) and in the absence of strong singularities (e.g., this is the case of homogeneous isotropic materials) the crack initiation is brutal, that is, a big crack appears after a positive time t _i > 0. Conversely, in the presence of a point x of strong singularity, a crack will depart from x at the initial time of loading and with zero velocity. We prove these facts for admissible cracks belonging to the large class of closed one-dimensional sets with a finite number of connected components. The main tool we employ to address the problem is a local minimality result for the functional where , k > 0 and f is a suitable Carathéodory function. We prove that if the uncracked configuration u of Ω relative to a boundary displacement ψ has at most uniformly weak singularities, then configurations (uΓ, Γ) with small enough are such that .     

一类Signorini问题的边界变分不等式

本文讨论了一类简化的Signorini问题。首先将原问题和一个边值问题建立联系，其次将原问题的解分解为不带不等边界条件的变分方程的解和一个变分不等式的解。然后利用边值问题的边界积分方程将变分不等式等价地化解为边界变分不等式。这样原求区域上的第一类椭圆变分不等式问题化解为求一个区域上的变分方程和一个边界变分不等式。最后说明了边界变分不等式解的存在唯一性。文末计算了柱面和半无限刚性基础的摩擦接触问题。结论表明文中方法具有较好的精度。     

On the Fundamental Solution of a Linearized Uehling–Uhlenbeck Equation

In this paper we describe the fundamental solution of the equation that is obtained linearizing the Uehling–Uhlenbeck equation around the steady state of Kolmogorov type f(k) = k ^−7/6. Detailed estimates on its asymptotics are obtained.     

Thin Walled Beams with Residual Stress

A one-dimensional variational problem for an anisotropic, partially inhomogeneous, residually stressed, rectangular thin-walled beam is derived, by Γ-convergence, from the three-dimensional theory of linear elasticity with residual stress.      

The laminar hole pressure for Newtonian fluids

For flows with wall turbulence the hole pressure, P _H , was shown empirically by Franklin and Wallace (J Fluid Mech, 42, 33–48, 1970) to depend solely on R ₊, the Reynolds number constructed from the friction velocity and the hole diameter b. Here this dependence is extended to the laminar regime by numerical simulation of a Newtonian fluid flowing in a plane channel (gap H) with a deep tap hole on one wall. Calculated hole pressures are in good agreement with experimental values, and for two hole sizes are well represented by: (P _H −P _HS )/τ _w = √(k ² + c ² R ₊²)−k, where the Stokes hole pressure P _HS /τ _w = s (b/H)³, k, c, s are fitted constants, and τ _w is the wall shear stress.     

Near-field tip vortex behind a swept wing model

The near-field flow structure of a tip vortex behind a sweptback and tapered NACA 0015 wing was investigated and compared with a rectangular wing at the same lift force and Re=1.81×10⁵. The tangential velocity decreased with the downstream distance while increased with the airfoil incidence. The core radius was about 3% of the root chord c _r, regardless of the downstream distance and α for α<8°. The core axial velocity was always wake-like. The core Γ_c and total Γ_o circulation of the tip vortex remained nearly constant up to x/c _r=3.5 and had a Γ_c/Γ_o ratio of 0.63. The total circulation of the tip vortex accounted for only about 40% of the bound root circulation Γ_b. For a rectangular wing, the axial flow exhibited islands of wake- and jet-like velocity distributions with Γ_c/Γ_o=0.75 and Γ_o/Γ_b=0.70. For the sweptback and tapered wing tested, the inner region of the tip vortex flow exhibited a self-similar behavior for x/c _r≥1.0. The lift force computed from the spanwise circulation distributions agreed well with the force-balance data. A large difference in the lift-induced drag was, however, observed between the wake integral method and the inviscid lifting-line theory.