An approach of dynamic mesh adaptation for simulating 3‐dimensional unsteady moving‐immersed‐boundary flows

In this paper, we present an approach of dynamic mesh adaptation for simulating complex 3‐dimensional incompressible moving‐boundary flows by immersed boundary methods. Tetrahedral meshes are adapted by a hierarchical refining/coarsening algorithm. Regular refinement is accomplished by dividing 1 tetrahedron into 8 subcells, and irregular refinement is only for eliminating the hanging points. Merging the 8 subcells obtained by regular refinement, the mesh is coarsened. With hierarchical refining/coarsening, mesh adaptivity can be achieved by adjusting the mesh only 1 time for each adaptation period. The level difference between 2 neighboring cells never exceeds 1, and the geometrical quality of mesh does not degrade as the level of adaptive mesh increases. A predictor‐corrector scheme is introduced to eliminate the phase lag between adapted mesh and unsteady solution. The error caused by each solution transferring from the old mesh to the new adapted one is small because most of the nodes on the 2 meshes are coincident. An immersed boundary method named local domain‐free discretization is employed to solve the flow equations. Several numerical experiments have been conducted for 3‐dimensional incompressible moving‐boundary flows. By using the present approach, the number of mesh nodes is reduced greatly while the accuracy of solution can be preserved.     

First principles investigations of Fe,Co, Ni in model honeycomb carbon networks

The behaviors of ferromagnetic transition metals of the first period: Fe, Co and Ni are examined within density functional theory calculations in two dimensional carbon extended networks using model structure LiC₆. Around geometry optimized structures, the energy-volume equations of states considering non magnetic and spin polarized configurations established ferromagnetic ground states with magnetizations –reduced with respect to the metals’– of 2 μ_B for FeC₆ and 1 μ_B for CoC₆ while no magnetic solution could be identified for NiC₆. In the D_6h point group of the P6/mmm space group l_m decomposition of the d states results with increasing energy into doublet state E_1g with d(x²-y²) and d(xy); singlet state A_1g d(z²) and doublet state E_2g d(xz) and d(yz) lying on E_F and responsible of the onset of magnetic moments. This was mirrored via molecular orbital approach with a construct of Fe embedded between two extended carbon networks thus validating the model structure proposed for TC₆ compounds. The 100% polarization in one spin channel allows proposing potential uses in spintronics applications.     

Two-dimensional vortex sheets for the nonisentropic Euler equations: Nonlinear stability

We show the short-time existence and nonlinear stability of vortex sheets for the nonisentropic compressible Euler equations in two spatial dimensions, based on the weakly linear stability result of Morando and Trebeschi (2008) [20]. The missing normal derivatives are compensated through the equations of the linearized vorticity and entropy when deriving higher-order energy estimates. The proof of the resolution for this nonlinear problem follows from certain a priori tame estimates on the effective linear problem in the usual Sobolev spaces and a suitable Nash–Moser iteration scheme.     

Decay rate for viscoelastic wave equation of variable coefficients with delay and dynamic boundary conditions

In this paper, we consider a viscoelastic wave equation of variable coefficients in the presence of past history with nonlinear damping and delay in the internal feedback and dynamic boundary conditions. Under suitable assumptions, we establish an explicit and general decay rate result without imposing restrictive assumption on the behavior of the relaxation function at infinity by Riemannian geometry method and Lyapunov functional method.     

A short note on the entrainment and exit boundary conditions

An attempt is made to find out the suitable entrainment and exit boundary conditions in laminar flow situations. Streamfunction vorticity formulation of the Navier–Stokes equations are solved by ADI method. Two‐dimensional laminar plane wall jet flow is used to test different forms of the boundary conditions. Results are compared with the experimental and similarity solution and the proper boundary condition is suggested. The Kind 1 boundary condition is recommended. It consists of zero first derivative condition for velocity variable and for streamfunction equation, mixed derivative at the entrainment and exit boundaries. Copyright © 2005 John Wiley & Sons, Ltd.     

ON A PAIR OF NON-ISOMETRIC ISOSPECTRAL DOMAINS WITH FRACTAL BOUNDARIES AND THE WEYL-BERRY CONJECTURE

ＯＮＡＰＡＩＲＯＦＮＯＮＩＳＯＭＥＴＲＩＣＩＳＯＳＰＥＣＴＲＡＬＤＯＭＡＩＮＳＷＩＴＨＦＲＡＣＴＡＬＢＯＵＮＤＡＲＩＥＳＡＮＤＴＨＥＷＥＹＬＢＥＲＲＹＣＯＮＪＥＣＴＵＲＥＳＬＥＥＭＡＮ，Ｂ．Ｄ．ＣＨＥＮＨＵＡＭａｎｕｓｃｒｉｐｔｒｅｃ...     

Efficient Collocation Schemes for Singular Boundary Value Problems

We discuss an error estimation procedure for the global error of collocation schemes applied to solve singular boundary value problems with a singularity of the first kind. This a posteriori estimate of the global error was proposed by Stetter in 1978 and is based on the idea of Defect Correction, originally due to Zadunaisky. Here, we present a new, carefully designed modification of this error estimate which not only results in less computational work but also appears to perform satisfactorily for singular problems. We give a full analytical justification for the asymptotical correctness of the error estimate when it is applied to a general nonlinear regular problem. For the singular case, we are presently only able to provide computational evidence for the full convergence order, the related analysis is still work in progress. This global estimate is the basis for a grid selection routine in which the grid is modified with the aim to equidistribute the global error. This procedure yields meshes suitable for an efficient numerical solution. Most importantly, we observe that the grid is refined in a way reflecting only the behavior of the solution and remains unaffected by the unsmooth direction field close to the singular point.     

CPL冷凝单管融化过程的双倒易边界元分析

1引言CPL（毛细泵两相流回路）能传输较高的热负荷，而且不需要循环泵、阀门等运动部件，重量轻可靠性好；在飞行器热控制方面有很好的应用前景[1]。CPL在恶劣的空间环境中运行，需要防止工作介质出现冷冻，为此，我们已经用解析[2]和数值的方法[3]进行了初步的分析，得到了一些可供工程设计参考应用的结果。另外，在恶劣的空间环境中，如果出现冻结，需要融化起动，热管的安全设计需要研究它的融化特性，以确保热管在空间能正常运行。双倒易边界元方法用Laplace基本解[4]，通过对一类偏微分方程两侧进行转化，将它全部转化为纯边界积分…