Breathers for a Relativistic Nonlinear Wave Equation

 A new class of one-dimensional relativistic nonlinear wave equations with a singular δ-type nonlinear term is considered. The sense of the equations is defined according to the least-action principle. The energy and momentum conservation is established. The main results are the existence of time-periodic finite-energy solutions, the existence of global solutions and soliton-type asymptotics for a class of finite-energy initial data. (Accepted May 28, 2002) Published online November 12, 2002 Communicated by G. FRIESECKE     

Optimization of a ribbon diode with magnetic insulation for increasing the current density in a high-current relativistic electron beam

The geometry of the ribbon diode of the U-2 accelerator is optimized to increase both the current density and the total current of the relativistic electron beam for its subsequent injection into the plasma of a multimirror GOL-3 trap. Beam simulation in the diode was performed using the POISSON-2 applied software modified on the basis of the results obtained using the theory of a planar diode in an inclined magnetic field. As a result of the optimization, the diode geometry and the magnetic field configuration were found that should provide a factor of 1.5–2 increase in the current density in experiments with a small angular divergence of electron velocities. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 3, pp. 25–35, May–June, 2009.     

On Asymptotic Stability of Kink for Relativistic Ginzburg–Landau Equations

We prove the asymptotic stability of kink for the nonlinear relativistic wave equations of the Ginzburg–Landau type in one space dimension: for any odd initial condition in a small neighborhood of the kink, the solution, asymptotically in time, is the sum of the kink and dispersive part described by the free Klein–Gordon equation. The remainder converges to zero in a global norm.     

The rheological behaviour of Ca(OH)2 suspensions

The rheological behaviour of Ca(OH)₂ suspensions is investigated, predominantly at a solid volume fraction of 0.25. The influence of standing without being subject to shear (“contact time”) is distinguished from that of being sheared (“shearing time”). The results are interpreted on the basis of the “elastic floc” model of energy dissipation during flow, with a view to the problem whether, in addition to an energy dissipation term related to the viscous drag experienced by particles moving within flocs, there should be an independent energy dissipation term related to fluid movement in the flocs when they change volume or shape. It appears that this additional energy dissipation term is not necessary, if the increase in viscous friction, experienced by two particles which are close together, is taken into account. Paper, presented at the First Conference of European Rheologists at Graz, April 14–16, 1982. A short version has been published in [18].     

Decay Structure for Symmetric Hyperbolic Systems with Non-Symmetric Relaxation and its Application

This paper is concerned with the decay structure for linear symmetric hyperbolic systems with relaxation. When the relaxation matrix is symmetric, the dissipative structure of the systems is completely characterized by the Kawashima–Shizuta stability condition formulated in Umeda et al. (Jpn J Appl Math 1:435–457, 1984) and Shizuta and Kawashima (Hokkaido Math J 14:249–275, 1985) and we obtain the asymptotic stability result together with the explicit time-decay rate under that stability condition. However, some physical models which satisfy the stability condition have non-symmetric relaxation term (for example, the Timoshenko system and the Euler–Maxwell system). Moreover, it had been already known that the dissipative structure of such systems is weaker than the standard type and is of the regularity-loss type (see Duan in J Hyperbolic Differ Equ 8:375–413, 2011; Ide et al. in Math Models Meth Appl Sci 18:647–667, 2008; Ide and Kawashima in Math Models Meth Appl Sci 18:1001–1025, 2008; Ueda et al. in SIAM J Math Anal 2012; Ueda and Kawashima in Methods Appl Anal 2012). Therefore our purpose in this paper is to formulate a new structural condition which includes the Kawashima–Shizuta condition, and to analyze the weak dissipative structure for general systems with non-symmetric relaxation.     

A Quasi-Realizable Cubic Low-Reynolds Eddy-Viscosity Turbulence Model with a New Dissipation Rate Equation

A non-linear relationship of the Reynolds stresses in function of the strain rate and vorticity tensors, with terms up to third order, is developed. Anisotropies in the normal stresses, influence from streamline curvature or rotation of the reference frame, and swirl effects are accounted for. The relationship is linked to ak–ε model with a modified transport equation for the dissipation rate. A new low-Reynolds source term is introduced and a model parameter is written in terms of dimensionless rate-of-strain and vorticity. The model is checked on different realizability constraints. It is shown that practically all constraints are fulfilled. The model is numerically tested on a fully developed channel and pipe flow, both stationary and rotating. The plane jet–round jet anomaly is addressed. Finally, the model is applied to the flow over a backward-facing step. Results are compared with a linear low-Reynolds k–ε model and the shear stress transport model. This revised version was published online in July 2006 with corrections to the Cover Date.     

