On Aerodynamic Forces for Viscous Compressible Flow

The paper is aimed at reviewing and adding some new results to our recent work on a force theory for viscous compressible flows around a finite body. It has been proposed to analyze aerodynamic forces directly in terms of fluid elements of nonzero vorticity and density gradient. Let ρ denote the density, u the velocity, and ω the vorticity. It is demonstrated that for largely separated flows about bluff bodies, there are two major source elements: R _e(x) =−?u ²∇ρ·∇ϕ and V _e(x) =−u×ω·∇ϕ, where ϕ is an acyclic potential, generated by the solid body moving with unit velocity in the negative direction of the force considered. In particular, under mild conditions, the (unique) choice of ϕ enforces that the elements R _e(x) and V _e(x) decay rapidly away from the body. Four kinds of finite body are considered: a circular cylinder, a sphere, a hemi sphere-cylinder, and a delta wing of elliptic section—all in transonic-to-supersonic regimes. From an extensive numerical study carried out for these bodies, it is found that these two elements contribute to 95% or more of the total drag or lift for all the cases under consideration. Moreover, R _e(x) due to density gradient becomes progressively important relative to V _e(x) due to vorticity as the Mach number increases. The present method of force analysis enables effective analysis and assessment of relative importance of aerodynamics forces, contributed from individual flow structures. The analysis could therefore be very much useful in view of the rapid growth in numerical fluid dynamics; detailed (either local or global) flow information is often available. The paper is dedicated to Sir James Lighthill in honor of his great contributions to aeronautics on the occasion of the publication of his collected works. Received 3 January 1997 and accepted 11 April 1997     

Entire Solutions with Merging Fronts to Reaction–Diffusion Equations

We deal with a reaction–diffusion equation u _t = u _xx + f(u) which has two stable constant equilibria, u = 0, 1 and a monotone increasing traveling front solution u = φ(x + ct) (c > 0) connecting those equilibria. Suppose that u = a (0 < a < 1) is an unstable equilibrium and that the equation allows monotone increasing traveling front solutions u = ψ₁(x + c ₁ t) (c ₁ < 0) and ψ₂(x + c ₂ t) (c ₂ > 0) connecting u = 0 with u = a and u = a with u = 1, respectively. We call by an entire solution a classical solution which is defined for all . We prove that there exists an entire solution such that for t≈ − ∞ it behaves as two fronts ψ₁(x + c ₁ t) and ψ₂(x + c ₂ t) on the left and right x-axes, respectively, while it converges to φ(x + ct) as t→∞. In addition, if c > − c ₁, we show the existence of an entire solution which behaves as ψ₁( − x + c ₁ t) in and φ(x + ct) in for t≈ − ∞.     

Global Dynamics of Blow-up Profiles in One-dimensional Reaction Diffusion Equations

We consider reaction diffusion equations of the prototype form u _t = u _xx + λ u + |u| ^p-1 u on the interval 0 < x < π, with p > 1 and λ > m ². We study the global blow-up dynamics in the m-dimensional fast unstable manifold of the trivial equilibrium u ≡ 0. In particular, sign-changing solutions are included. Specifically, we find initial conditions such that the blow-up profile u(t, x) at blow-up time t = T possesses m + 1 intervals of strict monotonicity with prescribed extremal values u ₁, . . . ,u _m . Since u _k = ± ∞ at blow-up time t = T, for some k, this exhausts the dimensional possibilities of trajectories in the m-dimensional fast unstable manifold. Alternatively, we can prescribe the locations x = x ₁, . . . ,x _m of the extrema, at blow-up time, up to a one-dimensional constraint. The proofs are based on an elementary Brouwer degree argument for maps which encode the shapes of solution profiles via their extremal values and extremal locations, respectively. Even in the linear case, such an “interpolation of shape” was not known to us. Our blow-up result generalizes earlier work by Chen and Matano (1989), J. Diff. Eq. 78, 160–190, and Merle (1992), Commun. Pure Appl. Math. 45(3), 263–300 on multi-point blow-up for positive solutions, which were not constrained to possess global extensions for all negative times. Our results are based on continuity of the blow-up time, as proved by Merle (1992), Commun. Pure Appl. Math. 45(3), 263–300, and Quittner (2003), Houston J. Math. 29(3), 757–799, and on a refined variant of Merle’s continuity of the blow-up profile, as addressed in the companion paper Matano and Fiedler (2007) (in preparation). Dedicated to Palo Brunovsky on the occasion of his birthday.     

