基于CFD/CSD的舵面控制气动弹性数值模拟方法

在动态网格上通过耦合求解流动控制方程和结构动力学方程, 发展了一种舵面控制下飞行器运动响应过程中气动弹性数值模拟研究方法.流动控制方程采用N-S方程, 结构动力学采用线性模态叠加方法, 其中流动控制方程空间离散采用基于非结构网格的有限体积方法, 对流通量采用计算HLLC格式, 非定常时间离散采用基于LU-SGS的双时间步长方法.模拟中, 气动运动和结构变形在双时间步长方法推进过程中采用改进松耦合方法, 气动网格与结构网格之间信息交换采用无限平板样条法实现, 飞行器的运动和变形采用基于重叠网格和Delaunay图映射变形网格相结合的方法进行处理.采用多个考核算例对发展的数值方法进行考核验证, 结果表明该方法可以高效精确模拟舵面开环控制下飞行器运动响应过程中的气动弹性特性.      

ρ-VOF:一种可清晰模拟自由界面的两相流方法

文章通过对EFM(effective field modeling)模型进行简化, 消除了原模型的非守恒性项和非双曲性特性项, 发展了一种基于密度的气液两相流模拟方法: ρ-VOF方法.利用体积分数信息对控制单元内的自由界面进行重构, 得到了控制单元内流体的空间分布, 并采用AUSM+-up格式获得考虑气液流体接触间断信息的对流通量.新方法可统一处理激波间断和接触间断的相互作用, 保持自由界面的尖锐性, 并且其计算量与自由界面的空间复杂度无关.最后, 数值模拟了液体激波管气液激波管和气体激波跨二维液滴传播等问题, 并与文献结果进行对比, 验证了本方法在气液两相流模拟中的准确性.      

基于非结构变形网格的间断装配法原理

在激波捕捉法计算得到的流场基础上采用辨识算法得到初始间断位置, 从ALE方程出发, 考虑离散几何守恒律, 采用变形网格和网格重构技术解决计算过程中间断运动和变形, 新旧网格之间流场采用高精度信息传递方法保持时间精度, 建立了基于非结构动网格技术的间断装配方法.通过激波管问题的二维模拟, 模拟了初始间断分解为激波和接触间断激波遇到固壁反射后与接触间断相交的非定常流动过程, 对这种新方法的基本原理进行了介绍.      

移动网格与Level Set耦合方法在气-液两相流数值模拟中的应用

本文采用自适应移动网格与Level Set函数相耦合的方法来实现气-液两相流的数值模拟与计算.作为自适应网格方法的一种,移动网格方法主要是为了解决发展方程的计算问题而设计的方法.文中给出了移动网格的生成方程,并针对方程的非线性,给出了一种半隐式的离散方法用于进行求解.本文将移动网格方法与Level Set方法相耦合,将控制流体运动的Navier-Stokes方程以及追踪相界面的Level Set方程转换到曲线坐标下,应用一套曲线坐标方程组来同时描述气、液两相流的运动规律,成功实现了对气-液两相流问题的数值模拟.通过对顶盖驱动流的计算以及对液滴沉降现象的模拟计算,验证了本文方法的可靠性.本文对常重力与微重力下两气泡融合的发展规律进行了数值模拟,通过分析对比,得到了重力对两气泡融合变形的影响规律.     

自然对流换热的高精度非结构网格有限体积法

用非结构网格有限体积法求解自然对流换热时,传统的对流项离散格式难以兼顾数值精度与计算效率,我们发展了一种耦合高精度格式的延迟修正方法,用于对流项的离散.高Re数下方腔驱动流数值计算验证了该方法具有较高的计算精度和较好的稳定性.Boussinesq流体的自然对流换热数值模拟,表明该方法能有效克服高Ra数时数值计算发散,可准确捕捉自然对流换热问题中不同偏心率下的等温线和流线分布特征.     

梯度增量level set/虚拟网格法模拟波浪结构物相互作用

针对波浪结构物相互作用的数值模拟问题,本文结合虚拟网格法和梯度增量level set (GALS)法发展了一个基于直角网格的多相流计算方法,其有能力模拟多个动物体与强非线性自由表面的相互作用。该模型采用时间半隐式有限差分法在固定直角网格上求解不可压缩两相流Navier-Stokes方程。采用虚拟网格浸入边界法通过在物体内部布置虚拟网格以考虑浸入边界对流场的影响。GALS方法捕捉强非线性自由表面如波浪翻卷和射流,并通过同时演化level set (LS)方程和其梯度向量方程以提高经典LS方法的精度。基于提出的数值方法,模拟了一系列典型算例,包括交错布置两圆柱相向运动,横向激励下的液舱共振晃荡和圆柱常速入水问题。本文结果与数值和试验结果吻合较好,从不同方面验证了本文计算方法的精度和可靠性。     

