In this paper, we consider numerical approximations of the Cahn–Hilliard type phase-field crystal model and construct a fully discrete finite element scheme for it. The scheme is the combination of the finite element method for spatial discretization and an invariant energy quadratization method for time marching. It is not only linear and second-order time-accurate, but also unconditionally energy-stable. We prove the unconditional energy stability rigorously and further carry out various numerical examples to demonstrate the stability and the accuracy of the developed scheme numerically.
This work introduces a general strategy to develop well-balanced high-order Discontinuous Galerkin (DG) numerical schemes for systems of balance laws. The essence of our approach is a local projection step that guarantees the exactly well-balanced character of the resulting numerical method for smooth stationary solutions. The strategy can be adapted to some well-known different time marching DG discretisations. Particularly, in this article, Runge–Kutta DG and ADER DG methods are studied. Additionally, a limiting procedure based on a modified WENO approach is described to deal with the spurious oscillations generated in the presence of non-smooth solutions, keeping the well-balanced properties of the scheme intact. The resulting numerical method is then exactly well-balanced and high-order in space and time for smooth solutions. Finally, some numerical results are depicted using different systems of balance laws to show the performance of the introduced numerical strategy.
In this paper, an explicit time marching procedure for solving the non-hydrostatic shallow water equation (SWE) problems is developed. The procedure includes a process of prediction and several iterations of correction. In these processes, it is essential to accurately calculate the spatial derives of the physical quantities such as the temporal water depth, the average velocities in the horizontal and vertical directions, and the dynamic pressure at the bottom. The weighted-least-squares (WLS) meshless method is employed to calculate these spatial derivatives. Though the non-hydrostatic shallow water equations are two dimensional, on the focus of presenting this new time marching approach, we just use one dimensional benchmark problems to validate and demonstrate the stability and accuracy of the present model. Good agreements are found in the comparing the present numerical results with analytic solutions, experiment data, or other numerical results.
Decarbonization and sustainability efforts challenge gas turbine engineers to come up with creative strategies for reduction of emissions and efficiency increase over the whole operating range. Burner staging at part loads presents a flexible solution to achieve these goals through selective burner deactivation. Shutting off burners could also be required for combustion of increased H2 content at some conditions. Burner staging will create circumferential unevenness with patterns of hot and cold streaks that could excite blade rows through the entire turbine. This paper presents a parametric method for annular combustor staging patterns profile generation intended for use for forced response predictions from a limited number of combustor CFD calculations while keeping the key phenomenological features. Two cases with burner staging turbine inlet temperature distributions are considered and compared to a base case with uniform temperature distribution. The unsteady aerodynamic forcing was obtained from full wheel time marching unsteady computational fluid dynamics calculations. The results show that the hot streaks generate important and noticeable excitation sources. Additionally, the results show that the pattern generator could be used extensively before and after the unsteady calculations phase to minimize the excitation levels and the computational load.
AbstractThis work presents the development of an explicit/implicit particle finite element method (PFEM) for the 2D modeling of metal cutting processes. The purpose is to study the efficiency of implicit and explicit time integration schemes in terms of precision, accuracy and computing time. The formulation for implicit and explicit time marching schemes is developed, and a detailed study on the explicit solution steps is presented. The PFEM remeshing procedures for insertion and removal of particles have been improved to model the multiple scales of time and/or space of the solution. The detection and treatment of the rigid tool contact are presented for both, implicit and explicit schemes. The performance of explicit/implicit integration is studied with a set of different two-dimensional orthogonal cutting tests of AISI 4340 steel at cutting speeds ranging from 1 m/s up to 30 m/s. It was shown that if the correct selection of the time integration scheme is made, the computing time can decrease up to 40 times. It allows us to affirm that the computing time of the PFEM simulations can be excessive due to the used time marching scheme independently of the meshing process. As a practical result, a set of recommendations to select the time integration schemes for a given cutting speed are given. This is intended to minimize one of the negative constraints pointed out by the industry when using metal cutting simulators.