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     

Asymptotic Stability of Riemann Solutions for a Class of Multidimensional Systems of Conservation Laws with Viscosity

We prove the asymptotic stability of two-state nonplanar Riemann solutions for a class of multidimensional hyperbolic systems of conservation laws when the initial data are perturbed and viscosity is added. The class considered here is those systems whose flux functions in different directions share a common complete system of Riemann invariants, the level surfaces of which are hyperplanes. In particular, we obtain the uniqueness of the self-similar L^∞ entropy solution of the two-state nonplanar Riemann problem. The asymptotic stability to which the main result refers is in the sense of the convergence as t→∞ in L_loc¹ of the space of directions ξ = x/t. That is, the solution u(t, x) of the perturbed problem satisfies u(t, tξ)→R(ξ) as t→∞, in L_loc¹(ℝⁿ), where R(ξ) is the self-similar entropy solution of the corresponding two-state nonplanar Riemann problem.     

Parameter Dependence of Stable Manifolds for Delay Equations with Polynomial Dichotomies

We establish the existence of Lipschitz stable invariant manifolds for semiflows generated by a delay equation x′ = L(t)x _t + f (t, x _t , λ), assuming that the linear equation x′ = L(t)x _t admits a polynomial dichotomy and that f is a sufficiently small Lipschitz perturbation. Moreover, we show that the stable invariant manifolds are Lipschitz in the parameter λ. We also consider the general case of nonuniform polynomial dichotomies.     

Uniform Convergence for Approximate Traveling Waves in Linear Reaction–Diffusion–Hyperbolic Systems

In this paper we study linear reaction–hyperbolic systems of the form , (i = 1, 2, ..., n) for x > 0, t > 0 coupled to a diffusion equation for p ₀ = p ₀(x, y, θ, t) with “near-equilibrium” initial and boundary data. This problem arises in a model of transport of neurofilaments in axons. The matrix (k _ij ) is assumed to have a unique null vector with positive components summed to 1 and the v _j are arbitrary velocities such that . We prove that as the solution converges to a traveling wave with velocity v and a spreading front, and that the convergence rate in the uniform norm is , for any small positive α.     

Fourier Representation of Intermittent Flow in a Porous Medium

The field measurements and numerical results for intermittent flow regime in a sandy soil show that the time distributions of the soil water flux q(z, t), and the soil water content θ(z, t)at various depths are periodic in nature, where t is time and z is the depth (i.e., at the surface z = 0 and at depths z = − 5, − 10, − 15 cm, etc). The period of q(z, t) and θ(z, t) variations are generally determined by the sum of the duration of pulse and the duration between the initiation of two consecutive pulses of water at the soil surface. Fourier series models have been given for q(z, t) and θ(z, t) variations. The predicted Fourier results for these variations have been compared with the experimentally verified numerical results—designated as observed values. The results show that the amplitudes of these variations were damped exponentially with depth, and the phase shift increased linearly with depth.     

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

Systems of the form

On thermodynamic potentials in linear thermoelasticity

The four thermodynamic potentials, the internal energy u=u(ε_ij,s), the Helmholtz free energy f=f(ε_ij,T), the Gibbs energy g=g(σ_ij,T) and the enthalpy h=h(σ_ij,s) are derived, independently of each other, by using the Duhamel–Neumann extension of Hooke's law and an assumed linear dependence of the specific heat on temperature. A systematic procedure is then presented to express all thermodynamic potentials in terms of four possible pairs of independent state variables. This procedure circumvents a tedious transition from one potential to another, based on the formal change of variables, and inversions of the stress–strain and entropy–temperature relations. The general results are applied to uniaxial loading paths under isothermal, adiabatic, constant stress, and constant strain conditions. An interplay of adiabatic and isothermal elastic constants in the expressions for exchanged heat along certain thermodynamic paths is indicated.     

The Decay Rate and Higher Approximation of Mild Solutions to the Navier�CStokes Equations

In this paper, we consider v(t) = u(t) − e ^tΔ u ₀, where u(t) is the mild solution of the Navier–Stokes equations with the initial data u₀ ? L²(\mathbb Rⁿ)?Lⁿ(\mathbb Rⁿ){u_0\in L^2({\mathbb R}^n)\cap L^n({\mathbb R}^n)} . We shall show that the L ² norm of D ^β v(t) decays like t^{-\frac |b|-1 2-\frac n4}{t^{-\frac {|\beta|-1} {2}-\frac n4}} for |β| ≥ 0. Moreover, we will find the asymptotic profile u ₁(t) such that the L ² norm of D ^β (v(t) − u ₁(t)) decays faster for 3 ≤ n ≤ 5 and |β| ≥ 0. Besides, higher-order asymptotics of v(t) are deduced under some assumptions.     

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.     

Estimating the viscoelastic moduli of complex fluids using the generalized Stokes–Einstein equation

