PDF of distance function for level‐set flamelet library modelling

The difference between a presumed distribution of flamelet position and a numerically simulated distribution of distance function (a signed distance to flamelet) is investigated. It is shown that even if the distribution of flamelet position is symmetrical and close to Gaussian, the distribution of distance function away from the mean flame position is skewed towards the mean position and the mean of the distance function is also different from the distance to the mean position. The difference depends on the distance to the mean flame and the flame wrinkling amplitude. An extension method for the variance of the distance function and an upwind scheme for solving the re‐initialization equation are presented. A numerical simulation of a premixed turbulent flame is compared to experimental data. Copyright © 2003 John Wiley & Sons, Ltd.     

凹坑形非光滑表面微观摩擦磨损试验回归分析

采用试验优化设计理论,利用L16(215)二水平正交表设计试验,对具有凹坑形非光滑表面形态的试件进行微观摩擦磨损试验研究,运用多元线性回归正交设计处理试验数据,得到具有凹坑形非光滑表面形态试件的体积磨损率和摩擦系数与各试验因素关系(如非光滑单元体大小、距离、速度和载荷)的回归方程.结果表明:试验因素对体积磨损率的影响程度由大到小依次为载荷、速度、单元体大小和距离;体积磨损率随着距离、速度及载荷的增加而增大,随着非光滑单元体尺寸的增加而减小,但距离和速度的交互影响很小;试验因素对摩擦系数的影响程度由大到小依次为载荷、速度、距离和单元体大小;摩擦系数随着距离、速度及载荷的增加而减小,随着非光滑单元体尺寸的增加而增大,但距离、速度及载荷的交互作用也很小.     

初始压力和狭缝宽度对毫米量级狭缝内爆轰起爆距离的影响

为获得狭缝内爆轰起爆距离的影响因素,分别在高度为1.0 mm,宽度为10、20、30 mm的狭缝爆轰管内,对不同初始压力下(p₀=5.0~50.0 kPa)等当量比的乙烯/氧气预混气体进行了单次爆轰性能实验研究。根据烟迹法、高速摄影图片判定起爆位置,分析了初始压力和狭缝宽度对爆轰起爆距离的影响。结果表明:(1)p₀=21.0~30.0 kPa时,起爆距离随着狭缝宽度的增大而逐渐缩短;(2)p₀=35.0~42.5 kPa时,起爆距离随着狭缝宽度的增大先缩短后增大,在p₀=45.0~50.0 kPa时起爆距离随着狭缝宽度的增大基本保持不变;(3)3种狭缝宽度下,量纲一起爆距离随量纲一初始压力的变化曲线差异较大。     

A stochastic model for elasticity tensors with uncertain material symmetries

In this paper, we consider the probabilistic modeling of media exhibiting uncertainties on material symmetries. More specifically, we address both the construction of a stochastic model and the definition of a methodology allowing the numerical simulation (and consequently, the inverse experimental identification) of random elasticity tensors whose mean distance (in a sense to be defined) to a given class of material symmetry is specified. Following the eigensystem characterization of the material symmetries, the proposed approach relies on the probabilistic model derived in Mignolet and Soize (2008), allowing the variance of selected eigenvalues of the elasticity tensor to be partially prescribed. In this context, a new methodology (regarding in particular the parametrization of the model) is defined and illustrated in the case of transversely isotropic materials. The efficiency of the approach is demonstrated by computing the mean distance of the random elasticity tensor to a given material symmetry class, the distance and projection onto the space of transversely isotropic tensors being defined by considering the Riemmanian metric and the Euclidean projection, respectively. It is shown that the methodology allows the above distance to be (partially) reduced as the overall level of statistical fluctuations increases, no matter the initial distance of the mean model used in the simulations. A comparison between this approach and the initial nonparametric approach introduced in Soize (2008) is finally provided.     

Uncertainty analysis of deflagration-to-detonation run-up distance

Three methods were adopted to estimate the deflagration-to-detonation run-up distance in a smooth tube, which comprise (1) the measurement of the propagation speed of pressure or combustion waves compared with the C–J detonation speed, (2) the time of the onset of detonation by the emission of visible light and the trajectory of pressure wave or combustion waves, and (3) the trajectory intersection with the presence of retonation wave. A nonstationary cross-correlation technique was applied to evaluate the uncertainty in estimating the run-up distance. Evaluation of the pressure wave (pressure wave speed or the pressure wave trajectory) appears to be more suitable to determine the deflagration-to-detonation run-up distance. PACS 47.40.-x Communicated by Z. Jiang     

