Cascading multi-segment rupture process of the 2023 Turkish earthquake doublet on a complex fault system revealed by teleseismic P wave back projection method

In this study,the vertical components of broadband teleseismic P wave data recorded by China Earthquake Network are used to image the rupture processes of the February 6th,2023 Turkish earthquake doublet via back projection analysis.Data in two frequency bands(0.5-2 Hz and 1-3 Hz)are used in the imaging processes.The results show that the rupture of the first event extends about 200 km to the northeast and about 150 km to the southwest,lasting~90 s in total.The southwestern rupture is triggered by the northeastern rupture,demonstrating a sequential bidirectional unilateral rupture pattern.The rupture of the second event extends approximately 80 km in both northeast and west directions,lasting~35 s in total and demonstrates a typical bilateral rupture feature.The cascading ruptures on both sides also reflect the occurrence of selective rupture behaviors on bifurcated faults.In addition,we observe super-shear ruptures on certain fault sections with relatively straight fault structures and sparse aftershocks.

Projection and Contraction Methods for Nonlinear Complementarity Problem

We applied the projection and contraction method to nonlinear complementarity problem (NCP). Moveover, we proposed an inexact implicit method for (NCP) and proved the convergence.

Static output feedback stabilization for second-order singular systems using model reduction methods

In this paper,the static output feedback stabilization for large-scale unstable second-order singular systems is investigated.First,the upper bound of all unstable eigenvalues of second-order singular systems is derived.Then,by using the argument principle,a computable stability criterion is proposed to check the stability of secondorder singular systems.Furthermore,by applying model reduction methods to original systems,a static output feedback design algorithm for stabilizing second-order singular systems is presented.A simulation example is provided to illustrate the effectiveness of the design algorithm.

A class of fully third-order accurate projection methods for solving the incompressible Navier-Stokes equations

In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We then derive a sufficient condition for the continuous projection equations to be temporally third-order accurate approximations of the original Navier-Stokes equations by means of the localtruncation-error-analysis technique. The continuous projection equations are discretized temporally and spatially to third-order accuracy on the staggered grids, resulting in a fully third-order discrete projection scheme. The possibility to design higher-order projection methods is thus demonstrated in the present paper. A heuristic stability analysis is performed on this projection method showing the probability of its being stable. The stability of the present scheme is further verified through numerical tests. The third-order accuracy of the present projection method is validated by several numerical test cases.

PROJECTION METHODS AND APPROXIMATIONS FOR ORDINARY DIFFERENTIAL EQUATIONS

A formulation of a differential equation as projection and fixed point pi-Mem alloivs approximations using general piecnvise functions. We prone existence and uniqueness of the up proximate solution* convergence in the L2 norm and nodal supercnnvergence. These results generalize those obtained earlier by Hulme for continuous piecevjise polynomials and by Delfour-Dubeau for discontinuous pieceuiise polynomials. A duality relationship for the two types of approximations is also given.

Velocity Projection with Upwind Scheme Based on the Discontinuous Galerkin Methods for the Two Phase Flow Problem

The upwind scheme is very important in the numerical approximation of some problems such as the convection dominated problem, the two-phase flow problem, and so on. For the fractional flow formulation of the two-phase flow problem, the Penalty Discontinuous Galerkin (PDG) methods combined with the upwind scheme are usually used to solve the phase pressure equation. In this case, unless the upwind scheme is taken into consideration in the velocity reconstruction, the local mass balance cannot hold exactly. In this paper, we present a scheme of velocity reconstruction in some H(div) spaces with considering the upwind scheme totally. Furthermore, the different ways to calculate the nonlinear coefficients may have distinct and significant effects, which have been investigated by some authors. We propose a new algorithm to obtain a more effective and stable approximation of the coefficients under the consideration of the upwind scheme.

General Convergence Analysis for Three-step Projection Methods and Applications to Variational Problems

First a general model for a three-step projection method is introduced, and second it has been applied to the approximation solvability of a system of nonlinear variational inequality problems in a Hilbert space setting. Let H be a real Hilbert space and K be a nonempty closed convex subset of H. For arbitrarily chosen initial points x0, y0, z0 ∈ K, compute sequences xn, yn, zn such thatT ： K→ H is a nonlinear mapping onto K. At last three-step models are applied to some variational inequality problems.

