共查询到19条相似文献,搜索用时 250 毫秒
1.
2.
3.
4.
5.
针对二维颗粒流体系统Euler/Lagrange模型的有限元模拟,建立了三角网格生成的自适应算法。该算法能够根据颗粒分布与颗粒大小自适应地调整网格的疏密程度,使其网格密度在系统边界附近及颗粒边缘附近较大,而在其它地方较小。与此同时,网格的光滑化也提高了网格质量, 从而为颗粒流体系统介观尺度的有限元模拟奠定了基础。 相似文献
6.
本文针对三维柱形区域提出了定常/非定常Stokes方程基于一致分裂格式的维数分裂算法(DSA).文章推导了三维定常/非定常Stokes方程维数分裂方法的数值迭代格式.新算法的优势在于一系列的二维问题能够并行执行,而且数值计算中避免了三维网格的生成.大量的数值结果表明新算法既能获得最优收敛阶,而且能获得比采用四面体元求解更精确的逼近解.最后,通过采用并行求解新算法能够得到比较好的加速比和并行效率. 相似文献
7.
本文针对三维柱形区域提出了定常/非定常Stokes方程基于一致分裂格式的维数分裂算法(DSA).文章推导了三维定常/非定常Stokes方程维数分裂方法的数值迭代格式.新算法的优势在于一系列的二维问题能够并行执行,而且数值计算中避免了三维网格的生成.大量的数值结果表明新算法既能获得最优收敛阶,而且能获得比采用四面体元求解更精确的逼近解.最后,通过采用并行求解新算法能够得到比较好的加速比和并行效率. 相似文献
8.
本文从变域变分原理出发,提出了一种新的有限元方法,并将所提算法用于求解形状最优化问题。我们给出了最优化问题解的存在性及其离散形式。讨论了数值解的收敛性,给出了数值解与精确解之间的误差估计。最后给出了算法的具体实施步骤和一个具有代表性的数值算例,数值结果表明该算法是正确、高效的。 相似文献
9.
10.
11.
Hanquan Wang Yang Xiang 《International journal for numerical methods in engineering》2013,94(6):573-597
In this paper, we present an adaptive level set method for motion of high codimensional objects (e.g., curves in three dimensions). This method uses only two (or a few fixed) levels of meshes. A uniform coarse mesh is defined over the whole computational domain. Any coarse mesh cell that contains the moving object is further divided into a uniform fine mesh. The coarse‐to‐fine ratios in the mesh refinement can be adjusted to achieve optimal efficiency. Refinement and coarsening (removing the fine mesh within a coarse grid cell) are performed dynamically during the evolution. In this adaptive method, the computation is localized mostly near the moving objects; thus, the computational cost is significantly reduced compared with the uniform mesh over the whole domain with the same resolution. In this method, the level set equations can be solved on these uniform meshes of different levels directly using standard high‐order numerical methods. This method is examined by numerical examples of moving curves and applications to dislocation dynamics simulations. This two‐level adaptive method also provides a basis for using locally varying time stepping to further reduce the computational cost. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
12.
D. Guégan O. Allain A. Dervieux F. Alauzet 《International journal for numerical methods in engineering》2010,84(11):1376-1406
This paper discusses the contribution of mesh adaptation to high‐order convergence of unsteady multi‐fluid flow simulations on complex geometries. The mesh adaptation relies on a metric‐based method controlling the L p‐norm of the interpolation error and on a mesh generation algorithm based on an anisotropic Delaunay kernel. The mesh‐adaptive time advancing is achieved, thanks to a transient fixed‐point algorithm to predict the solution evolution coupled with a metric intersection in the time procedure. In the time direction, we enforce the equidistribution of the error, i.e. the error minimization in L ∞ norm. This adaptive approach is applied to an incompressible Navier–Stokes model combined with a level set formulation discretized on triangular and tetrahedral meshes. Applications to interface flows under gravity are performed to evaluate the performance of this method for this class of discontinuous flows. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
13.
14.
C. S. Jog R. B. Haber M. P. Bendse 《International journal for numerical methods in engineering》1994,37(8):1323-1350
Significant performance improvements can be obtained if the topology of an elastic structure is allowed to vary in shape optimization problems. We study the optimal shape design of a two-dimensional elastic continuum for minimum compliance subject to a constraint on the total volume of material. The macroscopic version of this problem is not well-posed if no restrictions are placed on the structure topoiogy; relaxation of the optimization problem via quasiconvexification or homogenization methods is required. The effect of relaxation is to introduce a perforated microstructure that must be optimized simultaneously with the macroscopic distribution of material. A combined analytical-computational approach is proposed to solve the relaxed optimization problem. Both stress and displacement analysis methods are presented. Since rank-2 layered composites are known to achieve optimal energy bounds, we restrict the design space to this class of microstructures whose effective properties can easily be determined in explicit form. We develop a series of reduced problems by sequentially interchanging extremization operators and analytically optimizing the microstructural design fields. This results in optimization problems involving the distribution of an adaptive material that continuously optimizes its microstructure in response to the current state of stress or strain. A further reduced problem, involving only the response field, can be obtained in the stress-based approach, but the requisite interchange of extremization operators is not valid in the case of the displacement-based model. Finite element optimization procedures based on the reduced displacement formulation are developed and numerical solutions are presented. Care must be taken in selecting the discrete function spaces for the design density and displacement response, since the reduced problem is a two-field, mixed variational problem. An improper choice for the solution space leads to instabilities in the optimal design similar to those encountered in mixed formulations of the Stokes problem. 相似文献
15.
A. Dutta C. V. Ramakrishnan 《International journal for numerical methods in engineering》1998,41(6):977-999
Design sensitivities for structures under transient dynamic loads with constraints on displacements and stresses are sensitive to proper space and time discretization. Accuracy within acceptable error limit is feasible when an appropriate time increment coupled with an optimal mesh is used. In this paper, we handle this problem by systematically achieving an adaptive mesh for a reasonably fine but constant time step. Design sensitivities calculated for a good number of examples demonstrate the behaviour of this integrated approach. Comparison is made in terms of total computational time between time-marching scheme and modal superposition method in the context of design sensitivity calculation. Optimal meshes are also obtained corresponding to adaptive time stepping and accurate values of design sensitivities are computed using the optimal mesh and the values of the time increment obtained adaptively. © 1998 John Wiley & Sons, Ltd. 相似文献
16.
G. Bugeda J. Oliver 《International journal for numerical methods in engineering》1993,36(18):3161-3185
The proposed methodology is based on the use of the adaptive mesh refinement (AMR ) techniques in the context of 2D shape optimization problems analysed by the finite element method. A suitable and very general technique for the parametrization of the optimization problem, using B-splines to define the boundary, is first presented. Then mesh generation, using the advancing frontal method, the error estimator and the mesh refinement criterion are studied in the context of shape optimization problems In particular, the analytical sensitivity analysis of the different items ruling the problem (B-splines. finite element mesh, structural behaviour and error estimator) is studied in detail. The sensitivities of the finite element mesh and error estimator permit their projection from one design to the next one leading to an a priori knowledge of the finite element error distribution on the new design without the necessity of any additional structural analysis. With this information the mesh refinement criterion permits one to build up a finite element mesh on the new design with a specified and controlled level of error. The robustness and reliability of the proposed methodology is checked by means of several examples. 相似文献
17.
R. Radovitzky M. Ortiz 《International journal for numerical methods in engineering》1998,43(4):607-619
A fully Lagrangian finite element method for the analysis of Newtonian flows is developed. The approach furnishes, in effect, a Lagrangian implementation of the compressible Navier–Stokes equations. As the flow proceeds, the mesh is maintained undistorted through continuous and adaptive remeshing of the fluid mass. The principal advantage of the present approach lies in the treatment of boundary conditions at material surfaces such as free boundaries, fluid/fluid or fluid/solid interfaces. In contrast to Eulerian approaches, boundary conditions are enforced at material surfaces ab initio and therefore require no special attention. Consistent tangents are obtained for Lagrangian implicit analysis of a Newtonian fluid flow which may exhibit compressibility effects. The accuracy of the approach is assessed by comparison of the solution for a sloshing problem with existing numerical results and its versatility demonstrated through a simulation of wave breaking. The finite element mesh is maintained undistorted throughout the computation by recourse to frequent and adaptive remeshing © 1998 John Wiley & Sons, Ltd. 相似文献
18.
Fredrik Bengzon Mats G. Larson 《International journal for numerical methods in engineering》2010,84(12):1451-1465
In this paper we consider finite element simulation of the mechanical response of an elastic solid immersed into a viscous incompressible fluid flow. For simplicity, we assume that the mechanics of the solid is governed by linear elasticity and the motion of the fluid by the Stokes equation. For this one‐way coupled multiphysics problem we derive an a posteriori error estimate using duality techniques. Based on the estimate we propose an adaptive algorithm that automatically constructs a suitable mesh for the fluid and solid computational domains given a specific goal quantity for the elastic problem. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
19.
We consider optimal shape design problems for polymer spin packs which are widely used in the production of synthetic fibers and nonwoven materials. The design goal is the minimization of the residence time of the polymer, which can be achieved by adjusting the wall shear stress along the boundary. Depending on the specific industrial setting we construct two tailored algorithms. First, we consider the design in three spatial dimensions based on a PDE constrained shape optimization problem. Here, the constraint is given by the Stokes flow. Second, we change the design goal and want to construct shapes in two spatial dimensions which allow for a lower bound on the wall shear stress. This can be incorporated as an additional state constraint. By relaxing this condition and employing the method of mapping we can pull-back the problem onto a fixed reference domain. We get an elegant formulation of this state constrained optimization problem, in which geometric constraints on the boundary can also be included. After discretization we end up with a large-scale NLP which can be handled by existing solvers. Finally, we present numerical results underlining the feasibility of our approach. 相似文献