Adaptive Refinement Strategy Research Articles

Fatigue crack growth in functionally graded materials using an adaptive phase field method with cycle jump scheme

Functionally graded materials provide versatility in adjusting the volume fractions of constituent materials to meet specific design requirements. However, this customization often introduces mode-mixity at the crack tip, posing challenges in predicting fracture under cyclic loading with discrete approaches and computationally expensive with conventional phase-field fracture models. To address these issues, this paper introduces an adaptive phase-field fracture formulation with cycle jump scheme to elegantly predict fatigue crack nucleation and growth in functionally graded materials. Within this framework, the effective properties at a point are estimated using the Mori–Tanaka homogenization scheme, while the crack growth due to cyclic load is captured by incorporating an additional fatigue degradation parameter. Moreover, the computational efficiency of the proposed framework is improved through an adaptive mesh refinement and explicit cycle jump scheme. The adaptive refinement scheme utilizes an error indicator derived from both the displacement solution and phase-field variable. The adaptive refinement scheme is integrated with efficient quadtree decomposition, which generates a hierarchical mesh structure. Hanging nodes resulting from the quadtree decomposition are efficiently handled using a polygonal finite element method. The proposed framework is validated against experimental and numerical results reported in the literature. Furthermore, we investigate the fatigue crack growth resistance across a broad range of material gradation directions, gaining valuable insights and identifying functionally graded materials with high fatigue resistance.

Lagrangian relaxation for continuous‐time optimal control of coupled hydrothermal power systems including storage capacity and a cascade of hydropower systems with time delays

AbstractThis work considers a short‐term, continuous time setting characterized by a coupled power supply system controlled exclusively by a single provider and comprising a cascade of hydropower systems (dams), fossil fuel power stations, and a storage capacity modeled by a single large battery. Cascaded hydropower generators introduce time‐delay effects in the state dynamics, which are modeled with differential equations, making it impossible to use classical dynamic programming. We address this issue by introducing a novel Lagrangian relaxation technique over continuous‐time constraints, constructing a nearly optimal policy efficiently. This approach yields a convex, nonsmooth optimization dual problem to recover the optimal Lagrangian multipliers, which is numerically solved using a limited memory bundle method. At each step of the dual optimization, we need to solve an optimization subproblem. Given the current values of the Lagrangian multipliers, the time delays are no longer active, and we can solve a corresponding nonlinear Hamilton–Jacobi–Bellman (HJB) Partial Differential Equation (PDE) for the optimization subproblem. The HJB PDE solver provides both the current value of the dual function and its subgradient, and is trivially parallelizable over the state space for each time step. To handle the infinite‐dimensional nature of the Lagrange multipliers, we design an adaptive refinement strategy to control the duality gap. Furthermore, we use a penalization technique for the constructed admissible primal solution to smooth the controls while achieving a sufficiently small duality gap. Numerical results based on the Uruguayan power system demonstrate the efficiency of the proposed mathematical models and numerical approach.

A phase field approach to the fracture simulation of quasi-brittle structures with initial state

To solve the simulation problem of crack propagation in solid structures with initial stress, this paper proposes a phase field approach that considers the initial state of the quasi-brittle structure. By using the initial equilibrium iteration method, the estimated initial state of the structure, including initial stress, is transformed into internal stress loads of the structure. Thus, the problem of non-convergence caused by inaccurate initial state estimation of the structure in engineering calculations was solved, and an improved subproblem staggered iteration method was adopted to simulate crack propagation of the structure under external loads. The universality and computational accuracy of this method for different phase field models were verified by using pre made notched thin plates for the expansion of mode I and II cracks in different initial states. Among them, the calculation error using the PF_CZM model is within 1%, and the calculation error of the AT2 model increases with the initial stress, but is less than 7%. At the same time, the crack growth problem of cracked structures under different uniform initial stresses is accurately simulated, which reflects the applicability of this method to actual engineering. In addition, the study also verified the universal applicability of the adaptive refinement strategy proposed based on the AT2 model for other phase field models.

Numerical study of bubble rise in plunging breaking waves

