Interlaminar Crack Propagation Analysis of ENF Specimens Based on the Cohesive Zone Model北大核心CSCD

Based on the classical laminated plate theory and the cohesive zone model, a theoretical model for general delamination cracked laminates was established for crack propagation of pure mode Ⅱ ENF specimens. Compared with the conventional beam theory, the proposed model fully considered the softening process of the cohesive zone and introduced the nonlinear behavior of ENF specimens before failure. The predicted failure load is smaller than that under the beam theory and closer to the experimental data in literatures. Compared with the beam theory with only fracture toughness considered, the proposed model can simultaneously analyze the influences of the interface strength, the fracture toughness and the initial interface stiffness on the load-displacement curves in ENF tests. The results show that, the interface strength mainly affects the mechanical behavior of specimens before failure, but has no influence on crack propagation. The fracture toughness is the main parameter affecting crack propagation, and the initial interface stiffness only affects the linear elastic loading stage. The cohesive zone length increases with the fracture toughness and decreases with the interface strength. The effect of the interface strength on the cohesive zone length is more obvious than that of the fracture toughness. When the adhesive zone tip reaches the half length of the specimen, the adhesive zone length will decrease to a certain extent. Copyright ©2022 Applied Mathematics and Mechanics. All rights reserved.     

Uniqueness to Some Inverse Source Problems for the Wave Equation in Unbounded Domains

This paper is concerned with inverse acoustic source problems in an unbounded domain with dynamical boundary surface data of Dirichlet kind.The measurement data are taken at a surface far away from the source support.We prove uniqueness in recovering source terms of the form f(x)g(t)and f(x1,x2,t)h(x3),where g(t)and h(x3)are given and x=(x1,x2,x3)is the spatial variable in three dimensions.Without these a priori information,we prove that the boundary data of a family of solutions can be used to recover general source terms depending on both time and spatial variables.For moving point sources radiating periodic signals,the data recorded at four receivers are prove sufficient to uniquely recover the orbit function.Simultaneous determination of embedded obstacles and source terms was verified in an inhomogeneous background medium using the observation data of infinite time period.Our approach depends heavily on the Laplace transform.     

Mesoscopic Numerical Study on Flow Boiling Heat Transfer Performance in Channels With Multiple Rectangular Heaters北大核心CSCD

The flow boiling phenomenon in a channel with multiple rectangular heaters under a constant wall temperature was numerically studied with the lattice Boltzmann method. The effects of spacings between heaters, heater lengths and heater surface wettabilities on the bubble morphology, the bubble area and the heat flux on the heater surface, were studied. The results show that, the bubble growth rate increases with the spacing between heaters. The larger the bubble area is, the earlier the nucleated bubbles will leave the heater surface. The corresponding boiling heat transfer performance increases by 12% with the spacing between heaters growing from 250 lattices to 1 000 lattices. On the other hand, the longer the heater length is, the earlier the bubble will nucleate and leave the heater surface, and the better the boiling heat transfer performance will be. The boiling heat transfer performance increases by 13% with the heater length rising from 16 lattices to 22 lattices. In addition, the bubble nucleates later on the hydrophilic surface than on the hydrophobic surface. Compared with the hydrophilic surface, the hydrophobic surface retains residual bubbles after the leaving of bubbles from the heater. The average heat flux and the bubble area of the hydrophilic surface are less than those of the hydrophobic surface. With the contact angle changing from 77° to 120°, the heat transfer performance increases by 26%. Finally, the orthogonal test results indicate that, the wettability of the heat exchanger surface has the greatest influence on the flow boiling heat transfer performance, while the heater length has the least influence. © 2022 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Simulation Study on Dam Break Flow Based on the B-Spline Material Point Method北大核心CSCD