Global Classical Solutions of the Relativistic Vlasov–Darwin System with Small Cauchy Data: The Generalized Variables Approach

We show that a smooth, small enough Cauchy datum launches a unique classical solution of the relativistic Vlasov–Darwin (RVD) system globally in time. A similar result is claimed in Seehafer (Commun Math Sci 6:749–769, 2008) following the work in Pallard (Int Mat Res Not 57191:1–31, 2006). Our proof does not require estimates derived from the conservation of the total energy, nor those previously given on the transversal component of the electric field. These estimates are crucial in the references cited above. Instead, we exploit the formulation of the RVD system in terms of the generalized space and momentum variables. By doing so, we produce a simple a priori estimate on the transversal component of the electric field. We widen the functional space required for the Cauchy datum to extend the solution globally in time, and we improve decay estimates given in Seehafer (2008) on the electromagnetic field and its space derivatives. Our method extends the constructive proof presented in Rein (Handbook of differential equations: evolutionary equations, vol 3. Elsevier, Amsterdam, 2007) to solve the Cauchy problem for the Vlasov–Poisson system with a small initial datum.     

Fractional conservation laws in optimal control theory

Using the recent formulation of Noether’s theorem for the problems of the calculus of variations with fractional derivatives, the Lagrange multiplier technique, and the fractional Euler–Lagrange equations, we prove a Noether-like theorem to the more general context of the fractional optimal control. As a corollary, it follows that in the fractional case the autonomous Hamiltonian does not define anymore a conservation law. Instead, it is proved that the fractional conservation law adds to the Hamiltonian a new term which depends on the fractional-order of differentiation, the generalized momentum and the fractional derivative of the state variable. Partially presented at FDA ’06—2nd IFAC Workshop on Fractional Differentiation and its Applications, 19–21 July 2006, Porto, Portugal.     

A new approach for the ellipsoidal statistical model

In this paper we aim to introduce a systematic way to derive relaxation terms for the Boltzmann equation based on the minimization problem for the entropy under moments constraints (Levermore in J. Stat. Phys. 83:1021–1065, 1996; Schneider in M2AN 38:541–561, 2004). In particular the moment constraints and corresponding coefficients are linked with the eigenfunctions and eigenvalues of the linearized collision operator through the Chapman–Enskog expansion. Then we deduce from this expansion a single relaxation term of BGK type. Here we stop the moments constraints at order two in the velocity v and recover the ellipsoidal statistical model (Holway in Rarefied Gas Dynamics, vol I, pp 193–215, 1966).      

Navier–Stokes Equations with Navier Boundary Conditions for an Oceanic Model

We consider the Navier–Stokes equations in a thin domain of which the top and bottom surfaces are not flat. The velocity fields are subject to the Navier conditions on those boundaries and the periodicity condition on the other sides of the domain. This toy model arises from studies of climate and oceanic flows. We show that the strong solutions exist for all time provided the initial data belong to a “large” set in the Sobolev space H ¹. Furthermore we show, for both the autonomous and the nonautonomous problems, the existence of a global attractor for the class of all strong solutions. This attractor is proved to be also the global attractor for the Leray–Hopf weak solutions of the Navier–Stokes equations. One issue that arises here is a nontrivial contribution due to the boundary terms. We show how the boundary conditions imposed on the velocity fields affect the estimates of the Stokes operator and the (nonlinear) inertial term in the Navier–Stokes equations. This results in a new estimate of the trilinear term, which in turn permits a short and simple proof of the existence of strong solutions for all time.     

Existence of positive solutions of systems of Volterra nonlinear difference equations

In this paper, a classification scheme for eventually positive solutions of a class of two-dimensional Volterra nonlinear difference equations is given in terms of asymptotic magnitudes. Some necessary as well as sufficient conditions for the existence of such solutions are provided without any monotonicity conditions on the nonlinear term. Published in Neliniini Kolyvannya, Vol. 9, No. 1, pp. 37–47, January–March, 2006.     

FD-method for a nonlinear eigenvalue problem with discontinuous eigenfunctions

An algorithm for the solution of a nonlinear eigenvalue problem with discontinuous eigenfunctions is developed. The numerical technique is based on a perturbation of the coefficients of a differential equation combined with the Adomian decomposition method for the nonlinear term of the equation. The proposed approach provides an exponential convergence rate dependent on the index of the trial eigenvalue and on the transmission coefficient. Numerical examples support the theory. Published in Neliniini Kolyvannya, Vol. 10, No. 1, pp. 126–143, January–March, 2007.     

