An Analytical Solution for Three-Dimensional Diffraction of Plane P-Waves by a Hemispherical Alluvial Valley with Saturated Soil Deposits

An analytical solution for the three-dimensional scattering and diffraction of plane P-waves by a hemispherical alluvial valley with saturated soil deposits is developed by employing Fourier-Bessel series expansion technique. Unlike previous studies, in which the saturated soil deposits were simulated with the single-phase elastic theory, in this paper, they are simulated with Biot's dynamic theory for saturated porous media, and the half space is assumed as a single-phase elastic medium. The effects of the dimensionless frequency, the incidence angle of P-wave and the porosity of soil deposits on the surface displacement magnifications of the hemispherical alluvial valley are investigated. Numerical results show that the existence of a saturated hemispherical alluvial valley has much influence on the surface displacement magnifications. It is more reasonable to simulate soil deposits with Biot's dynamic theory when evaluating the displacement responses of a hemispherical alluvial valley with an incidence of P-waves.

A three-dimensional analytical solution of the microheater temperature before bubble nucleation

A three-dimensional analytical solution of the microheater temperature based on heat diffusion equation is developed and compared with experimental results. Dimensionless parameters are introduced to analyze the temperature rise time and the distribution under steady state. To study the microheater temperatures before bubble nucleation, a set of working fluids and microheaters are considered. It is shown that the dimensionless time xi(-)(0) required for the temperature rise from room to 95% of the steady state temperature is about 75, not dependent on working fluids and microheaters. Heat transfer to the surrounding liquid is mainly caused by conduction, not by convection and radiation mechanisms. The microheater length affects the surface temperature uniformity, while its width influences the steady temperatures significantly, yielding the transition from heterogeneous to homogeneous nucleation mechanism from square microheaters to narrow line microheaters. 

Analytical solution of a multidimensional Langevin equation at high friction limits and probability passing over a two-dimensional saddle

The analytical solution of a multidimensional Langevin equation at the overdamping limit is obtained and the probability of particles passing over a two-dimensional saddle point is discussed. These results may break a path for studying further the fusion in superheavy elements synthesis.

Numerical and Analytical Study on the Pull-in Instability of Micro-Structure under Electrostatic Loading

The one-mode analysis method on the pull-in instability of micro-structure under electrostatic loading is presented. Taylor series are used to expand the electrostatic loading term in the one-mode analysis method, which makes analytical solution available. The one-mode analysis is the combination of Galerkin method and Cardan solution of cubic equation. The one-mode analysis offers a direct computation method on the pull-in voltage and displacement. In low axial loading range, it shows little difference with the established multi-mode analysis on predicting the pull-in voltages for three different structures (cantilever, clamped-clamped beams and the plate with four edges simply-supported) studied here. For numerical multi-mode analysis, we also show that using the structural symmetry to select the symmetric mode can greatly reduce both the computation effort and the numerical fluctuation.

Seismic Responses of a Hemispherical Alluvial Valley to SV Waves: A Three-Dimensional Analytical Approximation

An analytical solution to the three-dimensional scattering and diffraction of plane SV-waves by a saturated hemispherical alluvial valley in elastic half-space is obtained by using Fourier-Bessel series expansion technique. The hemispherical alluvial valley with saturated soil deposits is simulated with Biot's dynamic theory for saturated porous media. The following conclusions based on numerical results can be drawn: (1) there are a significant differences in the seismic response simulation between the previous single-phase models and the present two-phase model; (2) the normalized displacements on the free surface of the alluvial valley depend mainly on the incident wave angles, the dimensionless frequency of the incident SV waves and the porosity of sediments; (3) with the increase of the incident angle, the displacement distributions become more complicated; and the displacements on the free surface of the alluvial valley increase as the porosity of sediments increases.

An asymptotic solution of hydrodynamic equations at the mid-plane of a plume in the earths mantle

