Second-order two-scale analysis and numerical algorithms for the hyperbolic–parabolic equations with rapidly oscillating coefficients

We study the hyperbolic–parabolic equations with rapidly oscillating coefficients. The formal second-order two-scale asymptotic expansion solutions are constructed by the multiscale asymptotic analysis. In addition, we theoretically explain the importance of the second-order two-scale solution by the error analysis in the pointwise sense. The associated explicit convergence rates are also obtained. Then a second-order two-scale numerical method based on the Newmark scheme is presented to solve the equations. Finally, some numerical examples are used to verify the effectiveness and efficiency of the multiscale numerical algorithm we proposed.

Statistical second-order two-scale analysis and computation for heat conduction problem with radiation boundary condition in porous materials

This paper discusses a statistical second-order two-scale（SSOTS） analysis and computation for a heat conduction problem with a radiation boundary condition in random porous materials.Firstly,the microscopic configuration for the structure with random distribution is briefly characterized.Secondly,the SSOTS formulae for computing the heat transfer problem are derived successively by means of the construction way for each cell.Then,the statistical prediction algorithm based on the proposed two-scale model is described in detail.Finally,some numerical experiments are proposed,which show that the SSOTS method developed in this paper is effective for predicting the heat transfer performance of porous materials and demonstrating its significant applications in actual engineering computation.

Second-order two-scale computational method for ageing linear viscoelastic problem in composite materials with periodic structure

The correspondence principle is an important mathematical technique to compute the non-ageing linear viscoelastic problem as it allows to take advantage of the computational methods originally developed for the elastic case. However, the correspon- dence principle becomes invalid when the materials exhibit ageing. To deal with this problem, a second-order two-scale （SOTS） computational method in the time domain is presented to predict the ageing linear viscoelastic performance of composite materials with a periodic structure. First, in the time domain, the SOTS formulation for calcu- lating the effective relaxation modulus and displacement approximate solutions of the ageing viscoelastic problem is formally derived. Error estimates of the displacement ap- proximate solutions for SOTS method are then given. Numerical results obtained by the SOTS method are shown and compared with those by the finite element method in a very fine mesh. Both the analytical and numerical results show that the SOTS computational method is feasible and efficient to predict the ageing linear viscoelastic performance of composite materials with a periodic structure.

Stochastic Second-Order Two-Scale Method for Predicting the Mechanical Properties of Composite Materials with Random Interpenetrating Phase

In this paper,a stochastic second-order two-scale(SSOTS)method is proposed for predicting the non-deterministic mechanical properties of composites with random interpenetrating phase.Firstly,based on random morphology description functions(RMDF),the randomness of the material properties of the constituents as well as the correlation among these random properties are fully characterized through the topologies of the constituents.Then,by virtue of multiscale asymptotic analysis,the random effective quantities such as stiffness parameters and strength parameters along with their numerical computation formulae are derived by a SSOTS strategy combined with the Monte-Carlo method.Finally,the SSOTS method developed in this paper shows an excellent computational accuracy,and therefore present an important advance towards computationally efficient multiscale modeling frameworks considering microstructure uncertainties.

Dynamic thermo-mechanical coupled response of random particulate composites:A statistical two-scale method

This paper focuses on the dynamic thermo-mechanical coupled response of random particulate composite materials. Both the inertia term and coupling term are considered in the dynamic coupled problem. The formulation of the problem by a statistical second-order two-scale （SSOTS） analysis method and the algorithm procedure based on the finite-element difference method are presented. Numerical results of coupled cases are compared with those of uncoupled cases. It shows that the coupling effects on temperature, thermal flux, displacement, and stresses are very distinct, and the micro- characteristics of particles affect the coupling effect of the random composites. Furthermore, the coupling effect causes a lag in the variations of temperature, thermal flux, displacement, and stresses.

Second-order two-scale method for bending behavior analysis of composite plate with 3-D periodic configuration and its approximation

This paper considers the bending behaviors of composite plate with 3-D periodic configuration.A second-order two-scale(SOTS)computational method is designed by means of construction way.First,by 3-D elastic composite plate model,the cell functions which are defined on the reference cell are constructed.Then the effective homogenization parameters of composites are calculated,and the homogenized plate problem on original domain is defined.Based on the Reissner-Mindlin deformation pattern,the homogenization solution is obtained.And then the SOTS’s approximate solution is obtained by the cell functions and the homogenization solution.Second,the approximation of the SOTS’s solution in energy norm is analyzed and the residual of SOTS’s solution for 3-D original in the pointwise sense is investigated.Finally,the procedure of SOTS’s method is given.A set of numerical results are demonstrated for predicting the effective parameters and the displacement and strains of composite plate.It shows that SOTS’s method can capture the 3-D local behaviors caused by3-D micro-structures well.

Heat transfer of nanofluids considering nanoparticle migration and second-order slip velocity

