基于有限元方法的极小曲面造型

讨论极小曲面方程的求解。极小曲面方程是一个高度非线性的二阶椭圆偏微分方程，求解十分困难。该文基于有限元方法，使用一个简单而有效的线性化策略，将问题转化为一系列线性问题，从而大大简化了求解过程。数值结果表明该方法简单有效，能产生合理的结果。     

推广Dirichlet 方法用于B 样条极小曲面设计

为弥补当前NURBS系统无法有效设计工程所急需的B样条极小曲面的缺陷,将构造Bézier极小曲面的Dirichlet方法成功地推广到了B样条极小曲面设计.提出了插值控制网格边界的B样条曲面模型,运用B样条基函数的求导公式及求值割角算法,将计算极小曲面内部控制顶点的问题转化为一个线性方程组的求解,从而避免了强非线性问题所导致的困惑,极大地提高了运算效率.最后,用大量实例对理论和算法进行了验证.     

三重周期极小曲面夹芯结构的隔声性能研究

极小曲面夹芯结构具有高比刚度、良好吸能等优点,在实际工程中展现了广泛的应用前景,而其隔声性能并未得到充分的研究.基于此,利用数值分析和试验方法研究了 D型和I-WP型三重周期极小曲面夹芯结构的隔声性能.首先,基于声振耦合理论,建立极小曲面夹芯结构有限元模型;然后使用声阻抗管法验证了有限元模型的精度;最后,利用验证后的有限元模型系统分析了夹芯厚度、面板厚度、弹性模量、极小曲面常数和周期常数等重要参数对极小曲面夹芯结构隔声性能的影响.研究发现I-WP型极小曲面夹芯结构在4400Hz-5000Hz频段,隔声量可达60dB.本研究有望为高速列车等领域轻质结构被动隔声提供新的解决方案.     

旋转曲面变换PSO 算法解非线性最优控制问题

针对利用粒子群优化算法进行多极值点函数优化时,存在陷入局部极小点和搜寻效率低的问题.提出旋转曲面变换方法,将被优化函数映射到一个同胚曲面上.它将当前局部极小点变换为全局最大点,并保持被优化函数值在当前局部极小点以下部分的形状不变,从而克服陷入局部极小点的问题.最后将其用于解一个非线性系统的最优控制问题,实验结果证明了该方法的可行性和有效性.     

插值边界的四边网格离散极小曲面建模方法

如何实现极小曲面的快速三维建模,是几何设计与计算领域中的难点和热点问题.给定一条封闭的边界离散折线,本文研究如何构造以其为边界的四边网格离散极小曲面.首先从曲面的内蕴微分几何度量出发,给出了离散四边网格极小曲面的数学定义;然后利用保长度边界投影、四边网格生成、径向基函数插值映射和非线性优化技术,提出了由给定边界离散折线快速构造离散四边网格极小曲面的一般技术框架.最后通过若干建模实例验证了本文方法的有效性.该方法可实现四边网格极小曲面的高质量建模,在建筑几何领域具有一定的应用价值.     

任意边界下的极小曲面造型问题

讨论了任意边界下的极小曲面造型问题,提出了一个用B-样条函数做任意有界区域上极小曲面造型的新方法.基本思想是:对B-样条函数加权,权函数为节点到区域边界的距离函数,使用加权的B-样条函数空间作为有限元子空间,从而可以得到极小曲面在任意边界上的B-样条曲面近似解.结果表明,新方法得到的极小曲面具有良好的光顺性,且计算精度高.     

三周期极小曲面结构设计及应用综述EI北大核心CSCD

三周期极小曲面是一种隐函数表示的代数曲面,具有许多优良的性质.增材制造技术的快速发展极大地增强了复杂几何与拓扑结构的制造能力,三周期极小曲面作为几何建模工具越来越受到关注.首先,从数学表达、性质、应用以及几何设计方法等方面对三周期极小曲面的研究现状进行介绍;其次,对三周期极小曲面在力学、传热传质、组织工程和声学等方面的性能及相关应用进行总结,从几何建模方法的角度对三周期极小曲面进行分类梳理,将现有工作分为规则单元法、参数单元法、区域拼接法以及整体优化法4类,对各方法特点进行分析;最后,结合实际应用对该领域面临的挑战进行总结,并展望未来工作与发展趋势.     