空孔与运动裂纹相互作用的动焦散线实验研究

为研究预制裂纹不同偏移距离时运动裂纹与空孔的相互作用规律,采用动态焦散线实验系统,将预制裂纹的偏移距离设定为唯一变量,对含空孔的有机玻璃(PMMA)试件进行冲击三点弯实验。研究表明,存在两个临界距离:（6 mm (2 R )、9 mm (3 R )）,在该偏移距离下,裂纹扩展轨迹、动态断裂特性发生显著改变:(1) 预制裂纹偏移距离不大于3 mm时,裂纹贯穿空孔,发生二次起裂,且二次起裂的速度与应力强度因子显著大于一次起裂,无偏移时裂纹轨迹的分形维数为最小值;(2) 偏移距离增大至6 mm时,裂纹不再贯穿空孔,空孔对裂纹先吸引后排斥,裂纹速度与应力强度因子先减小后增大,裂纹轨迹的分形维数达到最大值;(3) 偏移距离大于6 mm时,空孔对裂纹的吸引作用逐渐减小,大于9 mm后,空孔对裂纹的吸引不再显著,裂纹起裂后即向落锤加载方向扩展直至贯穿试件。     

A Priori Assessment of Prediction Confidence for Data-Driven Turbulence Modeling

Although Reynolds-Averaged Navier–Stokes (RANS) equations are still the dominant tool for engineering design and analysis applications involving turbulent flows, standard RANS models are known to be unreliable in many flows of engineering relevance, including flows with separation, strong pressure gradients or mean flow curvature. With increasing amounts of 3-dimensional experimental data and high fidelity simulation data from Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS), data-driven turbulence modeling has become a promising approach to increase the predictive capability of RANS simulations. However, the prediction performance of data-driven models inevitably depends on the choices of training flows. This work aims to identify a quantitative measure for a priori estimation of prediction confidence in data-driven turbulence modeling. This measure represents the distance in feature space between the training flows and the flow to be predicted. Specifically, the Mahalanobis distance and the kernel density estimation (KDE) technique are used as metrics to quantify the distance between flow data sets in feature space. To examine the relationship between these two extrapolation metrics and the machine learning model prediction performance, the flow over periodic hills at Re = 10595 is used as test set and seven flows with different configurations are individually used as training sets. The results show that the prediction error of the Reynolds stress anisotropy is positively correlated with Mahalanobis distance and KDE distance, demonstrating that both extrapolation metrics can be used to estimate the prediction confidence a priori. A quantitative comparison using correlation coefficients shows that the Mahalanobis distance is less accurate in estimating the prediction confidence than KDE distance. The extrapolation metrics introduced in this work and the corresponding analysis provide an approach to aid in the choice of data source and to assess the prediction performance for data-driven turbulence modeling.     

并列仿生鱼自主游动的数值模拟研究

采用自适应浸没边界方法,数值模拟了两条并列仿生鱼在不同间距下的同相位和反相位自主游动。通过与相同游动参数下单条仿生鱼自主游动的比较,可以得到结论：（1）并列自主游动的速度小于单条鱼游动的速度,并且随着间距减小而减小。相同间距下,反相位并列游动时的速度大于同相位;（2）并列游动时的阻力均大于单条鱼,且随着间距减小而增大;（3）反相位游动时功率消耗随间距减小而增大,但同相位游动则相反;（4）当间距大于0.5个身长时,推进效率都略有增大。所以,在综合考虑了游动速度和推进效率两方面的因素后认为,鱼群中并列相邻的两条鱼应当以反相位摆动,且侧向间距保持在0.5个身长以上,这与自然界中观察的结果是一致的。     

Interaction between two spheres placed in tandem arrangement in steady and pulsating flow