The dam break flow poses a common free surface flow problem in hydraulic engineering, and its ac⁃ curate simulation is of great engineering significance. The B⁃spline material point method (BSMPM), as an im⁃ proved algorithm of the material point method (MPM), has optimized accuracy and convergence in material point calculations and unique algorithmic advantages in free surface flow problems. Based on the BSMPM, a weakly compressible BSMPM (WC⁃BSMPM) was developed through introduction of an artificial equation of state. The simulation of the dam break flow problem was carried out, with the effects of the order of the B⁃ spline interpolation basis function on the simulation results analyzed. The results show that, the simulated fluid wavefront position, the wavefront velocity and the elevation variation at a given position are basically consistent with the existing experimental results. As the order of the basis function increases, the computation time will lengthen for about 1.5 times. However, the computation times of the BSMPM of different orders will uniformly increase approximately linearly with the background grid size. The validity of the WC⁃BSMPM simulation of the dam break flow problem was verified. The research provides a new idea and method for the simulation of dam break flow problems. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Influence of Multiple Micro Cracks on the Damage Behavior of a Macro-Crack Tip北大核心CSCD

The solution of an infinite plane containing a macro crack and a cluster of micro cracks under uniaxial tensile load was presented based on Muskhelishvili’s complex function method and the stepwise recursive method. The stress field and stress intensity factor K were obtained. Combined with the damage mechanics, damage parameter D of the macro-crack tip and the micro-crack tip under uniaxial tension was redefined, and the influence of different damage zone forms on the damage of the crack tip was analyzed. The results show that, both the chain-distribution and the reverse-chain-distribution micro cracks have an amplifying effect on the macro crack growth, and the damage parameter increases with the decrease of the inclination angle of the micro crack and the reduction of the distance between the macro crack and the micro cracks. For a relatively small inclination angle of the micro crack, the damage parameters of the macro crack and the micro crack heightens, and the damage parameter of the macro crack increases with the micro-crack length. For evenly distributed micro cracks in the continuous damage zone, the micro cracks have an amplifying effect on the macro-crack growth, and the damage parameter of the macro crack increases with the micro-crack number. © 2022 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Numerical Simulation of Transient Non-Isothermal Viscoelastic Couette Flows Based on the SPH Method北大核心CSCD

Based on the smoothed particle hydrodynamics (SPH) method, the transient non-isothermal viscoelastic flows were numerically simulated. First, the viscoelastic Couette flow based on the Oldroyd-B model under isothermal condition was simulated. Then, the simulation was extended to the non-isothermal case, in which the Reynolds exponential model was adopted to evaluate the dependence of the viscosity and the relaxation time on the temperature. The accuracy and effectiveness of the SPH method for simulating transient non-isothermal viscoelastic flows were verified through comparison with the finite volume method and evaluation of numerical convergence. The different flow characteristics of the non-isothermal flow compared with those of the isothermal flow were discussed. The effects of the temperature dependence coefficient and the Péclet number on the flow physics were analyzed. The numerical results demonstrate that, the SPH method can accurately and effectively simulate transient non-isothermal viscoelastic flow problems. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Study on Hydrodynamics Characteristics of a Single Bubble in Viscoelastic Fluid at Low Weissenberg Numbers北大核心CSCD

The VOF method was used to numerically study the upward motion of a single bubble in viscoelastic fluid, and the Oldroyd-B model was applied to describe the fluid viscoelastic property. At low Weissenberg numbers (Wi≤1), the effects of the viscous force, the relaxation time, the surface tension and the viscosity ratio on the rising motion of the bubble were studied. The results show that, under relatively large viscous and elastic forces (such as Ga= 2, Wi≥0.5 and β = 0.2), the bubble exhibits the phenomenon of “a pointed rear end”, and this phenomenon intensifies with the increase of the elasticity and the decrease of the surface tension. Otherwise, under a relatively weak elasticity (such as Wi= 0.1), the phenomenon of “a pointed rear end” disappears, and the bubble bears a hat-like shape. For a large surface tension (such as Eo = 1), the bubble bears a longitudinally elongated ellipse-like shape without distinct tail features. The effect of the surface tension on the bubble in viscoelastic fluid is like that in viscous fluid. The bubble has 2 types of rising motions, namely, “continuous acceleration” to a stable velocity and “acceleration-deceleration-reacceleration” to a stable velocity, and the bubble rising velocity in viscoelastic fluid is higher than that in pure viscous fluid. The elastic stress around the bubble is influenced by the viscosity and the relaxation time of the fluid, and with the decrease of the fluid viscosity or/ and the increase of the relaxation time, the incidence of the elastic stress becomes wide. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

