The initial-value problem for elastodynamics of incompressible bodies

We prove short-time well-posedness of the Cauchy problem for incompressible strongly elliptic hyperelastic materials. Our method consists in:

1. Reformulating the classical equations in order to solve for the pressure gradient (The pressure is the Lagrange multiplier corresponding to the constraint of incompressibility.) This formulation uses both spatial and material variables.
2. Solving the reformulated equations by using techniques which are common for symmetric hyperbolic systems. These are:

1. Using energy estimates to bound the growth of various Sobolev norms of solutions.
2. Finding the solution as the limit of a sequence of solutions of linearized problems.

Our equations differ from hyperbolic systems, however, in that the pressure gradient is a spatially non-local function of the position and velocity variables.     

Dynamics as a mechanism preventing the formation of finer and finer microstructure

We study the dynamics of pattern formation in the one-dimensional partial differential equation $$u_u - (W'(u_x ))_x - u_{xxt} + u = 0{\text{ (}}u = u(x,t),{\text{ }}x \in (0,1),{\text{ }}t > 0)$$ proposed recently by Ball, Holmes, James, Pego & Swart [BHJPS] as a mathematical “cartoon” for the dynamic formation of microstructures observed in various crystalline solids. Here W is a double-well potential like 1/4((u _x )² ?1)². What makes this equation interesting and unusual is that it possesses as a Lyapunov function a free energy (consisting of kinetic energy plus a nonconvex “elastic” energy, but no interfacial energy contribution) which does not attain a minimum but favours the formation of finer and finer phase mixtures: $$E[u,u_t ] = \int\limits_0^1 {(\frac{{u_t^2 }}{2} + W(u_x ) + \frac{{u^2 }}{2})dx.}$$ Our analysis of the dynamics confirms the following surprising and striking difference between statics and dynamics, conjectured in [BHJPS] on the basis of numerical simulations of Swart & Holmes [SH]:

?While minimizing the above energy predicts infinitely fine patterns (mathematically: weak but not strong convergence of all minimizing sequences (u _nv_n) of E[u,v] in the Sobolev space W ^{1 p}(0, 1)×L²(0,1)), solutions to the evolution equation of ball et al. typically develop patterns of small but finite length scale (mathematically: strong convergence in W ₁ ^p(0,1)×L²(0,1) of all solutions (u(t),u_t(t)) with low initial energy as time t → ∞).

Moreover, in order to understand the finer details of why the dynamics fails to mimic the behaviour of minimizing sequences and how solutions select their limiting pattern, we present a detailed analysis of the evolution of a restricted class of initial data — those where the strain field u _x has a transition layer structure; our analysis includes proofs that

?at low energy, the number of phases is in fact exactly preserved, that is, there is no nucleation or coarsening

?transition layers lock in and steepen exponentially fast, converging to discontinuous stationary sharp interfaces as time t → ∞

?the limiting patterns — while not minimizing energy globally — are ‘relative minimizers’ in the weak sense of the calculus of variations, that is, minimizers among all patterns which share the same strain interface positions.

     

Marangoni convection at alcohol aqueous solutions-air interfaces

For aqueousn-heptanol solutions and in a nearly two-dimensional flow, two strikingfeatures have been detected:

1. a shift of the minimum of the surface tension
2. a discrepancy between the observed Marangoni flow velocities and the expected ones from static surface tension values.

A qualitative explanation is given.     

The approach of solutions of nonlinear diffusion equations to travelling front solutions

The paper is concerned with the asymptotic behavior as t → ∞ of solutions u(x, t) of the equation u_t—u_xx—∞;(u)=O, x∈(—∞, ∞) , in the case ∞(0)=∞(1)=0, ∞′(0)<0, ∞′(1)<0. Commonly, a travelling front solution u=U(x-ct), U(-∞)=0, U(∞)=1, exists. The following types of global stability results for fronts and various combinations of them will be given.

1. Let u(x, 0)=u ₀(x) satisfy 0≦u ₀≦1. Let \(a\_ = \mathop {\lim \sup u0}\limits_{x \to - \infty } {\text{(}}x{\text{), }}\mathop {\lim \inf u0}\limits_{x \to \infty } {\text{(}}x{\text{)}}\) . Then u approaches a translate of U uniformly in x and exponentially in time, if a_? is not too far from 0, and a₊ not too far from 1.
2. Suppose \(\int\limits_{\text{0}}^{\text{1}} {f{\text{(}}u{\text{)}}du} > {\text{0}}\) . If a _? and a ₊ are not too far from 0, but u₀ exceeds a certain threshold level for a sufficiently large x-interval, then u approaches a pair of diverging travelling fronts.
3. Under certain circumstances, u approaches a “stacked” combination of wave fronts, with differing ranges.

     

Guidance properties of nonlinear planar waveguides

We discuss the propagation of electromagnetic waves through a stratified dielectric. The ability of such a device to support guided waves depends upon the way in which the refractive index varies across the layers. In the present discussion, we show how nonlinear effects and appropriate stratification can be used to obtain any one of the following behaviours:

1. guidance occurs only at low power.
2. guidance occurs only at high power.
3. guidance occurs at all powers.
4. there is no guidance.

The situation (i) is obtained by using materials with a defocusing dielectric response, whereas the situation (ii) is obtained for suitable configurations of self-focusing materials. The situations (iii) and (iv) can be obtained by using either defocusing or self-focusing materials. By seeking solutions of a particular form, we reduce the problem to the study of solutions in the Sobolev space H ¹(?) of a second-order differential equation. The discussion of defocusing nonlinearities is based in the study of the global behaviour of the branch of solutions bifurcating from a simple eigenvalue. For self-focusing nonlinearities we use a variational approach.     

Flow of a dipolar fluid between parallel plates

The flow of a dipolar fluid between two parallel plates with and without heat transfer is studied. The following cases are discussed:

1. Isothermal flow due to the relative motion of the plates,
2. Isothermal flow due to a constant pressure gradient with the plates at rest,
3. Nonisothermal flow with linearly varying plate temperatures.

Case (ii) is of particular interest to the experimentalists as it shows the effect of the material constants even when there are no externally applied dipolar tractions on the plates.     

Stress concentrations around horizontal circular tunnels

The problem of axial variation of stress concentrations at the periphery and normal to the axis of a circular tunnel is solved by means of the three-dimensional photoelasticity technique, under the following conditions:

1. The center lines of two horizontal tunnels of equal diameter (2r) are separated by a distanceK and include an angle α.
2. K and α assume values of 0, 3r, 7/2r, 4r and 30 deg, 60 deg, 90 deg, respectively.
3. The tunnels are located in a uniform, uniaxial stress field normal to the axes of the tunnels.

     

A logrithmic bound on the location of the poles of the scattering matrix

It is shown that the complex poles z of the scattering matrix satisfy the inequality: Im z≧a+b log ¦z¦, b>0, in three instances of classical scattering in three space dimensions described by the wave equation u_{t t}?c²Δu+qu=0.

1. c and q smooth with c=1 and q=0 for ¦x¦>p, all rays going to infinity, and the energy form positive definite.
2. c=1 and q=0 outside of a convex body on which u=0.
3. c=1, q bounded and measurable, q=0 for ¦x¦>p, and the energy form not necessarily positive definite.

     

Determination of stress-intensity factors from photoelastic data with applications to surface-flaw problems

ATaylor-series correction to the maximum inplane shear stress was studied as a means of extending the data zone in photoelastic determination of stress-intensity factors beyond the singular region of a two-degree-of-freedom analysis. Convergence properties were obtained by comparing with several complete two-dimensional solutions. Experiments were performed on two kinds of three-dimensional problems, plates containing surface flows in both bending and extension. Results were analyzed by both a two-degree-of-freedom and aTaylor-series correction method (TSCM). Results were compared to theories of F. W. Smith and A. S. Kobayashi and R. C. Shah. It was concluded that:

1. The TSCM program converges rapidly to accurateK _I values and will accommodate the scatter inherent in experimental data if the series is properly truncated.
2. The TSCM program is essentially equivalent to the two-parameter representation when only the crack-surface effects dominate.
3. When effects other than crack surfaces are important, TSCM requires more terms but still predictsK _I with reasonable accuracy.

     

Die Berechnung des Blasenwachstums beim Verdampfen von Flüssigkeiten an Heizflächen als numerische Lösung der Erhaltungsgleichungen

A numerical procedure on the basis of the Marker and Cell-method [1] was developed in order to solve the conservation equations for mass, momentum and energy for the case of bubble growth on a heating surface. This procedure was used to calculate steam bubble growth on a horizontal stainless steel heating surface under saturated pool boiling conditions at a system pressure of 1 bar and different superheatings. The essential results obtained are:

-Good agreement was found between calculations and experiments concerning bubble growth rates, bubble shape and temperature field in the liquid surrounding the bubble.

-During its growth the bubble penetrates the temperature boundary layer formed in the liquid on the heating surface, simultaneously liquid is displaced aside.

-The microlayer evaporation fraction of the total bubble growth increases with growth time from 20 % to 50%.

     

Oscillation modes of laval nozzle flow

Experimental investigations of Laval nozzle flow show for relatively low supply to exit pressure ratios, which correspond to shock wave positions close to the nozzle throat, three different, oscillatory instabilities.

1. Shock pattern oscillations where the root of a λ-like shock front remains nearly in constant position, but where the proportion between the normal part and the oblique part of the shock changes periodically.
2. Shock wave and separation bubble oscillations where the motion of the shock wave is accompanied by displacements of the separation bubble.
3. Flow rate oscillations where the shock waves leave periodically through the nozzle throat in upstream direction.

     

A fast and accurate graphical method for designing an air conditioning cooling coil

The aim of this paper is to develop a fast and simple accurate graphical method for designing the required cooling coil for an air conditioning system in which both the sensible and latent heat are transferred. The method has the following advantages:

1. Direct solution utilizing only Psychrometer charts.
2. Solution of combined heat and mass transfer problems occuring at pressures other than that of atmosphere.
3. Direct determination of boundary temperature at which dehumidification begins for that type of problem where a portion of the surface is in a dry condition.
4. Less effort required than trial and error method in determination of air condition leaving a counterflow coil of a given area.

The calculated values of the method showed a good agreement with the experimental results. The average deviation for the total heat is about + 15 % and for sensible heat is + 9 %, which are of positive nature and on the safe side for practical design purposes. The method is also useful for similar practical application.     

Experimental determination of viscosity of rocks

An important parameter involved in the viscoelastic deformation of structural materials is the coefficient of “solid” viscosity. Determination of this parameter is necessary, if it is to be used in structural design. This paper deals with pertinent analytical considerations concerning solid viscosity and describes the procedures followed in the determination of parameters for structural and true viscosity of a Queenston limestone. The following three techniques were used:

1. Relaxation technique
2. Uniaxial compressive loading
3. Cantilever-beam loading

The results obtained are in close mutual agreement except for (2) above, where experimental conditions were different from those in (1) and (3). A quasi-periodic behavior of strain is indicated. It has been shown that the solid viscous parameter is a transient property and may depend on such factors as applied load, time, grain size, grain-packing in a material, and the direction of testing. It has been concluded that coefficients of true and structural solid viscosity of materials can be determined for a given set of conditions.     

Effects of artificial cracks and poisson's ratio upon photoelastic stress-intensity determination

A series of stress-freezing photoelastic experiments were performed with multiple replications upon edge-cracked strips for three types of “cracks” in current use:

1. Rectangular slots 0.152 mm wide,
2. 1.59-mm-wide slots terminating in a 30-deg vee notch of approximately 0.025-mm root radius, and
3. Natural cracks (approximately 0.0025-mm root radius).

Stress-intensity results were compared with the Gross-Srawley analysis; in addition (1) was compared with Savin's solution. It was concluded that (2) and (3) yield the same results but (1) was slightly higher. Both (2) and (3) were about 12 percent higher than the Gross-Srawley results. This is shown to be related to a Poisson's ratio effect.     

Hamilton principle for an inviscid compressible fluid model in Eulerian coordinates

It is shown that the well-known variational principles for the ideal compressible fluid model in Eulerian coordinates have the following deficiencies:

1. They are not related to the corresponding variational principles in Lagrangian coordinates;
2. The variation procedure in these variational problems does not lead to the equations of motion themselves in the Euler form; rather it leads to relations which correspond to definite classes of solutions of the Euler equations. Here allowance for the equations of the constraints imposed by the adiabaticity and continuity conditions limits the region of application of these variational principles to only potential flows;
3. More general results, involving flows other than potential, are achieved by artificial selection of certain additional constraint conditions imposed on the quantities being varied, and in this case additional clarification is required to ascertain whether any inviscid compressible fluid flow is the extremum of the corresponding variational problem.

A new formulation of the Hamilton principle for the inviscid compressible fluid in Eulerian coordinates is suggested which is free from these deficiencies.     

Some pending problems in polymer melt rheology,as seen from the point of view of Doi's slip-link model

An exposition is given of results, as obtained with the aid of Doi's sliplink model, being considered as the most simple version of the famous “reptation model”. It turns out that this model which exhibits three distinct phases of relaxation (an extremely fast phase, an equilibration phase and a slow disengagement phase) is capable of explaining several peculiar features of polymer melt rheology:

1. The molecular mass dependence of the breadth of the rubber plateau in the storage modulus, of the zero-shear viscosity and of the normal-stress coefficients.
2. The molecular mass independence of the equilibrium (shear and tensile) compliances for monodisperse polymers (semi-quantitative prediction).
3. The seemingly contradictory sensitivity of these compliances for the breadth of the molecular mass distribution.
4. The critical value of the shear stress at which melt fracture occurs in capillary flow.
5. An equilibration phase in tensile experiments on unvolcanized rubber.

In this evaluation optical (flow birefringence) measurements are preferentially used.     

Variable isochromatic/isopachic-fringe visibility in photoholoelasticity

This paper presents a new “hybrid” method whereby the ratio of the isochromatic-fringe visibility/isopachic-fringe visibility may be easily and continuously varied. This simple procedure merely combines a conventional polariscope with a holographic system. A variable beam splitter permits an incoherent superposition of the reconstruction of a doubly exposed hologram with real-time isochromatics, either dark or light field. By varying the ratios of the above two, in the image plane, numerous interesting results may be obtained including:

1. Isochromatics only, without errors in position
2. Isochromatics-isopachic fringes identical to those obtained through classical interferometry
3. Isochromatic-isopachic fringes whereby the amplitude modulation between the two may be minimized
4. Continuously variable isopachic/isochromatic-fringe visibility.

     

On the synthesis of acoustic sources with controllable near fields

In this paper we present a strategy for the synthesis of acoustic sources with controllable near fields in free space and constant depth homogeneous ocean environments. We first present the theoretical results at the basis of our discussion and then, to illustrate our findings we focus on the following three particular examples:

• 1.acoustic source approximating a prescribed field pattern in a given bounded sub-region of its near field.
• 2.acoustic source approximating different prescribed field patterns in given disjoint bounded near field sub-regions.
• 3.acoustic source approximating a prescribed back-propagating field in a given bounded near field sub-region while maintaining a very low far field signature.

For each of these three examples, we discuss the optimization scheme used to approximate their solutions and support our claims through relevant numerical simulations.     

C 1 Solutions for Semi-Implicit Systems of Differential Equations

The present note is a continuation of the author??s effort to study the existence of continuously differentiable solutions to the semi-implicit system of differential equations (1) $$f(x^{\prime}(t)) = g(t, x(t))$$ (2) $$\quad x(0) = x_0,$$ where

${\quad\Omega_g \subseteq \mathbb{R} \times\mathbb{R}^n}$ is an open set containing (0, x ₀) and ${g:\Omega_g \rightarrow\mathbb{R}^n}$ is a continuous function,

${\quad\Omega_f \subseteq \mathbb{R}^n}$ is an open set and ${f:\Omega_f\rightarrow\mathbb{R}^n}$ is a continuous function.

The transformation of (1)?C(2) into a solvable explicit system of differential equations is trivial if f is locally injective around an element ${\gamma\in \Omega_f\cap f^{-1}(g(0,x_0))}$ . In this paper, we study (1)?C(2) when such a translation is not possible because of the inherent multivalued nature of f ^?1.     

Some characteristics of supersonic flow of an electrically conductive gas in an MHD-channel

We study supersonic flows of an electrically conductive gas in crossed electric and magnetic fields [1] in the presence of shock waves. It is shown that three steady flow regimes can exist, and that these are defined by the electrical conductivity of the gas as a function of temperature and density.

1. The normal regime is characterized by a tendency for the shock to move toward the channel entrance on increase of the static pressure at the channel exit. The steady regime of this type exists and is stable.
2. The anomalous regime (formally constructed) is characterized by a tendency for the shock to move toward the exit on increase of the static pressure at the channel exit. This regime is unstable and the flow in the MHD-channel may be either entirely supersonic or entirely subsonic.
3. The limiting (boundary) regime is intermediate between the normal and anomalous regimes and is characterized by the fact that the stationary position of the shock wave and its amplitude are not uniquely defined. Steady flow in this case is not unique.

This study involves formal construction both of the solution to the steady-state problem and the corresponding nonsteady-state problem [4]. The establishment of a steady regime in the solution of the unsteady problem, is at the same time, a verification of its stability.