During the occurrence of plunging wave breaking, a substantial number of multi-scale bubbles are generated. These submerged bubbles persist for extended periods and contribute to the distinct acoustic and optical characteristics of the wake. In this study, we utilize high-fidelity simulations, combined with the adaptive refinement strategy to accurately track bubbles of multi-scales during the entire rising stage. Unlike previous studies, our emphasis is specifically on investigating the process of bubble rising during plunging wave breaking. Comprehensive statistical analyses are performed and characteristics of bubbles across various scales are also provided. Our findings reveal that most bubbles are concentrated in small scales, while larger bubbles rapidly ascend to the surface or undergo fragmentation into smaller bubbles through breaking cascades eventually. A distinct stratification of bubble size distribution along the depth direction is observed. Bubble velocity distributions are also important characteristics that are frequently neglected in studies of plunging wave breaking. Bubbles primarily spread along the spanwise direction, with a uniform distribution of velocity in this dimension. The velocity distribution of bubbles displays asymmetric tails that extend to higher velocities, and within this high-velocity regime, a power law behavior is observed, similar to the size distributions. Ultimately, the flow field is left with only a few small bubbles, moving at an exceedingly low speed. Furthermore, dynamical evolution of bubble rise in plunging wave breaking is described in detail and we analyze the intricate interactions between bubbles and turbulent flows. We observe that vortices are predominantly generated in close proximity to the bubbles, and bubble motion plays a crucial role in initiating turbulent flows. Simultaneously, these vortices contribute to the fragmentation of large-scale bubbles, transforming them into smaller counterparts due to turbulent fluctuations.

Dynamic crack growth in orthotropic brittle materials using an adaptive phase-field modeling with variable-node elements

In this paper, crack growth in orthotropic materials subjected to dynamic loading is numerically studied using an adaptive phase-field method. The study starts with a coarse structured mesh and the adaptive refinement strategy based on a user-defined threshold on the phase-field variable is proposed for computational efficiency, and variable-node elements are employed to treat the hanging nodes as a result of local adaptive refinement. The Hughes-Hilbert-Taylor (HHT) time integration scheme is adopted for temporal discretization. The directionality of orthotropic materials is represented by a penalized second-order structural matrix, which is incorporated in the crack face energy density. Through numerical examples, the influence of the material orientation on the dynamic crack growth in orthotropic materials is studied and the reliability of the proposed framework is validated.

Low-cycle fatigue crack growth in brittle materials: Adaptive phase-field modeling with variable-node elements

This work presents an adaptive phase field framework for the simulation of crack nucleation and propagation in brittle materials subject to low-cycle loading. We introduce a fatigue history strain parameter within the phase-field framework to capture the fatigue effect. The resulting coupled differential equations are solved using a staggered iteration scheme. In order to improve the computational efficiency of the phase-field model, an adaptive refinement scheme is introduced. This scheme utilizes an artificial threshold that takes into consideration of both phase-field variables and cumulative fatigue history field variables to determine the elements to be refined. The handling of the hanging nodes resulting from mesh refinement is accomplished through the application of the variable-node element technique, offering a flexible and efficient mesh integration technique. We validate our proposed numerical scheme through five representative numerical examples: single-edge cracked specimen, compact tension specimen, double-edge cracked panel, plate with a hole and plate with multi-cracks and hole. The results demonstrate that the adaptive mesh refinement scheme significantly reduces computational costs without compromising the accuracy of the numerical predictions.

Non‐parametric measure approximations for constrained multi‐objective optimisation under uncertainty

AbstractIn this article, we propose non‐parametric estimations of robustness and reliability measures approximation error, employed in the context of constrained multi‐objective optimisation under uncertainty (OUU). These approximations with tunable accuracy permit to capture the Pareto front in a parsimonious way, and can be exploited within an adaptive refinement strategy. First, we illustrate an efficient approach for obtaining joint representations of the robustness and reliability measures, allowing sharper discrimination of Pareto‐optimal designs. A specific surrogate model of these objectives and constraints is then proposed to accelerate the optimisation process. Secondly, we propose an adaptive refinement strategy, using these tunable accuracy approximations to drive the computational effort towards the computation of the optimal area. To this extent, an adapted Pareto dominance rule and Pareto optimal probability computation are formulated. The performance of the proposed strategy is assessed on several analytical test‐cases against classical approaches. We also illustrate the method on an engineering application, performing shape OUU of an Organic Rankine Cycle turbine.

Towards an error indicator-based [formula omitted]-adaptive refinement scheme in kinematic upper-bound limit analysis with the presence of seepage forces

