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.

A frame-invariant scheme for the geometrically exact beam using rotation vector parametrization

While frame-invariant solutions for arbitrarily large rotational deformations have been reported through the orthogonal matrix parametrization, derivation of such solutions purely through a rotation vector parametrization, which uses only three parameters and provides a parsimonious storage of rotations, is novel and constitutes the subject of this paper. In particular, we employ interpolations of relative rotations and a new rotation vector update for a strain-objective finite element formulation in the material framework. We show that the update provides either the desired rotation vector or its complement. This rules out an additive interpolation of total rotation vectors at the nodes. Hence, interpolations of relative rotation vectors are used. Through numerical examples, we show that combining the proposed update with interpolations of relative rotations yields frame-invariant and path-independent numerical solutions. Advantages of the present approach vis-a-vis the updated Lagrangian formulation are also analyzed.

Spatial Variability of the Depth of Weathered and Engineering Bedrock using Multichannel Analysis of Surface Wave Method

In this paper an attempt has been made to evaluate the spatial variability of the depth of weathered and engineering bedrock in Bangalore, south India using Multichannel Analysis of Surface Wave (MASW) survey. One-dimensional MASW survey has been carried out at 58 locations and shear-wave velocities are measured. Using velocity profiles, the depth of weathered rock and engineering rock surface levels has been determined. Based on the literature, shear-wave velocity of 330 ± 30 m/s for weathered rock or soft rock and 760 ± 60 m/s for engineering rock or hard rock has been considered. Depths corresponding to these velocity ranges are evaluated with respect to ground contour levels and top surface levels have been mapped with an interpolation technique using natural neighborhood. The depth of weathered rock varies from 1 m to about 21 m. In 58 testing locations, only 42 locations reached the depths which have a shear-wave velocity of more than 760 ± 60 m/s. The depth of engineering rock is evaluated from these data and it varies from 1 m to about 50 m. Further, these rock depths have been compared with a subsurface profile obtained from a two-dimensional (2-D) MASW survey at 20 locations and a few selected available bore logs from the deep geotechnical boreholes.

Study of carotid arterial plaque stress for symptomatic and asymptomatic patients

Stroke is one of the leading causes of death in the world, resulting mostly from the sudden ruptures of atherosclerosis carotid plaques. Until now, the exact plaque rupture mechanism has not been fully understood, and also the plaque rupture risk stratification. The advanced multi-spectral magnetic resonance imaging (MRI) has allowed the plaque components to be visualized in-vivo and reconstructed by computational modeling. In the study, plaque stress analysis using fully coupled fluid structure interaction was applied to 20 patients (12 symptomatic and 8 asymptomatic) reconstructed from in-vivo MRI, followed by a detailed biomechanics analysis, and morphological feature study. The locally extreme stress conditions can be found in the fibrous cap region, 85% at the plaque shoulder based on the present study cases. Local maximum stress values predicted in the plaque region were found to be significantly higher in symptomatic patients than that in asymptomatic patients (200±43. kPa vs. 127±37. kPa, p=0.001). Plaque stress level, defined by excluding 5% highest stress nodes in the fibrous cap region based on the accumulative histogram of stress experienced on the computational nodes in the fibrous cap, was also significantly higher in symptomatic patients than that in asymptomatic patients (154±32. kPa vs. 111±23. kPa, p<0.05). Although there was no significant difference in lipid core size between the two patient groups, symptomatic group normally had a larger lipid core and a significantly thinner fibrous cap based on the reconstructed plaques using 3D interpolation from stacks of 2D contours. Plaques with a higher stenosis were more likely to have extreme stress conditions upstream of plaque throat. The combined analyses of plaque MR image and plaque stress will advance our understanding of plaque rupture, and provide a useful tool on assessing plaque rupture risk.

A Doubly curved Quadrilateral finite-element for the analysis of laminated anisotropic thin shells of revoution

The details of development of the stiffness matrix for a doubly curved quadrilateral element suited for static and dynamic analysis of laminated anisotropic thin shells of revolution are reported. Expressing the assumed displacement state over the middle surface of the shell as products of one-dimensional first order Hermite polynomials, it is possible to ensure that the displacement state for the assembled set of such elements, is geometrically admissible. Monotonic convergence of total potential energy is therefore possible as the modelling is successively refined. Systematic evaluation of performance of the element is conducted, considering various examples for which analytical or other solutions are available.