插值边界曲线的NURBS 近似极小曲面设计

探索性地设计了一个插值给定边界曲线的NURBS 近似极小曲面算法,弥补了当前NURBS 系统无法有效地设计工程所急需的一般NURBS 极小曲面的缺陷.运用NURBS 曲面的节点插入、Hybrid 多项式逼近等多种技术,将NURBS 曲面转化为相对简单的分片Bézier 曲面求解,并运用各子曲面片的控制顶点优化、整体曲面不断更新的迭代方法,成功地得到高精度的近似分片Bézier 极小曲面.最后,可以按用户的各种要求选择运用相应不同的迭代逼近算法,求取插值给定边界曲线的近似NURBS 极小曲面.     

一种隐式曲面交互调整的新方法

提出了一种对隐式曲面形状进行交互调整的新方法,为隐式曲面的调整提供了两种交互工具,分别是对曲面上点的位置调整和法向调整.该方法以调整后的位置和法向为新曲面的插值条件建立目标函数,极小化该目标函数求解曲面参数的变化量,从而确定新的隐式曲面.文中采用拟牛顿法和序列二次规划法(SQP)求解该非线性优化问题.在调整过程中用粒子的方法对隐式曲面进行绘制,实现了对隐式曲面形状的实时交互调整.最后用实例说明了新方法的有效性.     

散乱分布数据曲面重构的光顺-有限元方法

提出了一种基于散乱分布的数据点重构三维曲面的有限元方法.根据最佳逼近与数据光顺理论建立正定的目标泛函,采用有限元最佳拟合使泛函极小化,求得最优解.通过八节点等参数有限元插值计算,重新构造出三维曲面.这种光顺-有限元方法有效地抑制了输入数据上误差噪声的影响,与有限元拟合方法相比,所需的输入数据点少,重构的曲面逼近精度高、光顺性好.数值实验表明,该方法简单,便于应用.     

Local Discontinuous Galerkin Method for Surface Diffusion and Willmore Flow of Graphs

In this paper, we develop a local discontinuous Galerkin (LDG) finite element method for surface diffusion and Willmore flow of graphs. We prove L ² stability for the equation of surface diffusion of graphs and energy stability for the equation of Willmore flow of graphs. We provide numerical simulation results for different types of solutions of these two types of the equations to illustrate the accuracy and capability of the LDG method.     

Construction of minimal subdivision surface with a given boundary

The fascinating characters of minimal surface make it to be widely used in shape design. While the flexibility and high quality of subdivision surface make it a powerful mathematical tool for shape representation. In this paper, we construct minimal subdivision surfaces with given boundaries using the mean curvature flow, a second order geometric partial differential equation. This equation is solved by a finite element method where the finite element space is spanned by the limit functions of an extended Loop’s subdivision scheme proposed by Biermann et al. Using this extended Loop’s subdivision scheme we can treat a surface with boundary, thereby construct the perfect minimal subdivision surfaces with any topology of the control mesh and any shaped boundaries.     

An energy-preserving algorithm for nonlinear Hamiltonian wave equations with Neumann boundary conditions

In this paper, a novel energy-preserving numerical scheme for nonlinear Hamiltonian wave equations with Neumann boundary conditions is proposed and analyzed based on the blend of spatial discretization by finite element method (FEM) and time discretization by Average Vector Field (AVF) approach. We first use the finite element discretization in space, which leads to a system of Hamiltonian ODEs whose Hamiltonian can be thought of as the semi-discrete energy of the original continuous system. The stability of the semi-discrete finite element scheme is analyzed. We then apply the AVF approach to the Hamiltonian ODEs to yield a new and efficient fully discrete scheme, which can preserve exactly (machine precision) the semi-discrete energy. The blend of FEM and AVF approach derives a new and efficient numerical scheme for nonlinear Hamiltonian wave equations. The numerical results on a single-soliton problem and a sine-Gordon equation are presented to demonstrate the remarkable energy-preserving property of the proposed numerical scheme.     

