Search

The Proton Exchange Membrane Fuel Cell (PEMFC) performances are strongly impacted by the compression of the Gas Diffusion Layer (GDL). Despite its fibrous microstructure, this material is usually considered as a continuous medium and characterized with uniform loading. However, the GDL is subjected to a heterogeneous compression onto rib/channel patterns in the fuel cell assembly. In the present study, a complex behavior of the GDL response is experimentally revealed when the material is loaded with a rib/channel pattern, compared to uniform compression. The tests are simulated by finite element modeling using a classical strain-dependent elastic law, using parameters fitted from uniform compression experiments. It is shown that the numerical results do not reproduce the effect of pattern observed experimentally. Hypotheses to interpret these results involve mechanisms at the fiber microscale including fiber fracture, cross-link breakage and fiber rearrangement, which are exacerbated by larger material deformation caused by the heterogeneous loading. • Unexpected impact of rib/channel patterns on GDL compression. • Experimental investigation of the effect of rib/channel width on GDL deformation. • Numerical study with classical model unable to reproduce experimental results. • Interpretations based on mechanisms at the fiber microscale. [ABSTRACT FROM AUTHOR]

Paper, a porous-fibrous network, is made up of hydrophilic fibres, which are notably susceptible to deformations due to variations in moisture content. The response of paper sheets to full or partial wetting, based on the water-induced fibre swelling, has mostly been predicted in the past using highly idealised two-dimensional network geometries. The purpose of the present paper is to establish how precisely the two-stage process of determining effective hygro-expansivity under uniform wetting by a two-dimensional model followed by bending analysis using a continuum model predicts the curl of a fully three-dimensional network. In the present multi-scale modelling work, the moisture-induced homogenised sheet scale curl, derived from hygro-elastic deformations using non-uniform moisture profiles, is analysed based on an idealised non-woven three-dimensional fibrous network finite element model. First, the effective paper sheet swelling properties, derived computationally for the uniform wetting condition, are compared with analytically homogenised hygro-elastic properties. Whereas the analytically predicted elastic modulus has an inaccuracy of only ∼ 5%, the hygro-expansivity coefficient is over-predicted by roughly a third for a network with straight fibres. The difference is smaller if waviness of the fibres and wrap-around in the fibre bonds are taken into account. Fibre alignment turns out to have little influence. Subsequently, paper sheet curl predictions made by the analytical model are compared with that computed from the three-dimensional network model. The inaccuracy of the fully analytical model, based on the two-dimensional network description, relative to the fully detailed three-dimensional computational model with straight fibres is on the order of 40%–50%. The fidelity of the simple model is primarily constrained by the inaccuracy of the effective hygro-expansion coefficient. If the expansivity extracted from the three-dimensional simulations under uniform wetting is inserted, the inaccuracy drops to less than 10%. The continuum bending description hence performs adequately, provided it is based on a reliable estimation of the uniform expansion behaviour of the material. • Hygroscopic curl of paper sheets may be predicted using homogenised models. • We assess the accuracy of such predictions. • As a reference simulations of an idealised three-dimensional fibre network are used. • The predictions are accurate provided an accurate hygro-expansivity is used. • Two-dimensional network models may not be sufficient for this. [ABSTRACT FROM AUTHOR]

The stability, robustness, accuracy, and efficiency of space-time finite element methods crucially depend on the choice of approximation spaces for test and trial functions. This is especially true for high-order, mixed finite element methods which often must satisfy an inf-sup condition in order to ensure stability. With this in mind, the primary objective of this paper and a companion paper is to provide a wide range of explicitly stated, conforming, finite element spaces in four dimensions. In this paper, we construct explicit high-order conforming finite elements on 4-cubes (tesseracts); our construction uses tools from the recently developed 'Finite Element Exterior Calculus'. With a focus on practical implementation, we provide details including Piola-type transformations, and explicit expressions for the volumetric, facet, face, edge, and vertex degrees of freedom. In addition, we establish important theoretical properties, such as the exactness of the finite element sequences, and the unisolvence of the degrees of freedom. [ABSTRACT FROM AUTHOR]

This paper delves into the numerical analysis of the extended Fisher-Kolmogorov (EFK) equation within open bounded convex domains R ⊂ R d , where d ≤ 3. Two distinct finite element schemes are introduced, namely the semi-discrete and fully-discrete finite element approximations. The existence and uniqueness of solutions are established for both the semi-discrete and fully-discrete finite element approximations. Error bounds are investigated across different scenarios, including comparisons between the semi-discrete and exact solutions, as well as between the semi-discrete and fully-discrete solutions, along with the fully-discrete and exact solutions. An effective algorithm has been proposed to solve the nonlinear system resulting from the fully-discrete finite element approximation at each time step. The paper also provides computed numerical error results and showcases a variety of numerical experiments to further illustrate and support the findings. [ABSTRACT FROM AUTHOR]

A novel higher-order multi-scale computational approach for efficient nonlocal gradient elasticity simulation of composite materials is proposed in this paper. Based on the use of decoupling technique, raw fourth-order nonlocal gradient equations with periodically oscillatory coefficients are decoupled as the new second-order multi-scale equations for reducing the continuity requirements when implementing multi-scale simulation. Next, a novel macro-micro coupled multi-scale computational model with higher-order corrected terms is rigorously devised for high-resolution multi-scale and nonlocal analysis of composite materials via multi-scale asymptotic analysis. Then, the local error analysis of higher-order multi-scale solutions in the point-wise sense illustrates that establishing higher-order asymptotic corrected terms plays a vital role in investigating the micromechanical mechanisms. Moreover, a rigorous global error estimation with an explicit rate of higher-order multi-scale solutions is first derived in the energy norm sense. In addition, an efficient multi-scale numerical algorithm is presented to effectively simulate nonlocal gradient elasticity problems of composite materials based on finite element method (FEM). And we also derive the convergence estimation of the proposed numerical algorithm. Finally, numerical experiments are conducted to gauge the efficiency and accuracy of the proposed higher-order multi-scale method. This paper offers a unified higher-order two-scale computational framework for enabling nonlocal gradient elasticity behavior simulation of composite materials. [ABSTRACT FROM AUTHOR]