Effects of tempered structure on the near threshold fatigue crack growth behavior of a coarse grained high strength steel

Near threshold fatigue crack growth behavior of a high strength steel under different temper levels was investigated. It is found that the observed variations in ΔKth could predominantly be attributed to roughness induced crack closure. The closure-free component of the threshold stress intensity range, ΔKeff,th showed a systematic variation with monotonic yield strength.

Instability of laminated composite thin-walled open-section beams

Instability of thin-walled open-section laminated composite beams is studied using the finite element method. A two-noded, 8 df per node thin-walled open-section laminated composite beam finite element has been used. The displacements of the element reference axis are expressed in terms of one-dimensional first order Hermite interpolation polynomials, and line member assumptions are invoked in formulation of the elastic stiffness matrix and geometric stiffness matrix. The nonlinear expressions for the strains occurring in thin-walled open-section beams, when subjected to axial, flexural and torsional loads, are incorporated in a general instability analysis. Several problems for which continuum solutions (exact/approximate) are possible have been solved in order to evaluate the performance of finite element. Next its applicability is demonstrated by predicting the buckling loads for the following problems of laminated composites: (i) two layer (45°/−45°) composite Z section cantilever beam and (ii) three layer (0°/45°/0°) composite Z section cantilever beam.

A simple Bayesian decision-theoretic design for dose-finding trials

A flexible and simple Bayesian decision-theoretic design for dose-finding trials is proposed in this paper. In order to reduce the computational burden, we adopt a working model with conjugate priors, which is flexible to fit all monotonic dose-toxicity curves and produces analytic posterior distributions. We also discuss how to use a proper utility function to reflect the interest of the trial. Patients are allocated based on not only the utility function but also the chosen dose selection rule. The most popular dose selection rule is the one-step-look-ahead (OSLA), which selects the best-so-far dose. A more complicated rule, such as the two-step-look-ahead, is theoretically more efficient than the OSLA only when the required distributional assumptions are met, which is, however, often not the case in practice. We carried out extensive simulation studies to evaluate these two dose selection rules and found that OSLA was often more efficient than two-step-look-ahead under the proposed Bayesian structure. Moreover, our simulation results show that the proposed Bayesian method's performance is superior to several popular Bayesian methods and that the negative impact of prior misspecification can be managed in the design stage.

A H-1-NMR STUDY OF BILIRUBIN IX-ALPHA SOLUBILIZATION BY CHOLATE MICELLES - APPLICATION OF NUCLEAR OVERHAUSER EFFECTS

The solubilization of bilirubin IX-Alpha in aqueous solution by sodium cholate micelles has been examined by 270 MHz 1H-NMR spectroscopy. Incorporation of bilirubin into the micelles is accompanied by specific shifts of bilirubin vinyl and bridgehead protons and the C18 and C19 methyl groups of the steroid. The observed chemical shifts show a monotonic concentration dependence suggesting that changes in aggregation size are continuous. Nuclear Overhauser effects (NOE) have been shown to be a useful probe or micellization. A 4:1 cholate/bilirubin mixture has been investigated by difference NOE spectroscopy. The observation of intermolecular nuclear Overhauser effects between peripheral protons of bilirubin and cholate are diagnostic of spatially proximate groups. Inter-cholate nuclear Overhauser effects increase in magnitude upon bilirubin incorporation suggesting closer packing of steroid molecules on solubilization of the pigment. Intramolecular nuclear Overhauser effects observed for solubilized bilirubin are consistent with a compact intramolecularly hydrogen-bonded conformation resembling that determined for bilirubin in the solid state.

Continuation of limit cycles near saddle homoclinic points using splines in phase space

