An icing physics study by using lifetime-based molecular tagging thermometry technique

An experimental investigation was conducted to quantify the unsteady heat transfer and phase changing process within small icing water droplets in order to elucidate underlying physics to improve our understanding of the important micro-physical process of icing phenomena. A novel, lifetime-based molecular tagging thermometry (MTT) technique was developed and implemented to achieve temporally-and-spatially resolved temperature distribution measurements to reveal the time evolution of the unsteady heat transfer and dynamic phase changing process within micro-sized water droplets in the course of icing process. It was found that, after a water droplet impinged onto a frozen cold surface, the liquid water at the bottom of the droplet would be frozen and turned to solid ice rapidly, while the upper portion of the droplet was still in liquid state. As the time goes by, the interface between the liquid phase water and solid phase ice was found to move upward continuously with more and more liquid water within the droplet turned to solid ice. Interestingly, the averaged temperature of the remaining liquid water within the small icing droplet was found to increase, rather than decrease, continuously in the course of icing process. The temperature increase of the remaining liquid water is believed to be due to the heat release of the latent heat during solidification process. The volume expansion of the water droplet during the icing process was found to be mainly upward to cause droplet height growth rather than radial to enlarge the contact area of the droplet on the test plate. As a result, the spherical-cap-shaped water droplet was found to turn to a prolate-spheroid-shaped ice crystal with cusp-like top at the end of the icing process. The required freezing time for the water droplets to turn to ice crystals completely was found to depend on the surface temperature of the test plate strongly, which would decrease exponentially as the surface temperature of the frozen cold test plate decreases.     

An Integral Method to Calculate Water Saturation in a Centrifuge Experiment to Determine Capillary Pressure

The centrifuge method is commonly used to determine the capillary pressure of a porous medium, and the original approximating method for data analysis developed by Hassler and Brunner is still being used. Its limitations are, however, not well understood. Application to analyze experiments where one of the assumptions was obviously violated had been given in the literature. While the result appeared to be quite reasonable, it was not clear how close was it to reality. One of the objectives of this paper is to review the assumptions that is required to develop this method, so that the experimental condition in which it is applicable can be established. The other objective is to derive a completely different solution technique to this problem. There is no need to assume that the ratio of the inlet radius to the exit radius of the core to the center of the centrifuge be close to 1. With the freedom from this limitation it is, therefore, possible to construct machines at lower cost and to improve on the data quality by allowing longer cores.     

The effect of variations in the creep exponent on the buckling of circular cylindrical shells

The senior author solved the problem of axially symmetrical creep buckling of thin circular cylindrical shells subjected to uniform axial compression. In that analysis the constitutive equation was a power law, and the exponent was taken to be equal to three. The purpose of this work was to extend the solution to a range of values of the creep exponent, n. To cope with the increasing algebraic complexity, a digital computer was employed in two ways: to generate the set of equations symbolically, and then to solve these equations. The computer programs were used to generate numerical solutions for the cases in which n was equal to 3, 5, 7 and 9. Two simple extrapolation techniques were then employed to obtain approximate solutions to the critical time problem for values of n up to 29.     

Sloshing and energy dissipation in an egg: SPH simulations and experiments

Tall structures, such as towers and bridges, can oscillate at excessive magnitudes when subjected to wind and earthquake loads. Liquid sloshing absorbers can be used to suppress these excessive oscillations by tuning the frequency of the sloshing to the critical frequency of the structure. Sloshing absorbers are simple structures consisting of a partially full container of liquid with a free surface. Tuning ensures that significant amounts of harmful energy can be extracted from the structure to the sloshing liquid. However, there needs to be a rapid means of dissipating this energy to avoid its returning back to the structure (then back to the liquid periodically).A hen׳s egg seems to have evolved to efficiently dissipate energy to protect its embryo using sloshing of its liquid content. Hence, the potential to implement the egg׳s unique properties as a sloshing absorber for structural control, is the main focus of this study. Numerical simulations, using Smoothed Particle Hydrodynamics (SPH), and experimental comparisons are presented in this paper. One objective is to demonstrate the ability of SPH to simulate complex free surface behaviour in three dimensions. Such a tool is then useful to identify different dissipation modes. Effects of fill volume and viscosity on the rate of dissipation, are also investigated.     

