A numerical method for separation of stresses in photo-orthotropic elasticity

A numerical method is suggested for separation of stresses in photo-orthotropic elasticity using the numerical solution of compatibility equation for orthotropic case. The compatibility equation is written in terms of a stress parameter S analogous to the sum of principal stresses in two-dimensional isotropic case. The solution of this equation provides a relation between the normal stresses. The photoelastic data give the shear stress and another relation between the two normal stresses. The accuracy of the numerical method and its application to practical problems are illustrated with examples.

An absolute digital positioning numerical control system

A module containing all the functional components required for a digital absolute positioning process of one axis of a machine tool has been designed and constructed. Circuit realization makes use of integrated circuit elements.

Numerical study of heat transfer in laminar film boiling by the finite-difference method

The finite-difference form of the basic conservation equations in laminar film boiling have been solved by the false-transient method. By a judicious choice of the coordinate system the vapour-liquid interface is fitted to the grid system. Central differencing is used for diffusion terms, upwind differencing for convection terms, and explicit differencing for transient terms. Since an explicit method is used the time step used in the false-transient method is constrained by numerical instability. In the present problem the limits on the time step are imposed by conditions in the vapour region. On the other hand the rate of convergence of finite-difference equations is dependent on the conditions in the liquid region. The rate of convergence was accelerated by using the over-relaxation technique in the liquid region. The results obtained compare well with previous work and experimental data available in the literature.

Vibrations of Multifiber Composite Shells-Some Numerical Results

A study of vibrations of multifiber composite shells is presented. Special attention is paid to the effect of composition of different fibers on the frequency spectrum of a freely vibrating cylindrical shell. The numerical results indicate clustering of frequency spectrum of a freely vibrating cylindrical composite shell as compared with the isotropic shell, and the spectrum varies considerably with the composition of the constituent materials.

Numerical integration over arbitrary polygonal domains based on Schwarz-Christoffel conformal mapping

This paper presents a new numerical integration technique oil arbitrary polygonal domains. The polygonal domain is mapped conformally to the unit disk using Schwarz-Christoffel mapping and a midpoint quadrature rule defined oil this unit disk is used. This method eliminates the need for a two-level isoparametric mapping Usually required. Moreover, the positivity of the Jacobian is guaranteed. Numerical results presented for a few benchmark problems in the context of polygonal finite elements show that the proposed method yields accurate results.

Stationary crack tip fields in elastic-plastic solids: an overview of recent numerical simulations

In this paper, an overview of some recent numerical simulations of stationary crack tip fields in elastic-plastic solids is presented. First, asymptotic analyses carried out within the framework of 2D plane strain or plane stress conditions in both pressure insensitive and pressure sensitive plastic solids are reviewed. This is followed by discussion of salient results obtained from recent computational studies. These pertain to 3D characteristics of elastic-plastic near-front fields under mixed mode loading, mechanics of fracture and simulation of near-tip shear banding process of amorphous alloys and influence of crack tip constraint on the structure of near-tip fields in ductile single crystals. These results serve to illustrate several important features associated with stress and strain distributions near the crack tip and provide the foundation for understanding the operative failure mechanisms. The paper concludes by highlighting some of the future prospects for this field of study.

Singular value decomposition based computationally efficient algorithm for rapid dynamic near-infrared diffuse optical tomography

Purpose: A computationally efficient algorithm (linear iterative type) based on singular value decomposition (SVD) of the Jacobian has been developed that can be used in rapid dynamic near-infrared (NIR) diffuse optical tomography. Methods: Numerical and experimental studies have been conducted to prove the computational efficacy of this SVD-based algorithm over conventional optical image reconstruction algorithms. Results: These studies indicate that the performance of linear iterative algorithms in terms of contrast recovery (quantitation of optical images) is better compared to nonlinear iterative (conventional) algorithms, provided the initial guess is close to the actual solution. The nonlinear algorithms can provide better quality images compared to the linear iterative type algorithms. Moreover, the analytical and numerical equivalence of the SVD-based algorithm to linear iterative algorithms was also established as a part of this work. It is also demonstrated that the SVD-based image reconstruction typically requires O(NN2) operations per iteration, as contrasted with linear and nonlinear iterative methods that, respectively, requir O(NN3) and O(NN6) operations, with ``NN'' being the number of unknown parameters in the optical image reconstruction procedure. Conclusions: This SVD-based computationally efficient algorithm can make the integration of image reconstruction procedure with the data acquisition feasible, in turn making the rapid dynamic NIR tomography viable in the clinic to continuously monitor hemodynamic changes in the tissue pathophysiology.

