Numerical modelling of shocks in solids with elastic-plastic conditions

Results are presented for a range of one- and two-dimensional shock-wave problems in elastic-plastic and hydrodynamic metals. These problems were solved numerically using the Flux-Corrected Transport (FCT) technique which achieves high resolution without non-physical oscillations, especially at shock fronts, and has not been used before in elastic-plastic solids. The two-dimensional problems were solved using both operator- and non-operator-split techniques to highlight the significant differences between these techniques when solving shock-wave problems in elastic-plastic solids. Comparisons of the elastic-plastic solutions with the hydrodynamic solutions are made and illustrate the importance of including elastic-plastic conditions when modelling the behaviour of solids. Also, the errors in these solutions that are due to the initial conditions are discussed in detail.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Plastic size effect analysis of lamellar composites using a discrete dislocation plasticity approach

Plastic size effect analysis of lamellar composites consisting of elastic and elastic-plastic layers is performed using a discrete dislocation plasticity approach, which is based on applying periodic homogenization to the superposition method for discrete dislocation plasticity. In this approach, the decomposition of displacements into macro and perturbed components circumvents the calculation of superposing displacement fields induced by dislocations in an infinitely homogeneous medium, resulting in two periodic boundary value problems specialized for the analysis of representative volume elements. The present approach is verified by analyzing a model lamellar composite that includes edge dislocations fixed at interfaces. The plastic size effects due to dislocation pile-ups at interfaces are also analyzed. The analysis shows that, strain hardening in elastic-plastic layers arises depending on two factors, namely the thickness and stiffness of elastic layers; and the gap between slip planes in adjacent elastic-plastic layers. In the case where the thickness of elastic layers is several dozen nm, strain hardening in elastic-plastic layers is restrained as the gap of the slip planes decreases. This particular effect is attributed to the long range stress due to pile-ups in adjacent elastic-plastic layers.     

Navier-Stokes equations with imposed pressure and velocity fluxes

Boundary value problems for Stokes and Navier-Stokes equations with non-standard boundary conditions are studied. Included is the case where the pressure or its normal derivative is given on some part of the boundary or the pressure is given up to a constant but given velocity flux. First, a variational formulation is introduced which is shown to be equivalent to the Stokes equations with the non-standard boundary conditions under consideration. The existence and uniqueness of the solution of the variational problem are studied. Secondly, most of the results obtained for the Stokes equations are extended to the case of the Navier-Stokes equations. The final section is devoted to numerical experiments, flows in pipes and physiological flows.     

A Method Based on the Radon Transform for Three-Dimensional Elastodynamic Problems of Moving Loads

An integral transform procedure is developed to obtain fundamental elastodynamic three-dimensional (3D) solutions for loads moving steadily over the surface of a half-space. These solutions are exact, and results are presented over the entire speed range (i.e., for subsonic, transonic and supersonic source speeds). Especially, the results obtained here for the tangential loads (one of these loads is along the direction of motion and the other is orthogonal to the direction of motion) are quite new in the literature. The solution technique is based on the use of the Radon transform and elements of distribution theory. It also fully exploits as auxiliary solutions the ones for the corresponding plane-strain and anti-plane shear problems (the latter problems are 2D and uncoupled from each other). In particular, it should be noticed that the plane-strain problem here is completely analogous to the original 3D problem, not only with respect to the field equations but also with respect to the boundary conditions. This makes the present technique more advantageous than other techniques, which require first the determination of a fictitious auxiliary plane-strain problem through the solution of an integral equation. Our approach becomes particularly simple when there is no angular dependence in the boundary conditions (i.e., when axially symmetric problems regarding their boundary conditions are considered). On the contrary, no such constraint needs to be fulfilled as regards the material response (and, therefore, the governing equations of the problem) and/or also possible loss of axisymmetry due to the generation of shock (Mach-type) waves in the medium. Therefore, the solution technique can easily deal with general 3D problems having a largely arbitrary radial dependence in the boundary conditions and involving: (i) material anisotropy in static and dynamic situations, and (ii) asymmetry caused by changes in the nature of governing PDEs due to the existence of different velocity regimes (super-Rayleigh, transonic, supersonic) in dynamic situations. This revised version was published online in July 2006 with corrections to the Cover Date.     

On pressure boundary conditions for thermoconvective problems