This paper presents a new computational strategy for kinematic upper bound limit analysis in the presence of seepage forces with an improved mesh refinement scheme. In particular, the original adaptive refinement scheme is enhanced with a simple but efficient error-indicator of the nodal plastic dissipation for high-order elements. The proposed refinement criteria facilitate a precise quantitative evaluation of errors in the nodal plastic energy dissipation across each element and ensure that the rate of total refined element number gradually decreases with adaptive step. Adhering to two-dimensional steady state seepage condition, numerical details regarding the calculation of total water head distributions for the seepage field are provided. In a similar manner as treating the unit weight of the soil, the effects of seepage forces are incorporated as body forces in the upper bound formulation. Detailed numerical procedure of the proposed error indicator-based h-adaptive refinement scheme incorporating the inclusion of seepage forces is addressed and implemented in the in-house code. Several benchmark problems, including both homogeneous and non-homogeneous soils, are analyzed to evaluate the excellent performance of the proposed error indicator-based h-adaptive refinement scheme in upper-bound limit analysis and assess the influence of seepage forces on the stability of geostructures through comparison with that for dry condition.

Numerical Modeling of Pumping-Induced Earth Fissures Using Coupled Quasi-Static Material Point Method

Pumping-induced earth fissuring is a typical fluid–solid coupling problem involving fracture evolution. Despite recent advances in the numerical modeling of earth fissures, a holistic approach capable of predicting their evolution and understanding their behavior from fracture mechanisms is absent. This paper presents a coupled numerical method based on the framework of the material point method (MPM) and fracture mechanics for the process-oriented simulation of pumping-induced tensile earth fissures. In this method, fissure behavior is driven by the real-time stress field affected by soil consolidation. After an earth fissure initiates in the region of maximum tensile stress, it propagates in the direction of the maximum circumferential stress around the fissure tip. The great scale difference between a fissure tip and the entire model is addressed with an adaptive refinement strategy. Through the simulation of a physical model experiment reproducing pumping-induced earth fissures with the presence of a bedrock ridge, the complete fissuring process, including initiation, propagation, and arrest, is presented. The simulated fracture time (24 h after drainage), location (directly above the bedrock ridge), final length (165 mm), and morphology evolution (from straight to zigzag) of the main fissure basically were in agreement with the experimental results. The stress intensity factors (SIFs) around fissure tips, which are important fracture parameters that have not been discussed in existing earth fissure analyses, also are captured during the fissure growth. The proposed method quantifies the inhibition of geostatic stress fields in earth fissure development and correlate fissure morphology with the stress field at a fissure tip.

Enhanced Floating Isogeometric Analysis

The numerical simulation of additive manufacturing techniques promises the acceleration of costly experimental procedures to identify suitable process parameters. We recently proposed Floating Isogeometric Analysis (FLIGA), a new computational solid mechanics approach, which is mesh distortion-free in one characteristic spatial direction. FLIGA emanates from Isogeometric Analysis and its key novel aspect is the concept of deformation-dependent “floating” of individual B-spline basis functions along one parametric axis of the mesh. Our previous work showed that FLIGA not only overcomes the problem of mesh distortion associated to this direction, but is also ideally compatible with material point integration and enjoys a stability similar to that of conventional Lagrangian mesh-based methods. These features make the method applicable to the simulation of large deformation problems with history-dependent constitutive behavior, such as additive manufacturing based on polymer extrusion. In this work, we enhance the first version of FLIGA by (i) a novel quadrature scheme which further improves the robustness against mesh distortion, (ii) a procedure to automatically regulate floating of the basis functions (as opposed to the manual procedure of the first version), and (iii) an adaptive refinement strategy. We demonstrate the performance of enhanced FLIGA on relevant numerical examples including a selection of viscoelastic extrusion problems.

An adaptive finite element method for Riesz fractional partial integro-differential equations

AbstractThe Riesz fractional derivative has been employed to describe the spatial derivative in a variety of mathematical models. In this work, the accuracy of the finite element method (FEM) approximations to Riesz fractional derivative was enhanced by using adaptive refinement. This was accomplished by deducing the Riesz derivatives of the FEM bases to work on non-uniform meshes. We utilized these derivatives to recover the gradient in a space fractional partial integro-differential equation in the Riesz sense. The recovered gradient was used as an a posteriori error estimator to control the adaptive refinement scheme. The stability and the error estimate for the proposed scheme are introduced. The results of some numerical examples that we carried out illustrate the improvement in the performance of the adaptive technique.

Adaptive discretization refinement for discrete models of coupled mechanics and mass transport in concrete

