Two kinds of contact problems in three-dimensional icosahedral quasicrystals

Two kinds of contact problems, i.e., the frictional contact problem and the adhesive contact problem, in three-dimensional （3D） icosahedral quasicrystals are dis- cussed by a complex variable function method. For the frictional contact problem, the contact stress exhibits power singularities at the edge of the contact zone. For the adhe- sive contact problem, the contact stress exhibits oscillatory singularities at the edge of the contact zone. The numerical examples show that for the two kinds of contact problems, the contact stress exhibits singularities, and reaches the maximum value at the edge of the contact zone. The phonon-phason coupling constant has a significant effect on the contact stress intensity, while has little impact on the contact stress distribution regu- lation. The results are consistent with those of the classical elastic materials when the phonon-phason coupling constant is 0. For the adhesive contact problem, the indentation force has positive correlation with the contact displacement, but the phonon-phason cou- pling constant impact is barely perceptible. The validity of the conclusions is verified.

Three-dimensional forward modeling for magnetotelluric sounding by finite element method

A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.

Three-dimensional land FD-CSEM forward modeling using edge finite-element method

A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.

An Efficient Three-Dimensional Coupled Normal Mode Model and Its Application to Internal Solitary Wave Problems

We present an efficient three-dimensional coupled-mode model based on the Fourier synthesis technique. In principle, this model is a one-way model, and hence provides satisfactory accuracy for problems where the forward scattering dominates. At the same time, this model provides an efficiency gain of an order of magnitude or more over two-way coupled-mode models. This model can be applied to three-dimensional range-dependent problems with a slowly varying bathymetry or internal waves. A numerical example of the latter is demonstrated in this work. Comparisons of both accuracy and efficiency between the present model and a benchmark model are also provided.

THE CAPTURE RANGES AND THE DANGEROUS RANGES FOR THE THREE-DIMENSIONAL TWO-AIRCRAFT COMBAT PROBLEM

In this paper, we use differential game theory to study the three-dimensional two-aircraft air-to-air combat problem. We give the ways to determine the Capture Ranges (CR) and the Dangerous Ranges (DR) for these two aircraft according to the target entry directions, barrier and isochronic lines respectively. The simulations are given by referring to two sets of real aircraft parameters. After discussing the simulation results, we have obtained some conclusions that match the real air-to-air combat situation quite well.

Exact analytical solution to three-dimensional phase change heat transfer problems in biological tissues subject to freezing

Analytically solving a three-dimensional （3-D） bioheat transfer problem with phase change during a freezing process is extremely difficult but theoretically important. The moving heat source model and the Green function method are introduced to deal with the cryopreservation process of in vitro biomaterials. Exact solutions for the 3-D temperature transients of tissues under various boundary conditions, such as totally convective cooling, totally fixed temperature cooling and a hybrid between them on tissue surfaces, are obtained. Furthermore, the cryosurgical process in living tissues subject to freezing by a single or multiple cryoprobes is also analytically solved. A closed-form analytical solution to the bioheat phase change process is derived by considering contributions from blood perfusion heat transfer, metabolic heat generation, and heat sink of a cryoprobe. The present method is expected to have significant value for analytically solving complex bioheat transfer problems with phase change.

The application of three-dimensional textbook in English teaching

The paper discusses the application of three-dimensional textbook in English teaching. The three-dimensional textbook plays an incomparable function than the traditional textbook and thus has an extensive application in the college English education. It demonstrates the application through the usage of the three-dimensional textbook in the real teaching and shows that excessive usage of sound and pictures or drawings may distract the attention of students. It raises several suggestions at the end.

Three-dimensional frequency-domain full waveform inversion based on the nearly-analytic discrete method

The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the time-domain.This study applies the NAD method to solving three-dimensional(3D)acoustic wave equations in the frequency-domain.This forward modeling approach is then used as the“engine”for implementing 3D frequency-domain full waveform inversion(FWI).In the numerical modeling experiments,synthetic examples are first given to show the superiority of the NAD method in forward modeling compared with traditional finite difference methods.Synthetic 3D frequency-domain FWI experiments are then carried out to examine the effectiveness of the proposed methods.The inversion results show that the NAD method is more suitable than traditional methods,in terms of computational cost and stability,for 3D frequency-domain FWI,and represents an effective approach for inversion of subsurface model structures.

Interactions of penny-shaped cracks in three-dimensional solids

The interaction of arbitrarily distributed penny-shaped cracks in three-dimensional solids is analyzed in this paper. Using oblate spheroidal coordinates and displacement functions, an analytic method is devel- oped in which the opening and the sliding displacements on each crack surface are taken as the basic unknown functions. The basic unknown functions can be expanded in series of Legendre polynomials with unknown coefficients. Based on superposition technique, a set of governing equations for the unknown coefficients are formulated from the traction free conditions on each crack surface. The boundary collocation procedure and the average method for crack-surface tractions are used for solving the governing equations. The solution can be obtained for quite closely located cracks. Numerical examples are given for several crack problems. By comparing the present results with other existing results, one can conclude that the present method provides a direct and efficient approach to deal with three-dimensional solids containing multiple cracks.