Prediction of Flow and Heat Transfer through Stationary and Rotating Ribbed Ducts Using a Non-linear <Emphasis Type="Italic">k</Emphasis>−<Emphasis Type="Italic">ε</Emphasis> Model

The present paper deals with the prediction of three-dimensional fluid flow and heat transfer in rib-roughened ducts of square cross-section, which are either stationary, or rotate in orthogonal mode. The main objective is to assess how a recently developed variant of a cubic non-linear k−ε model (proposed by Craft et al. Flow Turbul Combust 63:59–80, 1999) can predict three-dimensional flow and heat transfer characteristics through stationary and rotating ribbed ducts. The present paper discusses turbulent air flow and heat transfer through two different configurations, namely: (I) a stationary square duct with “in-line” normal and (II) a square duct with normal ribs in a “staggered” arrangement under stationary and rotating conditions, with the axis of rotation normal to the flow direction and parallel to the ribs. In this paper the flow and thermal predictions of the linear k−ε model (EVM) are also included, as a set of baseline predictions. The mean flow predictions show that both linear and non-linear k−ε models can successfully reproduce most of the measured data for stream-wise and cross-stream velocity components. Moreover, the non-linear model is able to produce better results for the turbulent stresses. The heat transfer predictions show that both EVM and NLEVM2, the more recent variant of the non-linear k−ε, with the algebraic length-scale correction term, overestimate the measured Nusselt numbers for both geometries examined. While the EVM with the differential length-scale correction term underestimates heat transfer levels, the Nusselt number predictions with the NLEVM2 and the ‘NYP’ term are in close agreements with the measured data. Comparisons with our earlier work, Iacovides and Raisee (Int J Heat Fluid Flow, 20:320–328, 1999), show that the NLEVM2 thermal predictions are of similar quality to those of a second-moment closure.     

Global Attractor for a Wave Equation with Nonlinear Localized Boundary Damping and a Source Term of Critical Exponent

The paper addresses long-term behavior of solutions to a damped wave equation with a critical source term. The dissipative frictional feedback is restricted to a subset of the boundary of the domain. This paper derives inverse observability estimates which extend the results of Chueshov et al. (Disc Contin Dyn Syst 20:459–509, 2008) to systems with boundary dissipation. In particular, we show that a hyperbolic flow under a critical source and geometrically constrained boundary damping approaches a smooth finite-dimensional global attractor. A similar result for subcritical sources was given in Chueshov et al. (Commun Part Diff Eq 29:1847–1876, 2004). However, the criticality of the source term in conjunction with geometrically restricted dissipation constitutes the major new difficulty of the problem. To overcome this issue we develop a special version of Carleman’s estimates and apply them in the context of abstract results on dissipative dynamical systems. In contrast with the localized interior damping (Chueshov et al. Disc Contin Dyn Syst 20:459–509, 2008), the analysis of a boundary feedback requires a more careful treatment of the trace terms and special tangential estimates based on microlocal analysis.      

Solutions of the Dirac–Fock Equations and the Energy of the Electron-Positron Field

We consider atoms with closed shells, i.e. the electron number N is 2, 8, 10,..., and weak electron-electron interaction. Then there exists a unique solution γ of the Dirac–Fock equations with the additional property that γ is the orthogonal projector onto the first N positive eigenvalues of the Dirac–Fock operator . Moreover, γ minimizes the energy of the relativistic electron-positron field in Hartree–Fock approximation, if the splitting of into electron and positron subspace is chosen self-consistently, i.e. the projection onto the electron-subspace is given by the positive spectral projection of. For fixed electron-nucleus coupling constant g:=α Z we give quantitative estimates on the maximal value of the fine structure constant α for which the existence can be guaranteed.     

Polynomial Filtration Laws for Low Reynolds Number Flows Through Porous Media

