Development and elaboration of numerical method for simulating gas–liquid–solid three-phase flows based on particle method

A numerical method for simulating gas–liquid–solid three-phase flows based on the moving particle semi-implicit (MPS) approach was developed in this study. Computational instability often occurs in multiphase flow simulations if the deformations of the free surfaces between different phases are large, among other reasons. To avoid this instability, this paper proposes an improved coupling procedure between different phases in which the physical quantities of particles in different phases are calculated independently. We performed numerical tests on two illustrative problems: a dam-break problem and a solid-sphere impingement problem. The former problem is a gas–liquid two-phase problem, and the latter is a gas–liquid–solid three-phase problem. The computational results agree reasonably well with the experimental results. Thus, we confirmed that the proposed MPS method reproduces the interaction between different phases without inducing numerical instability.     

Advanced response surface method for mechanical reliability analysis

Based on the classical response surface method (RSM), a novel RSM using improved experimental points (EPs) is presented for reliability analysis. Two novel points are included in the presented method. One is the use of linear interpolation, from which the total EPs for determining the RS are selected to be closer to the actual failure surface; the other is the application of sequential linear interpolation to control the distance between the surrounding EPs and the center EP, by which the presented method can ensure that the RS fits the actual failure surface in the region of maximum likelihood as the center EPs converge to the actual most probable point (MPP). Since the fitting precision of the RS to the actual failure surface in the vicinity of the MPP, which has significant contribution to the probability of the failure surface being exceeded, is increased by the presented method, the precision of the failure probability calculated by RS is increased as well. Numerical examples illustrate the accuracy and efficiency of the presented method.     

A reformulated synthetic turbulence generation method for a zonal RANS–LES method and its application to zero-pressure gradient boundary layers

A synthetic turbulence generation (STG) method for subsonic and supersonic flows at low and moderate Reynolds numbers to provide inflow distributions of zonal Reynolds-averaged Navier–Stokes (RANS) – large-eddy simulation (LES) methods is presented. The STG method splits the LES inflow region into three planes where a local velocity signal is decomposed from the turbulent flow properties of the upstream RANS solution. Based on the wall-normal position and the local flow Reynolds number, specific length and velocity scales with different vorticity content are imposed at the inlet plane of the boundary layer. The quality of the STG method for incompressible and compressible zero-pressure gradient boundary layers is shown by comparing the zonal RANS–LES data with pure LES, pure RANS, and direct numerical simulation (DNS) solutions. The distributions of the time and spanwise wall-shear stress, Reynolds stress distributions, and two point correlations of the zonal RANS–LES simulations are smooth in the transition region and in good agreement with the pure LES and reference DNS findings. The STG approach reduces the RANS-to-LES transition length to less than four boundary-layer thicknesses.     

The moiré method applied to three-dimensional elastic problems

The purpose of this paper is to develop the technique of moiré for three-dimensional strain and stress analysis. The classical problem of contact of a sphere with a semi-infinite space is studied, and experimental results are compared with those of Hertz's solution. Results show that the method is accurate and has a high sensitivity. The results prove also that the strain distribution in a very small region with a steep gradient can be detected quite precisely.Paper was presented at the 1964 SESA Spring Meeting held in Salt Lake City, Utah, on May 6–8.     

A multibody-based dynamic simulation method for electrostatic actuators

A numerical simulation method is developed to analyze the dynamic responses of electrostatic actuators, which are electromechanically-coupled systems. The developed method can be used to determine the dynamic responses of cantilever-type switches, which are an example of typical MEMS (Micro-Electro-Mechanical System) devices driven by an electrostatic force. We propose the approach that adopts a point charge to deal with electric field effects between electrodes. This approach may be considered as a lumped parameter model for the electrostatic interactions. An advantage of this model may be the easy incorporation of the electrostatic effects between electrodes into a multibody dynamics analysis algorithm. The resulting equations contain the variables for position, velocity, and electric charge to describe the motion of the masses and the charges on the electrodes in a system. By solving these equations simultaneously, the dynamic response of an electrostatically-driven system can be correctly simulated. In order to realize this approach, we implement the procedures into RecurDyn, the multibody dynamics software developed by the authors. The developed numerical simulation tool was evaluated by applying it to cantilever-type electrostatic switches in many different driving conditions. The results suggest that the developed tool may be useful for predicting behaviors of electrostatic actuators in testing as well as in design.     

