Material point method simulation of hydro-mechanical behaviour in twophase porous geomaterials: A state-of-the-art review

The material point method(MPM)has been gaining increasing popularity as an appropriate approach to the solution of coupled hydro-mechanical problems involving large deformation.In this paper,we survey the current state-of-the-art in the MPM simulation of hydro-mechanical behaviour in two-phase porous geomaterials.The review covers the recent advances and developments in the MPM and their extensions to capture the coupled hydro-mechanical problems involving large deformations.The focus of this review is aiming at providing a clear picture of what has or has not been developed or implemented for simulating two-phase coupled large deformation problems,which will provide some direct reference for both practitioners and researchers.

On the “Onion Husk” Algorithm for Approximate Solution of the Traveling Salesman Problem

The paper describes some implementation aspects of an algorithm for approximate solution of the traveling salesman problem based on the construction of convex closed contours on the initial set of points (“cities”) and their subsequent combination into a closed path (the so-called contour algorithm or “onion husk” algorithm). A number of heuristics related to the different stages of the algorithm are considered, and various variants of the algorithm based on these heuristics are analyzed. Sets of randomly generated points of different sizes (from 4 to 90 and from 500 to 10,000) were used to test the algorithms. The numerical results obtained are compared with the results of two well-known combinatorial optimization algorithms, namely the algorithm based on the branch and bound method and the simulated annealing algorithm. .

Sign-Changing Solutions for Superlinear Kirchhoff Type Problem via the Nehari Method

In this paper, we consider a class of Kirchhoff type problem with superlinear nonlinearity. A sign-changing solution with exactly two nodal domains will be obtained by combining the Nehari method and an iterative technique.

Exploring the Application Effect of Flipped Classroom Combined with Problem-Based Learning Teaching Method in Clinical Skills Teaching of Standardized Training for Resident Doctors of Traditional Chinese Medicine