STABILIZATION EFFECT OF FRICTIONS FOR TRANSONIC SHOCKS IN STEADY COMPRESSIBLE EULER FLOWS PASSING THREE-DIMENSIONAL DUCTS

Transonic shocks play a pivotal role in designation of supersonic inlets and ramjets.For the three-dimensional steady non-isentropic compressible Euler system with frictions,we constructe a family of transonic shock solutions in rectilinear ducts with square cross-sections.In this article,we are devoted to proving rigorously that a large class of these transonic shock solutions are stable,under multidimensional small perturbations of the upcoming supersonic flows and back pressures at the exits of ducts in suitable function spaces.This manifests that frictions have a stabilization effect on transonic shocks in ducts,in consideration of previous works which shown that transonic shocks in purely steady Euler flows are not stable in such ducts.Except its implications to applications,because frictions lead to a stronger coupling between the elliptic and hyperbolic parts of the three-dimensional steady subsonic Euler system,we develop the framework established in previous works to study more complex and interesting Venttsel problems of nonlocal elliptic equations.

THE CHARACTERISTIC FINITE ELEMENT ALTERNATING-DIRECTION METHOD AND ANALYSIS FOR THREE-DIMENSIONAL NUMERICAL RESERVOIR SIMULATION

Petroleum science has made remarkable progress in organic geochcmistry and in the research into the theories of petroleum origin, its transport and accumulation. In estimating the oil-gas resources of a basin, the knowledge of its evolutionary history and especially the numerical computation of fluid flow and the history of its changes under heat is vital. The mathematical model can be described as a coupled system of nonlinear partial differentical equations with initial-boundary value problems. This thesis, from actual conditions such as the effect of fluid compressibility and the three-dimensional characteristic of large-scale science-engineering computation, we put forward a kind of characteristic finite element alternating-direction schemes and obtain optimal order estimates in L^2 norm for the error in the approximate assumption.

NONLINEAR THREE-DIMENSIONAL ANALYSIS OF COMPOSITELAMINATED PLATES

In this paper. an analytic method is. presented. for the research of nonlinear Three-dimensional problems of composite laminated plates. The perturbation method and the variational principle are used to satisfy the basic differential equations and the boundary conditions of the three-dimensional theory of elasticity. The nonlinear three-dimensional problems are studied .for composite anisotropic circular laminas and laminates subjected to transverse loading. The perturbation series solutions of high accuracy are obtained. A large number of results show that transverse normal stress and transverse shear stresses are very important in the nonlinear three-dimensional analysis of laminated plates.

THE REGULARIZED SOLUTION APPROXIMATION OF FORWARD/BACKWARD PROBLEMS FOR A FRACTIONAL PSEUDO-PARABOLIC EQUATION WITH RANDOM NOISE

This paper deals with the forward and backward problems for the nonlinear fractional pseudo-parabolic equation ut+(-Δ)^(s1)ut+β(-Δ)^(s2)u=F(u,x,t)subject o random Gaussian white noise for initial and final data.Under the suitable assumptions s1,s2andβ,we first show the ill-posedness of mild solutions for forward and backward problems in the sense of Hadamard,which are mainly driven by random noise.Moreover,we propose the Fourier truncation method for stabilizing the above ill-posed problems.We derive an error estimate between the exact solution and its regularized solution in an E‖·‖Hs22norm,and give some numerical examples illustrating the effect of above method.

The Forward and Inverse Problem Based on Magneto-Acoustic Tomography with Current Injection

The Magneto-acoustic Tomography with Current Injection (MAT-CI) is a new biological electrical impedance imaging technique that combines Electrical Impedance Tomography (EIT) with Ultrasonic Imaging (UI), which possesses the non-invasive and high-contrast of the EIT and the high-resolution of the UI. The MAT-CI is expected to acquire high quality image and embraces a wide application. Its principle is to put the conductive sample in the Static Magnetic Field(SMF) and inject a time-varying current, during which the SMF and the current interact and generate the Lorentz Force that inspire ultrasonic signal received by the ultrasonic transducers positioned around the sample. And then according to related reconstruction algorithm and ultrasonic signal, electrical conductivity image is obtained. In this paper, a forward problem mathematical model of the MAT-CI has been set up to deduce the theoretical equation of the electromagnetic field and solve the sound source distribution by Green’s function. Secondly, a sound field restoration by Wiener filtering and reconstruction of current density by time-rotating method have deduced the Laplace’s equation that caters to the current density to further acquire the electrical conductivity distribution image of the sample through iteration method. In the end, double-loop coils experiments have been conducted to verify its feasibility.