This paper designs, analyzes, and implements a finite element analysis for the simulation of the nonlinear advection-diffusion model. The key feature of the proposed work is to provide novel a priori error estimates in Bôchner norms L 2 (0 , T ; H 0 1 (Ω)) and L ∞ (0 , T ; H 0 1 (Ω)) simultaneously. To achieve the desired error estimates for both semidiscrete and fully discrete schemes, our idea is to introduce the Ritz projection operator that allows us to enforce the known estimates for P 1 and P 2 finite elements in the formulation of the desired results. Further, this paper establishes the stability, existence, and uniqueness in the precise normed spaces. Finally, we prove the fully discrete error estimates using the Backward Euler and Crank-Nicolson's schemes which are corroborated by numerical experiments to validate the developed theory on various domains for both schemes. To the best of author's knowledge, this is the first finite element simulation that provides the estimates in the L 2 Bôchner norm for any such types of time-dependent problems. [ABSTRACT FROM AUTHOR]

Grinding damage of hard and brittle materials has a direct impact on part strength and subsequent polishing efficiency. Damage suppression is one of the key research issues in hard and brittle material processing. Due to the uncertain damage evolution law, there is still a lack of optimization methods for the multi-step grinding process under given damage constraints. In this paper, the extended finite element method is used to establish a model between the prefabricated crack depth and final crack depth, which reveals the influence of various existing subsurface damages on crack propagation. It is found that only when the prefabricated crack depth is close to the indentation crack depth on the perfect surface, the final crack depth will be affected to some extent. To optimize the grinding process to achieve low damage, a semi-empirical model for predicting grinding damage is proposed which is suitable for process optimization and does not rely on some difficult-to-obtain parameters. Finally, based on the damage evolution law and prediction model, the core content of this paper is proposed, that is, a multi-step grinding process optimization method with variable grinding depth. This method can minimize the damage while ensuring the same processing efficiency, the damage depth is reduced by 25 %. The effectiveness of the optimization method and the correctness of the damage evolution law are proved by experiments. This method is suitable for engineering applications and has practical significance for understanding the damage evolution mechanism, which is helpful for the suppression of grinding damage and the optimization of the grinding process. • Reveal the evolution mechanism of various existing subsurface damage on crack propagation by XFEM. • Propose a semi-empirical model for grinding damage prediction suitable for process optimization. • Propose a multi-step grinding process optimization method for hard brittle materials with variable grinding depth. • Prove the effectiveness of the optimization method and the evolutionary model through experiments. [ABSTRACT FROM AUTHOR]

In this paper, further in-depth research on chip formation mechanism is conducted on the cutting of metallic glasses. The study of chips contributes to a better understanding of the machining nature of materials, which is of significant importance for improving product quality, reducing energy consumption, and resource utilization. This paper employs a combination of simulation and experimental methods, wihch not only analyzes the mechanism of chip formation in depth, but also investigates the time-varying characteristics of cutting forces. The results indicate that, when comparing simulation with experiments, the average error in cutting forces is 7.62%, and the chip morphology is also quite similar. Observations in the simulation show that as the cutting thickness increases from 7.96 μm to 23.87 μm, the free surface of the chip exhibits serrated features, and the main shear zone of the chip becomes clearer. Cutting forces exhibit periodic variations, with alternating peaks and valleys representing the formation of shear zones in the chip. In the experiment, observed chips displays a wrinkled free surface and a smooth non-free surface, with some chip edges exhibiting a burr phenomenon. By selecting appropriate machining parameters, the generation of secondary shear bands on the free surface was effectively reduced. By observing the chip morphology, it is found that the chip formation is accompanied by the material removal of plasticity and brittleness, which is the result of the combined action of heat, stress, and free volume. It should be considered to lower temperature to reduce the phenomenon of chip adhesion, providing guidance for further optimizing the cutting process of metallic glass. This paper can provide to readers further understanding of the chip formation mechanism in milling metal glass. • 2D finite element modeling was used to simulate the chip formation. • The shear zones appeared on the serrated chips. • The experimental results verified the accuracy of the simulation model. • The chips of BMG have smooth non-free surface and wrinkled free surface. • The serrated chips are the result of heat, stress, and free volume. [ABSTRACT FROM AUTHOR]