Numerical Study of Heat Transfer From Pin-Fin Heat Sink Using Steady and Pulsated Impinging Jets

This paper investigates numerically the heat transfer characteristics of confined slot jet impingement on a pin-fin heat sink. A variety of pin-fin heat sinks is investigated, and the resulting enhancement of heat transfer studied. The distribution of heat transfer coefficient on the top surface of the base plate and that along the fin height are examined. Both steady and pulsated jets are studied. It is observed that for a steady jet impingement on a pin-fin heat sink, the effective heat transfer coefficient increases with fin height, leading to a corresponding decrease in base plate temperature for the same heat flux. In the case of pulsated jets, the influence of pulse frequency and the Reynolds number is examined, and their effect on the effective heat transfer coefficient is studied.

A numerical study of projectile impact on thin aluminium plates

This article deals with a simulation-based Study of the impact of projectiles on thin aluminium plates using LS-DYNA by modelling plates with shell elements and projectiles with solid elements. In order to establish the required modelling criterion in terms of element size for aluminium plates, a convergence Study of residual velocity has been carried Out by varying mesh density in the impact zone. Using the preferred material and meshing criteria arrived at here, extremely good prediction of test residual velocities and ballistic limits given by Gupta et al. (2001) for thin aluminium plates has been obtained. The simulation-based pattern of failure with localized bulging and jagged edge of perforation is similar to the perforation with petalling seen in tests. A number Of simulation-based parametric studies have been carried out and results consistent with published test data have been obtained. Despite the robust correlation achieved against published experimental results, it would be prudent to conduct one's own experiments, for a final correlation via the present modelling procedure and analysis with the explicit LS-DYNTA 970 solver. Hence, a sophisticated ballistic impact testing facility and a high-speed camera have been used to conduct additional tests on grade 1100 aluminium plates of 1 mm thickness with projectiles Of four different nose shapes. Finally, using the developed numerical simulation procedure, an excellent correlation of residual velocity and failure modes with the corresponding test results has been obtained.

Performance Analysis of Reconstruction Techniques for Frequency-Domain Optical-Coherence Tomography

We address the issue of noise robustness of reconstruction techniques for frequency-domain optical-coherence tomography (FDOCT). We consider three reconstruction techniques: Fourier, iterative phase recovery, and cepstral techniques. We characterize the reconstructions in terms of their statistical bias and variance and obtain approximate analytical expressions under the assumption of small noise. We also perform Monte Carlo analyses and show that the experimental results are in agreement with the theoretical predictions. It turns out that the iterative and cepstral techniques yield reconstructions with a smaller bias than the Fourier method. The three techniques, however, have identical variance profiles, and their consistency increases linearly as a function of the signal-to-noise ratio.

A numerical study of crack tip constraint in ductile single crystals

In this work, the effect of crack tip constraint on near-tip stress and deformation fields in a ductile FCC single crystal is studied under mode I, plane strain conditions. To this end, modified boundary layer simulations within crystal plasticity framework are performed, neglecting elastic anisotropy. The first and second terms of the isotropic elastic crack tip field, which are governed by the stress intensity factor K and T-stress, are prescribed as remote boundary conditions and solutions pertaining to different levels of T-stress are generated. It is found that the near-tip deformation field, especially, the development of kink or slip shear bands, is sensitive to the constraint level. The stress distribution and the size and shape of the plastic zone near the crack tip are also strongly influenced by the level of T-stress, with progressive loss of crack tip constraint occurring as T-stress becomes more negative. A family of near-tip fields is obtained which are characterized by two terms (such as K and T or J and a constraint parameter Q) as in isotropic plastic solids.

A numerical study of the effect of loading conditions on tubular hydroforming

