比例延迟微分方程的连续有限元法

利用连续有限元法求解比例延迟微分方程,在一致网格下,给出比例延迟微分方程连续有限元解的整体收敛阶,数值实验验证了理论结果的正确性.     

Energy-Preserving $H^1$--Galerkin Schemes for the Hunter-Saxton Equation

We consider the numerical integration of the Hunter–Saxton equation, which models the propagation of weakly nonlinear orientation waves. For the equation, we present two weak forms and their Galerkin discretizations. The Galerkin schemes preserve the Hamiltonian of the equation and can be implemented with cheap H~1 elements. Numerical experiments confirm the effectiveness of the schemes.     

A class of discontinuous Petrov–Galerkin methods. II. Optimal test functions

We lay out a program for constructing discontinuous Petrov–Galerkin (DPG) schemes having test function spaces that are automatically computable to guarantee stability. Given a trial space, a DPG discretization using its optimal test space counterpart inherits stability from the well posedness of the undiscretized problem. Although the question of stable test space choice had attracted the attention of many previous authors, the novelty in our approach lies in the fact we identify a discontinuous Galerkin (DG) framework wherein test functions, arbitrarily close to the optimal ones, can be locally computed. The idea is presented abstractly and its feasibility illustrated through several theoretical and numerical examples. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

Short-cut distillation columns based on orthogonal polynomials of a discrete variable

Increasing competition in the process industries enforces theapplication of mathematical simulation techniques both in thedesign phase and in the operating phase of a plant. A basicapparatus for separation processes is the distillation column.Its rigorous (tray by tray) mathematical modelling results ina system of simultaneous nonlinear equations (algebraic in thesteady-state case, differential-algebraic in the dynamic case).For high (and realistic) numbers of trays and components, thesesystems may become rather large (thousands of equations). Inaddition, realistic plant models often include several distillationcolumns. As a consequence, the numerical solution of these modelsmay become difficult and time-consuming. This has led to attemptsto model the distillation columns less rigorously with the aimof achieving a considerable reduction in the number of equations.The name shortcut distillation columns is common for modelsof this type. The present paper uses a discrete weighted residualmethod for the development of short-cut models. It suggestsa Galerkin method based on orthogonal polynomials in a discretevariable: the tray number. It is a remarkable advantage of thistechnique that even very coarse models satisfy all global balancesexactly.     

A weak Galerkin method for diffraction gratings

This paper concerns a weak Galerkin method (WGM) for the diffraction of a time-harmonic incident wave impinging upon a one-dimensional periodic grating structure. The existence and uniqueness of the weak Galerkin solution to the grating problem are established using a variational approach. The convergence rate of the proposed WGM is systematically analyzed. Numerical simulations are presented to verify the efficiency of the WGM for solving grating problems.     

Coupled fluid–solid interaction under shock wave loading

This paper considers the treatment of fluid–solid interaction problems under shock wave loading, where the solid experiences large bulk Lagrangian displacements. This work addresses the issues associated with using a level set as a generalized interface for fluid–solid coupling where the fluid–solid interface is embedded in an unstructured fluid grid. We outline the formulation used for the edge‐based unstructured‐grid Euler solver. The identification of the fluid–solid interface on the unstructured fluid mesh uses a super‐sampled ??² projection technique, which in conjunction with a Lagrangian interface position, permits fast identification of the interface and the concomitant imposition of boundary conditions. The use of a narrow‐band approach for the identification of the wetted interface is presented with the details of the construction of interface conditions. A series of two and three‐dimensional shock‐body computations are presented to demonstrate the validity of the current approach on problems with static and dynamic interfaces, including projectile/shock interaction simulations. Copyright © 2010 John Wiley & Sons, Ltd.     

Discontinuous Galerkin methods for the Navier–Stokes equations using solenoidal approximations

An interior penalty method and a compact discontinuous Galerkin method are proposed and compared for the solution of the steady incompressible Navier–Stokes equations. Both compact formulations can be easily applied using high‐order piecewise divergence‐free approximations, leading to two uncoupled problems: one associated with velocity and hybrid pressure, and the other one only concerned with the computation of pressures in the elements interior. Numerical examples compare the efficiency and the accuracy of both proposed methods. Copyright © 2009 John Wiley & Sons, Ltd.     

Nonlinear Normal Modes of a Parametrically Excited Cantilever Beam

We investigate theoretically thenonlinear normal modes of a vertical cantilever beam excited by aprincipal parametric resonance. We apply directly the method ofmultiple scales to the governing nonlinear nonautonomousintegral-partial-differential equation and associated boundary conditions.In the absence of damping, it is shown that the system has nonlinear normal modes, as defined by Rosenberg, even in the presence of the parametric excitation.We calculate the spatial correction to the linear mode shapedue to the effects of the inertia and curvature nonlinearities andthe parametric excitation. We compare the result obtained withthe direct approach with that obtained using a single-mode Galerkindiscretization.The deviation between the two predictions increases as the oscillationamplitude increases.     

Are upwind techniques in multi-phase flow models necessary?

Two alternatives of primary variables are compared for two-phase flow in heterogeneous media by solving fully established benchmarks. The first combination utilizes pressure of the wetting fluid and saturation of the non-wetting fluid as primary variables, while the second employs capillary pressure of the wetting fluid and pressure of the non-wetting fluid. While the standard Galerkin finite element method (SGFEM) is known to fail in the physical reproduction of two-phase flow in heterogeneous media (unless employing a fully upwind correction), the second scheme with capillary pressure as a primary variable without applying an upwind technique produces correct physical fluid behaviour in heterogeneous media, as observed from experiments.