Application for Superconvergence of Finite Element Approximations for the Elliptic Problem by Global and Local L²-Projection Methods

Numerical experiments are given to verify the theoretical results for superconvergence of the elliptic problem by global and local L2-Projection methods.

SOLVING SYMMETRIC MONOTONE LINEAR VARIATIONAL INEQUALITIES BY SOME MODIFIED LEVITIN-POLYAK PROJECTION METHODS

Abstract Some modified Levitin Polyak projection methods are proposed in this paper for solving monotone linear variational inequalityx∈Ω,(x′-x) T(Hx+c)≤0,\ x′∈Ω.It is pointed out that there are similar methods for solving a general linear variational inequality.

New algorithm for solving 3D incompressible viscous equations based on projection method

A new algorithm based on the projection method with the implicit finite difference technique was established to calculate the velocity fields and pressure.The calculation region can be divided into different regions according to Reynolds number.In the far-wall region,the thermal melt flow was calculated as Newtonian flow.In the near-wall region,the thermal melt flow was calculated as non-Newtonian flow.It was proved that the new algorithm based on the projection method with the implicit technique was correct through nonparametric statistics method and experiment.The simulation results show that the new algorithm based on the projection method with the implicit technique calculates more quickly than the solution algorithm-volume of fluid method using the explicit difference method.

基于Projection VOF方法的熔模铸造蜡模充型过程数值模拟

针对熔模铸造蜡模充型过程数值模拟软件自主开发,建立了流动场和温度场数学模型,采用Cross-WLF粘度模型描述蜡料的流变性能,基于Projection VOF方法进行蜡料流动过程速度场、压力场以及自由表面的求解。通过试验测量获取了K512型蜡料的热物理性能和流变性能的相关数据,拟合得到粘度模型参数。分别在蜡料常粘度和变粘度模型两种条件下对扳机铸件蜡模充型过程进行模拟,结果表明,所建立蜡模充型过程数学模型能较好地描述蜡料的动态流变特性和传热行为。

Self-adaptive one-dimensional nonlinear finite element method based on element energy projection method

The element energy projection （EEP） method for computation of super- convergent resulting in a one-dimensional finite element method （FEM） is successfully used to self-adaptive FEM analysis of various linear problems, based on which this paper presents a substantial extension of the whole set of technology to nonlinear problems. The main idea behind the technology transfer from linear analysis to nonlinear analysis is to use Newton＇s method to linearize nonlinear problems into a series of linear problems so that the EEP formulation and the corresponding adaptive strategy can be directly used without the need for specific super-convergence formulation for nonlinear FEM. As a re- sult, a unified and general self-adaptive algorithm for nonlinear FEM analysis is formed. The proposed algorithm is found to be able to produce satisfactory finite element results with accuracy satisfying the user-preset error tolerances by maximum norm anywhere on the mesh. Taking the nonlinear ordinary differential equation （ODE） of second-order as the model problem, this paper describes the related fundamental idea, the imple- mentation strategy, and the computational algorithm. Representative numerical exam- ples are given to show the efficiency, stability, versatility, and reliability of the proposed approach.

Low side lobe pattern synthesis using projection method with genetic algorithm for truncated cone conformal phased arrays

A hybrid method for synthesizing antenna＇s three dimensional （3D） pattern is proposed to obtain the low sidelobe feature of truncated cone conformal phased arrays. In this method, the elements of truncated cone conformal phased arrays are projected to the tangent plane in one generatrix of the truncated cone. Then two dimensional （2D） Chebyshev amplitude distribution optimization is respectively used in two mutual vertical directions of the tangent plane. According to the location of the elements, the excitation current amplitude distribution of each element on the conformal structure is derived reversely, then the excitation current amplitude is further optimized by using the genetic algorithm （GA）. A truncated cone problem with 8x8 elements on it, and a 3D pattern desired side lobe level （SLL） up to 35 dB, is studied. By using the hybrid method, the optimal goal is accomplished with acceptable CPU time, which indicates that this hybrid method for the low sidelobe synthesis is feasible.