水陆两栖飞机静水面高速滑行性能数值计算与试验分析

针对水陆两栖飞机静水面高速滑行过程的运动响应大、流场强非线性等问题，提出了一种基于传统动网格技术的"状态预估——精确计算"的数值模拟方法:通过求解Reynolds平均N-S方程结合运动方程来模拟飞机静水面滑行时的流场特征和运动特性，数值模拟方法为隐式有限体积法，湍流模型采用k-ω（SST Menter）结合壁函数进行处理，自由液面捕捉采用VOF方法；数值计算时，首先采用粗网格对简化后的飞机在不同航速下的姿态和升沉进行快速预估，再将飞机置于预估状态下进行精确网格划分，最后进行精确数值计算分析.为了验证数值模拟结果的正确性，在物理水池中进行了静水拖曳试验，将数值计算结果与试验结果进行对比分析可得:数值计算与水池试验的流场特征吻合，且阻力、姿态和升沉的计算精度达到90%，验证了数值模拟方法的可行性.      

基于广义交替数值通量的局部间断Galerkin方法求解二维波动方程

波的传播往往在复杂的地质结构中进行,如何有效地求解非均匀介质中的波动方程一直是研究的热点.本文将局部间断Galekin(local discontinuous Galerkin, LDG)方法引入到数值求解波动方程中.首先引入辅助变量,将二阶波动方程写成一阶偏微分方程组,然后对相应的线性化波动方程和伴随方程构造间断Galerkin格式;为了保证离散格式满足能量守恒,在单元边界上选取广义交替数值通量,理论证明该方法满足能量守恒性.在时间离散上,采用指数积分因子方法,为了提高计算效率,应用Krylov子空间方法近似指数矩阵与向量的乘积.数值实验中给出了带有精确解的算例,验证了LDG方法的数值精度和能量守恒性;此外,也考虑了非均匀介质和复杂计算区域的计算,结果表明LDG方法适合模拟具有复杂结构和多尺度结构介质中的传播.     

复杂地形下带自由表面水流的分离涡模拟

将分离涡模型(DES),即一种RANS和LES的混合模型,应用于带自由表面的地表水流运动,建立一套数值仿真模型.模型基于有限体积法,水平面内采用非结构计算网格,垂向为结构化网格,对流项离散格式采用二阶TVD格式,并行基于Open MP语言库.算例表明DES模型有助于揭示复杂地形条件下带自由表面水流的大涡拟序结构.     

Stokes流动中液滴变形的数值模拟

数值模拟研究Stokes流动中的液滴变形问题。采用边界积分方法计算液滴表面的运动,建立在网格能量最小化基础上的自适应表面离散及网格局部重建方法以准确的分辨出变形液滴的形状,并将数值模拟结果与实验结果进行了比较,模拟结果表明,本方法能够较好的模拟液滴变形的动力学特性。     

视频分镜头边界检测的一种新方法

定义了用于渐变边界检测的累积帧差,提出了一种基于高斯模型的镜头边界检测新方法,对突变边界和渐变边界分别采用不同参数的高斯模型,实现了阈值的自适应选择,克服了现有的双阈值比较法中下边界阈值难以自适应设定的困难.初步实验结果表明该方法能准确地检测突变和渐变边界.     

The wavelet response as a multiscale NDT method

We analyze interfaces by using reflected waves in the framework of the wavelet transform. First, we introduce the wavelet transform as an efficient method to detect and characterize a discontinuity in the acoustical impedance profile of a material. Synthetic examples are shown for both an isolated reflector and multiscale clusters of nearby defects. In the second part of the paper we present the wavelet response method as a natural extension of the wavelet transform when the velocity profile to be analyzed can only be remotely probed by propagating wavelets through the medium (instead of being directly convolved as in the wavelet transform). The wavelet response is constituted by the reflections of the incident wavelets on the discontinuities and we show that both transforms are equivalent when multiple scattering is neglected. We end this paper by experimentally applying the wavelet response in an acoustic tank to characterize planar reflectors with finite thicknesses.     

A numerical study of three-dimensional liquid sloshing in tanks