The interaction among two spheres in tandem formation are studied for a Reynolds number of 300 using both steady and pulsating inflow conditions. The purpose is to further investigate the force characteristics as well as the shedding patterns of the two spheres as the separation distance is changed from 1.5 to 12 sphere diameters. The method used for the simulations is the volume of solid (VOS) method, an approach based on the volume of fluid (VOF) method. Comparisons with other computational methods have shown VOS to accurately resolve the flow field around solid spheres. The results show that the separation distance plays a significant role in changing the flow patterns and shedding frequencies at moderate separation distances, whereas effect on drag is observed even at a separation distance of 12 diameters.     

基于偏转距离近似模型的动能撞击小行星防御任务脉冲轨道优化研究

小行星撞击对地球上的生命存在重大潜在威胁,动能撞击是目前最易实现且成熟度最高的防御方案.动能撞击任务的一种轨道优化指标为最大化偏转距离(即小行星被偏转前后近地距的改变量),若用数值积分的方法精确计算偏转距离, 会导致优化效率较低.在动能撞击任务的设计初期, 可以对动力学模型及偏转距离的计算方法进行简化,以提升优化效率. 本文首先将高精度模型简化为二体模型,分析了两种经典偏转距离解析模型的适用条件,同时提出一种基于近地点时刻预估的偏转距离近似模型; 考虑运载约束,将化学推进变轨简化为脉冲推力变轨,建立了直接转移(两脉冲及三脉冲)和行星借力飞行转移(单次及两次借力)的动能撞击轨道优化模型,利用遗传算法求解了优化问题. 以偏转小行星Apophis为例, 相比于解析模型,验证了本文提出的近似模型可以同时提升最优性、降低求解复杂性. 优化结果表明,三脉冲直接转移方案与两脉冲直接转移方案的最优偏转效果基本一致,借力飞行转移方案相比于直接转移方案对偏转距离的提升效果并不明显.在动能撞击任务的前期设计中, 可以基于二体模型进行防御效果的快速评估,虽然对计算偏转距离存在一定误差, 但对防御窗口的优化结果影响不大. 进一步,数值求解偏转距离时, 可通过引入主要引力摄动项(金星、地球、木星)修正二体模型,使其与高精度模型之间的求解误差在1%以下.      

UNIQUENESS OF SOLUTION OF FIELD POINT OF SINGULAR SOURCE OUTSIDE-REGION-DISTRIBUTION METHOD

IntroductionAmethod,inwhichaproblemissolvedbysingularsource(s)outside-region-distribution(s),wheretheboundaryconditionoftheproblemissatisfied,iscalled"SingularSourceOutside-Region-DistributionMethod"(SSORDM).Usnallynon-singularintegralequationisobtainedbysingularsource(s)outSide-region-distribunonconhnuously.Asanearlyuse,singUlarsourcedistribuhoncanbefoundinthemethodofsource-sinkdistributioninfluidmechanicsasisstatedin[I],however,noproblemwasreducedtoanintegralequationatthattime.Asforsol…     

Uniquenes of solution of field point of singular source outside-region-distribution method

The uniqueness of solution of field point, inside a convex region due to singular source(s) with kernel function decreasing with distance increasing, outsidergeion-distribution(s) such that the boundary condition expressed by the response of the source(s) is satisfied, is proved by using the condition of kernel function decreasing with distance increasing and an integral inequality. Examples of part ofthese singular sources such as Kelvin’s point force, Point-Ring-Couple (PRC) etc.are given. The proof of uniqueness of solution of field point in a twisted shaft of revolution due to PRC distribution is given as an example of application.     

一种新型危险品仓库结构设计及其安全距离

安全距离是危险品仓库建设和研究中重点关注的问题之一。为减小危险品仓库的安全距离,结合现行规范和危险品仓库建设现状,针对一种由浅埋式库房主体、顶部堆土和钢筋混凝土分配板组成的新型危险品仓库形式开展了3组缩尺模型野外爆炸试验,记录了各组试验的爆炸过程,统计了冲击波超压峰值和爆炸破片的飞散范围,给出了爆炸冲击波的安全距离,分析了分配板、库房强度等因素对冲击波传播和破片飞散的影响。研究表明,这种新型危险品仓库可实现定向泄爆,有效限制库房两侧及后方爆炸冲击波的传播和爆炸破片的飞散,使库房两侧及后方的安全距离最大减小77%;与覆土库相比,库房后方的安全距离可减小约50%。钢筋混凝土分配板是新型危险品仓库的重要组成部分,同无分配板库房相比,最大可使后方安全距离减小30%。与波纹钢库房主体相比,强度较高的钢筋混凝土库房主体可使库房后方的安全距离最大减小38%。     