Fast prediction of the fatigue behavior of short-fiber-reinforced thermoplastics based on heat build-up measurements: application to heterogeneous cases

Short-fiber-reinforced thermoplastics components for structural applications are usually very complex parts as stiffeners, ribs and thickness variations are used to compensate the quite low material intrinsic stiffness. These complex geometries induce complex local mechanical fields but also complex microstructures due to the injection process. Accounting for these two aspects is crucial for the design in regard to fatigue of these parts, especially for automotive industry. The aim of this paper is to challenge an energetic approach, defined to evaluate quickly the fatigue lifetime, on three different heterogeneous cases: a classic dog-bone sample with a skin-core microstructure and two structural samples representative of the thickness variations observed for industrial components. First, a method to evaluate dissipated energy fields from thermal measurements is described and is applied to the three samples in order to relate the cyclic loading amplitude to the fields of cyclic dissipated energy. Then, a local analysis is detailed in order to link the energy dissipated at the failure location to the fatigue lifetime and to predict the fatigue curve from the thermomechanical response of one single sample. The predictions obtained for the three cases are compared successfully to the Wöhler curves obtained with classic fatigue tests. Finally, a discussion is proposed to compare results for the three samples in terms of dissipation fields and fatigue lifetime. This comparison illustrates that, if the approach is leading to a very relevant diagnosis on each case, the dissipated energy field is not giving a straightforward access to the lifetime cartography as the relation between fatigue failure and dissipated energy seems to be dependent on the local mechanical and microstructural state.     

Determination of crack-speed history and tip locations for cracks moving with nonuniform velocity

This paper describes the use of Rayleigh waves as observed with the reflected-light method of caustics to obtain information pertaining to crack speed and tip location. The basic idea is to use the Doppler-shift effect as applied to elastic-wave emission from a running crack to geometrically reconstruct crack-tip locations. In addition to the usual information pertaining to the stress-intensity factor at the time when the photograph was taken, the process yields a complete history of crack-propagation velocities even along curved-crack paths.     

Modelling of the dyeing process of packed cotton threads using reactive dyes

The Method of Volume Averaging is used to model the process of dyeing textile threads on bobbins. This analysis allows one to upscale the relevant information at the micro-scale, composed of the textile fibres of the thread in contact with the dyeing bath fluid, to the macro-scale, consisting of the bobbins of threads inside the equipment. The final mathematical model consists of two equations, one for the fluid phase external to the thread and the other for the fluid phase internal to the thread. In order to solve the partial differential equations obtained in the mathematical model, the authors developed a computation code using the Method of Finite Volumes. This code utilized a system of generalized coordinates to facilitate application of the boundary conditions to different bobbin geometries. The numerical results for the kinetics of dyeing packed cotton threads with reactive dyes are compared to the experimental results obtained in Brazilian textile industries, leading to good agreement between theory and experiment. This demonstrates that the model developed in this paper is able to predict the operational conditions to be used in the textile industries, minimizing the consumption of dyes and other auxiliaries necessary for the dyeing process.     

Numerical solutions of Navier-Stokes equations with an integrated compartment method (ICM)

The most common numerical methods that are used by physical scientists to approximate partial differential equations employ finite differences and/or finite elements. In addition, compartment analyses have been adopted by ecological system analysts to simulate the evolution of processes governed by differential equations without spatial derivatives. An integrated compartment method (ICM) is proposed to combine the merits of these three numerical techniques. The basic procedures of the ICM are first to discretize the region of interest into compartments, then to apply three integral theorems of vectors to transform the volume integral to the surface integral, and finally to use interpolation to relate the interfacial values in terms of compartment values to close the system. These procedures are applied to the Navier-Stokes equations to yield the computational algorithm from which computer programs can be coded. The computer code is designed to solve one-, two-, or three-dimensional problems as desired. The program is applied to two simple cases: wake formation behind an obstacle in a channel and circulatory motion of a body of fluid in the square cavity. These preliminary applications have shown promising results.     

Interpolation-based reduced-order modelling for steady transonic flows via manifold learning