ON THE MONOTONE CONVERGENCE OF THE PROJECTED ITERATION METHODS FOR LINEAR COMPLEMENTARITY PROBLEMS

Under suitable conditions,the monotone convergence about the projected iteration method for solving linear complementarity problem is proved and the influence of the involved parameter matrix on the convergence rate of this method is investigated.

Program evaluation and its application to equipment based on super-efficiency DEA and gray relation projection method

For the gray attributes of the equipment program and its difficulty to carry out the quantitative assessment of the equipment program information, the gray relation projection method is simply reviewed. Combining the super-data envelopment analysis（DEA） model and the gray system theory, a new super-DEA for measuring the weight is proposed, and a gray relation projection model is established to rank the equipment programs. Finally, this approach is used to evaluate the equipment program. The results are verified valid and can provide a new way for evaluating the equipment program.

Projected Runge-Kutta methods for constrained Hamiltonian systems

Projected Runge-Kutta （R-K） methods for constrained Hamiltonian systems are proposed. Dynamic equations of the systems, which are index-3 differential-algebraic equations （DAEs） in the Heisenberg form, are established under the framework of Lagrangian multipliers. R-K methods combined with the technique of projections are then used to solve the DAEs. The basic idea of projections is to eliminate the constraint violations at the position, velocity, and acceleration levels, and to preserve the total energy of constrained Hamiltonian systems by correcting variables of the position, velocity, acceleration, and energy. Numerical results confirm the validity and show the high precision of the proposed method in preserving three levels of constraints and total energy compared with results reported in the literature.

Superconvergence of nonconforming finite element penalty scheme for Stokes problem using L^2 projection method

A modified penalty scheme is discussed for solving the Stokes problem with the Crouzeix-Raviart type nonconforming linear triangular finite element. By the L^2 projection method, the superconvergence results for the velocity and pressure are obtained with a penalty parameter larger than that of the classical penalty scheme. The numerical experiments are carried out to confirm the theoretical results.

A Second-Order Method for the Electromagnetic Scattering from a Large Cavity

In this paper,we study the electromagnetic scattering from a two dimen- sional large rectangular open cavity embedded in an infinite ground plane,which is modelled by Helmholtz equations.By introducing nonlocal transparent boundary con- ditions,the problem in the open cavity is reduced to a bounded domain problem.A hypersingular integral operator and a weakly singular integral operator are involved in the TM and TE cases,respectively.A new second-order Toeplitz type approximation and a second-order finite difference scheme are proposed for approximating the hyper- singular integral operator on the aperture and the Helmholtz in the cavity,respectively. The existence and uniqueness of the numerical solution in the TE case are established for arbitrary wavenumbers.A fast algorithm for the second-order approximation is pro- posed for solving the cavity model with layered media.Numerical results show the second-order accuracy and efficiency of the fast algorithm.More important is that the algorithm is easy to implement as a preconditioner for cavity models with more general media.

A stable implicit nodal integration-based particle finite element method(N-PFEM)for modelling saturated soil dynamics

In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.

SELF-ADAPTIVE STRATEGY FOR ONE-DIMENSIONAL FINITE ELEMENT METHOD BASED ON ELEMENT ENERGY PROJECTION METHOD

Based on the newly-developed element energy projection （EEP） method for computation of super-convergent results in one-dimensional finite element method （FEM）, the task of self-adaptive FEM analysis was converted into the task of adaptive piecewise polynomial interpolation. As a result, a satisfactory FEM mesh can be obtained, and further FEM analysis on this mesh would immediately produce an FEM solution which usually satisfies the user specified error tolerance. Even though the error tolerance was not completely satisfied, one or two steps of further local refinements would be sufficient. This strategy was found to be very simple, rapid, cheap and efficient. Taking the elliptical ordinary differential equation of second order as the model problem, the fundamental idea, implementation strategy and detailed algorithm are described. Representative numerical examples are given to show the effectiveness and reliability of the proposed approach.