Search results

The purpose of this paper is twofold. Firstly, to present a detailed account of the generalized Lagrangian formulation of Hui and Zhao, in which the stream function ζ and Lagrangian distance λ, are used as independent variables, and secondly to assess and compare the performance of various flux limiters in this formulation with their corresponding performance in the Eulerian formulation. The generalized Lagrangian formulation is obtained by a transformation from the cartesian co‐ordinates (x, y) to the Lagrangian co‐ordinates (λ, ζ). In this manner, the number of independent variables for steady, 3‐D flow is reduced from four to three, placing this formulation on the same footing as the Eulerian formulation even for steady flows (as opposed to the conventional Lagrangian formulation which apparently still requires four independent variables even for steady flows). The generalized Lagrangian formulation with the Godunov scheme (using flux limiters) appears to have distinct advantages over the corresponding Eulerian formulation, particularly with respect to accuracy. Furthermore, the method requires no grid generation.

The purpose of this paper is the development of a new density-based (DB) semi-Lagrangian method to speed up the conventional pressure-based (PB) semi-Lagrangian methods.

The semi-Lagrangian-based solvers are typically PB, i.e. semi-Lagrangian pressure-based (SLPB) solvers, where a Poisson equation is solved for obtaining the pressure field and ensuring a divergence-free flow field. As an elliptic-type equation, the Poisson equation often relies on an iterative solution, so it can create a challenge of parallel computing and a bottleneck of computing speed. This study proposes a new DB semi-Lagrangian method, i.e. the semi-Lagrangian artificial compressibility (SLAC), which replaces the Poisson equation by a hyperbolic continuity equation with an added artificial compressibility (AC) term, so a time-marching solution is possible. Without the Poisson equation, the proposed SLAC solver is faster, particularly for the cases with more computational cells, and better suited for parallel computing.

The study compares the accuracy and the computing speeds of both SLPB and SLAC solvers for the lid-driven cavity flow and the step-flow problems. It shows that the proposed SLAC solver is able to achieve the same results as the SLPB, whereas with a 3.03 times speed up before using the OpenMP parallelization and a 3.35 times speed up for the large grid number case (512 × 512) after the parallelization. The speed up can be improved further for larger cases because of increasing the condition number of the coefficient matrixes of the Poisson equation.

This paper proposes a method of avoiding solving the Poisson equation, a typical computing bottleneck for semi-Lagrangian-based fluid solvers by converting the conventional PB solver (SLPB) to the DB solver (SLAC) through the addition of the AC term. The method simplifies and facilitates the parallelization process of semi-Lagrangian-based fluid solvers for modern HPC infrastructures, such as OpenMP and GPU computing.

The arbitrary Lagrangian—Eulerian (ALE) formulation, which is already well established in the hydrodynamics and fluid‐structure interaction fields, is extended to materials with memory, namely, non‐ linear path‐dependent materials. Previous attempts to treat non‐ linear solid mechanics with the ALE description have, in common, the implicit interpolation technique employed. Obviously, this implies a numerical burden which may be uneconomical and may induce to give up this formulation, particularly in fast‐transient dynamics where explicit algorithms are usually employed. Here, several applications are presented to show that if adequate stress updating techniques are implemented, the ALE formulation could be much more competitive than classical Lagrangian computations when large deformations are present. Moreover, if the ALE technique is interpreted as a simple interpolation enrichment, adequate—in opposition to distorted or locally coarse—meshes are employed. Notice also that impossible computations (or at least very involved numerically) with a Lagrangian code are easily implementable in an ALE analysis. Finally, it is important to observe that the numerical examples shown range from a purely academic test to real engineering simulations. They show the effective applicability of this formulation to non‐linear solid mechanics and, in particular, to impact, coining or forming analysis.

This paper, which is concerned with fluid‐structure interaction analysis, is a sequel to our earlier paper which gave an introduction to the numerical treatment of such systems. The paper is divided into five main sections. In the first two, a state‐of‐the‐art review on near‐field and far‐field fluid structure interaction is presented. In attempting to highlight where current research should be directed, only the most widely used computer codes are reviewed in the third section. Conclusions are presented in the fourth section.

The purpose of this study is to design numerical simulations of bubbly flow around a cylinder to better understand the characteristics of flow around a rigid obstacle.

The bubbly flow around a circular cylinder was numerically simulated using a semi-Lagrangian–Lagrangian method composed of a vortex-in-cell method for the liquid phase and a Lagrangian description of the gas phase. Additionally, a penalization method was applied to account for the cylinder inside the flow. The slip condition of the bubbles on the cylinder’s surface was enforced, and the outflow conditions were applied to the liquid flow at the far field.

The simulation clarified the characteristics of a bubbly flow around a circular cylinder. The bubbles were shown to move around and separate from both sides of the cylinder, because of entrainment by the liquid shear layers. Once the bubbly flow fully developed, the bubbles distributed into groups and were dispersed downstream of the cylinder. A three-dimensional vortex structure of various scales was also shown to form downstream, whereas a quasi-stable two-dimensional vortex structure was observed upstream. Overall, the proposed method captured the characteristics of a bubbly flow around a cylinder well.

A semi-Lagrangian–Lagrangian approach was applied to simulate a bubbly flow around a circular cylinder. The simulations provided the detail features of these flow phenomena.

This paper presents a novel framework for solving elliptic partial differential equations (PDEs) over irregularly spaced meshes on bounded domains.

Second‐generation wavelet construction gives rise to a powerful generalization of the traditional hierarchical basis (HB) finite element method (FEM). A framework based on piecewise polynomial Lagrangian multiwavelets is used to generate customized multiresolution bases that have not only HB properties but also additional qualities.