A Problem on the Spherical Maximal Functions

The classical Lebesgue’s theorem shows that for any locally integrable functionf defined on R~n where B(x;ε)are n-dimensional balls with center x and radius ε. A extendedproblem of the above is,does     

A Numerical Calculation Algorithm of the Inclination Function Based on the Representation of SO(3) Group

In the determination of the Earth gravity field in satellite geodesy, the inclination functions represent the projection of data observed along the orbital plane of a satellite orbit into the sphere in the terrestial reference frame. The inclination functions in this work is studied from a group theoretical perspective. The inclination functions are proved to generate a representation of the SO(3) group. An orthogonal relation of the inclination functions is derived and some recurrence relations...     

Asymptotic of the Solutions of the Initial Boundary Value Problem for the Diffusion Equations for Semiconductors（Ⅱ）

The paper deal with the asymptotic behavior of the solutions to the initial boundary value problem for unipolar drift diffusion equations for semiconductors. Under the proper assumptions on doping profile and initial value, we prove that the smooth solutions to these evolutionary problems tend to the unique stationary solution exponentially as time tends to infinity.     

连续综合控制系统的状态反馈广义H_(2)控制北大核心CSCD

基于Lyapunov稳定性理论、矩阵分析法、线性矩阵不等式等方法,对同时带有控制输入和干扰输入的奇异摄动时变时滞不确定控制系统进行广义H_(2)控制研究.设计一个记忆状态广义H_(2)控制器,给出具体设计方法的判定定理.并对时滞依赖和时滞独立两种情形下采用新的引理,推出保守性相对更小的稳定性判据.对所得结论进行线性化处理,用数值样例验证了该文所得结论的有效性和可行性.指出在零到奇异摄动上界的整个区间范围内,闭环系统渐近稳定,扩大了广义H_(2)稳定空间,缩小了L_(2)-L_(∞)的性能指标.通过与相关文献进行稳定态指标对比,展示出该文所得方法具有一定的优越性和较小的保守性,并且适用于标准和非标准情形.     

单个碳粒燃烧的活化能近似逼近解

本文把活化能近似逼近方法应用于单个碳粒在静止高温气体中的燃烧问题。对非均相反应2C+O₂→2CO和C+CO₂→2CO及处于平衡态下的气相反应2CO+O2?2CO₂都加以了考虑。研究发现碳粒的燃烧过程主要在扩散极限阶段完成,但在加热阶段与扩散极限阶段之间有四个短暂但相当复杂的阶段。本文还举例对活化能近似逼近解法与精确的数值解法进行了比较,比较结果令人满意。     

Superconductivity and Ferromagnetism in Nonadiabatic Systems with Magnetic Impurity: A Step Beyond the Migdal Theorem

The possibility of the ferromagnetic ordering of a paramagnetic impurity in nonadiabatic superconducting systems is investigated. The effect of the relative shift of the Fermi surface by the internal magnetic field, the exchange interaction of the impurity ions through the conductivity electrons, and the spin–orbit interaction of the nonmagnetic impurity are taken into account. The problem is solved in the linear approximation with respect to the nonadiabaticity by taking the vertex and crossing diagrams corresponding to the electron–phonon and the electron–impurity interactions into account. We obtain basic equations of the superconductivity theory for nonadiabatic systems with the ferromagnetic ordering of the impurity spins and show that the nonadiabaticity alters the superconducting transition temperature T _c and the critical impurity concentration. The behavior of the magnetic-ordering temperature T _C as a function of the impurity concentration c in the superconductive state of the nonadiabatic system is also investigated. We obtain the phase diagram (T,c) and show that the nonadiabaticity effects lead to the enlargement of the domain where the superconductivity and the ferromagnetism exist simultaneously.     

双线性Fourier乘子交换子的有界性与紧性

