共查询到20条相似文献，搜索用时 0 毫秒
1.
This paper presents a contribution to level‐set reinitialization in the context of discontinuous Galerkin finite element methods. We focus on high‐order polynomials for the discretization and level set geometries, which are comparable to the element size. In contrast to hyperbolic and geometric reinitialization techniques, our method relies on solving a nonlinear elliptic PDE iteratively. We critically compare two different variants of the algorithm experimentally in numerical studies. The results demonstrate that the method is stable for nontrivial test cases and shows high‐order accuracy. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
2.
The reinitialization, which is required to regularize the level set function, can be computationally expensive and hence is a determining factor for the overall efficiency of a level set method. However, it often has a significantly adverse impact on the accuracy of the level set solution. This short note is meant to shed light on the efficiency and accuracy issues of the reinitialization process. Using just one clearly defined level set propagation test case with an analytical solution the solutions obtained using a recently proposed efficient lower‐order constrained reinitialization (CR) scheme and standard low‐ and high‐order reinitialization schemes are juxtaposed to evidence the superiority of the novel CR formulation. It is shown that maintaining the location of the zero level set during the reinitialization is crucial for the accuracy and that the displacement caused by standard high‐order reinitialization schemes clearly outweighs the benefit of the high‐order smoothing of the level set function. Finally, results of a three‐dimensional problem are concisely reported to demonstrate the general applicability of the CR scheme. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
3.
A hybrid scheme for reinitializing the level set function and its gradient within the frame work of the augmented level set method is presented. It is based on first dividing the domain into an interfacial region (i.e. nodes close to the interface) and its complement. Within the interfacial region, the level set and its gradient are updated explicitly through a modified version of Newton's method (Chopp, 2001, SIAM J. Sci. Comput. 23 230‐244) and is implemented here within the context of Hermite polynomials. In the region away from the interface, the solution pertains to a semi‐Lagrangian implementation of the reinitialization equations, which are solved based on Hermite polynomials and are time marched with a single step and a multipoint scheme. It is shown that for various exercises, the present method predicts the signed distance function and its gradient to 4th and 3rd order (in space), respectively with regards to the L1, L2, and L ∞ norms, provided the level set field is sufficiently smooth. A range of test cases are also considered from the literature, where the present method is compared with existing methods and shown to be generally more accurate. Moreover, the well‐known issue of volume loss due to reinitialization is addressed successfully with the current implementation, even for objects that are of the size of one grid cell, and whose local radius of curvature falls below the local grid size. For both time marching schemes, it is shown that the L2 and L ∞ errors decay to negligible levels, are smooth in space, and do not exhibit temporal oscillations. Finally the performance of the hybrid scheme is evaluated by applying it on various kinematic test cases. For solid body rotation problems (zero deformation flow field), the benefit stemming from hybrid reinitialization is marginal. When applied to kinematic cases involving severe deformation, such as the standard vortex flow, the reinitialization strategy helps maintain a smooth level set field, which prevents serious numerical errors from developing.Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
4.
马景槐 《应用数学和力学(英文版)》1997,(1)
I.Viscoplastic-DynamicsConstitutiveRelationThestructurematerialwouldarisesimultaneouslytheelasticity,viscosityandplasticityindynamichighstraincondition.Atpresent,theovcrstressmodel,viscoplasticitymodel,quasi-linearconstitutiverelationtheoryandintrinsicvar… 相似文献
5.
马景槐 《应用数学和力学(英文版)》1997,18(1):61-68
This paper presents the optimal control variational principle for Perzyna model which is one of the main constitutive relation of viscoplasticity in dynamics. And it could also be transformed to solve the parametric quadratic programming problem. The FEM form of this problem and its implementation have also been discussed in the paper. 相似文献
6.
In this paper we describe and evaluate a geometric mass‐preserving redistancing procedure for the level set function on general structured grids. The proposed algorithm is adapted from a recent finite element‐based method and preserves the mass by means of a localized mass correction. A salient feature of the scheme is the absence of adjustable parameters. The algorithm is tested in two and three spatial dimensions and compared with the widely used partial differential equation (PDE)‐based redistancing method using structured Cartesian grids. Through the use of quantitative error measures of interest in level set methods, we show that the overall performance of the proposed geometric procedure is better than PDE‐based reinitialization schemes, since it is more robust with comparable accuracy. We also show that the algorithm is well‐suited for the highly stretched curvilinear grids used in CFD simulations. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
7.
O. Pironneau 《International Journal of Computational Fluid Dynamics》2018,32(2-3):121-130
ABSTRACTEulerian variational formulations for deformable solids, with or without fluids around them, end up, after implicit time discretisation, as large non-linear systems for the velocities in the moving domains. Handling moving domains and moving boundaries requires careful meshing procedures; on the other hand, the detection of contact is particularly simple with a distance function. Then at every time step, a variational inequality can be used to update the velocities. This article gives new implementation details and two new complex simulations: a very soft bouncing ball in an axisymmetric flow and a disk hit by a club. 相似文献
8.
主要研究了扩展有限元法(extended finite element method, XFEM)在处理弱不连续问题时不同改进函数形式对XFEM数值求解精度的影响,阐述了各种改进函数影响XFEM求解精度的关键因素,指出校正的扩展有限元法(corrected-XFEM)能够提高数值求解精度的实质在于它拓展了改进结点域,即将常规扩展有限元法(standard-XFEM)的改进结点域增加一层作为corrected-XFEM的改进结点域,文中建议延拓corrected-XFEM的改进结点域,即在corrected-XFEM的改进结点域基础上再增加一层改进结点. 利用水平集函数表征材料内部的不连续界面,推导了XFEM求解的支配方程,给出了一种改进单元的数值积分方案以及改进单元处高精度应力的求解方法. 含夹杂问题的数值计算结果表明：建议的延拓corrected-XFEM改进结点域的方法能够明显提高XFEM的数值求解精度. 相似文献
9.
10.
Tetsuya Miyaoka 《International Journal of Computational Fluid Dynamics》2013,27(10):677-686
The purpose of this article is to present a technique to optimally control river flood using a drainage basin considering a moving boundary. The main theme of this article is to obtain outflow discharge from the drainage basin that maintains the water level at a downstream point and empties the drainage basin as soon as possible. The water flow phenomenon inside the drainage basin when a river flood occurs is considered. This phenomenon can be analysed by the finite element method considering a moving boundary. The optimal control theory can be implemented to obtain the optimal control discharge. The finite element analysis with a moving boundary is introduced in the optimal control theory. A new boundary condition on the downstream side of the river is proposed. This condition is formulated by the solitary wave condition based on the basic water level being capable of representing natural water surface. As a numerical study, optimal control of shallow water flow is carried out for the Tsurumi River and its drainage basin model. 相似文献
11.
In this paper we present a problem we have encountered using a stabilized finite element method on fixed grids for flows with interfaces modelled with the level set approach. We propose a solution based on enriching the pressure shape functions on the elements cut by the interface. The enrichment is used to enable the pressure gradient to be discontinuous at the interface, thus improving the ability to simulate the behaviour of fluids with different density under a gravitational force. The additional shape function used is local to each element and the corresponding degree of freedom can therefore be condensed prior to assembly, making the implementation quite simple on any existing finite element code. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
12.
‘Fluid buckling’ is a phenomenon characterized mainly by the existence of fluid toroidal oscillations during flow. It appears when a high viscosity fluid flows vertically against a flat surface and may occur in industrial applications, as in injection molding of a propergol in complex‐shaped cavities. These coiling or folding oscillations appear during the mold filling stage, leading to air entrapment. To understand and to model this free surface flow problem, a convected level set method is proposed. First, a sinus filter is applied to the distance function to get a smooth truncation far from the interface. Second, the reinitialization is embedded in the transport equation model, avoiding it as a separate step during calculation. In order to validate the method, numerical results are presented on classical interface capturing benchmarks. Finally, results are shown on two‐dimensional and three‐dimensional viscous jet buckling problems. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
13.
为了在三角形和四边形网格上采用水平集方法模拟轴对称爆轰波阵面与曲率相关的运动,假设爆轰波阵面的法向速度是曲率的线性函数,通过坐标变换得到了轴对称坐标系下的水平集方程。水平集方程的曲率无关项采用正格式离散,曲率项采用伽辽金等参有限元方法空间离散,时间离散采用半隐格式。算例表明,在轴对称三角形网格和四边形网格上,含曲率的水平集方程的离散格式为强一阶精度。给出了三角形和四边形混合网格上不光滑界面以曲率收缩的运动过程,收缩过程未出现不稳定现象。多个爆轰波阵面相互作用的算例说明本文的格式可有效地模拟曲率相关的爆轰波的相互作用问题 相似文献
14.
The vapor deposition chemical reaction processes, which are of extremely extensive applications, can be classified as a mathematical model by the following governing nonlinear partial differential equations containing velocity vector, temperature field, pressure field, and gas mass field. The mixed finite element (MFE) method is employed to study the system of equations for the vapor deposition chemical reaction processes. The semidiscrete and fully discrete MFE formulations are derived. And the existence and convergence (error estimate) of the semidiscrete and fully discrete MFE solutions are demonstrated. By employing MFE method to treat the system of equations for the vapor deposition chemical reaction processes, the numerical solutions of the velocity vector, the temperature field, the pressure field, and the gas mass field can be found out simultaneously. Thus, these researches are not only of important theoretical means, but also of extremely extensive applied vistas. 相似文献
15.
In this paper, we present a finite element model for free surface flows on fixed meshes. The main novelty of the approach, compared with typical fixed mesh finite element models for such flows, is that we take advantage of the particularities of free surface flow, instead of considering it a particular case of two‐phase flow. The fact that a given free surface implies a known boundary condition on the interface, allows us to solve the Navier–Stokes equations on the fluid domain uncoupled from the solution on the rest of the finite element mesh. This, together with the use of enhanced integration allows us to model low Froude number flows accurately, something that is not possible with typical two‐phase flow models applied to free surface flow. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
16.
A. Unger F. Trltzsch 《ZAMM - Journal of Applied Mathematics and Mechanics / Zeitschrift für Angewandte Mathematik und Mechanik》2001,81(7):447-456
We consider the problem of cooling milled steel profiles at a maximum rate subject to given bounds on the difference of temperatures in prescribed points of the steel profile. This leads to a nonlinear parabolic control problem with pointwise state constraints in a 2D domain. A method of instantaneous control is applied to set up a fast solution technique. 相似文献
17.
A method is developed for modeling fluid transport in domains that do not conform to the finite element mesh. One or more level set functions are used to describe the fluid domain. A background, non‐conformal mesh is decomposed into elements that conform to the level set interfaces. Enrichment takes place by adding nodes that lie on the interfaces. Unlike other enriched finite element methods, the proposed technique requires no changes to the underlying element assembly, element interpolation, or element quadrature. The complexity is entirely contained within the element decomposition routines. It is argued that the accuracy of the method is no less than that for eXtended Finite Element Methods (XFEM) with Heaviside enrichment. The accuracy is demonstrated using multiple numerical tests. In all cases, optimal rates of convergence are obtained for both volume and surface quantities. Jacobi preconditioning is shown to remove the ill‐conditioning that may result from the nearly degenerate conformal elements. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
18.
Masanori Shichitake 《International Journal of Computational Fluid Dynamics》2013,27(5):351-365
This paper presents an optimal control applied to water flow using the first and second order adjoint equations. The gradient of the performance function with respect to control variables is analytically obtained by the first order adjoint equation. It is not necessary to compute the Hessian matrix directly using the second order adjoint equation. Two numerical studies have been performed to show the adaptability of the present method. The performance of the second order adjoint method is compared with that of the weighted gradient method, Broyden–Fletcher–Goldfarb–Shanno method and Lanczos method. The precise forms of the adjoint equations and the gradient to use for the minimisation algorithm are derived. The computation by the Lanczos method is shown as superior to those of the other methods discussed in this paper. The message passing interface library is used for the communication of parallel computing. 相似文献
19.
The paper describes and compares the performance of two options for numerically representing the surface tension force in combination with the level set interface‐tracking method. In both models, the surface tension is represented as a body force, concentrated near the interface, but the technical implementation is different: the first model is based on a traditional level set approach in which the force is distributed in a band around the interface using a regularized delta function, whereas in the second, the force is partly distributed in a band around the interface and partly localized to the actual computational cells containing the interface. A comparative study, involving analysis of several two‐phase flows with moving interfaces, shows that in general the two surface tension models produce results of similar accuracy. However, in the particular case of merging and pinching‐off of interfaces, the traditional level set model of surface tension produces an error that results in non‐converging solutions for film‐like interfaces (i.e. ones involving large contact areas). In contrast, the second model, based on the localized representation of the surface tension force, displays consistent first‐order convergence. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
20.
Using a hybrid Lagrangian-Eulerian approach, a level set function–based immersed interface method (LS-IIM) is proposed for the interaction of a flexible body immersed in a fluid flow. The LS-IIM involves finite volume method for the fluid solver, Galerkin finite element method for the structural solver, and a block-iterative partitioned method–based fully implicit coupling between the two solvers. The novelty of the proposed method is a level set function–based direct implementation of fluid-solid interface boundary conditions in both the solvers. Another novelty is the computation of the level set function from a geometric method instead of differential equations commonly used in level set methods—the novel geometric as compared to the traditional method is found to be more accurate and less time-consuming. The LS-IIM is demonstrated as second-order accurate. Verification study is presented first separately for both the solvers and then together for four fluid-structure interaction (FSI) problems, with different levels of complexity including lid-driven flow, channel flow, and free-stream flow. Benchmark solutions are presented for two class of FSI problems: first, easy to set up and less time-consuming and, second, a reasonably challenging and complex FSI problem involving sharp edges and forced-motion of the flexible structure. The benchmark solutions are proposed at steady state for the first problem, after a verification study with two open-source solvers and, at periodic state, after a validation with published experimental results for the second problem. Our benchmark solutions may be useful for verification study in future. 相似文献