亚洲男人的天堂2018av,欧美草比,久久久久久免费视频精选,国色天香在线看免费,久久久久亚洲av成人片仓井空

This paper devises a novel lowest-order conforming virtual element method (VEM) for planar linear elasticity with the pure displacement/traction boundary condition. The main trick is to view a generic polygon $K$ as a new one $\widetilde{K}$ with additional vertices consisting of interior points on edges of $K$, so that the discrete admissible space is taken as the $V_1$ type virtual element space related to the partition $\{\widetilde{K}\}$ instead of $\{K\}$. The method is shown to be uniformly convergent with the optimal rates both in $H^1$ and $L^2$ norms with respect to the Lam\'{e} constant $\lambda$. Numerical tests are presented to illustrate the good performance of the proposed VEM and confirm the theoretical results.

We present a four-field Virtual Element discretization for the time-dependent resistive Magnetohydrodynamics equations in three space dimensions, focusing on the semi-discrete formulation. The proposed method employs general polyhedral meshes and guarantees velocity and magnetic fields that are divergence free up to machine precision. We provide a full convergence analysis under suitable regularity assumptions, which is validated by some numerical tests.

In this work we present a novel bulk-surface virtual element method (BSVEM) for the numerical approximation of elliptic bulk-surface partial differential equations (BSPDEs) in three space dimensions. The BSVEM is based on the discretisation of the bulk domain into polyhedral elements with arbitrarily many faces. The polyhedral approximation of the bulk induces a polygonal approximation of the surface. Firstly, we present a geometric error analysis of bulk-surface polyhedral meshes independent of the numerical method. Then, we show that BSVEM has optimal second-order convergence in space, provided the exact solution is $H^{2+3/4}$ in the bulk and $H^2$ on the surface, where the additional $\frac{3}{4}$ is due to the combined effect of surface curvature and polyhedral elements close to the boundary. We show that general polyhedra can be exploited to reduce the computational time of the matrix assembly. To demonstrate optimal convergence results, a numerical example is presented on the unit sphere.

The goal of this manuscript is to present a partitioned Model Order Reduction method that is based on a semi--implicit projection scheme to solve multiphysics problems. We implement a Reduced Order Method based on a Proper Orthogonal Decomposition, with the aim of addressing both time--dependent and time--dependent, parametrized Fluid--Structure Interaction problems, where the fluid is incompressible and the structure is thick and two dimensional.

We investigate the role of microstructural bridging on the fracture toughness of composite materials. To achieve this, a new computational framework is presented that integrates phase field fracture and cohesive zone models to simulate fibre breakage, matrix cracking and fibre-matrix debonding. The composite microstructure is represented by an embedded cell at the vicinity of the crack tip, whilst the rest of the sample is modelled as an anisotropic elastic solid. The model is first validated against experimental data of transverse matrix cracking from single-notched three-point bending tests. Then, the model is extended to predict the influence of grain bridging, brick-and-mortar microstructure and 3D fibre bridging on crack growth resistance. The results show that these microstructures are very efficient in enhancing the fracture toughness via fibre-matrix debonding, fibre breakage and crack deflection. In particular, the 3D fibre bridging effect can increase the energy dissipated at failure by more than three orders of magnitude, relative to that of the bulk matrix; well in excess of the predictions obtained from the rule of mixtures. These results shed light on microscopic bridging mechanisms and provide a virtual tool for developing high fracture toughness composites.

We introduced the least-squares ReLU neural network (LSNN) method for solving the linear advection-reaction problem with discontinuous solution and showed that the method outperforms mesh-based numerical methods in terms of the number of degrees of freedom. This paper studies the LSNN method for scalar nonlinear hyperbolic conservation law. The method is a discretization of an equivalent least-squares (LS) formulation in the set of neural network functions with the ReLU activation function. Evaluation of the LS functional is done by using numerical integration and conservative finite volume scheme. Numerical results of some test problems show that the method is capable of approximating the discontinuous interface of the underlying problem automatically through the free breaking lines of the ReLU neural network. Moreover, the method does not exhibit the common Gibbs phenomena along the discontinuous interface.