This paper presents a parametric reduced-order model (ROM) based on manifold learning (ML) for use in steady transonic aerodynamic applications. The main objective of this work is to derive an efficient ROM that exploits the low-dimensional nonlinear solution manifold to ensure an improved treatment of the nonlinearities involved in varying the inflow conditions to obtain an accurate prediction of shocks. The reduced-order representation of the data is derived using the Isomap ML method, which is applied to a set of sampled computational fluid dynamics (CFD) data. In order to develop a ROM that has the ability to predict approximate CFD solutions at untried parameter combinations, Isomap is coupled with an interpolation method to capture the variations in parameters like the angle of attack or the Mach number. Furthermore, an approximate local inverse mapping from the reduced-order representation to the full CFD solution space is introduced. The proposed ROM, called Isomap+I, is applied to the two-dimensional NACA 64A010 airfoil and to the 3D LANN wing. The results are compared to those obtained by proper orthogonal decomposition plus interpolation (POD+I) and to the full-order CFD model.     

生物力学与基因-献给周培源教授诞辰100周年

生物界包罗万象,其中有力的作用,所以有生物力学.自Galileo,Harvey, Boreli, Hooke, Euler,Young等创始以来,生物力学阐明了鸟飞鱼游,人体运动,血液循环,人工脏器等,对人世社会,有所贡献.生物力学的基础是质点力学,传统地用连续体力学的概念来简化.但近年做生物组织在应力的作用下改造的问题,引起了必须更改传统连续体力学的几个公理的问题.我们将仔细讨论这些公理,然后指出新公理存在的理由,是由于基因在细胞里的日常工作.基因不单主宰遗传,变异;并且忙着控制日常生活.不过,现在仅见其端倪.详细的情形,要等将来来阐发了.      

New accurate efficient modeling techniques for the vibration analysis of T-joint thin-walled box structures

We propose new accurate efficient modeling techniques for the vibration analysis of T-joint thin-walled box structures. The essence of the present techniques is to use beam elements to model thin-walled members of the joint, but the elements are based on an eight-degree-of-freedom (8-DOF) beam theory capable of handling warping and distortion. Two approaches are considered to model the interfacing joint region connected to three adjacent thin-walled box structures: the first one is to model the joint region with plate elements and the second one is to use a joint element derived to be consistent with nearby 8-DOF beam elements. The efficiency of the present techniques comes from the use of beam elements to model the box structures while the accuracy comes from the use of the higher-order beam theory accounting for warping and distortion. The procedures to match the dissimilar elements and to develop the joint element are also presented in this work. The effectiveness of the present approaches is demonstrated by numerical examples.     

Monte Carlo Simulation and Well Testing Applied in Evaluating Reservoir Properties in a Deforming Longwall Overburden

During longwall mining, the intact strata start to deform and fracture as the mining face progresses. Gob gas ventholes (GGVs) are drilled from the surface over a longwall panel before mining to capture methane from the fractured zone. Due to fracturing and bedding-plane separations, reservoir properties change extensively. This poses a major problem for venthole designers and methane control engineers and may become a safety and health concern for underground work force due to unexpected methane emissions: it is difficult to predict the location of major strata separations and their temporal magnitudes to best locate the ventholes. Measurements obtained at different times during longwall mining may not be helpful for this purpose as strata deformation is a dynamic process and the results from different tests may not be lumped together to analyze the data collectively. This article uses a combination of Monte Carlo (MC) simulation and well testing methods to analyze multiple data sets obtained from a GGV at different longwall face locations. The aim was to determine the magnitude of average strata separation before conducting well test analyses to determine the properties of a deformed reservoir. MC simulation was used to process cross-correlated and uncertain data distributions obtained from measurements to generate a set of normally distributed values for each data type. These values were further used to project the amount of strata separation to the timing of well test. Finally, well-test analyses were used to interpret test data and to evaluate reservoir properties.     

基于振动监测的静载荷滑动轴承接触摩擦故障诊断实验研究