Solving numerically hydrodynamical problems of incompressible fluids raises the question of handling first order derivatives (those of pressure) in a closed container and determining its boundary conditions. A way to avoid the first point is to derive a Poisson equation for pressure, although the problem of taking the right boundary conditions still remains. To remove this problem another formulation of the problem has been used consisting of projecting the master equations into the space of divergence‐free velocity fields, so pressure is eliminated from the equations. This technique raises the order of the differential equations and additional boundary conditions may be required. High‐order derivatives are sometimes troublesome, specially in cylindrical coordinates due to the singularity at the origin, so for these problems a low order formulation is very convenient. We research several pressure boundary conditions for the primitive variables formulation of thermoconvective problems. In particular we study the Marangoni instability of an infinite fluid layer and we show that the numerical results with a Chebyshev collocation method are highly correspondent to the exact ones. These ideas have been applied to linear stability analysis of the Bénard–Marangoni (BM) problem in cylindrical geometry and the results obtained have been very accurate. Copyright © 2002 John Wiley & Sons, Ltd.     

On pressure boundary conditions for the incompressible Navier-Stokes equations

The pressure is a somewhat mysterious quantity in incompressible flows. It is not a thermodynamic variable as there is no ‘equation of state’ for an incompressible fluid. It is in one sense a mathematical artefact—a Lagrange multiplier that constrains the velocity field to remain divergence-free; i.e., incompressible—yet its gradient is a relevant physical quantity: a force per unit volume. It propagates at infinite speed in order to keep the flow always and everywhere incompressible; i.e., it is always in equilibrium with a time-varying divergence-free velocity field. It is also often difficult and/or expensive to compute. While the pressure is perfectly well-defined (at least up to an arbitrary additive constant) by the governing equations describing the conservation of mass and momentum, it is (ironically) less so when more directly expressed in terms of a Poisson equation that is both derivable from the original conservation equations and used (or misused) to replace the mass conservation equation. This is because in this latter form it is also necessary to address directly the subject of pressure boundary conditions, whose proper specification is crucial (in many ways) and forms the basis of this work. Herein we show that the same principles of mass and momentum conservation, combined with a continuity argument, lead to the correct boundary conditions for the pressure Poisson equation: viz., a Neumann condition that is derived simply by applying the normal component of the momentum equation at the boundary. It usually follows, but is not so crucial, that the tangential momentum equation is also satisfied at the boundary.     

General forms of elastic-plastic matching equations for mode-Ⅲ cracks near crack line

Crack line analysis is an effective way to solve elastic-plastic crack problems. Application of the method does not need the traditional small-scale yielding conditions and can obtain sufficiently accurate solutions near the crack line. To address mode- Ⅲ crack problems under the perfect elastic-plastic condition, matching procedures of the crack line analysis method axe summarized and refined to give general forms and formulation steps of plastic field, elastic-plastic boundary, and elastic-plastic matching equations near the crack line. The research unifies mode-III crack problems under different conditions into a problem of determining four integral constants with four matching equations. An example is given to verify correctness, conciseness, and generality of the procedure.     

Effective Flux Boundary Conditions for Upscaling Porous Media Equations

We introduce a new algorithm for setting pressure boundary conditions in subgrid simulations of porous media flow. The algorithm approximates the flux in the boundary cell as the flux through a homogeneous inclusion in a homogeneous background, where the permeability of the inclusion is given by the cell permeability and the permeability of the background is given by the ambient effective permeability. With this approximation, the flux in the boundary cell scales with the cell permeability when that permeability is small, and saturates at a constant value when that permeability is large. The flux conditions provide Neumann boundary conditions for the subgrid pressure. We call these boundary conditions effective flux boundary conditions (EFBCs). We give solutions for the flux through ellipsoidal inclusions in two and three dimensions, assuming symmetric tensor permeabilities whose principal axes align with the axes of the ellipse. We then discuss the considerations involved in applying these equations to scale up problems in geological porous media. The key complications are heterogeneity, fluctuations at all length scales, and boundary conditions at finite scales.     

General forms of elastic-plastic matching equations for mode-Ⅲ cracks near crack line