Numerical solution of coupled nonlinear Schrödinger equation by Galerkin method

The coupled nonlinear Schrödinger equation models several interesting physical phenomena presents a model equation for optical fiber with linear birefringence. In this paper we derive a finite element scheme to solve this equation, we test this method for stability and accuracy, many numerical tests have been conducted. The scheme is quite accurate and describe the interaction picture clearly.     

On numerical stabilization in the solution of Saint-Venant equations using the finite element method

Solving the Saint-Venant equations by using numerical schemes like finite difference and finite element methods leads to some unwanted oscillations in the water surface elevation. The reason for these oscillations lies in the method used for the approximation of the nonlinear terms. One of the ways of smoothing these oscillations is by adding artificial viscosity into the scheme. In this paper, by using a suitable discretization, we first solve the one-dimensional Saint-Venant equations by a finite element method and eliminate the unwanted oscillations without using an artificial viscosity. Second, our main discussion is concentrated on numerical stabilization of the solution in detail. In fact, we first convert the systems resulting from the discretization to systems relating to just water surface elevation. Then, by using M-matrix properties, the stability of the solution is shown. Finally, two numerical examples of critical and subcritical flows are given to support our results.     

Locally Conservative Continuous Galerkin FEM for Pressure Equation in Two-Phase Flow Model in Subsurfaces

A typical two-phase model for subsurface flow couples the Darcy equation for pressure and a transport equation for saturation in a nonlinear manner. In this paper, we study a combined method consisting of continuous Galerkin finite element methods (CGFEMs) followed by a post-processing technique for Darcy equation and a nodal centered finite volume method (FVM) with upwind schemes for the saturation transport equation, in which the coupled nonlinear problem is solved in the framework of operator decomposition. The post-processing technique is applied to CGFEM solutions to obtain locally conservative fluxes which ensures accuracy and robustness of the FVM solver for the saturation transport equation. We applied both upwind scheme and upwind scheme with slope limiter for FVM on triangular meshes in order to eliminate the non-physical oscillations. Various numerical examples are presented to demonstrate the performance of the overall methodology.     

Finite Element Approximation to a Finite-Size Modified Poisson-Boltzmann Equation

The inclusion of steric effects is important when determining the electrostatic potential near a solute surface. We consider a modified form of the Poisson-Boltzmann equation, often called the Poisson-Bikerman equation, in order to model these effects. The modifications lead to bounded ionic concentration profiles and are consistent with the Poisson-Boltzmann equation in the limit of zero-size ions. Moreover, the modified equation fits well into existing finite element frameworks for the Poisson-Boltzmann equation. In this paper, we advocate a wider use of the modified equation and establish well-posedness of the weak problem along with convergence of an associated finite element formulation. We also examine several practical considerations such as conditioning of the linearized form of the nonlinear modified Poisson-Boltzmann equation, implications in numerical evaluation of the modified form, and utility of the modified equation in the context of the classical Poisson-Boltzmann equation.     

Numerical methods for calculating pressure distribution in gas bearings

In gas bearings, the pressure distribution is governed by a non-linear Reynolds equation. In order to solve this equation two numerical methods, the conservative difference scheme and the finite element method, are provided in this paper. They are superior to the finite difference method of Colemman [2]. Use of the finite element method is advocated because of its flexibility in solving the Reynolds equation.     

Mixed finite element method for analysis of viscoelastic fluid flow

In the present paper, numerical analysis of incompressible viscoelastic fluid flow is discussed using mixed finite element Galerkin method. Because Maxwellian viscoelasticity is assumed as the constitutive equation, stress components could not be eliminated from the governing equation system. Because of this, mixed finite element method is utilized to discretize the basic equations. For the solution procedures to solve discretized equation system, Newton-Raphson method for steady flow and perturbation method for unsteady flow is employed. As the numerical examples, comparison was made on the finite element computational results between by direct method and by mixed method. Effects of the viscoelasticity is analyzed for the flows at Reynold's numbers 30, 50 and 70.