For the 1D Poisson problem, we propose – for any given order of approximation – a compact closed‐form wavelet basis that block‐diagonalizes the stiffness matrix. With this wavelet choice, all coupling between the coarse scale and detail scales in the matrix is eliminated. In contrast, traditional higher‐order (n>1) HB do not exhibit this property. We also achieve full scale‐decoupling for the 2D Poisson problem on an irregular mesh. No traditional HB has this quality in 2D.

Similar techniques may be applied to scale‐decouple the multiresolution finite element (FE) matrices associated with more general elliptic PDEs.

By decoupling scales in the FE matrix, the wavelet formulation lends itself particularly well to adaptive refinement schemes.

The paper explains second‐generation wavelet construction in a Lagrangian FE context. For 1D higher‐order and 2D first‐order bases, we propose a particular choice of wavelet, customized to the Poisson problem. The approach generalizes to other elliptic PDE problems.

The purpose of this paper is to evaluate the health‐states of unit under test (UUT) in aerospace systems by means of unreliable test outcomes, and the evaluation results can provide a guide for engineers to carry out proper maintenance prior to total failure.

In this paper, the authors formulate the health‐state evaluation (HSE) problem with unreliable test outcomes based on Bayes rule, and develop the Lagrangian relaxation and adaptive genetic algorithm (LRAGA) to solve it. The solution scheme can be viewed as a two‐level coordinated solution framework for the HSE problem. At the top level, the Lagrange multipliers are updated by using AGA. At the bottom level, each of the sub‐problems is solved by using AGA.

The experimental results show that the HSE model appears promising and the LRAGA can obtain the higher quality solution and converge to it at a faster rate than conventional methods (i.e. Lagrangian relaxation (LR), genetic algorithm (GA), simulated annealing (SA) and Lagrangian relaxation and genetic algorithm (LRGA).

The proposed method for the HSE problem of large‐scale systems which include thousands of faults and tests needs to be verified further.

The HSE results for aerospace systems can help engineers to carry out a schedule for prompt maintenance prior to UUTs' failure, to avoid the consequences of total failure. It is important to improve aerospace systems' safety, reliability, maintainability, affordability, and reduce life cycle cost.

This paper constructs the HSE model with unreliable test outcomes based on the Bayes rule and proposes a method based on LRAGA to solve the HSE problem.

Large strain model can be formulated in terms of the Lagrangian or the Eulerian frame. In this paper, the Eulerian type large strain models are studied. Numerical examples on the Lagrangian and Eulerian types large strain models are investigated and compared. It is found that the differences in the choice of large strain model under large strain and rotation problems are noticeable but not significant if small load step is used for analysis. Furthermore, we have also found that unsymmetrical formulation instead of symmetrical formulation should be adopted for Eulerian type large strain models.

This paper aims to provide discussions of a numerical method for bubbly flows and the interaction between a vortex ring and a bubble plume.

Small bubbles are released into quiescent water from a cylinder tip. They rise under the buoyant force, forming a plume. A vortex ring is launched vertically upward into the bubble plume. The interactions between the vortex ring and the bubble plume are numerically simulated using a semi-Lagrangian–Lagrangian approach composed of a vortex-in-cell method for the fluid phase and a Lagrangian description of the gas phase.

A vortex ring can transport the bubbles surrounding it over a distance significantly depending on the correlative initial position between the bubbles and the core center. The motion of some bubbles is nearly periodic and gradually extinguishes with time. These bubble trajectories are similar to two-dimensional-helix shapes. The vortex is fragmented into multiple regions with high values of Q, the second invariant of velocity gradient tensor, settling at these regional centers. The entrained bubbles excite a growth rate of the vortex ring's azimuthal instability with a formation of the second- and third-harmonic oscillations of modes of 16 and 24, respectively.

A semi-Lagrangian–Lagrangian approach is applied to simulate the interactions between a vortex ring and a bubble plume. The simulations provide the detail features of the interactions.

The purpose of this paper is to develop a Vortex-In-Cell (VIC) method with the semi-Lagrangian scheme and apply it to the high-Re lid-driven cavity flow.

The VIC method is developed for simulating high Reynolds number incompressible flow. A semi-Lagrangian scheme is incorporated in the convection term to produce unconditional stability, which gets rid of the constraint of the convection Courant-Friedrichs-Lewy (CFL) condition; the adaptive time step is used to maintain the numerical stability of the diffusion term; and the velocity boundary condition is readily converted to the vorticity formulation to suit discontinuous boundary treatment. The VIC simulation results are compared with those produced by other gird methods reported in open literature studies.

The lid-driven cavity flow is simulated from Re = 100 to 100,000. Similar vortex birth mechanisms are exhibited though, but distinct flow characteristics are revealed. At Re = 100 to 7,500, the cavity flow is confirmed steady. At Re = 10,000, 15,000 and 20,000, the cavity flow is periodical with a primary vortex held spatially at the center. In particular, at Re = 100,000 highly turbulent characteristics is first revealed and an analogous primary vortex is formed but in motion rather than stationary, which is caused by the considerable flow separation at all the boundaries.

In the lid-driven cavity, the flow becomes extremely complex and highly turbulent at Re = 100,000, and the analogous primary vortex structure is observed. Boundary layer separation is observed at all walls, producing small vortices and causing the displacement of the analogous primary vortex. Such a finding original and has not yet been reported by other investigators. It may provide a basis for conducting in-depth studies of the lid-driven cavity flow.