Objective: To explore the application effect of flipped classroom combined with problem-based learning teaching method in clinical skills teaching of standardized training for resident doctors of traditional Chinese Medicine. Methods: The study used the experimental control method. The study lasted from September to November 2022. The subjects of this study were 49 students of standardized training for resident doctors of traditional Chinese Medicine from grades 2020, 2021 and 2022 of Dazhou integrated TCM & Western Medicine Hospital. They were randomly divided into experiment group (25) and control group (24). The experiment group adopted flipped classroom combined with problem-based learning teaching method, and the control group adopted traditional teaching method. The teaching content was 4 basic clinical skill projects, including four diagnoses of traditional Chinese Medicine, cardiopulmonary resuscitation, dressing change procedure, acupuncture and massage. The evaluation method was carried out by comparing the students’ performance and a self-designed questionnaire was used to investigate the students’ evaluation of the teaching method. Results: The test scores of total scores in the experimental group (90.12 ± 5.89) were all higher than those in the control group (81.47 ± 7.96) (t = 4.53, P P Conclusions: The teaching process of the flipped classroom combined with problem-based learning teaching method is conducive to improving the efficiency of classroom teaching, cultivating students’ self-learning ability, and enhancing students’ willingness to learn.

An inverse analysis of fluid flow through granular media using differentiable lattice Boltzmann method

This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.

Highly Accurate Golden Section Search Algorithms and Fictitious Time Integration Method for Solving Nonlinear Eigenvalue Problems

This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.

Optimal Shape Factor and Fictitious Radius in the MQ-RBF:Solving Ill-Posed Laplacian Problems

To solve the Laplacian problems,we adopt a meshless method with the multiquadric radial basis function(MQRBF)as a basis whose center is distributed inside a circle with a fictitious radius.A maximal projection technique is developed to identify the optimal shape factor and fictitious radius by minimizing a merit function.A sample function is interpolated by theMQ-RBF to provide a trial coefficient vector to compute the merit function.We can quickly determine the optimal values of the parameters within a preferred rage using the golden section search algorithm.The novel method provides the optimal values of parameters and,hence,an optimal MQ-RBF;the performance of the method is validated in numerical examples.Moreover,nonharmonic problems are transformed to the Poisson equation endowed with a homogeneous boundary condition;this can overcome the problem of these problems being ill-posed.The optimal MQ-RBF is extremely accurate.We further propose a novel optimal polynomial method to solve the nonharmonic problems,which achieves high precision up to an order of 10^(−11).

Numerical Exploration of Asymmetrical Impact Dynamics: Unveiling Nonlinearities in Collision Problems and Resilience of Reinforced Concrete Structures

This study provides a comprehensive analysis of collision and impact problems’ numerical solutions, focusing ongeometric, contact, and material nonlinearities, all essential in solving large deformation problems during a collision.The initial discussion revolves around the stress and strain of large deformation during a collision, followedby explanations of the fundamental finite element solution method for addressing such issues. The hourglassmode’s control methods, such as single-point reduced integration and contact-collision algorithms are detailedand implemented within the finite element framework. The paper further investigates the dynamic responseand failure modes of Reinforced Concrete (RC) members under asymmetrical impact using a 3D discrete modelin ABAQUS that treats steel bars and concrete connections as bond slips. The model’s validity was confirmedthrough comparisons with the node-sharing algorithm and system energy relations. Experimental parameterswere varied, including the rigid hammer’s mass and initial velocity, concrete strength, and longitudinal and stirrupreinforcement ratios. Findings indicated that increased hammer mass and velocity escalated RC member damage,while increased reinforcement ratios improved impact resistance. Contrarily, increased concrete strength did notsignificantly reduce lateral displacement when considering strain rate effects. The study also explores materialnonlinearity, examining different materials’ responses to collision-induced forces and stresses, demonstratedthrough an elastic rod impact case study. The paper proposes a damage criterion based on the residual axialload-bearing capacity for assessing damage under the asymmetrical impact, showing a correlation betweendamage degree hammer mass and initial velocity. The results, validated through comparison with theoreticaland analytical solutions, verify the ABAQUS program’s accuracy and reliability in analyzing impact problems,offering valuable insights into collision and impact problems’ nonlinearities and practical strategies for enhancingRC structures’ resilience under dynamic stress.

Analysis of the Boundary Stability of a Diffusion-Reaction System on a Nanolayer

In this paper, the focus is on the boundary stability of a nanolayer in diffusion-reaction systems, taking into account a nonlinear boundary control condition. The authors focus on demonstrating the boundary stability of a nanolayer using the Lyapunov function approach, while making certain regularity assumptions and imposing appropriate control conditions. In addition, the stability analysis is extended to more complex systems by studying the limit problem with interface conditions using the epi-convergence approach. The results obtained in this article are then tested numerically to validate the theoretical conclusions.

A Priori Error Analysis for NCVEM Discretization of Elliptic Optimal Control Problem

In this paper, we propose the nonconforming virtual element method (NCVEM) discretization for the pointwise control constraint optimal control problem governed by elliptic equations. Based on the NCVEM approximation of state equation and the variational discretization of control variables, we construct a virtual element discrete scheme. For the state, adjoint state and control variable, we obtain the corresponding prior estimate in H1 and L2 norms. Finally, some numerical experiments are carried out to support the theoretical results.

A Combination of Residual Distribution and the Active Flux Formulations or a New Class of Schemes That Can Combine Several Writings of the Same Hyperbolic Problem:Application to the 1D Euler Equations

We show how to combine in a natural way(i.e.,without any test nor switch)the conservative and non-conservative formulations of an hyperbolic system that has a conservative form.This is inspired from two different classes of schemes:the residual distribution one(Abgrall in Commun Appl Math Comput 2(3):341–368,2020),and the active flux formulations(Eyman and Roe in 49th AIAA Aerospace Science Meeting,2011;Eyman in active flux.PhD thesis,University of Michigan,2013;Helzel et al.in J Sci Comput 80(3):35–61,2019;Barsukow in J Sci Comput 86(1):paper No.3,34,2021;Roe in J Sci Comput 73:1094–1114,2017).The solution is globally continuous,and as in the active flux method,described by a combination of point values and average values.Unlike the“classical”active flux methods,the meaning of the point-wise and cell average degrees of freedom is different,and hence follow different forms of PDEs;it is a conservative version of the cell average,and a possibly non-conservative one for the points.This new class of scheme is proved to satisfy a Lax-Wendroff-like theorem.We also develop a method to perform nonlinear stability.We illustrate the behaviour on several benchmarks,some quite challenging.

A two-parameter multiple shooting method and its application to the natural vibrations of non-prismatic multi-segment beams

This paper presents an enhanced version of the standard shooting method that enables problems with two unknown parameters to be solved.A novel approach is applied to the analysis of the natural vibrations of Euler-Bernoulli beams.The proposed algorithm,named as two-parameter multiple shooting method,is a new powerful numerical tool for calculating the natural frequencies and modes of multi-segment prismatic and non-prismatic beams with different boundary conditions.The impact of the axial force and additional point masses is also taken into account.Due to the fact that the method is based directly on the fourth-order ordinary differential equation,the structures do not have to be divided into many small elements to obtain an accurate enough solution,even though the geometry is very complex.To verify the proposed method,three different examples are considered,i.e.,a three-segment non-prismatic beam,a prismatic column subject to non-uniformly distributed compressive loads,and a two-segment beam with an additional point mass.Numerical analyses are carried out with the software MATHEMATICA.The results are compared with the solutions computed by the commercial finite element program SOFiSTiK.Good agreement is achieved,which confirms the correctness and high effectiveness of the formulated algorithm.

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.

Hilbert’s First Problem and the New Progress of Infinity Theory

In the 19th century, Cantor created the infinite cardinal number theory based on the “1-1 correspondence” principle. The continuum hypothesis is proposed under this theoretical framework. In 1900, Hilbert made it the first problem in his famous speech on mathematical problems, which shows the importance of this question. We know that the infinitesimal problem triggered the second mathematical crisis in the 17-18th centuries. The Infinity problem is no less important than the infinitesimal problem. In the 21st century, Sergeyev introduced the Grossone method from the principle of “whole is greater than part”, and created another ruler for measuring infinite sets. The discussion in this paper shows that, compared with the cardinal number method, the Grossone method enables infinity calculation to achieve a leap from qualitative calculation to quantitative calculation. According to Grossone theory, there is neither the largest infinity and infinitesimal, nor the smallest infinity and infinitesimal. Hilbert’s first problem was caused by the immaturity of the infinity theory.

Solving Multi-Objective Linear Programming Problem by Statistical Averaging Method with the Help of Fuzzy Programming Method

A multi-objective linear programming problem is made from fuzzy linear programming problem. It is due the fact that it is used fuzzy programming method during the solution. The Multi objective linear programming problem can be converted into the single objective function by various methods as Chandra Sen’s method, weighted sum method, ranking function method, statistical averaging method. In this paper, Chandra Sen’s method and statistical averaging method both are used here for making single objective function from multi-objective function. Two multi-objective programming problems are solved to verify the result. One is numerical example and the other is real life example. Then the problems are solved by ordinary simplex method and fuzzy programming method. It can be seen that fuzzy programming method gives better optimal values than the ordinary simplex method.

On the Application of Adomian Decomposition Method to Special Equations in Physical Sciences

The current manuscript makes use of the prominent iterative procedure, called the Adomian Decomposition Method (ADM), to tackle some important special differential equations. The equations of curiosity in this study are the singular equations that arise in many physical science applications. Thus, through the application of the ADM, a generalized recursive scheme was successfully derived and further utilized to obtain closed-form solutions for the models under consideration. The method is, indeed, fascinating as respective exact analytical solutions are accurately acquired with only a small number of iterations.

Variational regularization method of solving the Cauchy problem for Laplace's equation: Innovation of the Grad–Shafranov(GS) reconstruction

The simplified linear model of Grad-Shafranov （GS） reconstruction can be reformulated into an inverse boundary value problem of Laplace＇s equation. Therefore, in this paper we focus on the method of solving the inverse boundary value problem of Laplace＇s equation. In the first place, the variational regularization method is used to deal with the ill- posedness of the Cauchy problem for Laplace＇s equation. Then, the ＇L-Curve＇ principle is suggested to be adopted in choosing the optimal regularization parameter. Finally, a numerical experiment is implemented with a section of Neumann and Dirichlet boundary conditions with observation errors. The results well converge to the exact solution of the problem, which proves the efficiency and robustness of the proposed method. When the order of observation error δ is 10-1, the order of the approximate result error can reach 10-3.

Analysis and management strategies for the low voltage problems of rural power grid based on the distribution network flow calculation method

For a long time, because of the lack of investment capital and enough attentions, the overall constructions of rural power grid were far behind than the urban power grid in Chongqing Jiangbei Power Company. The low voltage problems were highlighted in the rural power grid due to the characteristics of rural power grid. Using the distribution network flow calculation method, we evaluated the low voltage problems of the rural power grid which belongs to Chongqing Jiangbei Power Company. In addition, we collected the data of distribution transformers in electricity consumption peak period. Some practical management strategies were proposed by the analysis and evaluation of potential and appeared low voltage problems.

Modeling of Diffusion Transport through Oral Biofilms with the Inverse Problem Method

Aim The purpose of this study was to develop a mathe-matical model to quantitatively describe the passive trans-port of macromolecules within dental biofilms. Methodology Fluorescently labeled dextrans with different molecular mass （3 kD,10 kD,40 kD,70 kD,2 000 kD） were used as a series of diffusion probes. Streptococcus mutans,Streptococcus sanguinis,Actinomyces naeslundii and Fusobacterium nucleatum were used as inocula for biofilm formation. The diffusion processes of different probes through the in vitro biofilm were recorded with a confocal laser microscope. Results Mathematical function of biofilm penetration was constructed on the basis of the inverse problem method. Based on this function,not only the relationship between average concentration of steady-state and molecule weights can be analyzed,but also that between penetrative time and molecule weights. Conclusion This can be used to predict the effective concentration and the penetrative time of anti-biofilm medicines that can diffuse through oral biofilm. Further-more,an improved model for large molecule is proposed by considering the exchange time at the upper boundary of the dental biofilm.

A CONTINUATION HOMOTOPY METHOD FOR THE INVERSE PROBLEM OF OPERATOR IDENTIFICATION AND ITS APPLICATION

A new widly convergent method for solving the problem of operator identification is illustrated. Numerical simulations are carried out to test the feasibility and to study the general characteristics of the technique without the real measurement data. This technique is a direct application of the continuation homo-topy method for solving nonlinear systems of equations. It is found that this method does give excellent results in solving the inverse problem of the elliptic differential equations.