Darboux problem for a differential equation with fractional derivative

We find conditions for the unique solvability of the problem u _xy (x, y) = f(x, y, u(x, y), (D ₀ ^r u)(x, y)), u(x, 0) = u(0, y) = 0, x ∈ [0, a], y ∈ [0, b], where (D ₀ ^r u)(x, y) is the mixed Riemann-Liouville derivative of order r = (r ₁, r ₂), 0 < r ₁, r ₂ < 1, in the class of functions that have the continuous derivatives u _xy (x, y) and (D ₀ ^r u)(x, y). We propose a numerical method for solving this problem and prove the convergence of the method. __________ Translated from Neliniini Kolyvannya, Vol. 8, No. 4, pp. 456–467, October–December, 2005.     

A Note on the Existence of Solutions to the Oseen System in Lipschitz Domains

We study the boundary-value problem associated with the Oseen system in the exterior of m Lipschitz domains of an euclidean point space We show, among other things, that there are two positive constants and α depending on the Lipschitz character of Ω such that: (i) if the boundary datum a belongs to L^q(∂Ω), with q ∈ [2,+∞), then there exists a solution (u, p), with and u ∈ L^∞(Ω) if a ∈ L^∞(∂Ω), expressed by a simple layer potential plus a linear combination of regular explicit functions; as a consequence, u tends nontangentially to a almost everywhere on ∂Ω; (ii) if a ∈ W^1-1/q,q(∂Ω), with then ∇u, p ∈ L^q(Ω) and if a ∈ C^0,μ(∂Ω), with μ ∈ [0, α), then also, natural estimates holds.     

Gearhart–Prüss Theorem and Linear Stability for Riemann Solutions of Conservation Laws

We study the spectral and linear stability of Riemann solutions with multiple Lax shocks for systems of conservation laws. Using a self-similar change of variables, Riemann solutions become stationary solutions for the system u _t + (Df(u) − x I)u _x = 0. In the space of O((1 + |x|)^−η) functions, we show that if , then λ is either an eigenvalue or a resolvent point. Eigenvalues of the linearized system are zeros of the determinant of a transcendental matrix. On some vertical lines in the complex plane, called resonance lines, the determinant can be arbitrarily small but nonzero. A C ⁰ semigroup is constructed. Using the Gearhart–Prüss Theorem, we show that the solutions are O(e ^{γ t} ) if γ is greater than the real parts of the eigenvalues and the coordinates of resonance lines. We study examples where Riemann solutions have two or three Lax-shocks. Dedicated to Professor Pavol Brunovsky on his 70th birthday.     

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.     

Algebraic-Delay Differential Systems, State-Dependent Delay, and Temporal Order of Reactions

Systems of the form

Orientational anisotropy for Rouse eigenmodes during creep and recovery process

The Rouse model is a well established model for nonentangled polymer chains and its dynamic behavior under step strain has been fully analyzed in the literature. However, to the knowledge of the authors, no analysis has been made for the orientational anisotropy for the Rouse eigenmodes during the creep and creep recovery processes. For completeness of the analysis of the Rouse model, this anisotropy is calculated from the Rouse equation of motion. The calculation is simple and straightforward, but the result is intriguing in a sense that respective Rouse eigenmodes do not exhibit the single Voigt-type retardation. Instead, each Rouse eigenmode has a distribution in the retardation time. This behavior, reflecting the interplay among the Rouse eigenmodes of different orders under the constant stress condition, is quite different from the behavior under rate-controlled flow (where each eigenmode exhibits retardation/relaxation associated with a single characteristic time).List of abbreviations and symbols a Average segment size at equilibrium - A_p(t) Normalized orientational anisotropy for the p-th Rouse eigenmode defined by Eq. (14) - p-th Fourier component of the Brownian force (=x, y) - F_B(n,t) Brownian force acting on n-th segment at time t - G(t) Relaxation modulus - J(t) Creep compliance - J_R(t) Recoverable creep compliance - k_B Boltzmann constant - N Segment number per Rouse chain - Q_j(t) Orientational anisotropy of chain sections defined by Eq. (21) - r(n,t) Position of n-th segment of the chain at time t - S(n,t) Shear orientation function (S(n,t)=a^–2<u_x(n,t)u_y(n,t)>) - T Absolute temperature - u(n,t) Tangential vector of n-th segment at time t (u = r/n) - V(r(n,t)) Flow velocity of the frictional medium at the position r(n,t) - X_p(t), Y_p(t), and Z_p(t) x-, y-, and z-components of the amplitudes of p-th Rouse eigenmode at time t - Strain rate being uniform throughout the system - Segmental friction coefficient - ₀ Zero-shear viscosity - _p Numerical coefficients determined from Eq. (25) - Gaussian spring constant ( = 3k_BT/a²) - Number of Rouse chains per unit volume - (t) Shear stress of the system at time t - _steady Shear stress in the steadily flowing state - _R Longest viscoelastic relaxation time of the Rouse chain     

