Elasticity solutions for functionally graded annular plates subject to biharmonic loads

Based on England’s expansion formula for displacements, the elastic field in a transversely isotropic functionally graded annular plate subjected to biharmonic transverse forces on its top surface is investigated using the complex variables method. The material parameters are assumed to vary along the thickness direction in an arbitrary fashion. The problem is converted to determine the expressions of four analytic functions α (ζ), β (ζ), ? (ζ) and ψ (ζ) under certain boundary conditions. A series of simple and practical biharmonic loads are presented. The four analytic functions are constructed carefully in a biconnected annular region corresponding to the presented loads, which guarantee the single-valuedness of the mid-plane displacements of the plate. The unknown constants contained in the analytic functions can be determined from the boundary conditions that are similar to those in the plane elasticity as well as those in the classical plate theory. Numerical examples show that the material gradient index and boundary conditions have a significant influence on the elastic field.     

Miscible Displacements in Porous Media with Time-Dependent Injection Velocities

Miscible displacements in homogeneous porous media involving time-dependent injection velocities \(U(t)\) are analyzed. The displacements consist of periodic cycles that involve alternating stages of injection and production or of injection and soaking. Results of a linear stability analysis revealed that the growth rate of disturbances follows the overall trends of the velocity \(U(t)\) but with noticeable differences in periods of transition from production to injection. Furthermore, the growth rate of the time-dependent velocity was found to be smaller than that of a constant injection with a velocity equal to that corresponding to the minimum of \(U(t)\) while an overshoot was observed with respect to a displacement with a constant injection velocity equal to the maximum of \(U(t)\) . Nonlinear simulations revealed that the dynamics of the flow can be drastically changed from those of the corresponding constant injection velocity and the changes depend on the period of the cycles, the amplitude of the velocity and on whether the displacement is initiated through an injection or a soaking stage. The enhancement or attenuation of the instability in comparison with the constant injection velocity increases with the cycles’ period and the amplitude of the velocity in the injection stage, with the effects of the former being more prominent. It was also found that, beyond a certain critical cycle period, it is possible to observe instability and fingering in the case of time-dependent displacements that actually result in a net zero flow.     

A Dynamic Method for the Residual Stress Measurement During Polymer Crystallization

A new dynamic method based on bilayer system is proposed to characterize the residual stress formation during the crystallization of semi-crystalline polyethylene glycol 10000 (PEG10000). The resin is coated on a solid polymeric film to create a filmsubstrate compound. Its temperature field and dynamic deflection are monitored by synchronized optical and thermography cameras. The crystallization kinetics is first characterized from the former information. Then a simple dynamic model is proposed to relate the dynamic deflection with crystallization process. Residual stresses are established and in the range of 0–2.1 MPa. The generation of residual stresses is due to the edge constraints of the cantilever beam and to the increase of viscosity during solidification that allows the polymer to carry tensile loadings. The spherulite impingement is found to be important for this period from a microscopic view. Boundary condition should be well controlled to steer residual stresses. Such method is promising to measure residual stresses at the micro-scale for polymers to be spread on a flexible substrate and can mimic different mechanical situations of interest.     

Natural convections in conjugated single and double enclosures

The natural convection in single and double conjugated enclosures are numerically investigated. The single and double enclosures are formed by low conductance walls with finite thickness. The outside vertical surfaces of the conducting walls are of the third kind of boundary condition while the top and bottom outside surfaces are adiabatic. The problem studied is characterized by a dominant horizontal temperature gradient and the thermal boundary conditions at the cavity surfaces can not be specified in priori. Numerical results reveal the characteristics in such kind of enclosures and show the importance of the thermal boundary conditions on the natural convection in enclosures. It is also found that the natural convections in the conjugated double enclosures are basically the same, with a major difference in their fluid temperature levels.     

Combined heat and mass transfer by laminar natural convection from a vertical plate

Simultaneous heat and mass transfer in buoyancy-induced laminar boundary-layer flow along a vertical plate is studied for any ratio of the solutal buoyancy force to the thermal buoyancy force by using a new similarity transformation. The effects of the buoyancy ratio and Lewis number on the rates of heat and mass transfer are presented explicitly for most practical gaseous solutions (Pr=0.7, 0.21≤Sc≤2.1) and aqueous solutions (Pr=7, 140≤Sc≤1400). Very accurate correlations of the mass transfer and heat transfer rates are developed for the cases of single and combined buoyancy forces.     

The effects of micro-defects and crack on the mechanical properties of metal fiber sintered sheets