A numerical model NEWTANK (Numerical Wave TANK) has been developed to study three-dimensional (3-D) non-linear liquid sloshing with broken free surfaces. The numerical model solves the spatially averaged Navier–Stokes equations, which are constructed on a non-inertial reference frame having arbitrary six degree-of-freedom (DOF) of motions, for two-phase flows. The large-eddy-simulation (LES) approach is adopted to model the turbulence effect by using the Smagorinsky sub-grid scale (SGS) closure model. The two-step projection method is employed in the numerical solutions, aided by the Bi-CGSTAB technique to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate volume-of-fluid (VOF) method is used to track the distorted and broken free surface. Laboratory experiments are conducted for both 2-D and 3-D non-linear liquid sloshing in a rectangular tank. A linear analytical solution of 3-D liquid sloshing under the coupled surge and sway excitation is also developed in this study. The numerical model is first validated against the available analytical solution and experimental data for 2-D liquid sloshing of both inviscid and viscous fluids. The validation is further extended to 3-D liquid sloshing. The numerical results match with the analytical solution when the excitation amplitude is small. When the excitation amplitude is large where sloshing becomes highly non-linear, large discrepancies are developed between the numerical results and the analytical solutions, the former of which, however, agree well with the experimental data. Finally, as a demonstration, a violent liquid sloshing with broken free surfaces under six DOF excitations is simulated and discussed.     

Evaluation of the temporal structure of postural sway fluctuations based on a comprehensive set of analysis tools

The analysis of postural control has a long history. Traditionally, the amount of body sway is solely used as an index of postural stability. Although this leads to some extent to an effective evaluation of balance performance, the control mechanisms involved have not yet been fully understood. The concept of nonlinear dynamics suggests that variability in the motor output is not randomness but structure, providing the stimulus to reveal the functionality of postural sway. The present work evaluates sway dynamics by means of COP excursions in a quiet standing task versus a dual-task condition in three different test times (30, 60, 300 s). Besides the application of traditional methods–which estimate the overall size of sway–the temporal pattern of body sway was quantified via wavelet transform, multiscale entropy and fractal analysis. We found higher sensitivity of the structural parameters to modulations of postural control strategies and partly an improved evaluation of sway dynamics in longer recordings. It could be shown that postural control modifications take place on different timescales corresponding to the interplay of the sensory systems. A continued application of nonlinear analysis can help to better understand postural control mechanisms.     

Scattering by an E plane dielectric sheet in rectangular waveguide

A method for the computation of S-parameters associated with a rectangular waveguide with a rectangular dielectric obstacle is presented. The method uses the Method of Lines and discretization of the dielectric permittivity to treat the dielectric discontinuity along the discretization direction. In this way, the Method of Lines can be applied to analyze waveguide discontinuity with dielectric inhomogenousity along the discretization direction. The numerical results have been verified by both theoretical and experimental results.     

A resolution of some orthogonality and expansion-set problems in many-body scattering theory

A modification to the coupled channel equations for many-body scattering is introduced, based on the evaluation of the discontinuity equation for the matrix of transition operators. The resulting (modified) equations are shown to be in a form which resolves some problems associated with continuum states in an analysis of rearrangement collisions using the eigenstate expansion method.     

Suppression of reverberation distortions of a receiver signal using the water tank transfer function

The paper considers a method for suppressing the reverberant distortions of an underwater sound receiver signal during receiver calibration in a laboratory water tank. The method is based on using the water tank transfer function, which is a complex frequency-dependent coefficient that establishes, for the point of signal reception, the relation between the sound pressures in the reverberant sound field of the water tank and in the free sound field. The procedure for experimentally obtaining the water tank transfer function is considered. Examples of suppressing reverberant distortions during noise and pulse sound reception are presented.     

Centrifugal instability in rotating shallow water and the problem of the spiral structure in galaxies

A new instability predicte by theory to occur in rotating shallow water in which the rotation velocity has a discontinuity, in a regime where the flow velocity exceeds the characteristics velocity of the waves, has been found experimentally. The instability develops when the radial gradient of the angular velocity across the discontinuity is negative; such an instability is likely to be responsible for the formation of the spiral structure in galaxies which have a similar rotational velocity profile.     

A new method for obtaining exact analytical formulae for the roots of transcendental functions

A new method is proposed for the derivation of closed-form formulae for the zeros and poles of sectionally analytic functions in the complex plane. This method makes use of the solution of the simple discontinuity problem in the theory of analytic functions and requires the evaluation of real integrals only (for functions with discontinuity intervals along the real axis). Many transcendental equations of mathematical physics can be successfully solved by the present approach. An application to such an equation, the molecular field equation in the theory of ferromagnetism, is made and the corresponding analytical formulae are reported together with numerical results.     

Hysteresis magnetoresistance and micromagnetic modeling of Ni microbelts

Magnetization reversal process in a single nickel microbelt is investigated by measuring magnetoresistance (MR). A hysteresis MR curve is observed, which consists of a reversible bell-shaped curve under high magnetic field and an irreversible discontinuity under low magnetic field. A micromagnetic simulation is performed for understanding the underlying physics behind the MR curve. The hysteresis is ascribed to the competition of exchange energy and Zeeman energy, and the irreversible discontinuity is attributed to the formation of domain walls and vortices.