This paper summarizes the effective components of computer-aided design and computer-aided manufacturing (CAD-CAM) resin composites that contribute to achieving greater mechanical properties and further development. In silico multi-scale analysis, in silico nonlinear dynamic finite element analysis (FEA), and artificial intelligence (AI) were used to explore the effective components of CAD-CAM resin composites. The effects of the filler diameter and silane coupling ratio on the mechanical properties of CAD-CAM resin composites have been clarified through multi-scale analysis. The effects of the filler contents, and filler and monomer compositions have been investigated by AI algorithms. The fracture behavior of CAD-CAM composite crown was analyzed using in silico non-linear dynamic FEA. The longevity of CAD-CAM composite crown was assessed through step-stress accelerating life testing (SSALT). As the filler diameter decreases, there is an increase in elastic moduli and compressive strengths at the macroscale. At the nanoscale, a decrease in the filler diameter results in a decrease in the maximum value of the maximum principal strain. When the silane coupling ratio decreases, there is a decrease in the elastic modulus and compressive strength. According to the exhaustive search and feature importance analysis based on the AI algorithm, the combination of certain components was narrowed down to achieve a flexural strength of 269.5 MPa. The in silico non-linear FEA successfully detected the sign of the initial crack of the CAD-CAM composite molar crown. The SSALT revealed that CAD-CAM resin composite molar crowns containing nanofillers with a high fraction of resin matrix demonstrated great longevity. This paper summarized the effective components of CAD-CAM resin composites for their further development. The integration of in vitro and in silico approaches will expedite the advancement of CAD-CAM resin composites, offering benefits such as time efficiency and reduction of material waste for researchers and manufacturers. • The effective components of CAD-CAM resin composites (RC) have been summarized. • A smaller filler and a high silane coupling ratio enhance the mechanical properties of CAD-CAM RC. • Non-linear FEA has detected the sign of the initial crack in the CAD-CAM composite crown. • AI has been utilized to narrow down the effective components that achieve high flexural strength. • CAD-CAM composite crowns with a high fraction of resin matrix demonstrate remarkable longevity. [ABSTRACT FROM AUTHOR]

In order to solve the deficiencies of traditional board-level assembly CubeSat, such as low modularity, complex assembly process and poor scalability, a sandwich module assembly structure and method with no drive components is proposed, which can achieve quick connection and flexible assembly between modules. First, the proposed CubeSat standard module can be compatible with the standard interface of traditional CubeSat. The 4-in-1 interface of Mechanical-Electrical-Thermal-Data does not occupy the internal space of the standard module, and the electrical and thermal interfaces meet the requirements of CubeSat in-orbit work. Secondly, in order to improve the reliability of the connection between the modules, a plug-in connection scheme of the mechanical interface was designed. Finally, a genetic algorithm was used to optimize the interface size, so that it has better plug-in performance. The simulation results show that the performance of the male interface of titanium alloy material is the best, and the error between the finite element analysis and genetic algorithm results is less than 4.8 %, which is consistent with the actual tensile test results of the optimized prototype, and meets the mechanical requirements of the launch process and space orbit operation for CubeSat. The test results can guide production and provide a reference for subsequent engineering applications of CubeSat modular assembly. • In this paper, a standardized modular assembly method of CubeSat sandwich structure is designed. The assembly structure does not include driving parts, is simple in form, can be quickly connected, is flexible in assembly mode, and does not occupy the internal space of the standard module. • A standard 4-in-1 connection device was designed in this paper, including four interfaces: mechanical, electrical, thermal and data, which can achieve "plug and play". • This paper uses genetic algorithm to optimize the structure of the mechanical connection male joint, and obtains the optimal configuration of the male joint. This ensures the reliability, firmness and stability of mechanical connections, and meets the mechanical requirements of CubeSats during launch and in orbit. [ABSTRACT FROM AUTHOR]

This paper is devoted to a study of a history-dependent mixed hemivariational-variational inequality arising in contact mechanics. The contact problem concerns the deformation of a viscoelastic body with long memory, subject to a general friction law on one part of the boundary and a frictionless Signorini condition on another part of the boundary. The solution existence and uniqueness of the history-dependent mixed hemivariational-variational inequality based on a Lagrange multiplier approach are proved. Then, a fully discrete scheme is introduced and studied. The trapezoidal rule is used to approximate the integral in the history-dependent operator. For the spatial discretization, the linear finite elements are used to discretize the displacement field, and dual basis functions are used in the approximation of the Lagrange multiplier. Optimal order error estimates are derived for the displacement and the Lagrange multiplier under appropriate solution regularity assumptions. Numerical results are presented to illustrate the theoretical prediction of the convergence orders. [ABSTRACT FROM AUTHOR]

• An exact matrix-form solution is proposed for analyzing multiple non-circular tunnels in layered media. • Generalized complex variable theory, conformal mapping and fast Fourier transform are used in derivation. • The proposed solution avoids complicated coefficient determinations and iterative algorithm. • Effects of elastic properties, tunnels' position, equivalent layer and the layer thickness surrounding tunnels are discussed. This paper presents a matrix-form complex variable method that can predict the stress and displacement field around multiple tunnels in layered media. The proposed matrix-form method avoids complicated coefficient determinations and iterative algorithm, which effectively improves the computation efficiency. Firstly, the expressions for the stress and displacement field are given under the generalized complex variable theory. The conformal mapping technique and fast Fourier transform algorithm are used to formulate the matrix-form expressions of boundary conditions. Considering the boundary conditions of the tunnels and soil interfaces, the global matrix solution for multiple tunnels in layered media are obtained by assembling the matrix solutions of all single layers. The present theory is verified by comparison with the existing literature, finite element analysis and field data in engineering practices, and several numerical examples are given to reveal the effects of elastic properties, tunnels' position, equivalent layer and the layer thickness surrounding tunnels. [ABSTRACT FROM AUTHOR]

• A hybrid method by integrating the advantages of frequency-wavenumber method and finite element method. • The entire process simulation from the source to the local terrain has been achieved. • Seismic response of mountains due to point dislocation source. • The crustal velocity structure and source parameters influence the seismic response. This paper proposes a hybrid method to analyze seismic responses from the source to complex terrain. The hybrid method is based on the domain reduction idea, dividing the overall model into global domain (including the source and propagation path) and local domain (including complex terrain). For the propagation of seismic wave fields in the global domain, a frequency-wavenumber method based on the exact dynamic stiffness matrix is adopted. For the local domain, the finite element method combined with the artificial boundary substructure method is used to analyze the seismic response of complex terrain. The advantage of this hybrid method is that the frequency-wavenumber method can accurately and efficiently calculate the three-component ground motion of the global large-scale model, while the finite element method has a variety of elements and constitutive models, which can finely simulate the local complex terrain. We take the Zigong Xishan Park site as the research object, and verify the correctness of the hybrid method. Based on this, we use the hybrid method to simulate the seismic response of the site, focusing on the propagation process of seismic waves on the mountain surface, the spatial distribution of peak ground displacement, and the effects of crustal velocity structure and source parameters on the seismic response of the mountain. The results show that the propagation of seismic waves and the distribution of peak ground displacement are closely related to changes in terrain, and the crustal velocity structure and source parameters have significant effects on the seismic response of the mountain. [ABSTRACT FROM AUTHOR]

