Active vibration control of a flexible beam subjected to arbitrary, unmeasurable disturbance forces is investigated in this paper. The concept of independent modal space control is adopted. Both the feedforward and feedback control is implemented here to reduce the beam vibration. Because of the existence of the disturbance forces, the feedforward control is applied by employing the idea of force cancellation. A modal space disturbance force observer is then established in this paper to observe the disturbance modal forces for the feedforward control. For obtaining the feedforward and feedback control gains with the optimal sense, the nearly optimal control law is derived, where the modal disturbance forces are regarded as additional states. The vibration control performances and the asymptotic properties of the control law are discussed.

The aim of this paper is to develop a CFD solver that used to simulate the transient flow phenomena in a reflected shock tunnel. The transient flow phenomena in a shock tunnel include the reflected shock/boundary layer interaction and the starting process of nozzle flow that can affect the duration of test flow in actual conditions. To numerically simulate these transient flow features, a full-size, axisymmetric reflected shock tunnel model is used. The governing equations are a full Navier-Stokes equation, a species equation and a simplified polynomial correlation to simulate the real gas effects. The numerical code is developed based on the finite volume method coupled with the upwind Roe's scheme and the total variation diminishing (TVD) method. To increase the calculation efficiency, a multi-block and multi-mesh grid generation technique is employed in a huge computational domain. The present computational results have not only confirmed the theoretical characteristics of a shock tube, but have also qualitatively presented the phenomena of reflected shock/boundary layer interaction and the starting process of nozzle flow. This numerical code is a useful tool to demonstrate the actual flow phenomena and to assist the design of experiments.

Strong pseudo-steady Mach reflections in sulfur hexafluoride (SF6) are analyzed using the three-shock and local three-shock theories, where both the vibrationally-frozen (γ = 1.333) and -equilibrated (γ = 1.093) perfect-gas models are used to compare with existing experiments. The ranges of the incident shock Mach number and reflecting wedge angle studied are 1.49 ≤ Ms ≤ 5.95 and 10° ≤ θω ≤ 42°, respectively. It is found that predicted angles between the incident and reflected shocks from the local three-shock theory using the vibrationally-equilibrated fictitious perfect-gas model (i.e., γ = 1.093) agree closely with those, currently available in literature, measured experimentally; while these predicted angles obtained using the vibrationally-frozen perfect-gas model (i.e., γ = 1.333) differ significantly from the existing experiments. Taking the convex Mach stem curvature at the triple point into consideration, it is shown that both the triple point trajectory angle and the angle between the incident and reflected shocks of strong pseudo-steady Mach reflections in SF6 can be more accurately determined for wide ranges of Ms and θω from the three-shock theory using the vibrationally-equilibrated fictitious perfect-gas model than those without considering this effect.

Flow visualization was made for water to broaden our fundamental understanding of the physical process of the confined circular jet impingement. Transition (with laminar wall jet) and turbulent regimes are identified. A correlation of boundary layer thickness was developed in terms of relevant parameters.

Direct numerical simulations of 2D turbulent flows, freely decaying as well as forced, are performed to examine the mechanism of the enstrophy cascade and serve as a template of developing LES models. The stretching effect on the 2D vorticity gradients is emphasized on the analogy of the stretching effect on 3D vorticity. The enstrophy cascade rate, the Reynolds stresses and the associated eddy viscosity for 2D turbulence are correspondingly derived and investigated. Proposed herein is that the enstrophy cascade rate to be modeled in a large-eddy simulation can be and should be calculated using the only available large-eddy information, especially when the Reynolds number is not very large or when the flow is not stationary.

The simulation results suggest all Kolmogorov's, Kraichnan's, and Saffman's similarity spectra. The Kolmogorov's spectrum appears in front of forced wave numbers and creates a subrange of a zero enstrophy cascade rate and a constant energy cascade rate. The Saffman's spectrum is the dissipation spectrum at large wave numbers. Kraichnan's spectrum shows up at intermediate wave numbers when the Reynolds number is sufficiently high. When the Smagorinsky model is employed for a large eddy simulation, its inability of capturing the significant reverse cascade phenomenon as observed in the DNS data becomes a fatal defect. Nonetheless, if only the mean cascade rate is concerned, the required Smagorinsky constant is evaluated using the DNS data and compared with the theoretical prediction of the Kraichnan's spectrum.