Crack line analysis is an effective way to solve elastic-plastic crack problems.Application of the method does not need the traditional small-scale yielding conditions and can obtain sufficiently accurate solutions near the crack line. To address mode-Ⅲ crack problems under the perfect elastic-plastic condition,matching procedures of the crack line analysis method are summarized and refined to give general forms and formulation steps of plastic field,elastic-plastic boundary,and elastic-plastic matching equations near the crack line. The research unifies mode-Ⅲ crack problems under different conditions into a problem of determining four integral constants with four matching equations.An example is given to verify correctness,conciseness,and generality of the procedure.     

定态火焰在可燃预混气中产生的压力波

火焰在可燃预混气中传播时，在火焰面前方产生一道压力波。忽略点火及火焰的初期加速，仅考虑火焰达到稳定传播速度的情况。用Ｏｐｅｎｈｅｉｍ自相似解分析流场，得到相应的控制方程及定解条体；用自适应步长的四阶Ｒｕｎｇｅ－Ｋｕｔｔａ法对方程积分，讨论了流场压力波结构及弱激波近似声波解；认为火焰为间断面，能量释放在火焰面后瞬时完成。利用火焰面两侧的能量关系，得到了火焰位置、燃速及对应Ｃ－Ｊ条件的火焰位置、Ｃ－Ｊ燃速。     

A staggered grid‐based explicit jump immersed interface method for two‐dimensional Stokes flows

In this paper the explicit jump immersed interface method (EJIIM) is applied to stationary Stokes flows. The boundary value problem in a general, non‐grid aligned domain is reduced by the EJIIM to a sequence of problems in a rectangular domain, where staggered grid‐based finite differences for velocity and pressure variables are used. Each of these subproblems is solved by the fast Stokes solver, consisting of the pressure equation (known also as conjugate gradient Uzawa) method and a fast Fourier transform‐based Poisson solver. This results in an effective algorithm with second‐order convergence for the velocity and first order for the pressure. In contrast to the earlier versions of the EJIIM, the Dirichlét boundary value problem is solved very efficiently also in the case when the computational domain is not simply connected. Copyright © 2007 John Wiley & Sons, Ltd.     

State-space approach to two-temperature generalized thermoelasticity without energy dissipation of medium subjected to moving heat source

In this work,a model of two-temperature generalized thermoelasticity without energy dissipation for an elastic half-space with constant elastic parameters is constructed.The Laplace transform and state-space techniques are used to obtain the general solution for any set of boundary conditions.The general solutions are applied to a specific problem of a half-space subjected to a moving heat source with a constant velocity.The inverse Laplace transforms are computed numerically,and the comparisons are shown in figures to estimate the effects of the heat source velocity and the two-temperature parameter.     

Rolling without Slip on a Transversely Isotropic Thermoelastic Half-space

Rolling without slip by a rigid cylinder on a transversely isotropic, coupled thermoelastic half-space at constant subcritical speed is studied. The cylinder is of infinite length, surface heat convection is neglected, and a dynamic steady state of plane strain is treated. The unmixed problem of traction applied to a translating surface strip is addressed first. A robust asymptotic form of the exact transform solution, valid when Fourier heat conduction dominates any thermal relaxation effect, is extracted, and inverted analytically. Use of material characterization and identification of parameters that vanish in the isotropic limit or are invariant under an isothermal–thermoelastic transformation result in compact full-field solutions. These expressions are used to construct analytical solutions that satisfy the mixed boundary value problem and auxiliary conditions of rolling contact. For the hexagonal material zinc, calculations are made for contact zone width and temperature increases near onset of zone yield. Mathematics Subject Classifications (2000) 73B30, 73C25, 73C30, 73C35.     

Invariant solutions of the Westervelt model of nonlinear hydroacoustics without dissipation

We study three-dimensional Westervelt model of a nonlinear hydroacoustics without dissipation. We received all of its invariant submodels. We studied all invariant submodels described by the invariant solutions of rank 0 and 1. All invariant solutions of rank 0 and 1 are found either explicitly, or their search is reduced to the solution of the nonlinear integral equations. With a help of these invariant solutions we researched: (1) a propagation of the intensive acoustic waves (self-similar, axisymmetric, planar and one-dimensional) for which the acoustic pressure and a speed of its change, or the acoustic pressure and its derivative in the direction of one of the axes are specified at the initial moment of the time at a fixed point , (2) a spherically symmetric ultrasonic field for which the acoustic pressure and a speed of its change, or the acoustic pressure and its radial derivative are specified at the initial moment of the time at a fixed point. Solving of the boundary value problems describing these processes is reduced to the solving of nonlinear integral equations. We are established the existence and uniqueness of solutions of these boundary value problems under some additional conditions. Mechanical relevance of the obtained solutions is as follows: (1) these solutions describe nonlinear and diffraction effects in ultrasonic fields of a special kind, (2) these solutions can be used as a test solutions in the numerical calculations performed in studies of ultrasonic fields generated by powerful emitters. We found all the conservation laws of the first order for the Westerveld equation written in dimensionless variables.     