• This manuscript firstly establishes a half vehicle-pavement-continuous beam bridge (VPCB) coupled system analytical model. • It is innovatively that the bridge and vehicle initial conditions are considered into analytical model. • The pavement road surface local deterioration is simulated by increasing the local roughness level of the road surface. At present, the highway bridge field is in a period of vigorous development, and the application of continuous beam bridges is becoming more extensive, so the related safety research is particularly important. However, there are few studies on the vehicle-bridge coupled interaction of continuous beam bridges, and the influence of pavement is not considered. In order to ensure that the theoretical research can better serve the engineering practice, it is urgent to establish a more refined and more complex vehicle-bridge coupled model. In this paper, a vehicle-pavement-continuous beam bridge (VPCB) coupled system model is established. A half-vehicle model is adopted, the pavement is simulated by a continuous and uniform spring damper, and the bridge is a three-span continuous beam bridge. The motion equation of VPCB system is established and deduced using D 'Alembert principle and modal superposition method. The dynamic responses of VPCB system are calculated by Newmark- β method. Here, the continuous beam bridge mode functions are fitted based on the modal vector of the continuous beam calculated by the finite element model. The correctness of the method is verified by existing literatures. The effects of bridge pavement, vehicle speed, bridge second span length, the pavement equivalent spring stiffness coefficient k p and damping coefficient c p , road roughness on VPCB system are studied. The results show that the pavement can reduce the bridge dynamic amplification factor (DAF) and the vehicle acceleration to improve the ride comfort. The dynamic responses of the vehicle and the bridge contain each other's natural frequencies. The responses of VPCB system increase with the vehicle velocity, the vehicle responses are more affected than the bridge responses. The pavement equivalent stiffness and damping coefficient has little effect on the system. The responses of vehicle and bridge increase with the increase of the bridge second span length. When the vehicle travels each span of the three-span continuous beam bridge, the time history curve of the pavement deformation under the tire resembles a sine wave. When the local roughness level of the road surface increases, the system responses at the deteriorated position suddenly increases, and the responses of the subsequent position will also be affected. The larger the roughness grade, the more obvious the responses increase. The effect on the system responses of the whole pavement damage is more significant than the local pavement damage. [ABSTRACT FROM AUTHOR]

In this paper, a novel linearized second-order energy-stable fully discrete scheme for the nonlinear Benjamin-Bona-Mahony-Burgers (BBMB) equation is introduced, along with a superconvergence analysis of the conforming finite element method (FEM). Through skillful decomposition of the nonlinear term u ∇ ⋅ u , a new linearized second-order fully discrete scheme is developed. This scheme requires fewer iterations and exhibits higher computational efficiency compared to the nonlinear approach. Furthermore, it maintains energy stability, enabling direct numerical solution boundedness in the H 1 -norm, which represents an improvement over the L ∞ -norm boundedness reported in previous literature. Secondly, relying on this boundedness and the high-precision results of the bilinear element, we derive the unconditional error estimates for superclose and global superconvergence without any restrictions on the time step size Δ t and spatial mesh size h. Finally, numerical examples are presented to validate the theoretical analysis. [ABSTRACT FROM AUTHOR]

The aim of this paper is to study the superconvergent behavior of a two-grid finite element method (FEM) with 2-step backward differential formula (BDF2) for nonlinear dispersive wave equation. By introducing an auxiliary variable q = u t for the original variable u , the problem is transformed into a parabolic system. With the help of the combination technique of the interpolation and Ritz projection, the superclose and superconvergent estimates for the above two variables in H 1 -norm are obtained under the lower regularity of the exact solution compared with the standard Galerkin FEM. Finally, some numerical results are provided to show the correctness of the theoretical predictions and good performance of the proposed method. [ABSTRACT FROM AUTHOR]

This paper proposes and analyzes a mixed variational formulation for the Navier-Stokes equations with variable viscosity that depends nonlinearly on the velocity gradient. Differently from previous works in which augmented terms are added to the formulation, here we employ a technique that had been previously applied to the stationary Boussinesq problem and the Navier-Stokes equations with constant viscosity. Firstly, a modified pseudostress tensor is introduced involving the diffusive and convective terms and the pressure. Secondly, by using the incompressibility condition, the pressure is eliminated, and the gradient of velocity is incorporated as an auxiliary unknown to handle the aforementioned nonlinearity. As a consequence, a Banach spaces-based formulation is obtained, which can be written as a perturbed twofold saddle point operator equation. We address the continuous and discrete solvabilities of this problem by linearizing the perturbation and employing a fixed-point approach along with a particular case of a known abstract theory. Given an integer ℓ ≥ 0 , feasible choices of finite element subspaces include discontinuous piecewise polynomials of degree ≤ ℓ for each entry of the velocity gradient, Raviart-Thomas spaces of order ℓ for the pseudostress, and discontinuous piecewise polynomials of degree ≤ ℓ for the velocity as well. Finally, optimal a priori error estimates are derived, and several numerical results confirming in general the theoretical rates of convergence, and illustrating the good performance of the scheme, are reported. [ABSTRACT FROM AUTHOR]