In this study, the motion of incompressible viscous fluid in a two-dimensional domain is solved by the finite element method using the velocity-vorticity formulation. To demonstrate the model feasibility, first of all the steady Stokes flow in a square cavity are computed. The results of square cavity flow are comparable with the numerical solutions of Burggraff (1966, FDM) [1]. Then the unsteady Navier-Stokes flow are computed and compared with other models [1∼4]. The results reveal that finite element analysis is a very powerful approach in the realm of computational fluid mechanics.

Fatigue crack propagation data of a batch of AISI 4340 steel specimens are released in the present paper. The statistical nature of the data is specially emphasized, and a probabilistic fracture mechanics model is introduced to analyze the data. The stochastic differential equation associated with the model is then solved. The solution gives us the crack exceedance probability as well as the probability distribution of the random time to reach a specified crack size. These quantities are useful in the reliability assessment of structures made of the tested material. Comparing the analytical result with the experimental result, it is found that the proposed probabilistic fracture mechanics model can reasonably explain the experimental data. For those data that cannot be fitted well by the proposed model, methods of improvement are proposed in the present paper as well.

This study develops a kinematic/kinetic model to evaluate the effect of different handle angles on wrist loading during lifting task. An imaged-based motion analysis system is used to study the movement pattern, force and moment of the wrist joint among nine different handle angles. Six CCD cameras were used to record 3-D trajectories of limb-mounted markers based on the laboratory coordination system as defined by an 8-marker cube. Euler angles are used to describe the orientation of a distal segment reference frame relative to a proximal segment reference frame. Each segment of the upper extremities is regarded as a uniform rigid body with six degrees of freedom. The resultant loading of the wrist joints was determined using an inverse dynamic procedure.

This study indicates that tool handles can be designed or selected to reduce manual loading and the potential for injury during tool use. The mean curve of joint force and moment provided consultations and understandings of the wrist loading during lifting task. In this study, handles that kept the wrist joints in a dorsiflexed and radial deviated position, showed significant reduction in stresses around the surrounding soft tissue.

This study investigates experimentally the velocity of the center of mass with respect to the base of support while subjects step on slippery flooring. The moments of the joints of both legs are also investigated to gain further insights into the source of any correlation found in reacting to slippery perturbation. Twenty-two healthy subjects dressed with safety harness walked first without and then with slippery perturbation, guided by a metronome at 120 steps/min and 90 steps/min cadence. Data were collected from a motion analysis system and force plates. Subjects falling had distinguishably slower velocity of center of mass with respect to the lead stance foot compared to the subjects maintaining balance at contralateral toe off. Larger knee flexion moment and ankle plantar flexion moment in the perturbed leg were found among subjects regaining balance. Faster velocity of center of mass with respect to base of support is noteworthy in subjects regaining balance from slippery perturbation. The response of lower extremities, especially about the thigh and leg in the lead stance limb, were important to improve the velocity of the center of mass with respect to base of support.

The asymptotic fields in an elastic anisotropic composite wedge are considered for a wide range of boundary conditions. It is shown that the eigenfunctions for the near-field and far-field are dual as they are generated by the same set of eigenvalues in general. If the boundary conditions on the wedge faces are the same, an additional eigenfunction may appear in the far-field. Moreover the dual eigenfunctions are used to derive path-independent integrals that relate the near-field to the far-field.