假定T_σ是关于乘子σ的双线性Fourier乘子算子,其中σ满足如下Sobolev正则条件:对某个s∈(n,2n],有sup_(κ∈Z)‖σ_k‖W~s(R~(2m))∞.对于p_1,p_2,p∈(1,∞)且满足1/p=1/p_1+1/p_2和ω=(ω_1,ω_2)∈A_(p/t)(R~(2n)),建立了T_σ及其与函数b=(b_1,b_2)∈(BMO(R~n))~2生成的交换子T_(σ,b)由L~(p_1,λ)(ω_1)×L~(p_2,λ)(ω_2)到L~(p,λ)(v_w)的有界性;同时,在b_1,b_2∈CMO(R~n)(C_c~∞(R~n)在BMO拓扑下的闭包)的条件下,证明交换子T_(σ,b)是L~(p_1,λ)(ω_1)×L~(p_2,λ)(ω_2)到L~(p,λ)(v_w)的紧算子.为了得到主要结果,我们先后建立了几个双(次)线性极大函数在加多权Morrey空间上的有界性以及该空间中准紧集的判定.     

Complexity analysis of an interior point algorithm for the semidefinite optimization based on a kernel function with a double barrier term

In this paper, we establish the polynomial complexity of a primal-dual path-following interior point algorithm for solving semidefinite optimization(SDO) problems. The proposed algorithm is based on a new kernel function which differs from the existing kernel functions in which it has a double barrier term. With this function we define a new search direction and also a new proximity function for analyzing its complexity. We show that if q1 q2 1, the algorithm has O((q1 + 1) nq1+1/2(q1-q2)logn/ε)and O((q1 + 1)3q1-2q2+1/2(q1-q2)n~1/2 logn/ε) complexity results for large- and small-update methods, respectively.     

拟协调元空间的紧致性和拟协调元法的收敛性

本文首先讨论拟协调元空间的紧致性,把Rellich紧致定理推广到拟协调元空间序列,进而把广义Poincare、Friedrichs和Poincare-Friedrichs等不等式推广到拟协调元空间.然后讨论拟协调元法的收敛性和误差估计.本文证明了如果拟协调元空间具有逼近性和强连续性、满足单元秩条件且通过检验IPT,则近似解是收敛的.做为例子,我们证明了6参、9参、12参、15参、18参及21参拟协调元的收敛精度在L2,2(Ω)范数下分别是O(h_τ)、O(h_τ)、O(h_τ2)、O(h_τ2)、O(h_τ3)及O(h_τ4)量级.     

实际气体的幂律状态方程

为克服Onnes(昂内斯)气体状态方程参数多的缺点, 提出了描述实际气体的幂律状态方程.该模型仅包含两个参数,其中幂律函数的阶数可以为任意实数,刻画了实际气体偏离理想气体的程度.满足幂律状态方程的实际气体称为幂律气体.应用于描述氮气（N₂）和四氟甲烷（CF₄）两种实际气体的研究表明,与Onnes气体状态方程相比,幂律气体状态方程可以用较少的参数,准确地描述气体状态方程中压强和体积的幂律关系.此外,温度越低,幂律函数的阶数越小,反映了气体的实际状态越偏离理想气体.     

Coulomb expansion of a van der Waals C60 solid film

Scanning tunneling microscopy study revealed a van der Waals C60 solid film with 13% room-temperature lattice expansion on the GaAs(001) 2×4 surface. The mechanism involves fundamental Coulomb interaction due to charge transfer from the GaAs substrate. Theoretical calculation determines the charge transfer to be 1.76 electrons per C60 molecule. Oriented at its (110) crystallographic axis this film also distinguishes itself from those formed on all other semiconductor and metal substrates where only the low-energy (111) hexagonal packing of C60 molecules was developed. It is shown that this is due to the one-dimensional confinement effect of the anisotropic substrate, which may have the prospect of controlling crystal growth.     

Theoretical studies on solid heterogenous targets thickness in muon catalyzed fusion using transport theory and its comparison with available experimental values

Muon catalyzed fusion in solid heterogeneous H/D/T layers has been studied in recent years. When a negative muon enters the protium layer with a small admixture of tritium, the muonic atom is formed. molecular ion formation is highly dependent on the energy of the muonic tritium atom, μt, and this energy itself depends on the thickness of H/T and D₂ layers. In this paper, the transport theory is used to study the effect of the thickness of different layers (H/T and D₂) on the output kinetic energy of the μt atom. The calculated thickness by the used theoretical method are compared with the available experimental values and it is shown that they are in a good and acceptable agreement, and represents the reliability of used theoretical method.