This paper is contributed to explore how a Crank-Nicolson weak Galerkin finite element method (WG-FEM) addresses the singularly perturbed unsteady convection-diffusion equation with a nonlinear reaction term in 2D. The problem and some asymptotic behavior results are given for the exact solution and its derivatives with the parameter ε. These results are essential for proving the uniform convergence of the proposed WG-FEM. A tensor product Shishkin mesh featuring piece-wise discontinuous bilinear polynomials is employed to handle the layers for uniform convergence at different parameter values of ε. An error estimate ‖ u N − I N u ‖ W G is presented, where I N u denotes the vertex-edge-cell interpolation of the solution u , and ‖ ⋅ ‖ W G denotes the Weak Galerkin norm. The optimal uniform convergence order is demonstrated for semi-discrete and fully discrete schemes. Various numerical experiments are conducted to validate the optimal order of convergence demonstrated by the proposed method. [ABSTRACT FROM AUTHOR]

In this paper, we first investigate the well-posedness and regularity of mild solution to a stochastic active fluids model driven by the additive noise. A fully-discrete scheme is proposed for solving the given model, which is based on the finite element method for spatial discretization and the backward Euler method for temporal discretization. By overcoming the difficulty of error analysis caused by the discrete Laplacian operator, we obtain the convergence results of the developed scheme. Finally, some numerical examples are provided to validate the theoretical results and we also simulate the motion states of the active fluids. [ABSTRACT FROM AUTHOR]

Nonlocal problems involving fourth-order terms pose several difficulties such as numerical discretization and its related convergences analysis. In this paper, the well-posedness of the extended Fisher-Kolmogorov equation with nonlocal diffusion is first analyzed using the Faedo-Galerkin technique and the classical compactness arguments. Moreover, we adopt a BDF2 scheme for time discretization and a mixed Galerkin scheme for spatial discretization. Then, we derive the optimal order convergence rates of the fully discrete system. Finally, some numerical simulations and convergence results are provided to confirm the theoretical results and the accuracy of the proposed scheme. [ABSTRACT FROM AUTHOR]

This article introduces a unified mixed finite element framework based on a saddle-point formulation that applies to time-dependent fourth order linear and nonlinear problems with clamped, simply supported, and Cahn-Hilliard type boundary conditions. The classical mixed formulations lead to large matrix systems that demand huge storage and computational time making the schemes expensive, especially for the time-dependent problems. The proposed scheme circumvents this by employing biorthogonal basis functions that lead to sparse and positive-definite systems. The article discusses a mixed finite element method for the biharmonic problem and the time-dependent linear and nonlinear versions of the extended Fisher-Kolmogorov equations equipped with the aforementioned boundary conditions. The wellposedness of the scheme is discussed and a priori error estimates are presented for the semi-discrete and fully discrete finite element schemes. The numerical experiments validate the theoretical estimates derived in the paper. [ABSTRACT FROM AUTHOR]

In this paper, we consider a finite element approximation of an optimal control problem governed by a time fractional wave equation. We first prove the well-posedness and regularity of the optimal control problem. Furthermore, we discuss the fully discrete scheme and analyze its stability and a priori error estimate. Numerical experiments are given to illustrate the theoretical findings. • The weak solution to the time fractional wave equation is studied under low regular data. • The well-posedness and regularity of the optimal control problem are proven. • Numerical schemes exist the first-order and second-order accuracy in time and space respectively. [ABSTRACT FROM AUTHOR]

Positioning stages are significant in precision technologies under higher demands of operation and manufacturing. The performance of the selected guiding mechanisms is a decisive factor affecting the positioning accuracy of the stage. To restrict the parasitic motions along the degrees of constraint during motion transmission, a piezo-actuated compliant positioning stage with novel F-shaped linear guiding mechanisms (FLGM) is developed in this paper to enable a precision motion. With a simple and compact structure, a single FLGM can generate an approximate straight-line motion and enhance the transverse stiffness. Four FLGMs are arranged symmetrically in parallel to theoretically eliminate the parasitic motions, bringing in the merits of a low coupling rate and a wide bandwidth in the working direction. A stack piezoelectric actuator embedded in a forward bridge-type displacement amplifier is utilized to drive the guided stage. Analytical modeling of its kinematic, static, and dynamic characteristics is mathematically established by utilizing the pseudo-rigid-body method, compliance matrix method, and Lagrange's equation. Finite element analysis is conducted and a fabricated prototype is experimentally tested, validating the veracity of the analytical models. Experimental results show that the parasitic motion is constrained by the guiding mechanism effectively with a coupling rate below 0.95%, and a high natural frequency of 856.9 Hz is achieved simultaneously. • A novel F-shaped linear guiding mechanism with a simple and compact structure is developed. • The static and dynamic modeling of the designed positioning stage are conducted using the compliance matrix method and Lagrange's equation. • The proposed stage possesses a high natural frequency of 856.9 Hz and the ability to constrain parasitic motions, i.e., the cross-axis coupling rate is below 0.95%. [ABSTRACT FROM AUTHOR]

The contact quality between geometric features significantly affects the assembly accuracy and dynamic performance of mechanical products. However, accurate and efficient contact problem solving is currently a common challenge in the fields of tolerance analysis and tooth contact analysis (TCA). In this paper, a tolerance analysis method for gear mechanisms considering form deviations and elastic-plastic contact behaviors is proposed. Specifically, tolerance modeling of systematic and random deviations on tooth surfaces is represented by skin model shapes (SMS). Elastic-plastic contact behaviors including contact deformation, status, and stress distribution are modeled and calculated by a modified boundary element method (BEM). The accuracy and efficiency of this BEM algorithm are verified by comparing with the finite element method (FEM), indicating that the BEM can greatly improve the computational efficiency while maintaining high solution accuracy. A case study for a pair of spur gears demonstrates the significant effects of contact behaviors on transmission performance. This research provides an effective method to introduce mechanical contact behaviors in tolerance analysis of gear mechanisms. In addition, the contact simulation results can also be adopted in the analysis of vibration, reliability, and other aspects. • An accurate and efficient tolerance analysis method for gear mechanisms is proposed. • Real tooth surface topography and elastic-plastic contact behaviors are considered. • Systematic and random tooth surface deviations are represented by skin model shapes. • Elastic-plastic contact behaviors are modeled by a modified boundary element method. [ABSTRACT FROM AUTHOR]