采用静载荷滑动轴承试验台模拟轴颈-轴承从液体润滑状态逐步向干摩擦状态过渡时的接触摩擦故障的发生和发展过程,进而利用振动信号诊断滑动轴承的接触摩擦故障;通过对特征参数进行归一化处理得到无量纲特征参数,进而通过无量纲特征参数的适当数学组合得到无量纲诊断准则,并探讨了无量纲诊断准则的适用性.结果表明:利用所得到的无量纲诊断准则进行轴颈-轴承接触摩擦故障诊断时无须建立用于模式识别的标准模,可实现快捷方便的故障诊断;无量纲诊断准则对工况不敏感,而对故障更敏感,适用于可变工况下滑动轴承的故障诊断,且诊断成本较低.     

Systematic description of imperfect bifurcation behavior of symmetric systems

A systematic method is presented for describing experimental curves of force vs strain of a system with regular polygonal (dihedral group) symmetry subject to bifurcation behavior, with an aim toward overcoming the following problems : (1) it is difficult to judge whether the system is undergoing bifurcation or not ; (2) the perfect behavior of the system cannot be known due to the presence of initial imperfections ; (3) those curves are often qualitatively different from bifurcation diagrams predicted by mathematics. The tools employed are : the asymptotic theory for imperfect bifurcation, such as the Koiter law, and the stochastic theory of initial imperfections. The former theory is extended in this paper to the system with regular-polygonal symmetry to present asymptotic laws for recovering perfect curves with reference to the experimental ones. These laws are formulated for physically observable displacements, instead of the variables in the mathematical bifurcation diagrams, in order to make them readily applicable to the experimental curves. The stochastic theory is combined with an asymptotic law to develop a means to identify the multiplicity of the bifurcation point. The systematic method for describing the experimental curves developed in this manner is applied to the bifurcation analysis of regular-polygonal truss domes to testify its validity. Furthermore, this method is applied to the shear behavior of cylindrical sand specimens to show that they, in fact, are undergoing bifurcation, and, in turn, to demonstrate the importance of a viewpoint of bifurcation in the study of shear behavior of materials. The need of a dual viewpoint of bifurcation and plasticity in the study of constitutive relationship of materials is emphasized to conclude the paper.     

Dynamics of oblique detonations in ram accelerators

Time-accurate numerical simulations are used to study the dynamic development of oblique detonations on accelerating projectiles in ram accelerators. These simulations show that the oblique detonation can be stabilized on the projectile. The high pressure generated behind the detonation can result in accelerations up to 10⁶G and propel the projectile to velocities higher than 4.0 km/s. The detonation structure on the projectile is sensitive to the projectile geometry. A small change in the projectile shape is sufficient to alter the overall detonation structure and significantly affect the pressure distribution on the projectile. In order to maximize the thrust, an appropriate projectile shape has to be chosen to generate the detonation structure just behind the widest part of the projectile body. The projectile acceleration also has strong effects on the flow field and the detonation structure. During the acceleration, the location of the oblique detonation moves upstream from one reflected shock to another. However, one the detonation is stabilized behind the upstream shock, it remains at the new location until the transition to the next upstream shock occurs. In the simulations, the Non-Inertial-Source (NIS) technique was used to accurately represent of the projectile acceleration. Also, the Virtual-Cell-Embedding (VCE) method was employed to efficiently treat the complex projectile geometry on cartesian grids.     

Electro-osmotic pore pressures in soil due to an alternating electrical field

Pore pressure development in a soil specimen due to electro-osmosis under alternating current conditions is examined theoretically. Solutions to the governing equation are derived for one-dimensional flow with boundary conditions corresponding to an impervious (conventional no-flow boundary), a partially drained boundary, and a partially drained boundary with an intervening permeable zone between the boundary and the soil. These latter two boundary conditions can arise from details of pore pressure measuring systems at the specimen boundaries during laboratory experiments. An analysis of the solutions indicates that for a perfect no-flow boundary, excess pore pressures measured at an electrode consist of a steady state and rapidly-decaying transient response. The pore pressures exhibit a 45 degree phase shift relative to the applied electric current. The effect of the partially drained boundary is to reduce the peak to peak amplitude of the pore pressure and to increase the phase shift to as much as 90 degrees depending on the compressibility of the pore pressure measuring system. The effect of the impeded and partially drained boundary is to further reduce the amplitude of the pore pressures and to increase the phase shift to as much as 180 degrees depending on the relative permeability of the impeded boundary.     

