计算机图形学中3D打印研究进展

3D打印是一种加法加工制造技术,它改变了传统的减式加工技术,实现了许多传统加工制造技术无法完成的任务.3D打印不仅为给计算机图形学提出了新的研究内容与挑战,也有力地推动了计算机图形学相关技术的发展.文中首先介绍3D打印技术相对于传统制造技术的优势同时,简要介绍3D打印技术的分类、原理与过程.在此基础上,对3D打印在计算机图形学的研究进展进行综述,包括:结构强度分析、大物体打印、打印方向优化、自平衡优化、支撑结构设计与优化、降低打印成本、提高打印效率和质量等.最后对3D打印在计算机图形学领域的发展情况进行展望.     

一种求解运动曲面上对流扩散方程的三维水平集方法

给出了一种求解运动曲面上对流扩散方程的三维水平集算法. 水平集函数被用来表示曲面.曲面上的微分方程及其解通过水平集方法被延拓到包含曲面的一个小邻域中. 一种半隐式的Crank-Nicholson 格式被用来做时间推进, 中心差分和三阶加权实质无振荡(WENO) 格式被分别用来离散方程中的扩散项和对流项. 分析证明了它在标准的Courant-Friedrichs-Lewy (CFL) 条件下的稳定性. 数值算例显示了它能取得二阶精度.     

A regularity criterion for the solutions of 3D Navier-Stokes equations

In terms of two partial derivatives of any two components of velocity fields, we give a new criterion for the regularity of solutions of the Navier-Stokes equation in R³. More precisely, let u=(u¹,u²,u³) be a weak solution in (0,T)×R³. Then u becomes a classical solution if any two functions of _∂1u¹, _∂2u² and _∂3u³ belong to Lθ(0,T;Lr(R³)) provided with , .     

Improved 3D Surface Reconstruction via the Method of Fundamental Solutions

A new mathematical model of the modified bi-Helmholtz equation is proposed for the reconstruction of 3D implicit surfaces using the method of fundamental solutions. In the algorithm, we also show how to properly determine the parameter so that the spurious surface can be avoided. The main attraction of the proposed method is its simplicity. Four examples for the surface reconstruction are presented to validate the proposed numerical model.     

3D打印中的优化设计

3D打印技术的兴起在许多领域掀起了新的研究热点,《3D打印中的优化设计》应用课题针对几何处理领域中的许多相关问题做了系统性的研究,其中包括在3D打印中通过几何模型的结构优化以节省打印材料和打印时间的问题,三维几何模型保持特征的去噪和构建问题,以及三维几何模型的序贯重建和多余分支去除问题.相关问题中多涉及含有复杂约束的优化问题,以及压缩感知和稀疏优化问题.本文简要总结了上述问题的研究方法,并从几何优化的角度阐述了其中运筹学原理与方法的运用,表明了运筹学工具在相关领域研究中的重要性和有效性.     

三维PTT黏弹性液滴撞击固壁面问题的改进SPH模拟

基于光滑粒子动力学(smoothed particle hydrodynamics, SPH)方法,对三维Phan-Thien Tanner(PTT)黏弹性液滴撞击固壁面问题进行了数值模拟.为了有效地防止粒子穿透固壁,且缩减三维数值模拟所消耗的计算时间,提出了一种适合三维数值模拟的改进固壁边界处理方法.为了消除张力不稳定性问题,采用一种简化的人工应力技术.应用改进SPH方法对三维PTT黏弹性液滴撞击固壁面问题进行了数值模拟,精细地捕捉了液滴在不同时刻的自由面,讨论了PTT黏弹性液滴不同于Newton(牛顿)液滴的流动特征,分析了PTT拉伸参数对液滴宽度、高度和弹性收缩比等的影响.模拟结果表明,改进SPH方法能够有效而准确地描述三维PTT黏弹性液滴撞击固壁面问题的复杂流变特性和自由面变化特征.     

Numerical simulation of standing wave with 3D predictor-corrector finite difference method for potential flow equations

A three-dimensional (3D) predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa-tions with a free surface. The 3D irreg...     

3D Anisotropic Diffusion on GPUs by Closed-Form Local Tensor Computations

We present an efficient implementation of volumetric anisotropic image diffusion filters on modern programmable graphics processing units (GPUs), where the mathematics behind volumetric diffusion is effectively reduced to the diffusion in 2D images. We hereby avoid the computational bottleneck of a time consuming eigenvalue decomposition in $\mathbb{R}^3$. Instead, we use a projection of the Hessian matrix along the surface normal onto the tangent plane of the local isodensity surface and solve for the remaining two tangent space eigenvectors. We derive closed formulas to achieve this and prevent the GPU code from branching. We show that our most complex volumetric anisotropic diffusion filters gain a speed up of more than 600 compared to a CPU solution.     

三维泊松方程基于Richardson外推法的高阶紧致差分方法

基于Richardson外推法提出了数值求解三维泊松方程的高阶紧致差分方法.方法通过利用四阶和六阶紧致差分格式,分别在细网格和粗网格上求解,然后利用Richardson外推技术和算子插值方法,得到三维泊松方程在细网格上的六阶和八阶精度的数值解.数值实验结果验证了该方法的精确性和有效性.     