考虑中性层位置变化的微悬臂梁气体传感器静态模型的分析

微悬臂梁是一种高灵敏度的生化传感器。本文考虑吸附表面应力引起的中性层位置的变化,采用能量法建立了微悬臂梁在单层分子吸附稳定后的静态弯曲模型,并以表面吸附有水蒸汽分子的微悬臂梁为例,研究了微悬臂梁曲率半径随其厚度、杨氏模量及吸附分子间距的变化规律以及中性层位置变化对微悬臂梁传感器性能预测的影响,结果发现：1)微悬臂梁的曲率半径与其杨氏模量、厚度及吸附分子间距之间可以近似用一次、二次和八次函数关系表示;2）中性层变化导致的曲率半径计算误差,随着微悬臂梁厚度、杨氏模量的增加而减小,但影响较小,而吸附分子间距会对该相对误差产生明显影响;3）中性层位置变化会对微悬臂梁传感器灵敏度和表面应变预测产生明显的影响。     

装药比对中心炸药抛撒钨颗粒群影响的实验研究

在炸药抛撒颗粒的过程中,颗粒与炸药质量之比对抛撒颗粒的动能特性有重要影响。对中心为8710药柱,周围为钨颗粒的圆柱形装药进行了实验。实验中装药中心药柱的质量为200g,钨颗粒的质量与药柱的质量比(简称装药比)在2～7之间。实验结果表明,在有效作用距离内(2m～6m)颗粒数密度随装药比的增加呈正比增加;在近距离上(2m,3m),随装药比增大颗粒对纸靶的最大穿深明显下降,但在较远的距离上(4m,5m,6m),各装药颗粒对纸靶的最大穿深比较接近,随装药比增大而下降的规律不明显。颗粒对纸靶的平均穿深随装药比变化的规律与最大穿深相似。纸靶单位面积上的总穿深对装药比不敏感,随装药比的增大略成线性增大的趋势。     

Affordance-controlled bifurcations of action patterns in martial arts

Effects of participant-target distance and perceived handstriking efficiency on emergent behavior in the martial art of boxing were investigated, revealing affordance-controlled nonlinear dynamical effects (i.e. bifurcations) within the participant--target system. Results established the existence of critical values of scaled distances for emergence of first time excitations and annihilations of a diverse range of boxing actions i.e. on the appearance and dissolution of jabs, hooks and uppercuts. Reasons for the action diversity were twofold: (a) topological discontinuous changes (bifurcations) in the number of possible handstrikes, i.e. motor solutions to the hitting task; (b) fine modification of probabilities of emergence of striking patterns. Exploitation of a 'strikeability' affordance available in scaled distance-to-target information by boxers led to a diversity of emergent actions through a cascade of bifurcations in the task perceptual-motor work space. Data suggested that perceived efficiency (E) of an action changed as a function of scaled distance (D) and was correlated with the probability of occurrence of action patterns (P), exhibiting the following dependence P = P(E(D)). The implication is that probability of occurrence (P) depends on efficiency (E), which in turn depends on scaled distance (D) to the target. Accordingly, scaled distance-dependent perceived efficiency seems a viable candidate for a contextual (control) parameter to describe the nonlinear dynamics of striking actions in boxing.     

Identification of elastic property and loading fields from full-field displacement measurements

A method is introduced to identify simultaneously elastic properties and loading fields from a measured displacement field. Since the mechanical behavior of micro-electro-mechanical systems (MEMS) is governed by surface effects, this type of identification tool is thought to be of major interest. However, increasing the number of parameters to retrieve affects the redundancy necessary for an accurate identification. A finite-element formulation of a distance between measured and statically admissible (SA) displacement fields is shown to be equivalent to a standard least-squares distance to kinematically admissible (KA) fields if the used modeling is suitable. Any deviation from this equivalence is then the signature of a modeling error. Balancing the distance to KA and SA displacement fields allows one to retrieve unknown modeling parameters. This method is detailed on heterogeneous Euler–Bernoulli beams submitted to an unknown loading field and applied to experimental displacement fields of micro-cantilevers obtained with an electrostatic set-up. An elastic property field and a parameterized loading field are then identified, and the quality of the identification is assessed.     