Fracture produces new mesh fragments that introduce additional degrees of freedom in the system dynamics. Existing finite element method (FEM) based solutions suffer from an explosion in computational cost as the system matrix size increases. We solve this problem by presenting a graph-based FEM model for fracture simulation that is remeshing-free and easily scales to high-resolution meshes. Our algorithm models fracture on the graph induced in a volumetric mesh with tetrahedral elements. We relabel the edges of the graph using a computed damage variable to initialize and propagate fracture. We prove that non-linear, hyper-elastic strain energy is expressible entirely in terms of the edge lengths of the induced graph. This allows us to reformulate the system dynamics for the relabeled graph without changing the size of system dynamics matrix and thus prevents the computational cost from blowing up. The fractured surface has to be reconstructed explicitly only for visualization purposes. We simulate standard laboratory experiments from structural mechanics and compare the results with corresponding real-world experiments. We fracture objects made of a variety of brittle and ductile materials, and show that our technique offers stability and speed that is unmatched in current literature.

We consider a moving boundary problem with kinetic condition that describes the diffusion of solvent into rubber and study semi-discrete finite element approximations of the corresponding weak solutions. We report on both a priori and a posteriori error estimates for the mass concentration of the diffusants, and respectively, for the a priori unknown position of the moving boundary. Our working techniques include integral and energy-based estimates for a nonlinear parabolic problem posed in a transformed fixed domain combined with a suitable use of the interpolation-trace inequality to handle the interface terms. Numerical illustrations of our FEM approximations are within the experimental range and show good agreement with our theoretical investigation. This work is a preliminary investigation necessary before extending the current moving boundary modeling to account explicitly for the mechanics of hyperelastic rods to capture a directional swelling of the underlying elastomer.

In transient simulations of particulate Stokes flow, to accurately capture the interaction between the constituent particles and the confining wall, the discretization of the wall often needs to be locally refined in the region approached by the particles. Consequently, standard fast direct solvers lose their efficiency since the linear system changes at each time step. This manuscript presents a new computational approach that avoids this issue by pre-constructing a fast direct solver for the wall ahead of time, computing a low-rank factorization to capture the changes due to the refinement, and solving the problem on the refined discretization via a Woodbury formula. Numerical results illustrate the efficiency of the solver in accelerating particulate Stokes simulations.

This work proposes a novel model and numerical formulation for lubricated contact problems describing the mutual interaction between two deformable 3D solid bodies and an interposed fluid film. The solid bodies are consistently described based on nonlinear continuum mechanics allowing for finite deformations and arbitrary constitutive laws. The fluid film is modelled as a quasi-2D flow problem on the interface between the solids governed by the averaged Reynolds equation. The averaged Reynolds equation accounts for surface roughness utilizing spatially homogenized, effective fluid parameters and for cavitation through a positivity constraint imposed on the pressure field. In contrast to existing approaches, the proposed model accounts for the co-existence of frictional contact tractions and hydrodynamic fluid tractions at every local point on the contact surface of the interacting bodies and covers the entire range from boundary lubrication to mixed, elastohydrodynamic, and eventually to full film hydrodynamic lubrication in one unified modelling framework with smooth transition between these different regimes. Critically, the model relies on a recently proposed regularization scheme for the mechanical contact constraint combining the advantages of classical penalty and Lagrange multiplier approaches by expressing the mechanical contact pressure as a function of the effective gap between the solid bodies while at the same time limiting the minimal gap value occurring at the (theoretical) limit of infinitely high contact pressures. From a physical point of view, this approach can be considered as a model for the elastic deformation of surface asperities, with a bounded magnitude depending on the interacting solids' surface roughness. A consistent and accurate model behavior is demonstrated and validated by employing several challenging and practically relevant benchmark test cases.