The electrochemical discharge machining (ECDM) is an advanced and promising technology for fabricating micro structures. This method is effectively applicable for non-conductive hard and brittle materials micromachining, such as glass and ceramics. However, the ceramics high melting point and intense heat exchange effects during microchannel fabrication lead to the easy dissipation of thermal energy and the difficulty in processing. To address this issue and expanded the ECDM processing range. This paper dedicated to the investigation of the material removal mechanism and process characteristics for the alumina ceramic surfaces microchannels preparation. Firstly, this study explained the difference of the ECDM critical applied voltage and gas film generation mechanism by analyzing the volt-ampere characteristic curve of 20 wt% NaOH electrolyte with a tool immersed depth of 2 mm and 10 mm respectively. Then, a thermal removal finite element model (FEM) for ECDM was developed. The heat cumulative and intermittent cooling effects applied on the alumina ceramic surface dominated the materiel removal. Additionally, it can be found that the microchannel machining accuracy and bottom surface quality depended on the Gaussian heat source characteristics, gas film formation mechanism, heat accumulation effect, and temperatures spatial and temporal distribution. And a series of single factor experiments for evaluating the parametric studies such as the effect of discharge frequency, applied voltage, duty ratio and tool feed rate on the fabricated microchannel performance are executed. Finally, the complex three-dimensional structures of cross microchannel grid and zigzag microchannel array with variable depth were successfully fabricated on aluminum oxide ceramic surface. Notably, the results highlighted and enriched the ECDM stability, consistency and applicability for alumina ceramic surface microchannel fabrication. [ABSTRACT FROM AUTHOR]

In this paper we use the weighted discontinuous Galerkin finite element method to solve the dual-porosity-Brinkman coupling model. The dual-porosity model describes fluid flow in a dual-porosity medium composed of matrix and microfractures, while the Brinkman equation describes fluid flow in hydraulic fractures or conduits containing "debris filling" or "proppant filling". The two models are coupled by the four physically effective interface conditions, which are constructed based on the fundamental properties of the dual-porosity model and the Brinkman-Darcy model. We provide a weighted discontinuous discrete scheme for the model and derive the optimal error estimates. Through numerical experiments, we verify the properties of the model and the advantages of the numerical method, such as the optimal convergence rate of the numerical solution and the simulation of the flow characteristics around hydraulic fractures and conduits. [ABSTRACT FROM AUTHOR]

• The compressive buckling analysis of cracked nanosheets is performed by XFEM, considering small-scale effect. • The contribution of this work is the study of buckling of nanosheets in the presence of a crack. • The studied variables include small-scaling parameter, crack length, nanosheet dimensions, and boundary supports. • The results show that the cracking phenomenon and the increase of its length intensify the nonlocal effect. The process used to produce nanosheets may cause damage such as cracks in them. These cracks can adversely affect the buckling capacity of nanosheets. The main contribution of this research endeavor is to extend the previous studies on undamaged nanosheets to the cracked ones concerning the size-dependent effect. On the nanoscale, the assumption of material continuity, which is the fundamental assumption in classical continuum mechanics, is not valid. In this case, nonlocal non-classical theories are employed to analyze nanosheets. The present paper examines the nonlocal buckling of cracked graphene sheets under uniaxial and biaxial loads. To this end, the governing equations are developed based on the first-order shear deformation theory and using Eringen's nonlocal elasticity theory. Moreover, cracked nanosheets with various boundary conditions are analyzed via the extended finite element method. The variables under study are the small scaling parameter, crack length, nanosheet dimensional parameters, and support conditions. The results indicate that while the small scaling parameter reduces the critical buckling capacity of nanosheets, the presence of a crack in the nanosheet and the increase in its length intensify the nonlocal effect. In addition, as the aspect ratio of the nanosheet increases from one, the nonlocal effect decreases gradually. The small-scale effect increases by constraining the edges of the cracked nanosheets. [ABSTRACT FROM AUTHOR]

• A fully coupled three-dimensional THMC model (DEHydrate) is developed for simulating hydrate dissociation. • The numerical results obtained from the proposed model are in good agreement with those from other simulations. • Only a small fraction of gas produced from radial hydrate dissociation is collected from the wellbore. • The model provides an effective tool to predict the mechanical behavior during hydrate dissociation. Due to high energy density and huge reserves, many trial projects worldwide have been conducted to exploit natural gas hydrates (NGH). However, most of them did not meet the commercialization requirements because of submarine hazards, low gas production and sand production. Therefore, it is important to investigate the NGH dissociation process. In this paper, a fully coupled multiphase, strongly nonlinear three-dimensional (3D) thermal–hydro–mechanical–chemical model (DEHydrate) is developed to consider the multiphysics behaviours including solid-liquid-gas multiphase flow, heat transfer, NGH phase transition and solid deformation during hydrate dissociation. The model is solved using a fully implicit finite element method. The gas-liquid flow and heat transfer are simulated by low-order elements while the solid deformation is calculated by high- or low-order elements. The behaviours of the multiphase seepage and geomechanics are fully coupled and the physical properties of the reservoir are updated iteratively based on changes in pressure, temperature and displacement. The DEHydrate simulator provides an effective tool to predict the evolution of reservoir mechanical behavior during hydrate dissociation, including transient pressure, temperature, displacement and stress in the reservoir, as well as multiphase saturation, hydrate dissociation rate/mass and gas production rate/mass. During the NGH dissociation, the initial stage exhibits a rapid dissociation rate of up to 206.23 kg/(d·m) for NGH layer with a length of 600 m and a thickness of 22 m, which subsequently decreases to 11.64 kg/(d·m). Only a small portion of the gas generated during NGH dissociation is collected at the wellbore (approximately 35.6 % after 100 days), while the majority remains trapped in the reservoir. In addition, throughout NGH dissociation, the reservoir matrix continues to move toward the wellbore in the horizontal direction. While in the vertical direction, the reservoir matrix below the wellbore moves upwards, with decreasing displacement due to low-pressure diffusion, while the reservoir matrix above the wellbore moves downwards, with increasing displacement. The largest displacements occur near the wellbore during the early stage of NGH dissociation, while in the late stage, the largest displacements occur above the reservoir. [ABSTRACT FROM AUTHOR]