This work presents a rigorous three-dimensional finite element procedure to analyze belt transmission systems. The frictional contact behavior between the belt and the pulley, which accounts for the power loss of the system and the wear of the belt, is investigated in detail. In addition to adopting the transformation matrix to satisfy the geometric conditions on the contact surfaces, the proposed procedure also uses the modified elements with incremental Wilson displacement modes to improve the accuracy due to bending at the end zones of the contact area for the belt. To demonstrate the accuracy and feasibility of the proposed procedure, the analyses for flat and V belt drives are carried out. Excellent correlations between the calculated results and referenced theoretical/experimental solutions are found. The influences of friction coefficients on the deformation, normal and tangential contact forces on the contact surfaces are studied as well. Those will be helpful for the estimation of wear properties and operation efficiency for belt transmission systems.

The effects of blockage and swirl on the macro flow structures of the annular jet past a circular disc are experimentally studied through the time-averaged streamline patterns. In the blockage-effect regime, the flows present multiple modes, single bubble, dual rings, vortex breakdown, and triple rings, in different regimes of blockage ratio and swirl number. The topological models of the flow structures are proposed and discussed according to the measured flow fields to manifest the complex flow structures. The single bubble is a closed recirculation bubble with a stagnation point on the central axis. The dual-ring flow is an open-top recirculsation zone, in which a pair of counter-rotating vortex rings exists in the near wake. The fluids in the dual rings are expelled downstream through a central jet-like swirling flow. A vortex breakdown may occur in the central jet-like swirling flow if the exit swirl number exceeds critical values. When the vortex breakdown interacts with the dual rings, a complex triple-ring flow structure forms. Axial distributions of the local swirl number are presented and discussed. The local swirl number increases with the increase of the exit swirl number and attains the maximum in the dual-ring mode. At large exit swirl numbers where the vortex breakdown occurs, the local swirl number decreases drastically to a low value.

To study the effect of lubrication on the performance and dynamic loads of scroll compressors, a theoretical model is proposed in this paper. The model includes both thermal and dynamic analyses in which the lubrication theory is emphasized in particular. A computer simulation program is written accordingly, and the effects of different design parameters including oil flow rate, injection position, injection temperature and rotating speed are investigated. The output variables such as gas force, contact force, contact moment and thermal performance of the compressor are obtained. Some of the results are compared with experimentally measured data, and good agreements are found between the simulation results and the measured data. The applicability of the simulation program is thus verified.Analytical results from running this program can provide us valuable information on the design improvement of a scroll compressor.

Experiments for heat transfer characteristics of confined circular single jet impingement were conducted. The effect of jet Reynolds number, jet hole-to-plate spacing and heat flux levels on heat transfer characteristics of the heated target surface was examined and presented. The local heat transfer coefficient along the surface is measured and correlations of the stagnation point, local and average Nusselt number are developed and discussed. Finally, comparisons of the present data with existing results were also made.

Laminar heat transfer for large ranges of Reynolds numbers, rotational Reynolds numbers, and Prandtl numbers are studied numerically for incompressible fully developed flow in a circular straight pipe, which is rotating constantly about an axis perpendicular to its own axis under the constant wall temperature gradient condition. There exist four types of local Nusselt number distributions associated with the four different flow regimes for different parameters depending on the relative importance of different forces. Correlations of the averaged Nusselt number are also provided. When the Prandtl number is sufficiently large, the temperature distribution in the core is determined essentially by the secondary flow. Scaling analyses are provided for understanding the essential physics of the problem.

The flow-induced vibration of tubes in a rotated triangular array subject to cross flow is simulated numerically. In the study, the flow field around the tube bundle is computed by solving the continuity and Navier-Stokes equations with assumption of constant fluid properties, and the kε-model for turbulent Reynolds stress. With the flow field known, the fluid forces on the tube surfaces can be calculated, and then the displacement of each tube due to the fluid force can be evaluated. Iteration is needed to obtain the dynamic response of the tube structure in the fluid flow. The parameters in the study are inlet velocity of the cross flow and properties of the tube bundle including natural frequency, damping factor, and mass. Based on the tube response, the critical flow conditions of tube vibration are determined for varying mass damping. Once tube vibrations occur, it is shown that the vibrations of the tubes in the second and fourth tube rows are significant as compared to other tubes. The orbits of the tube vibration look like an ellipse with major axis in the cross-stream direction, implying large lift force on the tubes. The dominant frequency in the spectrum of lift coefficients of the tubes is the same as the natural frequency, and the corresponding amplitude is increased with increasing the inlet velocity. The calculated data predicted for the critical reduced velocity agrees well with the data by Kassera and Strohmeier [17].