We obtain the linear viscoelastic shear moduli of complex fluids from the time-dependent mean square displacement, <Δr ²(t)>, of thermally-driven colloidal spheres suspended in the fluid using a generalized Stokes–Einstein (GSE) equation. Different representations of the GSE equation can be used to obtain the viscoelastic spectrum, G˜(s), in the Laplace frequency domain, the complex shear modulus, G ^*(ω), in the Fourier frequency domain, and the stress relaxation modulus, G _r (t), in the time domain. Because trapezoid integration (s domain) or the Fast Fourier Transform (ω domain) of <Δr ²(t)> known only over a finite temporal interval can lead to errors which result in unphysical behavior of the moduli near the frequency extremes, we estimate the transforms algebraically by describing <Δr ²(t)> as a local power law. If the logarithmic slope of <Δr ²(t)> can be accurately determined, these estimates generally perform well at the frequency extremes. Received: 8 September 2000/Accepted: 9 March 2000     

On the free energy and variational theorem of elastic problem in thermal shock

In some investigations on variational principle for coupled thermoelastic problems, the free energy Φ(e_ij,θ) ,where the state variables are elastic strain e_ij and temperature increment θ, is expressed as Φ(e_ij,θ)=λ/2e_kke_ij=ue_k1e_k1-γe_kkθ-c/2 p θ²/T₀(0.1) This expression is employed only under the condition of |θ|≤T₀(absolute temperature of reference) But the value of temperature increment is great, even greater than T₀ in thermal shock. And the material properties (λ ,μ ,ν ,c , etc.) will not remain constant, they vary with θ. The expression of free energy for this condition.is derived in this paper. Equation (0.1) is its special case.Euler’s equations will be nonlinear while this expression of free energy has been introduced into variational theorem. In order to linearise, the time interval of thermal shock is divided into a number of time elements Δt_k, (Δt_k=t_k-t_k-1,k=1,2…,n), which are so small that the temperature increment θ_k within it is very small, too. Thus, the material properties may be defined by temperature field T_k-1=T(x₁,x₂,x₃,t_k-1) at instant t_k-1 , and the free energy Φ_k expressed by eg. (0.1) may be employed in element Δt_k.Hence the variational theorem will be expressed partly and approximately.     

Bifurcations of Heteroclinic Orbits

The search for traveling wave solutions of a semilinear diffusion partial differential equation can be reduced to the search for heteroclinic solutions of the ordinary differential equation ü − cu̇ + f(u) = 0, where c is a positive constant and f is a nonlinear function. A heteroclinic orbit is a solution u(t) such that u(t) → γ ₁ as t → −∞ and u(t) → γ ₂ as t → ∞ where γ ₁, γ ₂ are zeros of f. We study the existence of heteroclinic orbits under various assumptions on the nonlinear function f and their bifurcations as c is varied. Our arguments are geometric in nature and so we make only minimal smoothness assumptions. We only assume that f is continuous and that the equation has a unique solution to the initial value problem. Under these weaker smoothness conditions we reprove the classical result that for large c there is a unique positive heteroclinic orbit from 0 to 1 when f(0) = f(1) = 0 and f(u) > 0 for 0 < u < 1. When there are more zeros of f, there is the possibility of bifurcations of the heteroclinic orbit as c varies. We give a detailed analysis of the bifurcation of the heteroclinic orbits when f is zero at the five points −1 < −θ < 0 < θ < 1 and f is odd. The heteroclinic orbit that tends to 1 as t → ∞ starts at one of the three zeros, −θ, 0, θ as t → −∞. It hops back and forth among these three zeros an infinite number of times in a predictable sequence as c is varied.     

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≈ − ∞.     

On the Regularity Conditions of Suitable Weak Solutions of the 3D Navier–Stokes Equations

Let v and ω be the velocity and the vorticity of the a suitable weak solution of the 3D Navier–Stokes equations in a space-time domain containing z₀=(x₀, t₀)z_{0}=(x_{0}, t_{0}), and let Q_z0,r = B_x0,r ×(t₀ -r², t₀)Q_{z_{0},r}= B_{x_{0},r} \times (t_{0} -r^{2}, t_{0}) be a parabolic cylinder in the domain. We show that if either $\nu \times \frac{\omega}{|\omega|} \in L^{\gamma,\alpha}_{x,t}(Q_{z_{0},r})$\nu \times \frac{\omega}{|\omega|} \in L^{\gamma,\alpha}_{x,t}(Q_{z_{0},r}) with $\frac{3}{\gamma} + \frac{2}{\alpha} \leq 1, {\rm or} \omega \times \frac{\nu} {|\nu|} \in L^{\gamma,\alpha}_{x,t} (Q_{z_{0},r})$\frac{3}{\gamma} + \frac{2}{\alpha} \leq 1, {\rm or} \omega \times \frac{\nu} {|\nu|} \in L^{\gamma,\alpha}_{x,t} (Q_{z_{0},r}) with \frac3g + \frac2a ￡ 2\frac{3}{\gamma} + \frac{2}{\alpha} \leq 2, where L^{γ, α}_x,t denotes the Serrin type of class, then z₀ is a regular point for ν. This refines previous local regularity criteria for the suitable weak solutions.     