The effects of three types of defect (i.e., two micro defects—broken fibers and separation of fiber joints and one macro defect—crack) on the mechanical properties of porous metal fiber sintered sheets (MFSSs) are investigated by a combination of numerical simulation, analytical modeling, and experimental test. All simulations are based upon the previously developed micromechanics random beam model (Jin et al., 2013). Broken fibers are realized by removing cell edges (i.e., fibers between two joints) in an otherwise perfect model. Their induced decreases in the elastic moduli and strengths are found to be much lower than those of two dimensional (2D) foams and Kagome grids. For the defect in the form of separation of fiber joints, both analytical and numerical models are developed. The predicted linear decreases in the moduli and strengths (except for the compressive strength) with increasing number of separated fiber joints indicate that MFSSs be insensitive to the defect of joint separation. To explore the effect of crack, fracture toughness of MFSSs is measured and is found to be significantly higher than that of metal foams of the same relative density (i.e., volume fraction of the constituent solid material). The underlying ductile mechanism of MFSSs is further investigated by numerical simulations, showing that plastic deformation spreads all over the fibers in ligament rather than concentrates around crack tip. This study shows that MFSSs are superior in view of their resistance to the considered micro-defects and crack.     

Free convection heat transfer from an isothermal plate with arbitrary inclination

This paper presents a new formulation for the laminar free convection from an arbitrarily inclined isothermal plate to fluids of any Prandtl number between 0.001 and infinity. A novel inclination parameter is proposed such that all cases of the horizontal, inclined and vertical plates can be described by a single set of transformed equations. Moreover, the self-similar equations for the limiting cases of the horizontal and vertical plates are recovered from the transformed equations by setting=0 and=1, respectively. Heated upward-facing plates with positive and negative inclination angles are investigated. A very accurate correlation equation of the local Nusselt number is developed for arbitrary inclination angle and for 0.001 Pr .

---

Wärmeübertragung bei freier Konvektion an einer isothermen Platte mit beliebiger Neigung

Zusammenfasssung Diese Untersuchung stellt eine neue Formulierung der laminaren freien Konvektion von Flüssigkeiten mit einer Prandtl-Zahl zwischen 0,001 und unendlich an einer beliebig schräggestellten isothermen Platte dar. Ein neuer Neigungsparameter wird eingeführt, so daß alle Fälle der horizontalen, geneigten oder vertikalen Platte von einem einzigen Satz transformierter Gleichungen beschrieben werden können. Die unabhängigen Gleichungen für die beiden Fälle der horizontalen and vertikalen Platte wurden für=0 und=1 aus den transformierten Gleichungen wieder abgeleitet. Es wurden erwärmte aufwärtsgerichtete Platten mit positiven und negativen Neigungswinkeln untersucht. Eine sehr genaue Gleichung wurde für die lokale Nusselt-Zahl bei beliebigen Neigungswinkeln und für 0,001 Pr entwickelt.

---

Nomenclature C _p specific heat - f reduced stream function - g gravitational acceleration - Gr local Grashof number,g(T _w –T _w ) x³/v² - h local heat transfer coefficient - k thermal conductivity - n constant exponent - Nu local Nusselt number,hx/k - p pressure - Pr Prandtl number, v/ - Ra local Rayleigh number,g(T _w –T )J x³/v - T fluid temperature - T _w wall temperature - T temperature of ambient fluid - u velocity component in x-direction - v velocity component in y-direction - x coordinate parallel to the plate - y coordinate normal to the plate Greek symbols thermal diffusivity - thermal expansion coefficient - (Ra¦sin¦)^1/4/( Ra cos()^1/5 - pseudo-similarity variable, (y/) - dimensionless temperature, (T–T )/(T _w–T ) - ( Ra cos)^1/5+(Rasin)^1/4 - v kinematic viscosity - 1/[1 +(Ra cos)^1/5/( Ra¦sin)^1/4] - density of fluid - Pr/(1+Pr) - _w wall shear stress - angle of plate inclination measured from the horizontal - stream function - dimensionless dynamic pressure     

Free convection from a vertical plate with concentrated and distributed thermal sources

An analysis is developed for the laminar free convection from a vertical plate with uniformly distributed wall heat flux and a concentrated line thermal source embedded at the leading edge. We introduce a parameter=(1 +Q _L/Q_w)^–1=(1 + Ra_L/Ra_w)^–1 to describe the relative strength of the two thermal sources; and propose a unified buoyancy parameter=( Ra_L+ Ra_w)^1/5 with=1/(1 +Pr ^–1) to properly scale the dependent and independent variables. The variables are so defined that the resulting nonsimilar boundary-layer equations can describe exactly the buoyancy-induced flow from the dual sources with any relative strength to fluids of any Prandtl number from very small values to infinity. These nonsimilar equations are readily reducible to the self-similar equations of an adiabatic wall plume for=0, and to those of free convection from uniform flux plate for=1. Rigorous finite-difference solutions for fluids of Pr from 0.001 to are obtained over the entire range of from 0 to 1. The effects of both relative source strength and Prandtl number on the velocity profiles, temperature profiles, and the variations of wall temperature, are clearly illustrated.

---

Freie Konvektion an einer vertikalen Platte mit einer konzentrierten und einer gleichmäßig verteilten Wärmequelle

Zusammenfassung Für die freie Konvektion an einer vertikalen Platte mit einer gleichmäßig verteilten Wandwärmestromdichte und einer in der Vorderkante eingebetteten linienförmigen Wärmequelle wird eine Berechnungsmethode entwickelt. Zur Beschreibung der relativen Stärke der beiden Wärmequellen führen wir einen Parameter=(1 + Q_L/Q_w)^–1=(1 + Ra_L/Ra_w)^–1 ein und schlagen einen vereinheitlichten Auftriebsparameter=( Ra _L+ Ra _w)^1/5 mit=1/(1 +Pr ^–1 für die Skalierung der abhängigen und unabhängigen Variablen vor. Die Variablen werden so definiert, daß mit den sich ergebenden unabhängigen Grenzschichtgleichungen die von den beiden Wärmequellen beliebiger Stärke verursachte Auftriebsströmung von Fluiden beliebiger Prandtl-Zahl genau beschrieben werden kann. Diese unabhängigen Gleichungen können ohne weiteres auf die selbstähnlichen Gleichungen für den Fall einer lokalen Wärmezufuhr an einer sonst adiabatischen Wand für=0 und jenen der freien konvektion an einer Platte mit einheitlichem Wärmestrom für=1 zurückgeführt werden. Für Fluide mit der Prandtl-Zahl zwischen 0,001 und Unendlich werden nach der strengen finite Differenzen-Methode Lösungen im Bereich von zwischen 0 und 1 erhalten. Der jeweilige Einfluß der relativen Quellenstärke und der Prandtl-Zahl auf die Geschwindigkeits- und Temperaturprofile sowie die Veränderung der Wandtemperatur werden deutlich dargestellt.

---

Nomenclature C _f friction coefficient - C _p specific heat - f reduced stream function - g gravitational acceleration - k thermal conductivity - L width of the plate - Nu local Nusselt number - Pr Prandtl number - q _w wall heat flux - Q _L heat generated by the line source - Q _w heat released by the uniform-flux wall from 0 tox, q _w Lx - Ra _L local Rayleigh number, g T _L ^* x ³/( ) - Ra _w local Rayleigh number,g T _w ^* w ³/( ) - T fluid temperature - T temperature of ambient fluid - T _L ^* characteristic temperature of the line source,Q _{L/(C p} L) - T _w ^* characteristic temperature of the uniform flux wall, =q _w x/k=Q _w /(C _p L) - u velocity component in then-direction - U₀ dimensionless velocity,u/(/x) Ra _L ^2/5 - U ₁ dimensionless velocity,u/(/x) Ra _w ^2/5 - velocity component in they-direction - x coordinate parallel to the plate - y coordinate normal to the plate - thermal diffusivity - thermal expansion coefficient - pseudo-similarity variable,(y/x) - dimensionless temperature, (T–T )/(T _L ^* +T _w ^* ) - ₀ dimensionless temperature, (Ra_l)^1/5 (T–T )/T _L ^* - ₁ dimensionless temperature, (Ra_w)Ra_w)^1/5 (T–T )/T _w ^* - (Ra _L+Ra_w)^1/5 - kinematic viscosity - (1 +Ra _L/Ra_w)^–1=(1 +T _L ^* /T _w ^* )^–1=(1 + Q_L/Q_w)^–1 - density - Pr/(1 +Pr) - _w wall shear stress - stream function     

Non-elliptic deformation field near the tip of a mixed-mode crack in a compressible hyperelastic material

This paper gives an asymptotic analysis of the deformation field near the tip of an arbitrary mixed-mode crack in a compressible hyperelastic harmonic material which loses ellipticity at sufficiently large deformations. It is found that the near-tip deformation field is characterized by a localized non-elliptic deformation band issuing from the crack-tip and bounded by two curves of discontinuous deformation gradient. Explicit expression for the near-tip deformation field is obtained both inside and outside the localized deformation band. In particular, a simple relation is derived that determines the orientation of the deformation band in terms of two complex governing parameters of the near-tip fields inside and outside the deformation band, respectively.