An adaptive discretization refinement strategy for steady state discrete mesoscale models of coupled mechanics and mass transport in concrete is presented. Coupling is provided by two phenomena: the Biot’s theory of poromechanics and an effect of cracks on material permeability coefficient. The model kinematics is derived from rigid body motion of Voronoi cells obtained by tessellation of the domain. Starting with a coarse discretization, the density of Voronoi generator points is adaptively increased on the fly in regions where the maximum principal stress exceeds a chosen threshold. Purely elastic behavior is assumed in the coarse discretization, therefore no transfer of history/state variables is needed. Examples showing (i) computational time savings achieved via the adaptive technique and (ii) an agreement of the outputs from the fine and adaptive models during simulations of hydraulic fracturing and three-point bending combined with a fluid pressure loading are presented.

A continuous field adaptive mesh refinement algorithm for isogeometric topology optimization using PHT-Splines

This study proposes a continuous density-field based isogeometric topology optimization (ITO) using Polynomial splines over Hierarchical T-meshes (PHT-Splines). Taking the benefit of isogeometric methods, the design density values for the mesh are defined using spline basis functions as a continuous density function (CDF) through the control points. The material density at the control points will then appear as a facade over the domain, which is further employed to generate a very smooth topology and also to propose a marking criterion for an adaptive mesh refinement (AMR) strategy during the optimization process. During this adaptive meshing, the isogeometric domain mesh is refined during each adaptive step, properly tracking the density variation over the elements and sharpening the boundaries of the optimized topologies. By utilizing the rationale of the Geometry Independent Field approximaTion (GIFT) framework, any complex multi-patch geometry constructed in an industry-standard CAD package using Non-Uniform Rational B-Splines (NURBS) may be imported and discretized using PHT-Splines, including higher dimensions. This helps to maintain the geometrical accuracy through an initial coarse mesh and computational accuracy of the optimization iterations through the local refinement feature of the PHT-Splines. An adaptive first-neighbourhood smoothening algorithm based on the Shepard function is also proposed to obtain manufacturable topologies and to circumvent post-processing stages. The proposed method was observed to provide remarkably smooth topologies with an initial coarse mesh, thereby obtaining a computationally efficient domain to be solved. Compared to existing element density based methods, the proposed method demonstrates a substantial minimization in Degree’s-of-Freedom (DoF) requirement to arrive at clear and realistic solutions. Also, when compared to globally refined solutions, a 50%–90% reduction in DoF is achieved using the adaptive refinement strategy.

Adaptive phase-field modelling of fracture propagation in poroelastic media using the scaled boundary finite element method

A scaled boundary finite element-based phase field formulation is proposed to model two-dimensional fracture in saturated poroelastic media. The mechanical response of the poroelastic media is simulated following Biot’s theory, and the fracture surface evolution is modelled according to the phase field formulation. To avoid the application of fine uniform meshes that are constrained by the element size requirement when adopting phase field models, an adaptive refinement strategy based on quadtree meshes is adopted. The unique advantage of the scaled boundary finite element method is conducive to the application of quadtree adaptivity, as it can be directly formulated on quadtree meshes without the need for any special treatment of hanging nodes. Efficient computation is achieved by exploiting the unique patterns of the quadtree cells. An appropriate scaling is applied to the relevant matrices and vectors according the physical size of the cells in the mesh during the simulations. This avoids repetitive calculations of cells with the same configurations. The proposed model is validated using a benchmark with a known analytical solution. Numerical examples of hydraulic fractures driven by the injected fluid in cracks are modelled to illustrate the capabilities of the proposed model in handling crack propagation problems involving complex geometries.

Adaptive Nonintrusive Reconstruction of Solutions to High-Dimensional Parametric PDEs

Numerical methods for random parametric PDEs can greatly benefit from adaptive refinement schemes, in particular when functional approximations are computed as in stochastic Galerkin and stochastic collocations methods. This work is concerned with a nonintrusive generalization of the adaptive Galerkin finite element method with residual-based error estimation. It combines the nonintrusive character of a randomized least squares method with the a posteriori error analysis of stochastic Galerkin methods. The proposed approach uses the variational Monte Carlo method to obtain a quasi-optimal low-rank approximation of the Galerkin projection in a highly efficient hierarchical tensor format. We derive an adaptive refinement algorithm which is steered by a reliable error estimator. Opposite to stochastic Galerkin methods, the approach is easily applicable to a wide range of problems, enabling a fully automated adjustment of all discretization parameters. Benchmark examples with affine and (unbounded) lognormal coefficient fields illustrate the performance of the nonintrusive adaptive algorithm, showing the expected convergence rates of single-level strategies.

Cavity-Free Continuum Solvation: Implementation and Parametrization in a Multiwavelet Framework