The radiated noise from isotropic turbulence

The noise radiated from isotropic turbulence at low Mach numbers and high Reynolds numbers, as derived by Proudman (1952), was the first application of Lighthill's Theory of Aerodynamic Noise to a complete flow field. The theory presented by Proudman involves the assumption of the neglect of retarded-time differences and so replaces the second-order retarded-time and space covariance of Lighthill's stress tensor, T _ij , and in particular its second time derivative, by the equivalent simultaneous covariance. This assumption is a valid approximation in the derivation of the ² T _ij /t ² covariance at low Mach numbers, but is not justified when that covariance is reduced to the sum of products of the time derivatives of equivalent second-order velocity covariances as required when Gaussian statistics are assumed. When these assumptions are removed the changes to the analysis are substantial, but the change in the numerical result for the total acoustic power is small.This paper is based on an alternative analysis which does not neglect retarded times. It makes use of the Lighthill relationship, whereby the fourth-order T _ij retarded-time covariance is evaluated from the square of a similar second-order covariance, which is assumed known. In this derivation no statistical assumptions are involved. This result, using distributions for the second-order space-time velocity squared covariance based on the Direct Numerical Simulation (DNS) results of both Sarkar and Hussaini (1993) and Dubois (1993), is compared with a re-evaluation of Proudman's original model. These results are then compared with the sound power derived from a phenomenological model based on simple approximations to the retarded-time/space covariance of T _xx . Finally, the recent numerical solutions of Sarkar and Hussaini (1993) for the acoustic power are compared with the results obtained from the analytic solutions.     

Ellipsoidal Domains: Piecewise Nonuniform and Impotent Eigenstrain Fields

In association with multi-inhomogeneity problems, a special class of eigenstrains is discovered to give rise to disturbance stresses of interesting nature. Some previously unnoticed properties of Eshelby’s tensors prove useful in this accomplishment. Consider the set of nested similar ellipsoidal domains {Ω₁, Ω₂,⋯,Ω _N+1}, which are embedded in an infinite isotropic medium. Suppose that

Bifurcation From Stability to Instability for a Free Boundary Problem Arising in a Tumor Model

We consider a time-dependent free boundary problem with radially symmetric initial data: σ_t − Δσ + σ = 0 if and σ(r,0)=σ₀(r) in {r < R(0)} where R(0) is given. This is a model for tumor growth, with nutrient concentration (or tumor cells density) σ(r,t) and proliferation rate then there exists a unique stationary solution (σ_S(r), R_S), where R_S depends only on the number . We prove that there exists a number μ_*, such that if μ < μ_* . . . then the stationary solution is stable with respect to non-radially symmetric perturbations, whereas if μ > μ_* then the stationary solution is unstable.     

Lagrangian Time-Scales in Homogeneous Non-Gaussian Turbulence

Lagrangian time-scales in homogeneous non-Gaussian turbulence were studied using a one-dimensional Lagrangian Stochastic Model. The existence of two time-scales τ _L and T _L , one typical of the inertial subrange and the other which is an integral property, is outlined. Variations of the ratio T _L /τ _L in the plane skewness-flatness (S, F) are shown and a connection with the statistical constraint F ≥S ² + 1 is evidenced. The Lagrangian autocorrelation function ρ(t) of particle velocity was computed for some values of (S, F). It is shown that for small times, say t < T _L , the influence of non-Gaussianity is negligible and ρ(t) presents the same behaviour as in the Gaussian case regardless of variations in (S, F).As the time increases, departures from Gaussianity are observed and autocorrelation turns out to be always larger than in the Gaussiancase. This is supported by some considerations in terms of information entropy, which is shown to decrease with increasing departures from Gaussianity. Spectral analysis of Lagrangian velocity shows that non-Gaussianity is relevant only to large scales of the stochastic process and that the expected inertial subrange decay ω⁻² is attained by spectra of all simulations, except for one case in which the model probability density function is bimodal, due to the vicinity to the statistical limit. This revised version was published online in July 2006 with corrections to the Cover Date.     

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

We establish new properties of C ¹[−1, +∞)-solutions of the linear functional differential equation ẋ(t) = ax(t) + bx(qt) + hx(t−1) + cẋ(qt) + rẋ(t−1) in the neighborhood of the singular point t = +∞. __________ Translated from Neliniini Kolyvannya, Vol. 9, No. 2, pp. 170–177, April–June, 2006.     

Oscillation criteria for a second-order quasilinear neutral functional dynamic equation on time scales

Our aim is to establish some sufficient conditions for the oscillation of the second-order quasilinear neutral functional dynamic equation