We present a new algorithm for continuation of limit cycles of autonomous systems as a system parameter is varied. The algorithm works in phase space with an ordered set of points on the limit cycle, along with spline interpolation. Currently popular algorithms in bifurcation analysis packages compute time-domain approximations of limit cycles using either shooting or collocation. The present approach seems useful for continuation near saddle homoclinic points, where it encounters a corner while time-domain methods essentially encounter a discontinuity (a relatively short period of rapid variation). Other phase space-based algorithms use rescaled arclength in place of time, but subsequently resemble the time-domain methods. Compared to these, we introduce additional freedom through a variable stretching of arclength based on local curvature, through the use of an auxiliary index-based variable. Several numerical examples are presented. Comparisons with results from the popular package, MATCONT, are favorable close to saddle homoclinic points.

Estimation Of Properties Of Rare-Gas Solids From Compression Characteristics Of Closed-Shell Ions

We present a unified approach to repulsion in ionic and van der Waals solids based on a compressible-ion/atom model. Earlier studies have shown that repulsion in ionic crystals can be viewed as arising from the compression energy of ions, described by two parameters per ion. Here we obtain the compression parameters of the rare-gas atoms Ne. Ar. Kr and Xe by interpolation using the known parameters of related equi-electronic ions (e.g. Ar from S2-. Cl-, K- and Ca2-). These parameters fit the experimental zero-temperature interatomic distances and compressibilities of the rare-gas crystals satisfactorily. A hightemperature equation of state based on an Einstein model of thermal motions is used to calculate the thermal expansivities, compressibilities and their temperature derivatives for Ar. Kr and Xe. It is argued that an instability at higher temperatures represents the limit to which the solid can be superheated. beyond which sublimation must occur.

Finite element analysis of laminated anisotropic thin-walled open-section beams

A finite element analysis of thin-walled open-section laminated anisotropic beams is presented herein. A two-noded, 8 degrees of freedom per node thin-walled open-section laminated anisotropic beam finite element has been developed and used. The displacements of the element reference axes are expressed in terms of one-dimensional first order Hermite interpolation polynomials and line member assumptions are invoked in the formulation of the stiffness matrix. The problems of: 1. (a) an isotropic material Z section straight cantilever beam, and 2. (b) a single-layer (0°) composite Z section straight cantilever beam, for which continuum solutions (exact/approximate) are possible, have been solved in order to evaluate the performance of the finite element. Its applicability has been shown by solving the following problems: 3. (c) a two-layer (45°/−45°) composite Z section straight cantilever beam, 4. (d) a three-layer (0°/45°/0°) composite Z section straight cantilever beam.

A laminated anisotropic curved beam and shell stiffening finite element

The details of development of the stiffness matrix of a laminated anisotropic curved beam finite element are reported. It is a 16 dof element which makes use of 1-D first order Hermite interpolation polynomials for expressing it's assumed displacement state. The performance of the element is evaluated considering various examples for which analytical or other solutions are available.

Effect of Load on Dry Sliding Wear of Aluminum-Silicon Alloys

Aluminum-silicon alloy pins were slid against steel disks under nominally dry condition at a speed of 0.6 m s-1. Each pin was slid at a constant load for 5 min, the load being increased in suitable steps from 2 to 65 N. The results show the wear to increase almost monotonically with load, to be sensitive to the presence of silicon in the alloy, and to be insensitive to actual silicon content. The monotonic nature of wear rate-load characteristic suggests that one dominant wear mechanism prevails over the load range studied. Morphological studies of the pin surface and the debris support this contention and point to delamination as being the dominant mode of wear.

Stretching single DNA-molecules with temperature-stabilized optical tweezers

This thesis consists of two parts; in the first part we performed a single-molecule force extension measurement with 10kb long DNA-molecules from phage-λ to validate the calibration and single-molecule capability of our optical tweezers instrument. Fitting the worm-like chain interpolation formula to the data revealed that ca. 71% of the DNA tethers featured a contour length within ±15% of the expected value (3.38 µm). Only 25% of the found DNA had a persistence length between 30 and 60 nm. The correct value should be within 40 to 60 nm. In the second part we designed and built a precise temperature controller to remove thermal fluctuations that cause drifting of the optical trap. The controller uses feed-forward and PID (proportional-integral-derivative) feedback to achieve 1.58 mK precision and 0.3 K absolute accuracy. During a 5 min test run it reduced drifting of the trap from 1.4 nm/min in open-loop to 0.6 nm/min in closed-loop.