On the T-matrix for water-wave scattering problems

A rigid cylinder of infinite length is floating in the free surface of deep water. The cylinder is held fixed and a given time-harmonic wave of small amplitude is incident upon it. The corresponding linear two-dimensional boundary-value problem for a velocity potential φ is treated using the null-field method, and an expression for the T-matrix is obtained. (The T-matrix connects the diffraction potential away from the cylinder to the given incident potential.) Fundamental properties of the T-matrix are derived from considerations of energy and reciprocity. For regular wavetrains incident from the right or from the left, there are well-known relations between the corresponding reflection and transmission coefficients; these relations are recovered by specialising the equations satisfied by the T-matrix. Two extensions to water of constant finite depth are described: one uses multipole potentials whilst the other uses Havelock wavemaker functions; this second approach also leads to a new method for treating the problem of waves in a semi-infinite channel with an end-wall of arbitrary shape.     

Scattering of rayleigh lamb waves from a 2d-cavity in an elastic plate

The T-matrix method is applied to the problem of scattering of Rayleigh-Lamb modes from a twodimensional cavity in an elastic plate. A formal solution is obtained which is valid also for non-planar surfaces. Explicit expressions and numerical results are given for a plate with plane surfaces.     

A POD-based reduced-order finite difference extrapolating model with fully second-order accuracy for non-stationary Stokes equations

A proper orthogonal decomposition (POD) reduced-order finite difference (FD) extrapolating model is established for the channel flow with local expansion denoted by non-stationary Stokes equations. The POD-based reduced-order numerical model to produce the solutions on the time span [T₀, T] (T₀ ? T) are obtained by extrapolation and iteration from the very short time span [0, T₀] information. The guides to choose the number of POD basis and renew POD basis are provided, and an implementation for solving the POD-based reduced-order FD extrapolating model is given. Some numerical experiments are used to show that the POD-based reduced-order FD extrapolating model is feasible and efficient for simulating the channel flow with local expansion.     

Scattering of acoustic waves by elastic and viscoelastic obstacles immersed in a fluid

A scattering or T-matrix approach is presented for studying the scattering of acoustic waves by elastic and viscoelastic obstacles immersed in a fluid. A Kelvin-Voigt model is used to obtain the complex elastic moduli of the viscoelastic solid. The T-matris formulation is somewhat complicated because the wave equations and fields are quite different in the solid and fluid regions and are coupled by continuity conditions at the interface. We have obtained fairly extensive numerical results for prolate and oblate spheroids for a variety of scattering geometries. The backscattering, bistatic, absorption and extinction cross-section are presented as a function of the frequency of the incident wave.     

Adaptive numerical modeling using the hierarchical Fup basis functions and control volume isogeometric analysis

A novel adaptive algorithm that is based on new hierarchical Fup (HF) basis functions and a control volume formulation is presented. Because of its similarity to the concept of isogeometric analysis (IGA), we refer to it as control volume isogeometric analysis (CV-IGA). Among other interesting properties, the IGA introduced k-refinement as advanced version of hp-refinement, where every basis function of the nth order from one resolution level are replaced by a linear combination of more basis functions of the n+1th order at the next resolution level. However, k-refinement can be performed only on whole domain, while local adaptive k-refinement is not possible with classical B-spline basis functions. HF basis functions (infinitely differentiable splines) satisfy partition of unity, and they are linearly independent and locally refinable. Their main feature is execution of the adaptive local hp-refinement because any basis function of the nth order from one resolution level can be replaced by a linear combination of more basis functions of the n+1th order at the next resolution level providing spectral convergence order. The comparison between uniform vs hierarchical adaptive solutions is demonstrated, and it is shown that our adaptive algorithm returns the desired accuracy while strongly improving the efficiency and controlling the numerical error. In addition to the adaptive methodology, a stabilization procedure is applied for advection-dominated problems whose numerical solutions “suffer” from spurious oscillations. Stabilization is added only on lower resolution levels, while higher resolution levels ensure an accurate solution and produce a higher convergence order. Since the focus of this article is on developing HF basis functions and adaptive CV-IGA, verification is performed on the stationary one-dimensional boundary value problems.     

A method for measuring mode I crack tip constraint under static and dynamic loading conditions

A novel experimental technique for measuring crack tipT-stress, and hence in-plane crack tip constraint, in elastic materials has been developed. The method exploits optimal positioning of stacked strain gage rosette near a mode I crack tip such that the influence of dominant singular strains is negated in order to determineT-stress accurately. The method is demonstrated for quasi-static and low-velocity impact loading conditions and two values of crack length to plate width ratios (a/W). By coupling this new method with the Dally-Sanford single strain gage method for measuring the mode I stress intensity factorK _I , the crack tip biaxiality parameter is also measured experimentally. Complementary small strain, static and dynamic finite element simulations are carried out under plane stress conditions. Time histories ofK _I andT-stress are computed by regression analysis of the displacement and stress fields, respectively. The experimental results are in good agreement with those obtained from numerical simulations. Preliminary data for critical values ofK _I and β for dynamic experiments involving epoxy specimens are reported. Dynamic crack initiation toughness shows an increasing trend as β becomes more negative at higher impact velocities.     