From the partial differential equations of hydrodynamics governing the movements in the Earth's mantle of a Newtonian fluid with a pressure- and temperature-dependent viscosity, considering the bilateral symmetry of velocity and temperature distributions at the mid-plane of the plume, an analytical solution of the governing equations near the mid-plane of the plume was found by the method of asymptotic analysis. The vertical distribution of the upward velocity, viscosity and temperature at the mid-plane, and the temperature excess at the centre of the plume above the ambient mantle temperature were then calculated for two sets of Newtonian rheological parameters. The results obtained show that the temperature at the mid-plane and the temperature excess are nearly independent of the rheological parameters. The upward velocity at the mid-plane, however, is strongly dependent on the rheological parameters.

Analytical method of predicating the instabilities of a micro arch-shaped beam under electrostatic loading

An arch-shaped beam with different configurations under electrostatic loading experiences either the direct pull-in instability or the snap-through first and then the pull-in instability. When the pull-in instability occurs, the system collides with the electrode and adheres to it, which usually causes the system failure. When the snap-through instability occurs, the system experiences a discontinuous displacement to flip over without colliding with the electrode. The snap-through instability is an ideal actuation mechanism because of the following reasons: (1) after snap-through the system regains the stability and capability of withstanding further loading; (2) the system flips back when the loading is reduced, i.e. the system can be used repetitively; and (3) when approaching snap-through instability the system effective stiffness reduces toward zero, which leads to a fast flipping-over response. To differentiate these two types of instability responses for an arch-shaped beam is vital for the actuator design. For an arch-shaped beam under electrostatic loading, the nonlinear terms of the mid-plane stretching and the electrostatic loading make the analytical solution extremely difficult if not impossible and the related numerical solution is rather complex. Using the one mode expansion approximation and the truncation of the higher-order terms of the Taylor series, we present an analytical solution here. However, the one mode approximation and the truncation error of the Taylor series can cause serious error in the solution. Therefore, an error-compensating mechanism is also proposed. The analytical results are compared with both the experimental data and the numerical multi-mode analysis. The analytical method presented here offers a simple yet efficient solution approach by retaining good accuracy to analyze the instability of an arch-shaped beam under electrostatic loading.

Simple lattice Boltzmann model for traffic flows

A lattice Boltzmann model with 5-bit lattice for traffic flows is proposed. Using the Chapman-Enskog expansion and multi-scale technique, we obtain the higher-order moments of equilibrium distribution function. A simple traffic light problem is simulated by using the present lattice Boltzmann model, and the result agrees well with analytical solution.

A Study of Composite Beam with Shape Memory Alloy Arbitrarily Embedded under Thermal and Mechanical Loadings

The constitutive relations and kinematic assumptions on the composite beam with shape memory alloy (SMA) arbitrarily embedded are discussed and the results related to the different kinematic assumptions are compared. As the approach of mechanics of materials is to study the composite beam with the SMA layer embedded, the kinematic assumption is vital. In this paper, we systematically study the kinematic assumptions influence on the composite beam deflection and vibration characteristics. Based on the different kinematic assumptions, the equations of equilibrium/motion are different. Here three widely used kinematic assumptions are presented and the equations of equilibrium/motion are derived accordingly. As the three kinematic assumptions change from the simple to the complex one, the governing equations evolve from the linear to the nonlinear ones. For the nonlinear equations of equilibrium, the numerical solution is obtained by using Galerkin discretization method and Newton-Rhapson iteration method. The analysis on the numerical difficulty of using Galerkin method on the post-buckling analysis is presented. For the post-buckling analysis, finite element method is applied to avoid the difficulty due to the singularity occurred in Galerkin method. The natural frequencies of the composite beam with the nonlinear governing equation, which are obtained by directly linearizing the equations and locally linearizing the equations around each equilibrium, are compared. The influences of the SMA layer thickness and the shift from neutral axis on the deflection, buckling and post-buckling are also investigated. This paper presents a very general way to treat thermo-mechanical properties of the composite beam with SMA arbitrarily embedded. The governing equations for each kinematic assumption consist of a third order and a fourth order differential equation with a total of seven boundary conditions. Some previous studies on the SMA layer either ignore the thermal constraint effect or implicitly assume that the SMA is symmetrically embedded. The composite beam with the SMA layer asymmetrically embedded is studied here, in which symmetric embedding is a special case. Based on the different kinematic assumptions, the results are different depending on the deflection magnitude because of the nonlinear hardening effect due to the (large) deflection. And this difference is systematically compared for both the deflection and the natural frequencies. For simple kinematic assumption, the governing equations are linear and analytical solution is available. But as the deflection increases to the large magnitude, the simple kinematic assumption does not really reflect the structural deflection and the complex one must be used. During the systematic comparison of computational results due to the different kinematic assumptions, the application range of the simple kinematic assumption is also evaluated. Besides the equilibrium study of the composite laminate with SMA embedded, the buckling, post-buckling, free and forced vibrations of the composite beam with the different configurations are also studied and compared.