The heat transfer of a magnetohydrodynamics nanofluid inside an annulus considering the second-order slip condition and nanoparticle migration is theoret-ically investigated. A second-order slip condition, which appropriately represents the non-equilibrium region near the interface, is prescribed rather than the no-slip condition and the linear Navier slip condition. To impose different temperature gradients, the outer wall is subjected to q2, the inner wall is subjected to q1, and q1 〉 q2. A modified two-component four-equation non-homogeneous equilibrium model is employed for the nanofiuid, which have been reduced to two-point ordinary boundary value differential equations in the consideration of the thermally and hydrodynamically fully developed flow. The homotopy analysis method （HAM） is employed to solve the equations, and the h-curves are plotted to verify the accuracy and efficiency of the solutions. Moreover, the effects of the physical factors on the flow and heat transfer are discussed in detail, and the semi-analytical relation between NUB and NBT is obtained.

Two-scale analysis method for predicting heat transfer performance of composite materials with random grain distribution

A two-scale analysis (TSA) method for predicting the heat transfer performance of composite materials with the random distribution of same-scale grains is presented. First the representation of the materials with the random distribution is briefly described. Then the two-scale analysis formulation of heat transfer behavior of the materials with random grain distribution of small periodicity is formally derived by means of construction way for each cell. Finally the numerical result on the heat transfer parameters of composite materials is shown. The numerical result shows that TSA is effective to predict the heat transfer performance of composite materials with random grain distribution.

A Method for Calculation of Low-Frequency Slow Drift Motions Based on NURBS for Floating Bodies

Through a higher-order boundary element method based on NURBS （Non-uniform Rational B-splines）, the calculation of second-order low-frequency forces and slow drift motions is conducted for floating bodies. In the floating body＇s inner domain, an auxiliary equation is obtained by applying a Green function which satisfies the solid surface condition. Then, the auxiliary equation and the velocity potential equation are combined in the fluid domain to remove the solid angle coefficient and the singularity of the double layer potentials in the integral equation. Thus, a new velocity potential integral equation is obtained. The new equation is extended to the inner domain to reheve the irregular frequency effects; on the basis of the order analysis, the comparison is made about the contribution of all integral terms with the result in the second-order tow-frequency problem; the higher-order boundary element method based on NURBS is apphed to calculate the geometric position and velocity potentials; the slow drift motions are calculated by the spectrum analysis method. Removing the solid angle coefficient can apply NURBS technology to the hydrodynamic calculation of floating bodies with complex surfaces, and the extended boundary integral method can reduce the irregular frequency effects. Order analysis shows that free surface integral can be neglected, and the numerical results can also prove the correctness of order analysis. The results of second-order low-frequency forces and slow drift motions and the comparison with the results from references show that the application of the NURBS technology to the second-order low-frequency problem is of high efficiency and credible results.

The use of the node-based smoothed finite element method to estimate static and seismic bearing capacities of shallow strip footings

The node-based smoothed finite element method(NS-FEM)is shortly presented for calculations of the static and seismic bearing capacities of shallow strip footings.A series of computations has been performed to assess variations in seismic bearing capacity factors with both horizontal and vertical seismic accelerations.Numerical results obtained agree very well with those using the slip-line method,revealing that the magnitude of the seismic bearing capacity is highly dependent upon the combinations of various directions of both components of the seismic acceleration.An upward vertical seismic acceleration reduces the seismic bearing capacity compared to the downward vertical seismic acceleration in calculations.In addition,particular emphasis is placed on a separate estimation of the effects of soil and superstructure inertia on each seismic bearing capacity component.While the effect of inertia forces arising in the soil on the seismic bearing capacity is non-trivial,and the superstructure inertia is the major contributor to reductions in the seismic bearing capacity.Both tables and charts are given for practical application to the seismic design of the foundations.

A stable CS-FEM for the static and seismic stability of a single square tunnel in the soil where the shear strength increases linearly with depth

A numerical procedure using a stable cell-based smoothed finite element method(CS-FEM)is presented for estimation of stability of a square tunnel in the soil where the shear strength increases linearly with depth.The kinematically admissible displacement fields are approximated by uniform quadrilateral elements in conjunction with the strain smoothing technique,eliminating volumetric locking issues and the singularity associated with the MohreCoulomb model.First,a rich set of simulations was performed to compute the static stability of a square tunnel with different geometries and soil conditions.The presented results are in excellent agreement with the upper and lower bound solutions using the standard finite element method(FEM).The stability charts and tables are given for practical use in the tunnel design,along with a newly proposed formulation for predicting the undrained stability of a single square tunnel.Second,the seismic stability number was computed using the present numerical approach.Numerical results reveal that the seismic stability number reduces with an increasing value of the horizontal seismic acceleration(a_(h)),for both cases of the weightless soil and the soil with unit weight.Third,the link between the static and seismic stability numbers is described using corrective factors that represent reductions in the tunnel stability due to seismic loadings.It is shown from the numerical results that the corrective factor becomes larger as the unit weight of soil mass increases;however,the degree of the reduction in seismic stability number tends to reduce for the case of the homogeneous soil.Furthermore,this advanced numerical procedure is straightforward to extend to three-dimensional(3D)limit analysis and is readily applicable for the calculation of the stability of tunnels in highly anisotropic and heterogeneous soils which are often encountered in practice.