During the past several decades, nearly symmetrical tidal appearances, both the time of high water elevation and the tidal range, were noted along the western coast of Taiwan. Research corresponding to anomalous amplifications of tidal range have been well studied by Lin, et al. [1,2], but the particular symmetrical tidal appearance of high water time has not yet been clarified. In this paper, we firstly apply a harmonic analysis to tidal data compiling. Based on the results, the important symmetric appearances of high water time and tidal range can be located. Emphasis is given on the investigation of symmetrical distribution of high water time by superposing two progressive waves that interfere on a topographic shelf with specific width. A conceptual configuration of co-oscillating standing tides is proposed to demonstrate the symmetrical tidal range distributing in the space domain and oscillating in the time domain.

The fractal properties of premixed transient flames propagating downwards through a near-isotropic turbulent flow field in a fan-stirred cruciform burner were investigated. The long vertical section of the cruciform burner was used to provide a downward propagating premixed flame at 1 atm. The large horizontal vessel equipped with a pair of counter-rotating fans and perforated plates at each end was used to generate near-isotropic turbulence. Turbulent flame front images were obtained using high-speed laser sheet imaging for both methane-air and propane-air mixtures. The nondimensional turbulent intensity (u′/SL), Reynolds number based on the integral length scale, and turbulent Karlovitz number were varied from 1 to 10, from 698 to 6032, and from 0.05 to 1.43, respectively. Hundreds of runs for each experimental condition were carried out to obtain sufficient images of these turbulent transient flame fronts just in the central uniform region. These images were then processed to extract fractal dimension, inner and outer cutoffs using both the circle and the caliper methods. It was found that the mean fractal dimension is only 2.18, nearly independent of u′/SL, in support of recent Bunsen-flame results found by Gülder and his co-workers. This contradicts the findings of many previous studies in which the fractal dimension may approach asymptotically to a value of 2.33 when u′/SL > 3. The inner (εi) and outer (ε0) cutoffs are found to be nearly constant for all flames studied, where ε0 is an order of magnitude greater than εi and it is smaller than the integral length scale of unreacted turbulence. Finally, the present fractal characteristics cannot predict turbulent burning velocities correctly when the available fractal closure model was used, indicating a limit of the fractal analysis on prediction of turbulent burning velocities.

Exact electroelastic analysis of hollow and solid circular cylinders subjected to a uniform temperature change and electromechanical loads that do not vary axially are presented. Generalized plane strain and torsion of a piezoelectric cylindrical body is formulated in a state space framework. Uniformly valid exact solutions are obtained for hollow and solid circular cylinders subjected to extension, torsion, pressuring, shearing, electric loadings and a uniform temperature change.

An upwind finite-volume scheme is studied for solving the solutions of two dimensional Euler equations. It based on the MUSCL (Monotone Upstream Scheme for Conservation Laws) approach with the Roe approximate Riemann solver for the numerical flux evaluation. First, dissipation and dispersion relation, and group velocity of the scheme are derived to analyze the capability of the proposed scheme for capturing physical waves, such as acoustic, entropy, and vorticity waves. Then the scheme is greatly enhanced through a strategy on the numerical dissipation to effectively handle aeroacoustic computations. The numerical results indicate that the numerical dissipation strategy allows that the scheme simulates the continuous waves, such as sound and sine waves, at fourth-order accuracy and captures the discontinuous waves, such a shock wave, sharply as well as most of upwind schemes do. The tested problems include linear wave convection, propagation of a sine-wave packet, propagation of discontinuous and sine waves, shock and sine wave interaction, propagation of acoustic, vorticity, and density pulses in an uniform freestream, and two-dimensional traveling vortex in a low-speed freestream.