In this work, the mechanics of tubular hydroforming under various types of loading conditions is investigated. The main objective is to contrast the effects of prescribing fluid pressure or volume flow rate, in conjunction with axial displacement, on the stress and strain histories experienced by the tube and the process of bulging. To this end, axisymmetric finite element simulations of free hydroforming (without external die contact) of aluminium alloy tubes are carried out. Hill’s normally anisotropic yield theory along with material properties determined in a previous experimental study [A. Kulkarni, P. Biswas, R. Narasimhan, A. Luo, T. Stoughton, R. Mishra, A.K. Sachdev, An experimental and numerical study of necking initiation in aluminium alloy tubes during hydroforming, Int. J. Mech. Sci. 46 (2004) 1727–1746] are employed in the computations. It is found that while prescribed fluid pressure leads to highly non-proportional strain paths, specified fluid volume flow rate may result in almost proportional ones for the predominant portion of loading. The peak pressure increases with axial compression for the former, while the reverse trend applies under the latter. The implication of these results on failure by localized necking of the tube wall is addressed in a subsequent investigation.

Numerical Simulation of Groundwater Flow with Gradually Increasing Heterogeneity due to Clogging

Well injection replenishes depleting water levels in a well field. Observation well water levels some distance away from the injection well are the indicators of the success of a well injection program. Simulation of the observation well response, located a few tens of meters from the injection well, is likely to be affected by the effects of nonhomogeneous medium, inclined initial water table, and aquifer clogging. Existing algorithms, such as the U.S. Geological Survey groundwater flow software MODFLOW, are capable of handling the first two conditions, whereas time-dependent clogging effects are yet to be introduced in the groundwater flow models. Elsewhere, aquifer clogging is extensively researched in theory of filtration; scope for its application in a well field is a potential research problem. In the present paper, coupling of one such filtration theory to MODFLOW is introduced. Simulation of clogging effects during “Hansol” well recharge in the parts of western India is found to be encouraging.

Numerical studies on channel formation and growth during solidification: Effect of process parameters

In the present work, solidification of a hyper-eutectic ammonium chloride solution in a bottom-cooled cavity (i.e. with stable thermal gradient) is numerically studied. A Rayleigh number based criterion is developed, which determines the conditions favorable for freckles formation. This criterion, when expressed in terms of physical properties and process parameters, yields the condition for plume formation as a function of concentration, liquid fraction, permeability, growth rate of a mushy layer and thermophysical properties. Subsequently, numerical simulations are performed for cases with initial and boundary conditions favoring freckle formation. The effects of parameters, such as cooling rate and initial concentration, on the formation and growth of freckles are investigated. It was found that a high cooling rate produced larger and more defined channels which are retained for a longer durations. Similarly, a lower initial concentration of solute resulted in fewer but more pronounced channels. The number and size of channels are also found to be related to the mushy zone thickness. The trends predicted with regard to the variation of number of channels with time under different process conditions are in accordance with the experimental observations reported in the literature.

Single-resolution and multiresolution extended-Kalman-filter-based reconstruction approaches to optical refraction tomography

The problem of reconstruction of a refractive-index distribution (RID) in optical refraction tomography (ORT) with optical path-length difference (OPD) data is solved using two adaptive-estimation-based extended-Kalman-filter (EKF) approaches. First, a basic single-resolution EKF (SR-EKF) is applied to a state variable model describing the tomographic process, to estimate the RID of an optically transparent refracting object from noisy OPD data. The initialization of the biases and covariances corresponding to the state and measurement noise is discussed. The state and measurement noise biases and covariances are adaptively estimated. An EKF is then applied to the wavelet-transformed state variable model to yield a wavelet-based multiresolution EKF (MR-EKF) solution approach. To numerically validate the adaptive EKF approaches, we evaluate them with benchmark studies of standard stationary cases, where comparative results with commonly used efficient deterministic approaches can be obtained. Detailed reconstruction studies for the SR-EKF and two versions of the MR-EKF (with Haar and Daubechies-4 wavelets) compare well with those obtained from a typically used variant of the (deterministic) algebraic reconstruction technique, the average correction per projection method, thus establishing the capability of the EKF for ORT. To the best of our knowledge, the present work contains unique reconstruction studies encompassing the use of EKF for ORT in single-resolution and multiresolution formulations, and also in the use of adaptive estimation of the EKF's noise covariances. (C) 2010 Optical Society of America