Three-Dimensional Thermal-Stress Analysis of Semi-infinite Transversely Isotropic Composites

By making use of the direct integration method,an exact analysis of the general three-dimensional thermoelasticity problem is performed for the case of a transversely isotropic homogeneous half-space subject to local thermal and force loadings.The material plane of isotropy is assumed to be parallel to the limiting surface of the halfspace.By reducing the original thermoelasticity equations to the governing ones for individual stress-tensor components,the effect of material anisotropy in the stress field is analyzed with regard to the feasibility requirement,i.e.,the finiteness of the stress field at a distance from the disturbed area.As a result,the solution is constructed in the form of explicit analytical dependencies on the force and thermal loadings for various kinds of transversely isotropic materials and agrees with the basic principles of the continua mechanics.The solution can be efficiently used as a benchmark one for the direct computation of temperature and thermal stresses in transversely isotropic semi-infinite domains,as well as for the verification of solutions constructed by different means.

Three-dimensional numerical manifold method for heat conduction problems with a simplex integral on the boundary

The three-dimensional numerical manifold method(3D-NMM),which is based on the derivation of Galerkin's variation,is a powerful calculation tool that uses two cover systems.The 3D-NMM can be used to handle continue-discontinue problems and extend to THM coupling.In this study,we extended the 3D-NMM to simulate both steady-state and transient heat conduction problems.The modelling was carried out using the raster methods(RSM).For the system equation,a variational method was employed to drive the discrete equations,and the crucial boundary conditions were solved using the penalty method.To solve the boundary integral problem,the face integral of scalar fields and two-dimensional simplex integration were used to accurately describe the integral on polygonal boundaries.Several numerical examples were used to verify the results of 3D steady-state and transient heat-conduction problems.The numerical results indicated that the 3D-NMM is effective for handling 3D both steadystate and transient heat conduction problems with high solution accuracy.

Recursive super-convergence computation for multi-dimensional problems via one-dimensional element energy pro jection technique

This paper presents a strategy for computation of super-convergent solutions of multi-dimensional problems in the finite element method （FEM） by recursive application of the one-dimensional （1D） element energy projection （EEP） technique. The main idea is to conceptually treat multi-dimensional problems as generalized 1D problems, based on which the concepts of generalized 1D FEM and its consequent EEP formulae have been developed in a unified manner. Equipped with these concepts, multi-dimensional problems can be recursively discretized in one dimension at each step, until a fully discretized standard finite element （FE） model is reached. This conceptual dimension-by- dimension （D-by-D） discretization procedure is entirely equivalent to a full FE discretization. As a reverse D-by-D recovery procedure, by using the unified EEP formulae together with proper extraction of the generalized nodal solutions, super-convergent displacements and first derivatives for two-dimensional （2D） and three-dimensional （3D） problems can be obtained over the domain. Numerical examples of 3D Poisson＇s equation and elasticity problem are given to verify the feasibility and effectiveness of the proposed strategy.

NON-INTERIOR SMOOTHING ALGORITHM FOR FRICTIONAL CONTACT PROBLEMS

A new algorithm for solving the three-dimensional elastic contact problem with friction is presented. The algorithm is a non-interior smoothing algorithm based on an NCP-function. The parametric variational principle and parametric quadratic programming method were applied to the analysis of three-dimensional frictional contact problem. The solution of the contact problem was finally reduced to a linear complementarity problem, which was reformulated as a system of nonsmooth equations via an NCP-function. A smoothing approximation to the nonsmooth equations was given by the aggregate function. A Newton method was used to solve the resulting smoothing nonlinear equations. The algorithm presented is easy to understand and implement. The reliability and efficiency of this algorithm are demonstrated both by the numerical experiments of LCP in mathematical way and the examples of contact problems in mechanics.

基于改进正问题模型和模因算法的ECT图像重建

为提高电容层析成像的图像重建质量,提出了一种基于改进正问题模型和模因算法的图像重建算法,即Im-MA算法。通过对灵敏度矩阵进行平滑滤波和降维处理,改进正问题模型。使用L_(2)范数作为数据拟合的测度,并将全变差作为正则项,构建目标函数,将图像重建问题转化为最优化问题。采用麻雀搜索算法和天牛须搜索算法相结合的模因(memetic algorithm,MA)算法求解目标函数。仿真与实验结果表明:与常用的Landweber算法相比,Im-MA算法的重建误差更小,且重建图像与真实分布更加接近。Im-MA算法为解决ECT逆问题提供了一种有效的新方法。

一种求解包含问题的算法研究

一种实希尔伯特空间中求解包含问题的算法被提出,所提出的算法基于向前向后方法、压缩方法、惯性方法和无需搜索的自适应步长。算法的特点为迭代中多次使用惯性加速方法,且自适应步长随着迭代次数增加可能增大。在包含问题解集非空、一个映射极大单调、另一个映射单调且利普希茨连续的假设下,算法的强收敛性被证明。