A Semilinear Schrödinger Equation in the Presence of a Magnetic Field

We consider the semilinear stationary Schrödinger equation in a magnetic field: (–i+A)² u+V(x)u=g(x,|u|)u in ^N , where V is the scalar (or electric) potential and A is the vector (or magnetic) potential. We study the existence of nontrivial solutions both in the critical and in the subcritical case (respectively g(x,|u|)=|u|^{2 * –2} and |g(x,|u|)|c(1+|u| ^{p –2}), where 2<p<2^*). The results are obtained by variational methods. For g critical we use constrained minimization and for subcritical g we employ a minimax-type argument. In the latter case we also study the existence of infinitely many geometrically distinct solutions.     

Sufficient Conditions for the Existence of Bounded Solutions of Systems of Nonlinear Difference Equations

We obtain sufficient conditions for systems of nonlinear difference equations x(n + 1) = A(x(n))x(n) + f(n), n ∈ ℤ, where A(x) is a matrix function continuous on ℝ ^m , to have solutions in the space of bilateral number sequences. __________ Translated from Neliniini Kolyvannya, Vol. 8, No. 2, pp. 165–173, April–June, 2005.     

Vorticity and Regularity for Viscous Incompressible Flows under the Dirichlet Boundary Condition. Results and Related Open Problems

In reference [7] it is proved that the solution of the evolution Navier–Stokes equations in the whole of R ³ must be smooth if the direction of the vorticity is Lipschitz continuous with respect to the space variables. In reference [5] the authors improve the above result by showing that Lipschitz continuity may be replaced by 1/2-H?lder continuity. A central point in the proofs is to estimate the integral of the term (ω · ∇)u · ω, where u is the velocity and ω = ∇ × u is the vorticity. In reference [4] we extend the main estimates on the above integral term to solutions under the slip boundary condition in the half-space R ₊³. This allows an immediate extension to this problem of the 1/2-H?lder sufficient condition. The aim of these notes is to show that under the non-slip boundary condition the above integral term may be estimated as well in a similar, even simpler, way. Nevertheless, without further hypotheses, we are not able now to extend to the non slip (or adherence) boundary condition the 1/2-H?lder sufficient condition. This is not due to the “nonlinear" term (ω · ∇)u · ω but to a boundary integral which is due to the combination of viscosity and adherence to the boundary. On the other hand, by appealing to the properties of Green functions, we are able to consider here a regular, arbitrary open set Ω.      

Spatial and Dynamical Chaos Generated by Reaction–Diffusion Systems in Unbounded Domains

We consider in this article a nonlinear reaction–diffusion system with a transport term (L,∇ _x )u, where L is a given vector field, in an unbounded domain Ω. We prove that, under natural assumptions, this system possesses a locally compact attractor in the corresponding phase space. Since the dimension of this attractor is usually infinite, we study its Kolmogorov’s ɛ-entropy and obtain upper and lower bounds of this entropy. Moreover, we give a more detailed study of the spatio-temporal chaos generated by the spatially homogeneous RDS in . In order to describe this chaos, we introduce an extended (n + 1)-parametrical semigroup, generated on the attractor by 1-parametrical temporal dynamics and by n-parametrical group of spatial shifts ( = spatial dynamics). We prove that this extended semigroup has finite topological entropy, in contrast to the case of purely temporal or purely spatial dynamics, where the topological entropy is infinite. We also modify the concept of topological entropy in such a way that the modified one is finite and strictly positive, in particular for purely temporal and for purely spatial dynamics on the attractor. In order to clarify the nature of the spatial and temporal chaos on the attractor, we use (following Zelik, 2003, Comm. Pure. Appl. Math. 56(5), 584–637) another model dynamical system, which is an adaptation of Bernoulli shifts to the case of infinite entropy and construct homeomorphic embeddings of it into the spatial and temporal dynamics on . As a corollary of the obtained embeddings, we finally prove that every finite dimensional dynamics can be realized (up to a homeomorphism) by restricting the temporal dynamics to the appropriate invariant subset of .     