Exact analytic solutions for free convection boundary layers on a heated vertical plate with lateral mass flux embedded in a saturated porous medium

 The problem of the self-similar boundary flow of a “Darcy-Boussinesq fluid” on a vertical plate with temperature distribution T _w(x) = T _∞+A·x ^λ and lateral mass flux v _w(x) = a·x ^(λ−1)/2, embedded in a saturated porous medium is revisited. For the parameter values λ = 1,−1/3 and −1/2 exact analytic solutions are written down and the characteristics of the corresponding boundary layers are discussed as functions of the suction/ injection parameter in detail. The results are compared with the numerical findings of previous authors. Received on 8 March 1999     

Elastic stresses in bodies with elliptic,parabolic, and tunnel cracks

The paper gives explicit expressions of the elastic T-stress components T _I, T _II, and T _III for an elliptic crack in an unbounded body under uniform pressure and bending and expressions of all the T-stress components for parabolic and tunnel cracks under uniform loading. These formulas are derived by analyzing the asymptotic behavior of the stress components near the crack front using special harmonic functions. The dependence of the T-stresses on Poisson’s ratio, semiaxes and parametric angle of the elliptic crack is studied. The expressions of T _I, T _II, and T _III for a penny-shaped crack under arbitrary uniform pressure and bending follow as a special case from the respective expressions for an elliptic crack __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 8, pp. 57–70, August 2007.     

Stokes eigenfunctions and Galerkin projection of the disturbance equations in plane Poiseuille flow: a systematic analytical approach

In the present paper, the L ²-normalized Stokes eigenfunctions for plane Poiseuille flow, which form an orthonormal functional basis for the space of disturbances, are written in a general exponential form. Then, the evolution equations for the disturbances are Galerkin-projected on the considered basis functions, and all the terms of the resulting dynamical system are expressed systematically in analytical form. Finally, a numerical example is given in which the proposed basis functions are used for the simulation of the time evolution of the critical disturbance predicted by the energetic stability theory.     

A new formula for the C-matrix in the Somigliana identity

By making use of a convenient decomposition of the fundamental tractions, a new formula for the C-matrix in the Somigliana identity for a three- or two-dimensional elastic isotropic body is derived. This kind of formula is more advantageous for analytical and computational C-matrix evaluations than the currently well-known formula. A general closed analytical formula of the C-matrix for the case of any finite number of tangent planes to the boundary of the body at a non-smooth boundary point, presented in the final section of this paper, demonstrates the usefulness of the new formula.     

A finite element variational multiscale method for incompressible flows based on the construction of the projection basis functions

In this article, we present a finite element variational multiscale (VMS) method for incompressible flows based on the construction of projection basis functions and compare it with common VMS method, which is defined by a low‐order finite element space L_h on the same grid as X_h for the velocity deformation tensor and a stabilization parameter α. The best algorithmic feature of our method is to construct the projection basis functions at the element level with minimal additional cost to replace the global projection operator. Finally, we give some numerical simulations of the nonlinear flow problems to show good stability and accuracy properties of the method. Copyright © 2011 John Wiley & Sons, Ltd.     

Scattering of elastic waves by a periodic array of cracks in 3-D space

The problem of scattering of normal incident time harmonic plane elastic waves by a co-planar periodic array of cracks in 3-D space is investigated. The scattered waves consist of a superposition of an infinite number of wave modes [M, N]^T and [M, N]^L,M. N=0, 1, 2, , but only a finite number of them are propagating wave modes. The numerical calculation has been made for rectangular cracks and P wave incidence. The reflection coefficient of [O, O] order,R ₀ ³ , has been studied in detail for various wave numbers and parameters of the geometry for the problem. The reliability of the numerical calculation has been checked by an application of the balance of rates of energies. For an elongated rectangular crack,R ₀ ³ in the corresponding 2-D problem in [2] is recovered. The dynamic stress intensity factors around the crack edge have been obtained. The results as the wave number goes to zero have been compared with those in the correspoding static case. Good agreement is observed.     

On a Free Convection Problem Over a Vertical Flat Surface in a Porous Medium

