code2

自己写了一个二维Euler方程的间断有限元程序 上次发了一个三角形单元的程序 因为不是曲边单元 所以在圆柱后面容易形成涡 现在把程序改为曲边四边形单元了 没有涡出现 单元是8节点四边形单元 节点编号顺序是 1 5 2 6 3 7 4 8 也就是四个角点依次 是1 2 3 4 然后是边的中点编号 5 6 7 8. 时间推进采用 Runge-Kutta 方法 数值通量采用全局Lax-Friedrichs通量 仍然不能捕捉激波 因为没有做重构或者加人工粘性 等这个做出来了 再发一次。 程序没有进行优化 比如说内存的消耗没有优化 比如直边单元的边界积分仍然采 用了曲边的积分方法 增加了计算量 比如面积分、线积分都是采用的是Gauss- Legendre-Lobatto积分 积分精度会比一般的Gauss-Legendre积分精度低一阶 等 等问题。 二维的 纯属交流性质 就没有考虑这些问题 ^_^ 如果物面全部是直边 那么只要改变一个参数N 就可以获得不同的计算精度 且具 有谱精度 因为单元的节点是Gauss-Legendre-Lobatto积分点。 其实就是谱元法 （物面是曲边的情况我不清楚是不是也可以通过提高基函数的阶数 也就是增加N 来提高计算精度）(Wrote a two-dimensional Euler equations with discontinuous finite element program Last made ​ ​ a triangular element of the program, not curved edge unit is so easy to form a vortex in the cylinder behind the Program to curved edge quadrilateral element vortices appear Unit is the order of 8-node quadrilateral element node number is 15,263,748 which is the four corners of the points in turn Is 1234 and then the side of the midpoint of the number 5678. Time promote the use of Runge-Kutta method Numerical flux of the overall situation of Lax-Friedrichs, flux Still can not capture the shock wave did not do the reconstruction or artificial viscosity do it Zaifayici. The program is not optimized for example, memory consumption is not optimized such as straight-edge boundary integral of the unit is still mining Integral method to increase the amount of computation such as surface integral with a curved edge, the line integral using the Gauss- The Legendr)

说明：  本书详细介绍了格子波尔兹曼法的原理，边界条件，具体模型及数值方法，非常适合学习与应用。(This book introduce the theory of lattice boltzmann methods alone with it's boundary conditions,models and numerical methods. It's really a book for us to learn and use.)

速度曲线规划，在给定长度path 曲线上，根据各waypoint的速度约束，规划出平滑的速度曲线(speed curve planning.In a given length path curve, according to the speed of the waypoint constraints, planning a smooth speed curve.)

二维交错网格SIMPLE算法fortran计算程序，包括主程序和用户定义程序。(The two-dimensional to SIMPLE algorithm the fortran calculation procedures staggered grid, including the main program and the user-defined procedures.)

这是一个数组练习，程序的作用是查询一个一维数组（数值型数组）中每一个元素出现的次数，且元素出现的顺序不能改变(This is an array of exercises, the program' s role is to query a one-dimensional array (numeric array) in the number of occurrences of each element, and the element can not be changed in order of appearance)

使用惩罚函数法对机械设计中的问题进行优化处理(The use of penalty function method for mechanical design optimization problems deal with)

以前用VB6编写的一个计算器，支持常用的数学函数，支持弧度、角度转换，支持简单的函数作图。(A calculator writen by VB6, supports common mathematical functions, support radian, angle conversion, support simple function mapping.)

Fluent 中短语模拟水的受热蒸发相变的UDF程序代码(the UDF for phase change of water)

计算典型相关算法，求出最大投影向量,以及最大投影向量对应的相关系数 (calculate canonical correlations, find the largest correlation coefficients)

用MATLAB实现模拟立方体的重构,是计算机视觉中三维重构的重要的应用实验(Simulating cube reconstruction with MATLAB is an important application experiment of 3D reconstruction in computer vision.)