Forecasting bifurcation morphing: application to cantilever-based sensing

Two novel techniques are proposed to enhance the bifurcation morphing method as applied to cantilever-based sensors. First, nonlinear feedback excitations with added time delay are employed to minimize the sensitivity of the sensors to small variations in the unavoidable time delay. Second, a novel approach to forecast bifurcations is applied to the sensors. This approach significantly reduces the time required to obtain bifurcation diagrams. Both techniques are demonstrated experimentally in detecting mass variations of a test cantilever beam. This cantilever-based sensor operating based on the bifurcation morphing method is shown to be accurate, quick and robust when these techniques are utilized.     

Nonlinear and viscous effects on wave propagation in an elastic axisymmetric vessel

In this paper, a power series and Fourier series approach is used to solve the governing equations of motion in an elastic axisymmetric vessel with the assumption that the fluid is incompressible and Newtonian in a laminar flow. We obtain solutions for the wave speed and attenuation coefficient, analytically where possible, and show how these differ under a number of different conditions. Viscosity is found to reduce the wave speed from that predicted by linear wave theory and the nonlinear terms to increase the wave speed in comparison to the linear solution. For vessels with a wall stiffness in the arterial range, the reduction in the wave speed due to the viscous terms is approximately 10% and the increase due to the nonlinear terms is approximately 5%. This difference between the linear and nonlinear wave speeds was found to be largely constant irrespective of the number of terms considered in the power series for the velocity profile. The linear wave speed was found to vary weakly with stiffness, whilst the nonlinear wave speed was found to vary significantly with the stiffness, especially at low values of stiffness. The 10% variation in the wave speed due to the viscous terms was found to be constant with wall stiffness whilst the 5% variation due to the nonlinear terms was found to vary with wall stiffness. The importance of the number of terms considered in the power series is discussed showing that only a relatively small number is required in the viscous case to obtain accurate results.     

A minimal model for studying properties of the mode-coupling type instability in friction induced oscillations

A minimal two degree of freedom model is used to clarify from an intuitive perspective the physical mechanisms underlying the mode-coupling instability of self-excited friction induced oscillations. It is shown that simultaneous out-of-phase oscillations of friction force and displacement tangential to the friction force may lead to energy transfer from the frictional system to vibrational energy. Also it is shown that the friction force acts like a cross-coupling force linking motion normal to the contact surface to motion parallel to it and that a necessary condition for the onset of instability is that these friction-induced cross-coupling forces balance the corresponding structural cross-coupling forces of the system. Finally the origin and the role of phase shifts between oscillations normal and parallel to the contact surface is clarified with respect to the mode-coupling instability. It may be expected that the intuitive picture gained will be of considerable help for practical design purposes.     

SMA-induced snap-through of unsymmetric fiber-reinforced composite laminates

A theory is developed and experiments designed to study the concept of using shape memory alloy (SMA) wires to effect the snap-through of unsymmetric composite laminates. The concept is presented in the context of structural morphing, that is, a structure changing shape to adjust to changing conditions or to change operating characteristics. While the specific problem studied is a simplification, the overall concept is to potentially take advantage of structures which have multiple equilibrium configurations and expend power only to change the structure from one configuration to another rather than to continuously expend power to hold the structure in the changed configuration. The unsymmetric laminate could be the structure itself, or simply part of a structure. Specifically, a theory is presented which allows for the prediction of the moment levels needed to effect the snap-through event. The moment is generated by a force and support arrangement attached to the laminate. A heated SMA wire attached to the supports provides the force. The necessary SMA constitutive behavior and laminate mechanics are presented. To avoid dealing with the heat transfer aspects of the SMA wire, the theory is used to predict snap-through as a function of SMA wire temperature, which can be measured directly. The geometry and force level considerations of the experiment are discussed, and the results of testing four unsymmetric laminates are compared with predictions. Laminate strain levels vs. temperature and the snap-through temperatures are measured for the these laminates. Repeatability of the experimental results is generally good, and the predictions are in reasonable agreement with the measurements.