We present a multiwavelet-based implementation of a quantum/classical polarizable continuum model. The solvent model uses a diffuse solute-solvent boundary and a position-dependent permittivity, lifting the sharp-boundary assumption underlying many existing continuum solvation models. We are able to include both surface and volume polarization effects in the quantum/classical coupling, with guaranteed precision, due to the adaptive refinement strategies of our multiwavelet implementation. The model can account for complex solvent environments and does not need a posteriori corrections for volume polarization effects. We validate our results against a sharp-boundary continuum model and find a very good correlation of the polarization energies computed for the Minnesota solvation database.

Modeling and High-performance Trajectory Optimization of the Industrial Robot

Rapidity, accuracy and energy cost are the key performance indexes for of the work of industrial robot. In this paper, the shortest time and the minimum energy consumption are taken as the goal for the comprehensive optimization. The control vector parameterization method is used to reasonably refine the time grid of the optimization problem. Then, an efficient control strategy for non-uniform adaptive mesh refinement is proposed. According to different changing trends, time nodes are merged or inserted into the time grid. While reducing the number of discrete grids, it also ensures the accuracy of the solution. Taking the Manutec R3 industrial robot as the object, the trajectory optimization simulation calculation is carried out. According to the verification of the simulation computing, the proposed method can significantly improve the time grid refinement. Compared with the traditional CVP method, not only the calculation time is shortened, but also more accurate results can be obtained.

Adaptive interface-Mesh un-Refinement (AiMuR) based sharp-interface level-set-method for two-phase flow

Adaptive interface-Mesh un-Refinement (AiMuR) based Sharp-Interface Level-Set-Method (SI-LSM) is proposed for both uniform and non-uniform Cartesian-Grid. The AiMuR involves interface location based dynamic un-refinement (with merging of the four control volumes) of the Cartesian grid away from the interface. The un-refinement is proposed for the interface solver only. A detailed numerical methodology is presented for the AiMuR and ghost-fluid method based SI-LSM. Advantage of the novel as compared to the traditional SI-LSM is demonstrated with a detailed qualitative as well as quantitative performance study, involving the SI-LSMs on both coarse grid and fine grid, for three sufficiently different two-phase flow problems: dam break, breakup of a liquid jet and drop coalescence. A superior performance of AiMuR based SI-LSM is demonstrated - the AiMuR on a coarser non-uniform grid ( $$NU_{c}^{AiMuR}$$ ) is almost as accurate as the traditional SI-LSM on a uniform fine grid ( $$U_{f}$$ ) and takes a computational time almost same as that by the traditional SI-LSM on a uniform coarse grid ( $$U_{c}$$ ). The AMuR is different from the existing Adaptive Mesh Refinement (AMR) as the former involves only mesh un-refinement while the later involves both refinement and un-refinement of the mesh. Moreover, the proposed computational development is significant since the present adaptive un-refinement strategy is much simpler to implement as compared to that for the commonly used adaptive refinement strategies. The proposed numerical development can be extended to various other multi-physics, multi-disciplinary and multi-scale problems involving interfaces.

T-Splines for Isogeometric Analysis of the Large Deformation of Elastoplastic Kirchhoff–Love Shells

In this paper, we develop a T-spline-based isogeometric method for the large deformation of Kirchhoff–Love shells considering highly nonlinear elastoplastic materials. The adaptive refinement is implemented, and some relatively complex models are considered by utilizing the superiorities of T-splines. A classical finite strain plastic model combining von Mises yield criteria and the principle of maximum plastic dissipation is carefully explored in the derivation of discrete isogeometric formulations under the total Lagrangian framework. The Bézier extraction scheme is embedded into a unified framework converting T-spline or NURBS models into Bézier meshes for isogeometric analysis. An a posteriori error estimator is established and used to guide the local refinement of T-spline models. Both standard T-splines with T-junctions and unstructured T-splines with extraordinary points are investigated in the examples. The obtained results are compared with existing solutions and those of ABAQUS. The numerical results confirm that the adaptive refinement strategy with T-splines could improve the convergence behaviors when compared with the uniform refinement strategy.

Numerical comparison of three a posteriori error estimators for nonconforming finite element method

Abstract In this paper, we propose to compare three a posteriori error estimators namely equilibrated, star-based and residual based for the Poisson problem and the Stokes problem with lowest-order Crouzeix-Raviart finite element discretization. The numerical results are presented to compare the performance of the three estimators in an adaptive refinement strategy.