Ocean Waves Propagating over a Porous Seabed: Residual and Oscillatory Mechanisms

Two mechanisms for the wave-induced pore pressures in a porous seabed, i.e. oscillatory and residual excess pore pressures, have been observed in laboratory experiments and field measurements. Most previous investigations have focused on one of the mechanisms individually. In this paper, an analytical solution for the wave-induced residual pore pressure, which is not available yet, is derived, and compared with the existing experimental data. With the new solution, a parametric analysis is performed to clarify the applicable ranges of two mechanisms. Then, a simplified approximation for the prediction of wave-induced liquefaction potential is proposed for engineering practice.

A New Version of Smoothed Point Method to Simulate Linear Elastic Wave Propagation

By using the kernel function of the smoothed particle hydrodynamics (SPH) and modification of statistical volumes of the boundary points and their kernel functions, a new version of smoothed point method is established for simulating elastic waves in solid. With the simplicity of SPH kept, the method is easy to handle stress boundary conditions, especially for the transmitting boundary condition. A result improving by de-convolution is also proposed to achieve high accuracy under a relatively large smooth length. A numerical example is given and compared favorably with the analytical solution.

Numerical Investigation on Laser Transformation Hardening with Different Temporal Pulse Shapes

A 3-D numerical model for pulsed laser transformation hardening (LTH) is developed using the finite element method. In this model, laser spatial and temporal intensity distribution, temperature-dependent thermophysical properties of material, and multi-phase transformations are considered. The influence of laser temporal pulse shape on connectivity of hardened zone, maximum surface temperature of material and hardening depth is numerically investigated at different pulse energy levels. Results indicate that these hardening parameters are strongly dependent on the temporal pulse shape. For the rectangular temporal pulse shape, the temperature field obtained from this model is in excellent agreement with analytical solution, and the predicted hardening depth is favorably compared with experimental one. It should be pointed out that appropriate temporal pulse shape should be selected according to pulse energy level in order to achieve desirable hardening quality under certain laser spatial intensity distribution.

Variational Approach to the Thomas-Fermi Equation

By the semi-inverse method, a variational principle is obtained for the Thomas-Fermi equation, then the Ritz method is applied to solve an analytical solution, which is a much simpler and more efficient method.

Static Study of Cantilever Beam stiction under Electrostatic Force Influence

The model and analysis of the cantilever beam adhesion problem under the action of electrostatic force are given. Owing to the nonlinearity of electrostatic force, the analytical solution for this kind of problem is not available. In this paper, a systematic method of generating polynomials which are the exact beamsolutions of the loads with different distributions is provided. The polynomials are used to approximate the beam displacement due to electrostatic force. The equilibrium equation offers an answer to how the beam deforms but no information about the unstuck length. The derivative of the functional with respect to the unstuck length offers such information. But to compute the functional it is necessary to know the beam deformation. So the problem is iteratively solved until the results are converged. Galerkin and Newton-Raphson methods are used to solve this nonlinear problem. The effects of dielectric layer thickness and electrostatic voltage on the cantilever beamstiction are studied.The method provided in this paper exhibits good convergence. For the adhesion problem of cantilever beam without electrostatic voltage, the analytical solution is available and is also exactly matched by the computational results given by the method presented in this paper.