In this paper, we introduce and analyze a lowest-order locking-free weak Galerkin (WG) finite element scheme for the grad-div formulation of linear elasticity problems. The scheme uses linear functions in the interior of mesh elements and constants on edges (2D) or faces (3D), respectively, to approximate the displacement. An H (d i v) -conforming displacement reconstruction operator is employed to modify test functions in the right-hand side of the discrete form, in order to eliminate the dependence of the L a m e ´ parameter λ in error estimates, i.e., making the scheme locking-free. The method works without requiring λ ‖ ∇ ⋅ u ‖ 1 to be bounded. We prove optimal error estimates, independent of λ , in both the H 1 -norm and the L 2 -norm. Numerical experiments validate that the method is effective and locking-free. [ABSTRACT FROM AUTHOR]

The present paper investigates static deformation, free vibration, and patch loading response of one-dimensional (1D) hexagonal quasicrystal (QC) plates by the symplectic method. The Hamiltonian system is established by utilizing the variational principle, and then the symplectic analytical solutions of the extended displacement and traction vectors are obtained with simply supported boundary conditions. Numerical examples illustrate the effect of stacking sequence on the bending deformation of multilayered plates made up of QC and crystal materials. Moreover, the effects of the patch region, frequency, and phonon-phason coupling elastic coefficients on the corresponding problems are investigated. The symplectic analytical solutions presented in the article can be further used as benchmarks for numerical methods, such as the finite integral transformation and finite element methods. • The Hamiltonian system for dynamic equations of 1D hexagonal QC is established using the variational principle. • Analytic solutions of 1D hexagonal QC are rationally derived by the symplectic method. • The static deformation, free vibration, and patch loading response of 1D hexagonal QCs plates are investigated. • Some numerical tests are illustrated to justify the symplectic method. [ABSTRACT FROM AUTHOR]

• An accurate algorithm is presented for transcranial magneto-acoustic-electrical stimulation. • Unstructured mesh used in mesh discretization facilitates the pre-processing. • Our algorithm possesses higher calculation accuracy and faster convergence rate. • The developed model can achieve superior computational efficiency. • The proposed method is insensitive to mesh distortion. In this paper, a numerical method for simulating transcranial magneto-acoustic-electrical stimulation (TMAES) is presented. The smoothed finite element method (SFEM) is applied to calculate the magnetic field and the Rayleigh integral (RI) is used for solving the ultrasonic field. Three-node triangular and four-node tetrahedral elements are set as the background cells for constructing the smoothing domains (SDs). A gradient smoothing technique (GST) is applied over each SD to obtain the magnetic flux density. The discretized system equations are obtained by using the generalized smoothed Galerkin weakform. To solve the sound pressure, the ultrasonic transducer is discretized into a series of micro-elements. The acoustic pressure in the problem domain is then equal to the superposition of that generated by the individual micro-element. Finally, the current density can be calculated by the magneto-acoustic coupling effect. Numerical examples demonstrate that SFEM can improve accuracy, convergence rate, and efficiency to 10, 1.8, and 4 times those of FEM, respectively. In the calculation of coupled fields, SFEM-RI also improves the accuracy by 2.5 times. Additionally, the presented method possesses the property of insensitivity to mesh distortion. [ABSTRACT FROM AUTHOR]

The transient Stokes magnetohydrodynamic (Stokes MHD) equations are thoroughly investigated in this paper using the variational subgrid multiscale stabilized finite element approach. In order to get the optimal order of convergence, appropriate stabilization parameter expressions have been derived using Fourier analysis. The theoretical convergence analysis of the subgrid algebraic stabilized numerical scheme is also presented in this study. The numerical studies on multiply-driven cavity flow establish the credibility of the variational algebraic subgrid multiscale stabilized finite element approach. Further, the theoretically determined convergence rate is numerically validated. For different Hartmann, Reynolds, and magnetic force inclination angle values, the flow fields are traced. [ABSTRACT FROM AUTHOR]