On asymptotic properties of solutions of functional differential equations with linearly transformed argument

We establish new properties of solutions of the functional differential equation x′(t) = ax(t) + bx(t − r) + cx′(t − r) + px(qt) + hx′(qt) + f ₁(x(t), x(t − r), x′(t − r), x(qt), x′(qt)) in the neighborhood of the singular point t = +∞. __________ Translated from Neliniini Kolyvannya, Vol. 10, No. 1, pp. 144–160, January–March, 2007.     

Anisotropic Elastic Materials that Uncouple All Three Displacement Components, and Existence of One-Displacement Green's Function

It was shown in an earlier paper that, under a two-dimensional deformation, there are anisotropic elastic materials for which the antiplane displacement u ₃ and the inplane displacements u ₁, u ₂ are uncoupled but the antiplane stresses σ₃₁, σ₃₂ and the inplane stresses σ₁₁, σ₁₂, σ₂₂ remain coupled. The conditions for this to be possible were derived, but they have a complicated expression. In this paper new and simpler conditions are obtained, and a general anisotropic elastic material that satisfies the conditions is presented. For this material, and for certain monoclinic materials with the symmetry plane at x ₃ = 0, we show that the unnormalized Stroh eigenvectors a _k for k = 1, 2, 3 are all real. The matrix A =[a ₁, a ₂, a ₃] is a unit matrix when the material has a symmetry plane at x ₂ = 0. Thus any one of the u ₁, u ₂, u ₃ can be the only nonzero displacement, and the solution is a one-displacement field. Application to the Green's function due to a line of concentrated force f and a line dislocation with Burgers vector v in the infinite space, the half-space with a rigid boundary, and the infinite space with an elliptic rigid inclusion shows that one can indeed have a one-displacement field u ₁, u ₂ or u ₃. One can also have a two-displacement field polarized on a plane other than the (x ₁, x ₂)-plane. The material that uncouples u ₁, u ₂, u ₃ is not as restrictive as one might have thought. It can be triclinic, monoclinic, orthotropic, tetragonal, transversely isotropic, or cubic. However, it cannot be isotropic. This revised version was published online in August 2006 with corrections to the Cover Date.     

A Locally Quadratic Glimm Functional and Sharp Convergence Rate of the Glimm Scheme for Nonlinear Hyperbolic Systems

Consider the Cauchy problem for a strictly hyperbolic, N × N quasilinear system in one space dimension

Existence of Solutions with the Prescribed Flux of the Navier–Stokes System in an Infinite Cylinder

The existence and uniqueness of a solution to the nonstationary Navier–Stokes system having a prescribed flux in an infinite cylinder is proved. We assume that the initial data and the external forces do not depend on x₃ and find the solution (u, p) having the following form

A Sharp-Interface Limit for a Two-Well Problem in Geometrically Linear Elasticity

In the theory of solid-solid phase transitions the deformation of an elastic body is determined via a functional containing a nonconvex energy density and a singular perturbation. We study Frame indifference, within a linearized setting, requires that W depends only on the symmetric part of ∇u. The potential W is nonnegative and vanishes on two wells, i.e., for d = 2, on two lines in the space of matrices. We determine, for d = 2, the Gamma limit I₀ = Γ− lim _ɛ→0 I_ɛ. The limit I₀[u] is finite only for deformations u that fulfill W(∇u)=0 almost everywhere and have sharp interfaces where ∇u has jumps. For these u, I₀[u] equals the line integral over the interfaces of a surface energy density.     

Sharp Estimates of Bounded Solutions to Some Second-order Forced Dissipative Equations

Résumé A l’aide d’inégalités différentielles, on établit une estimation proche de l’optimalité pour la norme dans de l’unique solution bornée de u′′ + cu′ + Au = f(t) lorsque A = A ^* ≥ λ I est un opérateur borné ou non sur un espace de Hilbert réel H, V = D(A ^1/2) et λ, c sont des constantes positives, tandis que . By using differential inequalities, a close-to-optimal bound of the unique bounded solution of u′′ + cu′ + Au = f(t) is obtained whenever A = A ^* ≥ λ I is a bounded or unbounded linear operator on a real Hilbert space H, V = D(A ^1/2) and λ, c are positive constants, while .

---

     

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: