Phase-field elasticity model based on mechanical jump conditions

Computational models based on the phase-field method typically operate on a mesoscopic length scale and resolve structural changes of the material and furthermore provide valuable information about microstructure and mechanical property relations. An accurate calculation of the stresses and mechanical energy at the transition region is therefore indispensable. We derive a quantitative phase-field elasticity model based on force balance and Hadamard jump conditions at the interface. Comparing the simulated stress profiles calculated with Voigt/Taylor (Annalen der Physik 274(12):573, 1889), Reuss/Sachs (Z Angew Math Mech 9:49, 1929) and the proposed model with the theoretically predicted stress fields in a plate with a round inclusion under hydrostatic tension, we show the quantitative characteristics of the model. In order to validate the elastic contribution to the driving force for phase transition, we demonstrate the absence of excess energy, calculated by Durga et al. (Model Simul Mater Sci Eng 21(5):055018, 2013), in a one-dimensional equilibrium condition of serial and parallel material chains. To validate the driving force for systems with curved transition regions, we relate simulations to the Gibbs-Thompson equilibrium condition

The effect of surface tension in porous wave maker problems

Using a mixed-type Fourier transform of a general form in the case of water of infinite depth and the method of eigenfunction expansion in the case of water of finite depth, several boundary-value problems involving the propagation and scattering of time harmonic surface water waves by vertical porous walls have been fully investigated, taking into account the effect of surface tension also. Known results are recovered either directly or as particular cases of the general problems under consideration.

Finite Element Solution Of Surface-Tension Driven Flows In Laser Surface-Melting

Lasers are very efficient in heating localized regions and hence they find a wide application in surface treatment processes. The surface of a material can be selectively modified to give superior wear and corrosion resistance. In laser surface-melting and welding problems, the high temperature gradient prevailing in the free surface induces a surface-tension gradient which is the dominant driving force for convection (known as thermo-capillary or Marangoni convection). It has been reported that the surface-tension driven convection plays a dominant role in determining the melt pool shape. In most of the earlier works on laser-melting and related problems, the finite difference method (FDM) has been used to solve the Navier Stokes equations [1]. Since the Reynolds number is quite high in these cases, upwinding has been used. Though upwinding gives physically realistic solutions even on a coarse grid, the results are inaccurate. McLay and Carey have solved the thermo-capillary flow in welding problems by an implicit finite element method [2]. They used the conventional Galerkin finite element method (FEM) which requires that the pressure be interpolated by one order lower than velocity (mixed interpolation). This restricts the choice of elements to certain higher order elements which need numerical integration for evaluation of element matrices. The implicit algorithm yields a system of nonlinear, unsymmetric equations which are not positive definite. Computations would be possible only with large mainframe computers.Sluzalec [3] has modeled the pulsed laser-melting problem by an explicit method (FEM). He has used the six-node triangular element with mixed interpolation. Since he has considered the buoyancy induced flow only, the velocity values are small. In the present work, an equal order explicit FEM is used to compute the thermo-capillary flow in the laser surface-melting problem. As this method permits equal order interpolation, there is no restriction in the choice of elements. Even linear elements such as the three-node triangular elements can be used. As the governing equations are solved in a sequential manner, the computer memory requirement is less. The finite element formulation is discussed in this paper along with typical numerical results.

An absolute method for the determination of surface tension of liquids using pendent drop profiles

Drop formation at conical tips which is of relevance to metallurgists is investigated based on the principle of minimization of free energy using the variational approach. The dimensionless governing equations for drop profiles are computer solved using the fourth order Runge-Kutta method. For different cone angles, the theoretical plots of XT and ZT vs their ratio, are statistically analyzed, where XT and ZT are the dimensionless x and z coordinates of the drop profile at a plane at the conical tip, perpendicular to the axis of symmetry. Based on the mathematical description of these curves, an absolute method has been proposed for the determination of surface tension of liquids, which is shown to be preferable in comparison with the earlier pendent-drop profile methods.

Line tension of a two dimensional gas-liquid interface

In two dimensional (2D) gas-liquid systems, the reported simulation values of line tension are known to disagree with the existing theoretical estimates. We find that while the simulation erred in truncating the range of the interaction potential, and as a result grossly underestimated the actual value, the earlier theoretical calculation was also limited by several approximations. When both the simulation and the theory are improved, we find that the estimate of line tension is in better agreement with each other. The small value of surface tension suggests increased influence of noncircular clusters in 2D gas-liquid nucleation, as indeed observed in a recent simulation.

Interfacial tension model for catalytically driven nanorods

We present an analysis of the interfacial tension model for the movement of the catalytically driven nanorod. The model considers the convective reaction-diffusion equation for the production and diffusion of oxygen around the bimetallic nanorod. We solve the equation and find the concentration difference, which drives the nanorod. We use our expression to calculate the force on the nanorod and find that the result is within 20% of the results found earlier [ W. Paxton et al., J. Am. Chem. Soc. 128, 14881 (2006) ] by an approximate method. Unlike the earlier results, our results are valid from short to long lengths of the nanorod.

On latent heat of vaporization, surface tension, and temperature

A semitheoretical equation for latent heat of vaporization has been derived and tested. The average error in predicting the value at the normal boiling point in the case of about 90 compounds, which includes polar and nonpolar liquids, is about 1.8%. A relation between latent heat of vaporization and surface tension is also derived and is shown to lead to Watson's empirical relation which gives the change of latent heat of vaporization with temperature. This gives a physico-chemical justification for Watson's empirical relation and provides a rapid method of determining latent heats by measuring surface tension.

The effect of viscosity and surface tension on the Taylor instability of the surface of a cavitation bubble*

This paper is devoted to a consideration of the following problem: A spherical mass of fluid of density varrho1, viscosity μ1 and external radius R is surrounded by a fluid of density varrho2 and viscosity μ2.The fluids are immiscible and incompressible. The interface is accelerated radially by g1: to study the effect of viscosity and surface tension on the stability of the interface. By analyzing the problem in spherical harmonics the mathematical problem is reduced to one of solution of the characteristic determinant equation. The particular case of a cavity bubble, where the viscosity μ1 of the fluid inside the bubble is negligible in comparison with the viscosity μ2 of the fluid outside the bubble, is considered in some detail. It is shown that viscosity has a stabilizing role on the interface; and when g1 > T(n − 1) (n + 2)/R2(varrho2 − varrho1) the stabilizing role of both viscosity and surface tension is more pronounced than would result when either of them is taken individually.

On the investigation of singular stress field around radial cracks in annulii subjected to diametrical tension

Results of photoelastic investigation conducted on annulii containing a radial crack at inner edge and subjected to diametrical tension are reported. The cracks are oriented at 90°, 60° and 45° to the loading direction. The Stress-Intensity Factors (SIFs) were determined by analysing the crack-tip stress fields. Smith and Smith's method [Engng Fracture Mech.4, 357–366 (1972)] and a modified method developed earlier by the authors (to be published) were adopted in the evaluation of SIFs.

Stress intensity factor determination of radially cracked circular rings subjected to tension using photoelastic technique

A parametric study was carried out to determine the Stress Intensity Factor (SIF) in a cracked circular ring by using the photoelastic technique. The stress intensity factors for mode I deformation were determined by subjecting the specimens to the tensile loading from inner boundary and through the holes. The results of Non-Dimensional Stress Intensity Factor (NDSIF) variation with non-dimensional crack length for both methods of loading are compared with each other and with published results.

A lock-type cell for NMR studies at high hydrostatic pressures and low temperatures

Details of a simple and convenient high-pressure cell for continuous-wave, wide-line nuclear magnetic resonance investigation at high pressures and low temperatures are described. Experimental results obtained with the cell at 14*108 Pa and 77K for ammonium iodide are presented briefly.

Dynamics of freely suspended films with surface tension

We have studied the hydrodynamics of freely suspended membranes, liquid as well as crystalline, with surface tension. We find that nonlinear coupling to thermally excited undulations gives a singular contribution to the kinetic coefficients of these systems at low frequency and wavenumber. Our results differ in some important respects from those of Katz and Lebedev on this problem, and can be tested in mechanical impedance as well as time-correlation studies.