In this study, we use the method of homogenization to develop a filtration law in porous media that includes the effects of inertia at finite Reynolds numbers. The result is much different than the empirically observed quadratic Forchheimer equation. First, the correction to Darcy’s law is initially cubic (not quadratic) for isotropic media. This is consistent with several other authors (Mei and Auriault, J Fluid Mech 222:647–663, 1991; Wodié and Levy, CR Acad Sci Paris t.312:157–161, 1991; Couland et al. J Fluid Mech 190:393–407, 1988; Rojas and Koplik, Phys Rev 58:4776–4782, 1988) who have solved the Navier–Stokes equations analytically and numerically. Second, the resulting filtration model is an infinite series polynomial in velocity, instead of a single corrective term to Darcy’s law. Although the model is only valid up to the local Reynolds number, at the most, of order 1, the findings are important from a fundamental perspective because it shows that the often-used quadratic Forchheimer equation is not a universal law for laminar flow, but rather an empirical one that is useful in a limited range of velocities. Moreover, as stated by Mei and Auriault (J Fluid Mech 222:647–663, 1991) and Barree and Conway (SPE Annual technical conference and exhibition, 2004), even if the quadratic model were valid at moderate Reynolds numbers in the laminar flow regime, then the permeability extrapolated on a Forchheimer plot would not be the intrinsic Darcy permeability. A major contribution of this study is that the coefficients of the polynomial law can be derived a priori, by solving sequential Stokes problems. In each case, the solution to the Stokes problem is used to calculate a coefficient in the polynomial, and the velocity field is an input of the forcing function, F, to subsequent problems. While numerical solutions must be utilized to compute each coefficient in the polynomial, these problems are much simpler and robust than solving the full Navier–Stokes equations.     

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 Ω.      

Effects of Chemical Reaction and Double Dispersion on Non-Darcy Free Convection Heat and Mass Transfer

In this article, the effects of chemical reaction and double dispersion on non-Darcy free convection heat and mass transfer from semi-infinite, impermeable vertical wall in a fluid saturated porous medium are investigated. The Forchheimer extension (non-Darcy term) is considered in the flow equations, while the chemical reaction power–law term is considered in the concentration equation. The first order chemical reaction (n = 1) was used as an example of calculations. The Darcy and non-Darcy flow, temperature and concentration fields in this study are observed to be governed by complex interactions among dispersion and natural convection mechanisms. The governing set of partial differential equations were non-dimensionalized and reduced to a set of ordinary differential equations for which Runge–Kutta-based numerical technique were implemented. Numerical results for the detail of the velocity, temperature, and concentration profiles as well as heat transfer rates (Nusselt number) and mass transfer rates (Sherwood number) are presented in graphs.     

A Conceptual Study of Cavity Aeroacoustics Control Using Porous Media Inserts

In this study, an integrated flow simulation and aeroacoustics prediction methodology is applied to testing a sound control technique using porous inserts in an open cavity. Large eddy simulation (LES) combined with a three-dimensional Ffowcs Williams–Hawkings (FW–H) acoustic analogy is employed to predict the flow field, the acoustic sources and the sound radiation. The Darcy pressure – velocity law is applied to conceptually mimic the effect of porous media placed on the cavity floor and/or rear wall. Consequently, flow in the cavity could locally move in or out through these porous walls, depending on the local pressure differences. LES with “standard” subgrid-scale models for compressible flow is carried out to simulate the flow field covering the sound source and near fields, and the fully three-dimensional FW–H acoustic analogy is used to predict the sound field. The numerical results show that applying the conceptual porous media on cavity floor and/or rear wall could decrease the pressure fluctuations in the cavity and the sound pressure level in the far field. The amplitudes of the dominant oscillations (Rossiter modes) are suppressed and their frequencies are slightly modified. The dominant sound source is the transverse dipole term, which is significantly reduced due to the porous walls. As a result, the sound pressure in the far field is also suppressed. The preliminary study reveals that using porous-inserts is a promising technology for flow and sound radiation control.     

An Optimal Condition of Compactness for Elasticity Problems Involving one Directional Reinforcement

This paper deals with the homogenization of a homogeneous elastic medium reinforced by very stiff strips in dimension two. We give a general condition linked to the distribution and the stiffness of the strips, under which the nature of the elasticity problem is preserved in the homogenization process. This condition is sharper than the one used in Briane and Camar-Eddine (J. Math. Pures Appl. 88:483–505, 2007) and is shown to be optimal in the case where the strips are periodically arranged. Indeed, a fourth-order derivative term appears in the limit equation as soon as the condition is no more satisfied. In the periodic case the influence of oscillations in the medium surrounding the strips is also considered. The homogenization method is based both on a two-scale convergence for the strips and the use of suitable oscillating test functions. This allows us to obtain a distributional convergence of two of the three entries of the stress tensor contrary to the Γ-convergence approach of Briane and Camar-Eddine (J. Math. Pures Appl. 88:483–505, 2007).