The problem of the free convection boundary-layer flow over a semi-infinite vertical flat surface in a porous medium is considered, in which the surface temperature has a constant value T₁ at the leading edge, where T₁ is above the ambient temperature, and takes a value T₂ at a given distance L along the surface, varying linearly between these two values and remaining constant afterwards. Numerical solutions of the boundary-layer equations are obtained as well as solutions valid for both small and large distance along the surface. Results are presented for the three cases, when the temperature T₂ is greater, equal or less than the ambient temperature T_∞. In the first case, T₂ > T_∞, a boundary-layer flow develops along the surface starting with a flow associated with the temperature difference T₁ − T_∞ at the leading edge and approaching a flow associated with the temperature difference T₂ − T_∞ at large distances. In the second case, T₂ = T_∞, the convective flow set up on the initial part of the surface drives a wall jet in the region where the surface temperature is the same as ambient. In the final case, T₂ < T_∞, a singularity develops in the numerical solution at the point where the surface temperature becomes T_∞. The nature of this singularity is discussed.     

Instability and long-time evolution of cnoidal wave solutions for a Benney–Luke equation

We investigate numerically the stability of periodic traveling wave solutions (cnoidal waves) for a generalized Benney–Luke equation. By using a high-accurate Fourier spectral method, we find different kinds of evolution depending on the period of the perturbation. A cnoidal wave solution with period T is orbitally stable with regard to perturbations having the same period T, within certain range of wave velocities. This is a fact proved recently by Angulo and Quintero [Existence and orbital stability of cnoidal waves for a 1D boussinesq equation, International Journal of Mathematics and Mathematical Sciences (2007), in press, doi:10.1155/2007/52020] and our numerical experiments are consistent with their theory. In the present work we show numerically that cnoidal waves with period T become unstable when perturbed by small amplitude disturbances whose period is an integer multiple of T. Particularly, if the period of the perturbation is 2T, the evolution of the deviation of the solution from the orbit of the cnoidal wave is found to be approximately a time-periodic function. In other cases, the numerical experiments indicate a non-periodic behavior.     

Thermoeconomic design and analysis of constant cross-sectional area fins

A closed-form model for the second-law-based thermoeconomic optimization of constant cross-sectional area fins, is discussed with an example problem. In this approach, different monetary values are attached to the irreversible losses caused by the finite temperature difference heat transfer (T _{0Δ T} ) and pressure drop (T _{0Δ P} ) in the fin application. In addition, a simplified closed-form solution is presented for the case when the capital cost of the fin is negligible and only operational costs are considered. To illustrate the usefulness of the present analytical approach, the simplified cost optimized results are compared with the numerical results obtained from Poulikakos and Bejan's analysis, who have assumed same monetary values for T _{0Δ T} and T _{0Δ P} . Furthermore, the influence of important fin thermal, physical, geometrical and cost parameters on the optimum Reynolds numbers Re _Dopt and Re _Lopt are presented in algebraic forms, and also graphical results are shown for the case of pin and plate fin, as examples. Received on 26 January 1998     

A compact numerical algorithm for solving the time‐dependent mild slope equation

The mild slope equation has been widely used to describe combined wave refraction and diffraction. In this study, a new numerical algorithm is developed to solve the time‐dependent mild slope equation in a second‐order hyperbolic form. The numerical algorithm is based on a compact and explicit finite difference method that is second‐order accurate in both time and space. The algorithm has the similar structure to the leap‐frog method but is constructed on three time levels for the second‐order time derivative term. The numerical model has the capability of simulating transient wave motion by correctly predicting the speed of wave energy propagation, which is important for the real‐time forecast of the arrival time of storm waves generated in the far field. The model is validated against analytical solution for wave shoaling and experimental data for combined wave refraction and diffraction over a submerged elliptic shoal on a slope (Coastal Eng. 1982; 6 :255). Lastly, the realistic scale Homma's island (Geophys. Mag. 1950; 21 :199) is studied with the use of various wave periods of T = 720s, T = 120 s, and T = 24 s. These wave periods correspond to long, intermediate, and short waves for the given topography, respectively. Comparisons are made between numerical results and existing analytical solutions in terms of the wave amplification around the island, which serves as the indicator for the potential wave runup. Excellent agreements are obtained. The model runs on a PC (Pentium IV 1.8GHz) and the computer capacity allows the computation of a mesh system up to 3000 × 3000, which is equivalent to about 150 × 150 waves or a large area of 540km × 540km for a wave train with the period of T = 60 s. Copyright 2004 John Wiley & Sons, Ltd.     

Radial Flow in a Bounded Randomly Heterogeneous Aquifer