An Acceleration Technique for the Augmented IIM for 3D Elliptic Interface Problems

A new fast algorithm based on the augmented immersed interface method and a fast Poisson solver is proposed to solve three dimensional elliptic interface problems with a piecewise constant but discontinuous coefficient. In the new approach, an augmented variable along the interface, often the jump in the normal derivative along the interface is introduced so that a fast Poisson solver can be utilized. Thus, the solution of the Poisson equation depends on the augmented variable which should be chosen such that the original flux jump condition is satisfied. The discretization of the flux jump condition is done by a weighted least squares interpolation using the solution at the grid points, the jump conditions, and the governing PDEs in a neighborhood of control points on the interface. The interpolation scheme is the key to the success of the augmented IIM particularly. In this paper, the key new idea is to select interpolation points along the normal direction in line with the flux jump condition. Numerical experiments show that the method maintains second order accuracy of the solution and can reduce the CPU time by 20-50%. The number of the GMRES iterations is independent of the mesh size.     

A high order finite difference method with Richardsonextrapolation for 3D convection diffusion equation

In this paper, we extend the Sun and Zhang’s [24] work on high order finite difference method, which is based on the Richardson extrapolation technique and an operator interpolation scheme for the one and two dimensional steady convection diffusion equations to the three dimensional case. Firstly, we employ a fourth order compact difference scheme to get the fourth order accurate solution on the fine and the coarse grids. Then, we use the Richardson extrapolation technique by combining the two approximate solutions to get a sixth order accurate solution on coarse grid. Finally, we apply an operator interpolation scheme to achieve the sixth order accurate solution on the fine grid. During this process, we use alternating direction implicit (ADI) method to solve the resulting linear systems. Numerical experiments are conducted to verify the accuracy and effectiveness of the present method.     

H1-superconvergence of finite difference method based on Q1-element on quasi-uniform mesh for the 3D Poisson equation

In this paper, the finite difference (FD) method is considered for the 3D Poisson equation by using the Q₁-element on a quasi-uniform mesh. First, under the regularity assumption of , the H¹-superconvergence of the FD solution u_h based on the Q₁-element to the first-order interpolation function is obtained. Next, the H¹-superconvergence of the second-order interpolation postprocessing function based on the FD solution u_h to u is provided. Finally, numerical tests are presented to show the H¹-superconvergence result of the FD postprocessing function to u if .     

基于三维弹性问题的多极边界元法核函数分解

多极边界元法已经成功地应用于大规模工程计算中.得到并且证明了基于三维弹性问题的多极边界元法核函数分解的定理(定理1),完善了多击边界元法的数学理论.     

An $L^{\infty}$ Second Order Cartesian Method for 3D Anisotropic Interface Problems

A second order accurate method in the infinity norm is proposed for general three dimensional anisotropic elliptic interface problems in which the solution and its derivatives, the coefficients, and source terms all can have finite jumps across one or several arbitrary smooth interfaces. The method is based on the 2D finite element-finite difference (FE-FD) method but with substantial differences in method derivation, implementation, and convergence analysis. One of challenges is to derive 3D interface relations since there is no invariance anymore under coordinate system transforms for the partial differential equations and the jump conditions. A finite element discretization whose coefficient matrix is a symmetric semi-positive definite is used away from the interface; and the maximum preserving finite difference discretization whose coefficient matrix part is an M-matrix is constructed at irregular elements where the interface cuts through. We aim to get a sharp interface method that can have second order accuracy in the point-wise norm. We show the convergence analysis by splitting errors into several parts. Nontrivial numerical examples are presented to confirm the convergence analysis.     

A remark on the Beale-Kato-Majda criterion for the 3D MHD equations with zero kinematic viscosity

In this paper,we study the blow-up criterion of smooth solutions to the 3D magneto-hydrodynamic system in ˙ B 0 ∞,∞.We show that a smooth solution of the 3D MHD equations with zero kinematic viscosity in the whole space R 3 breaks down if and only if certain norm of the vorticity blows up at the same time.     

Block SSOR preconditionings for high-order 3D FE systems. III: Incomplete BSSOR preconditionings based on p-partitionings

This paper considers the construction of high quality IBSSOR preconditionings for p-adaptive hierarchical high-order 3D FE systems. Instead of the commonly-adopted domain decomposition approach to the selection of a block partitioning underlying the corresponding block SSOR preconditioning, we introduce the so-called p-partitionings based on subdividing the degrees of freedom in accordance with the levels of p-refinement. A theoretical analysis of IBSSOR preconditionings is presented. The results of numerical experiments with 3D FE systems arising from the hierarchical approximations of order p, 3 ≤ p ≤ 5, of the 3D equilibrium equations for linear elastic orthotropic materials show that the IBSSOR preconditionings suggested are more efficient than the IBSSOR preconditionings based on domain decomposition partitionings. Moreover, we demonstrate that the new preconditionings are more reliable when solving extremely ill-conditioned 3D FE systems. © 1997 John Wiley & Sons, Ltd.