The characteristics of the piezo-stack and the flexible mechanism are both influenced by temperature, which induces the performance change of piezoelectric actuators in low-temperature. This study investigates how temperature affects the dynamics model's stiffness, damping, driving force, and friction. An axial heat conduction model based on Fourier's heat conduction law was established, and finite element analysis was carried out to study the temperature distribution of the piezo-stack. Moreover, finite element analysis was carried out at various temperatures to investigate the impact of temperature on the flexible mechanism's deformation. Based on the aforementioned study, a dynamics model with variable stiffness and damping is proposed in this paper, and the simulation results are verified. In the experiment, the effects of temperature, voltage, and frequency on the output performance of the piezoelectric actuator were investigated separately using impedance and amplitude-frequency characteristic curves to determine the resonant frequency and damping ratio of the flexible mechanism at low temperatures. Theory and experiment show that high stiffness and low piezoelectric constants at low temperatures reduce the piezoelectric actuator's output displacement, friction reduction will enhance the driving displacement of the piezoelectric actuator, and the effect of damping on output displacement is insignificant. This study offers significant theoretical value and guidance for the use of piezoelectric actuators in low-temperature situations. • Temperature distribution across a piezo stack via mechanical series and electrical parallel connections. • Calculate stiffness, damping, and piezo constants at low temperatures, integrate them with a dynamics model for simulation. • Assess the influence of dynamics model parameters on piezo actuators through theoretical simulations and experiments. • Clarify the effect of each parameter of the dynamics model on the performance of piezoelectric actuators by means of experimental and theoretical simulations. [ABSTRACT FROM AUTHOR]

• Dynamic modelling of soft tissue deformation based on the dynamic model decomposition. • No prior knowledge of the system governing equation is required. • Construct a reduced-order model by extracting the essential features of the full model. This paper studies a new approach to dynamic modelling of soft tissue nonlinear deformation behavior based on dynamic model decomposition. This approach derives the deformation relationship between two continuous time points from correlated system snapshots. Subsequently, it constructs a reduced-order system by extracting most relevant information via thin-singular value decomposition to capture dynamic characteristics of the full-order deformation system. Unlike the existing methods which mainly rely on the establishment of complex constitutive equations for governing soft tissue deformation, the proposed method conducts tissue deformation directly from measured samples without involving constitutive governing equations. Simulation results show that the proposed method not only achieves real-time performance, but also sustains similar accuracy as the nonlinear finite element method. [ABSTRACT FROM AUTHOR]

In this paper, we consider the linearized backward Euler finite element scheme for the Navier-Stokes-Poisson-Nernst-Planck (NSPNP) equations. We aim to derive the unconditional error estimates of all variables in l ∞ (L 2) and l ∞ (H 1) norms. Compared with some previous results, we remove the mesh restriction. In addition, we derive two important physical properties for the proposed numerical scheme. Finally, three numerical examples are implemented to verify the theoretical analysis. [ABSTRACT FROM AUTHOR]

This paper is concerned with a class of nonlinear multi-term time fractional wave equations. Based on the Galerkin finite element method (FEM) in space and the L 2- 1 σ formula in time, combined with the corresponding fast algorithm, an efficient fully discrete scheme is constructed. The unconditional optimal convergence property in H 1 -norm of the proposed scheme is discussed by utilizing the time-space splitting technique. Furthermore, the optimal L 2 -norm convergence and the H 1 -norm superconvergence results of the method are obtained. Finally, several illustrative numerical experiments are performed to verify the theoretical results. Compared with the existing works, studying the convergence of FEMs for the nonlinear fractional wave equation is a great challenging task, and there are significant differences in research approaches. [ABSTRACT FROM AUTHOR]

In this paper, we present a new error analysis of a class of charge-conservative finite element methods for stationary inductionless magnetohydrodynamics (MHD) equations. The methods use the standard inf-sup stable Mini/Taylor–Hood pairs to discretize the velocity and pressure, and the Raviart–Thomas face element for solving the current density. Due to the strong coupling of the system and the pollution of the lower-order Raviart–Thomas face approximation in analysis, the existing analysis is not optimal. In terms of a mixed Poisson projection and the corresponding estimate in negative norms, we establish new and optimal error estimates. In particular, we prove that the method with the lowest-order Raviart–Thomas face element and Mini element provides the optimal accuracy for the velocity in L 2 -norm, and the method with the lowest-order Raviart–Thomas face element and P 2 − P 1 Taylor-Hood element supplies the optimal accuracy for the velocity in H 1 -norm and the pressure in L 2 -norm. Furthermore, we propose a simple recovery technique to obtain a new numerical current density of one order higher accuracy by re-solving a mixed Poisson equation. Numerical results are provided to verify the theoretical analysis. • New and optimal error analysis of charge-conservative finite element methods are established. • A simple recovery technique is proposed to obtain a new current density with higher accuracy. • Numerical examples are provided to verify the theoretical findings. [ABSTRACT FROM AUTHOR]

We present an application of stable finite element (FE) approximations of convection-diffusion initial boundary value problems (IBVPs) using a weighted least squares FE method, the automatic variationally stable finite element (AVS-FE) method [1]. The transient convection-diffusion problem leads to issues in classical FE methods as the differential operator can be considered a singular perturbation in both space and time. The stability property of the AVS-FE method, allows us significant flexibility in the construction of FE approximations in both space and time. Thus, in this paper, we take two distinct approaches to the FE discretization of the convection-diffusion problem: i) considering a space-time approach in which the temporal discretization is established using finite elements, and i i) a method of lines approach in which we employ the AVS-FE method in space whereas the temporal domain is discretized using the generalized- α method. We also consider another space-time technique in which the temporal direction is partitioned, thereby leading to finite space-time "slices" in an attempt to reduce the computational cost of the space-time discretizations. We present numerical verifications for these approaches, including numerical asymptotic convergence studies highlighting optimal convergence properties. Furthermore, in the spirit of the discontinuous Petrov-Galerkin (DPG) method by Demkowicz and Gopalakrishnan [2–6] , the AVS-FE method also leads to readily available a posteriori error estimates through a Riesz representer of the residual of the AVS-FE approximations. Hence, the norm of the resulting local restrictions of these estimates serves as error indicators in both space and time for which we present multiple numerical verifications in mesh adaptive strategies. [ABSTRACT FROM AUTHOR]