Photoviscoelasticity—An experimental method in viscoelastic stress analysis

In this paper, the photoviscoelasticity method of viscoelastic stress analysis has been discussed in detail.lt is shown that,in order to avoid the effects of shrinkage andaging in the test specimens, it is suggested that the specimens should be tempered for three days at a temperature of 60℃ before starting the experiments, and the temperature filtering arrangement is recommended in the experimental setups to keep the temperature absolutely constant. Besides the axi-symmetrical time-dependent stress state, the determination of the principle axes of the refraction tensor experimentally remains an unsufficiently solved problem. To avoid dynamic effect in the step wise loading, the time of measurement in every step should be limited in about one secend.     

Multiresolution method for bending of plates with complex shapes

A high-accuracy multiresolution method is proposed to solve mechanics problems subject to complex shapes or irregular domains. To realize this method, we design a new wavelet basis function, by which we construct a fifth-order numerical scheme for the approximation of multi-dimensional functions and their multiple integrals defined in complex domains. In the solution of differential equations, various derivatives of the unknown function are denoted as new functions. Then, the integral relations ...     

A Newton type iterative inverse method for heat-conduction problems

An inverse problem for identification of the coefficient in heat-conduction equation is considered. After reducing the problem to a nonlinear ill-posed operator equation, Newton type iterative methods are considered. The implicit iterative method is applied to the linearized Newton equation, and the key step in the process is that a new reasonable a posteriori stopping rule for the inner iteration is presented. Numerical experiments for the new method as well as for Tikhonov method and Bakushikskii method are given, and these results show the obvious advantages of the new method over the other ones.     

Hole-drilling method using grating rosette and Moiré interferometry

The hole-drilling method is one of the most wellknown methods for measuring residual stresses. To identify unknown plane stresses in a specimen, a circular hole is first drilled in the infinite plate under plane stress, then the strains resulting from the hole drilling is measured. The strains may be acquired from interpreting the Moire signature around the hole. In crossed grating Moire interferometry, the horizontal and vertical displacement fields （u and v） can be obtained to determinate two strain fields and one shearing strain field. In this paper, by means of Moire interferometry and three directions grating （grating rosette） developed by the authors, three displacement fields （u, v and s） are obtained to acquire three strain fields. As a practical application, the hole-drilling method is adopted to measure the relief strains for aluminum and fiber reinforced composite. It is a step by step method; in each step a single laminate or equivalent depth is drilled to find some relationships between the drilling depth and the residual strains relieved in the fiber reinforced composite materials.     

An improved proximal-based decomposition method for structured monotone variational inequalities

The proximal-based decomposition method was originally proposed by Chen and TebouUe(Math.Programming,1994,64:81-101 for solving convex minimization problems.This paper extends it to solving monotone variational inequalities associated with separable structures with the improvements that the restrictive assumptions on the involved parameters are much relaxed,and thus makes it practical to solve the subprob- lems easily.Without additional assumptions,global convergence of the new method is proved under the same mild assumptions on the problem's data as the original method.     

An ε-generalized gradient projection method for nonlinear minimax problems

In this paper, combining the techniques of ε-generalized gradient projection and Armjio’s line search, we present a new algorithm for the nonlinear minimax problems. At each iteration, the improved search direction is generated by an ε-generalized gradient projection explicit formula. Under some mild assumptions, the algorithm possesses global and strong convergence. Finally, some preliminary numerical results show that the proposed algorithm performs efficiently.     

An efficient differential quadrature method for fractional advection–diffusion equation

Nonlinear Dynamics - This article studies a direct numerical approach for fractional advection–diffusion equations (ADEs). Using a set of cubic trigonometric B-splines as test functions, a...     

Interval finite element method and its application on anti-slide stability analysis

The problem of interval correlation results in interval extension is discussed by the relationship of interval-valued functions and real-valued functions. The methods of reducing interval extension are given. Based on the ideas of the paper, the formulas of sub-interval perturbed finite element method based on the elements are given. The sub-interval amount is discussed and the approximate computation formula is given. At the same time, the computational precision is discussed and some measures of improving computational efficiency are given. Finally, based on sub-interval perturbed finite element method and anti-slide stability analysis method, the formula for computing the bounds of stability factor is given. It provides a basis for estimating and evaluating reasonably anti-slide stability of structures.