Flow to wells in nonuniform geologic formations is of central interest to hydrogeologists and petroleum engineers. There are, however, very few mathematical analyses of such flow. We present analytical expressions for leading statistical moments of vertically averaged hydraulic head and flux under steady-state flow to a well that pumps water from a bounded, randomly heterogeneous aquifer. Like in the widely used Thiem equation, we prescribe a constant pumping rate deterministically at the well and a constant head at a circular outer boundary of radius L. We model the natural logarithm Y = lnT of aquifer transmissivity T as a statistically homogeneous random field with a Gaussian spatial correlation function. Our solution is based on exact nonlocal moment equations for multidimensional steady state flow in bounded, randomly heterogeneous porous media. Perturbation of these nonlocal equations leads to a system of local recursive moment equations that we solve analytically to second order in the standard deviation of Y. In contrast to most stochastic analyses of flow, which require that log transmissivity be multivariate Gaussian, our solution is free of any distributional assumptions. It yields expected values of head and flux, and the variance–covariance of these quantities, as functions of distance from the well. It also yields an apparent transmissivity, T _a, defined as the negative ratio between expected flux and head gradient at any radial distance. The solution is supported by numerical Monte Carlo simulations, which demonstrate that it is applicable to strongly heterogeneous aquifers, characterized by large values of log transmissivity variance. The two-dimensional nature of our solution renders it useful for relatively thin aquifers in which vertical heterogeneity tends to be of minor concern relative to that in the horizontal plane. It also applies to thicker aquifers when information about their vertical heterogeneity is lacking, as is commonly the case when measurements of head and flow rate are done in wells that penetrate much of the aquifer thickness. Potential uses include the analysis of pumping tests and tracer test conducted in such wells, the statistical delineation of their respective capture zones, and the analysis of contaminant transport toward fully penetrating wells.     

Stability Properties of POD–Galerkin Approximations for the Compressible Navier–Stokes Equations

Fluid flows are very often governed by the dynamics of a mall number of coherent structures, i.e., fluid features which keep their individuality during the evolution of the flow. The purpose of this paper is to study a low order simulation of the Navier–Stokes equations on the basis of the evolution of such coherent structures. One way to extract some basis functions which can be interpreted as coherent structures from flow simulations is by Proper Orthogonal Decomposition (POD). Then, by means of a Galerkin projection, it is possible to find the system of ODEs which approximates the problem in the finite-dimensional space spanned by the POD basis functions. It is found that low order modeling of relatively complex flow simulations, such as laminar vortex shedding from an airfoil at incidence and turbulent vortex shedding from a square cylinder, provides good qualitative results compared with reference computations. In this respect, it is shown that the accuracy of numerical schemes based on simple Galerkin projection is insufficient and numerical stabilization is needed. To conclude, we approach the issue of the optimal selection of the norm, namely the H ¹ norm, used in POD for the compressible Navier–Stokes equations by several numerical tests. Received 21 April 1999 and accepted 18 November 1999     

A new analytical model for thermal stresses in multi-phase materials and lifetime prediction methods

Based on the fundamental equations of the mechanics of solid continuum, the paper employs an analytical model for determination of elastic thermal stresses in isotropic continuum represented by periodically distributed spherical particles with different distributions in an infinite matrix, imaginarily divided into identical cells with dimensions equal to inter-particle distances, containing a central spherical particle with or without a spherical envelope on the particle surface. Consequently, the multi-particle-（envelope）- matrix system, as a model system regarding the analytical modelling, is applicable to four types of multi-phase materials. As functions of the particle volume fraction v, the inter-particle distances dl, d2, d3 along three mutually per- pendicular axes, and the particle and envelope radii, R1 and R2, respectively, the thermal stresses within the cell, are originated during a cooling process as a consequence of the difference in thermal expansion coefficients of phases rep- resented by the matrix, envelope and particle. Analytical-（experimental）-computational lifetime prediction methods for multi-phase materials are proposed, which can be used in engineering with appropriate values of parameters of real multi-phase materials.     

Viscous dissipation effects on thermal entrance region heat transfer in pipes with uniform wall heat flux

The analytical solution to Graetz problem with uniform wall heat flux is extended by including the viscous dissipation effect in the analysis. The analytical solution obtained reduces to that of Siegel, Sparrow and Hallman neglecting viscous dissipation as a limiting case. The sample developing temperature profiles, wall and bulk temperature distributions and the local Nusselt number variations are presented to illustrate the viscous dissipation effects. It is found that the role of viscous dissipation on thermal entrance region heat transfer is completely different for heating and cooling at wall. In the case of cooling at wall, a critical value of Brinkman number, Br _c=−11/24, exists beyond which (−11/24<Br<0) the fluid bulk temperature will always be less than the uniform entrance temperature indicating the predominance of cooling effect over the viscous heating effect. On the other hand, with Br < Br _c the bulk temperature T _b will approach the wall temperature T _w at some downstream position and from there onward the bulk temperature T _b becomes less than the wall temperature T _w with T _w > B _b > T ₀ indicating overall heating effect for the fluid. The numerical results for the case of cooling at wall Br < 0 are believed to be of some interest in the design of the proposed artctic oil pipeline.