Fast computation of the wall distance in unsteady Eulerian fluid-structure computations

Highly nonlinear, turbulent, dynamic, fluid-structure interaction problems characterized by large structural displacements and deformations, as well as self-contact and topological changes, are encountered in many applications. For such problems, the Eulerian computational framework, which is often equipped with an embedded (or immersed) boundary method for computational fluid dynamics, is often the most appropriate framework. In many circumstances, it requires the computation of the time-dependent distance from each active mesh vertex of the embedding mesh to the nearest embedded discrete surface. Such circumstances include, for example, modeling turbulence using the Spalart-Allmaras or detached eddy simulation turbulence models and performing adaptive mesh refinement in order to track the boundary layer. Evaluating at each time step the distance to the wall is computationally prohibitive, particularly in the context of explicit-explicit fluid-structure time-integration schemes. Hence, this paper presents two complementary approaches for reducing this computational cost. The first one recognizes that many quantities depending on the wall distance are relatively insensitive to its inaccurate evaluation in the far field. Therefore, it simplifies a state-of-the-art algorithm for computing the wall distance accordingly. The second approach relies on an effective wall distance error estimator to update the evaluation of the wall distance function only when otherwise, a quantity of interest that depends on it would become tainted by an unacceptable level of error. The potential of combining both approaches for dramatically accelerating the computation of the wall distance is demonstrated with the Eulerian simulation of the inflation of a disk-gap-band parachute system in a supersonic airstream.     

动、静裂纹作用偏置效应的动焦散冲击实验

为了研究冲击荷载作用下脆性材料中运动裂纹与静止裂纹的相互作用,选取动态载荷下断裂行为与岩石材料类似且本身光学特性较好的有机玻璃（PMMA）作为实验材料,试件尺寸为220 mm×50 mm×5 mm,采用激光切割制作长度5 mm的预制裂纹和长度10 mm的静止裂纹,预制裂纹位于试件的底部边缘中心,静止裂纹的中心位于试件水平轴线。将静止裂纹偏置距离作为单一变量,采用数字激光动态焦散实验系统对含不同缺陷的PMMA进行三点弯曲实验,并结合几何分形理论研究不同偏置距离下运动裂纹的分形规律。实验结果表明:存在预制裂纹与静止裂纹的临界偏置距离（6 mm）,该条件下裂纹轨迹对应的分形维数值最大,裂纹轨迹的规则程度最低,裂纹破坏形态最复杂。当预制裂纹与静止裂纹的偏置距离在0～6 mm时,裂纹Ⅰ起裂后垂直向上扩展一段距离,与静止裂纹交汇,并停滞一段时间后发生二次起裂,直至贯穿试件,偏置距离和交汇点竖向坐标值呈近似线性函数关系。偏置距离的存在不会影响裂纹Ⅰ的起裂时间和应力强度因子,但会显著减小裂纹Ⅱ的动态应力强度因子,且停滞时长随偏置距离的增大而逐渐缩短。当偏置距离大于临界偏置距离时,运动裂纹不再与静止裂纹交汇而是呈拱状向试件上边缘扩展直至贯穿,裂纹的起偏时间、起偏位置也会出现明显的滞后现象。     

Intrinsic Geometry and Analysis of Diffusion Processes and L ∞-Variational Problems

The aim of this paper is twofold. First, we obtain a better understanding of the intrinsic distance of diffusion processes. Precisely, (a) for all n ≧ 1, the diffusion matrix A is weak upper semicontinuous on Ω if and only if the intrinsic differential and the local intrinsic distance structures coincide; (b) if n = 1, or if n ≧ 2 and A is weak upper semicontinuous on Ω, the intrinsic distance and differential structures always coincide; (c) if n ≧ 2 and A fails to be weak upper semicontinuous on Ω, the (non-)coincidence of the intrinsic distance and differential structures depend on the geometry of the non-weak-upper-semicontinuity set of A. Second, for an arbitrary diffusion matrix A, we show that the intrinsic distance completely determines the absolute minimizer of the corresponding L ^∞-variational problem, and then obtain the existence and uniqueness for given boundary data. We also give an example of a diffusion matrix A for which there is an absolute minimizer that is not of class C ¹. When A is continuous, we also obtain the linear approximation property of the absolute minimizer.