Constant velocity motion of a strip die on the boundary of a viscoelastic base

We consider a contact problem on the interaction of a rigid strip die with the boundary of a viscoelastic base. We assume that the die moves at a constant velocity on this boundary and is indented into it by a constant normal force. Friction in the die—surface contact region is neglected. The die base is corrugated in the direction perpendicular to the direction of motion. At the first stage, we determine the displacement of the base boundary due to the normal load applied to it. Then, at the second stage, we derive the integral equation of the contact problem for determining the contact pressure. At the third stage, we construct an approximate solution of this integral equation by using the modified Multhopp—Kalandiya method.     

On the critical problem of F. D. pressure treatment for laminar flows confined by permeable walls

In the present work the viscous (low Reynolds) flow in plane ducts confined by permeable walls has been studied. A simple model of the filtrating walls has been used, with the normal velocity component proportional to the pressure jump across the wall, resulting in a non-standard boundary value Navier-Stokes problem. A critical analysis of the appropriate boundary condition and pressure problem has led to the conclusions of employing a simple explicit finite volume approach, and of avoiding the use of higher order finite difference schemes. In this paper a special emphasis on the structure of the involved computational matrices has been given to illustrate the chosen algorithm. The latter yields a steady state solution that is second order accurate in space, and it has an accuracy in time of order ≤ Δt (the time step), due to the explicit treatment of the velocity boundary conditions along the membrane. The model has been tested to study the effects of the inlet/outlet conditions, Reynolds number and filtrating wall constant.     

On modeling the aggregate response of laminated thermoelastic media

Sufficient conditions are given that determine when solutions to general boundary value problems for inhomogeneous anisotropic laminated thermoelastic media are approximated, in the average, by solution of constant coefficient equations. These appropriate constant coefficients are computed and shown to be independent of the shape of the body and the loading on the body.     

A Direct theory of viscous fluid flow in channels

This paper deals with an Eulerian formulation of the theory of directed fluid sheets appropriate for incompressible, linear viscous fluid flow in channels with arbitrary shapes for their major boundaries which may be moving or fixed. Special cases of the theory are applied to a number of two-dimensional fluid flow problems and these solutions are in general discussed for unsteady flow. Specific applications include fluid flow in a channel whose boundaries are symmetric with respect to a middle plane in the channel, subjected to time-dependent pressure gradient at one end; and to lubrication problems in a general shaped channel when one of the channel walls is a fixed plane while the other is moving with a constant velocity. Flow of a viscous fluid with a free surface over a fixed boundary is also discussed.Dedicated to J. L. Ericksen on the occasion of his Sixtieth Birthday     

Dirichlet and Neumann boundary conditions for the pressure poisson equation of incompressible flow

In a recent paper Gresho and Sani showed that Dirichlet and Neumann boundary conditions for the pressure Poisson equation give the same solution. The purpose of this paper is to confirm this (for one case at least) by numerically solving the pressure equation with Dirichlet and Neumann boundary conditions for the inviscid stagnation point flow problem. The Dirichlet boundary condition is obtained by integrating the tangential component of the momentum equation along the boundary. The Neumann boundary condition is obtained by applying the normal component of the momentum equation at the boundary. In this work solutions for the Neumann problem exist only if a compatibility condition is satisfied. A consistent finite difference procedure which satisfies this condition on non-staggered grids is used for the solution of the pressure equation with Neumann conditions. Two test cases are computed. In the first case the velocity field is given from the analytical solution and the pressure is recovered from the solution of the associated Poisson equation. The computed results are identical for both Dirichlet and Neumann boundary conditions. However, the Dirichlet problem converges faster than the Neumann case. In the second test case the velocity field is computed from the momentum equations, which are solved iteratively with the pressure Poisson equation. In this case the Neumann problem converges faster than the Dirichlet problem.