气动力/热与结构多场耦合计算策略与方法研究

气动力/热与结构热响应多场耦合问题涉及物理因素众多且不可简单分割。本文对该耦合问题涉及的各物理场特点与耦合关系进行了研究,从基本物理特点和方程出发,研究了物理场间的耦合关系,并对一直沿用的多物理场耦合关系提出了新的观点与看法。以此为基础,研究获得了气动力/热与结构热响应多场耦合问题的数据流程,提出了针对该耦合问题特有的时间-空间耦合概念,并对其耦合策略与方法进行了深入探讨,为进一步开展相关计算研究工作奠定了理论基础。     

Chodorow型耦合腔慢波结构色散特性和耦合阻抗理论分析

本文建立了Chodorow型耦合腔慢波结构的解析模型, 利用并矢格林函数结合矩量法求解了场匹配方程, 给出了色散方程和耦合阻抗的计算式, 并数值计算出一个X波段Chodorow型慢波结构的高频特性. 结果表明, 本文方法的色散特性以及耦合阻抗与仿真软件HFSS计算的结果有很好的一致性, 且计算效率更高, 同时精度远高于等效电路法, 对工程设计有好的参考价值. 关键词： Chodorow型耦合腔慢波结构 色散特性 耦合阻抗 场匹配     

基于热超构材料的能量收集与热电转换特性

针对温差发电器的局限性,利用热超构材料的热场调控特性,提出了将温差发电器与二维扇形热超构材料能量收集结构进行集成,从而改善温差发电器的热电转换效率.基于有限元多物理场仿真软件COMSOL Multiphysics研究了不同材料对能量收集结构热场调控性能的影响,确定材料后对其进行热电性能仿真,仿真结果表明,能量收集结构可实现热流的有效调控,在同一仿真条件下能量收集中心的温度梯度相比自然材料提高了8倍.对不同尺寸温差发电器发电量进行研究,在此基础上综合考虑加工精度和测试难度,完成了能量收集结构3维建模及加工制造.搭建实验测试系统,使用热成像仪观测能量收集结构的温度分布,测试实验结果显示该能量收集结构可以有效调控热场,在相同冷热源条件下相比自然材料结构可以将温差发电器的工作效率提高3.2倍,对推动温差发电技术更加迅速地发展具有一定的现实意义.     

风力机翼型气动计算的降阶模型研究

将基于POD的降阶模型应用于风力机翼型的气动研究。首先应用CFD数值模拟得到一系列快照结果;应用基于本征正交分解(POD)的方法得到流场的一组基模态,认为对于所研究的问题,任一流场可以由这些基模态通过线性表达得到;对控制方程进行Galerkin投影,得到降阶模型,将离散求解N-S方程的问题转化为一组只有十几个自由度的常微分方程,从而减少计算时间,提高计算效率。并对二维翼型的绕流的定常和非定常问题进行了分析,计算结果表明,降阶模型可以较好地捕捉流动的特征,与直接CFD模拟相比计算精度相当,但大幅有效地提高了计算效率。     

类前缘防热层流场与热响应耦合计算研究

本文在以往对类前缘防热层热响应计算分析基础上,进一步研究实现了外流场高超音速NS方程数值计算表面气动加热与防热层结构热响应的耦合计算,这对于常用的非耦合计算方法来说是一进步,也为进一步开展外流场/结构热响应/热应力全耦合一体化计算研究和防热层表面吹气强化传热问题的场协同研究打下了基础。     

旋转柔性悬臂梁动力学的Bezier插值离散方法研究

在旋转柔性梁变形场描述中,引入Bezier插值离散方法.首先构建旋转运动悬臂梁物理模型,接着采用第二类Lagrange动力学方程和Bezier插值离散方法,在计入柔性梁横向弯曲变形引起的纵向缩短的情况下,推导了旋转柔性梁的刚柔耦合动力学方程,并编制旋转柔性梁的动力学仿真软件,然后通过仿真算例对系统的动力学问题进行研究.最后将仿真结果与有限元法、假设模态法进行分析比较,验证了提出的Bezier插值离散方法的正确性,并得出Bezier插值离散法的计算效率较高;计算精度符合工程实际需要,高速时计算精度大于假设模态法;Bezier插值离散方法在处理大柔性问题时比假设模态法合理.因此在多体系统动力学领域具有优良性能和应用价值的Bezier插值离散方法将具有推广价值.     

CFD在高速飞行器热气弹问题中的应用

采用气动力/热/结构耦合的方法对高速细长体飞行器结构热静气动弹性问题进行了研究.为保证耦合计算精度，达到准确预测热气动弹性特性的能力，气动力和气动热计算采用CFD数值模拟方法，热应力和热变形计算采用有限元方法并通过热考核试验验证.以该简单细长体飞行器模型为研究对象，对其热静气动弹性特性进行了计算与分析，计算结果表明:CFD/CSD耦合可准确模拟热气弹问题，且气动加热造成结构温升不均衡是结构变形的主导因素，力热耦合静气弹变形与单纯受力分析变形形式不同，对飞行器气动特性影响规律不同.准确预测飞行器热气动弹性特性对飞行器结构设计十分必要.      

涡轮动叶冷却结构设计方法Ⅲ:气热耦合计算

三维气热耦合数值模拟是涡轮叶片冷却结构详细设计中应用的热分析方法。三维气热耦合数值模拟的计算域实体模型建立借助于参数化方法,模型包括冷却结构尽可能多的细节、并保证足够高的模型精度。网格生成时,流体域采用结构化网格,并且需要在换热壁面上根据湍流模型的要求进行加密。将管网计算结果与三维数值模拟结果进行对比,发现管网计算的精度有限,而气热耦合计算方法能够捕捉更多流动换热现象,计算结果相对可信。对单个冷却结构方案进行三维数值计算的周期为5～15天,基本满足工程设计要求。     

基于电缆-无线耦合的随钻测井信号传输新方法及其数值模拟研究

信号传输是随钻测井的重要组成部分, 也是随钻测井进一步发展的瓶颈.针对这个问题, 研究了由Novatek^TM提出的基于钻杆内有线传输, 钻杆间电磁耦合传输的随钻信号高速传输方法.首先, 针对钻杆间的电磁耦合结构的对称性, 采用了高效的数值模式匹配法对该耦合结构进行了电磁建模; 然后, 为了优化钻杆间耦合结构的耦合性能, 利用上述的电磁建模方法对该耦合结构进行了仿真分析, 分析了该耦合结构各个参数对耦合性能的影响, 并得出了一些有意义的结论, 该结论能够指导随钻测井信号钻杆间耦合传输的优化设计.另外, 钻杆内传输线利用矩形传输线设计, 设计出了50 Ω 特性阻抗的矩形传输线, 并计算了相关损耗.最后, 通过仿真和实验验证了优化设计的钻杆间耦合传输结构具有非常好的传输性能, 并且仿真和实验结果符合良好. 关键词： 随钻信号传输 数值模式匹配 电磁耦合 矩形传输线     

双驱动电动换能器电磁场和声-结构耦合场分析

电动换能器的设计涉及电磁场和声-结构耦合问题,本文利用COMSOL多物理场分析软件提供的该两种分析模块进行组合,以电磁理论为基础推导出了适合COMSOL数值计算的公式,把两种耦合场的分析结合起来,利用电磁场分析过程计算所需要变量,代入声-结构耦合场进行分析,得到声场分布情况。文章也大致介绍了COMSOL的有限元建模和该两种模块的使用方法。本文对电动换能器结构如何优化不做讨论。     

A high-efficiency aerothermoelastic analysis method

In this paper,a high-efficiency aerothermoelastic analysis method based on unified hypersonic lifting surface theory is established.The method adopts a two-way coupling form that couples the structure,aerodynamic force,and aerodynamic thermo and heat conduction.The aerodynamic force is first calculated based on unified hypersonic lifting surface theory,and then the Eckert reference temperature method is used to solve the temperature field,where the transient heat conduction is solved using Fourier’s law,and the modal method is used for the aeroelastic correction.Finally,flutter is analyzed based on the p-k method.The aerothermoelastic behavior of a typical hypersonic low-aspect ratio wing is then analyzed,and the results indicate the following:(1)the combined effects of the aerodynamic load and thermal load both deform the wing,which would increase if the flexibility,size,and flight time of the hypersonic aircraft increase;(2)the effect of heat accumulation should be noted,and therefore,the trajectory parameters should be considered in the design of hypersonic flight vehicles to avoid hazardous conditions,such as flutter.     

A novel simulation theory and model system for multi-field coupling pipe-flow system

Due to the lack of a theoretical basis for multi-field coupling in many system-level models, a novel set of system-level basic equations for flow/heat transfer/combustion coupling is put forward. Then a finite volume model of quasi-1D transient flow field for multi-species compressible variable-cross-section pipe flow is established by discretising the basic equations on spatially staggered grids. Combining with the 2D axisymmetric model for pipe-wall temperature field and specific chemical reaction mechanisms, a finite volume model system is established; a set of specific calculation methods suitable for multi-field coupling system-level research is structured for various parameters in this model; specific modularisation simulation models can be further derived in accordance with specific structures of various typical components in a liquid propulsion system. This novel system can also be used to derive two sub-systems: a flow/heat transfer two-field coupling pipe-flow model system without chemical reaction and species diffusion; and a chemical equilibrium thermodynamic calculation-based multi-field coupling system. The applicability and accuracy of two sub-systems have been verified through a series of dynamic modelling and simulations in earlier studies. The validity of this system is verified in an air–hydrogen combustion sample system. The basic equations and the model system provide a unified universal theory and numerical system for modelling and simulation and even virtual testing of various pipeline systems.     

Nonlinear Energy Cascading in Turbulence during the Internal Reconnection Event at the Sino-United Spherical Tokamak

The characteristics of the energy transfer and nonlinear coupling among edge electromagnetic turbulence in thermal quench sub-period of the internal reconnection event(IRE) are studied at the sino-united spherical tokamak device using multiple Langmuir and magnetic probe arrays. The wavelet bispectral analysis and the modified Kim method are applied to investigate linear growth/damping and nonlinear energy transfer rates, along with multi-field turbulence interactions. The results show a multi-field nonlinear energy transfer from electrostatic to magnetic turbulence that results in two-mode coupling in magnetic turbulence, which may play a crucial role to trigger the IRE.     

A convergence study of a new partitioned fluid–structure interaction algorithm based on fictitious mass and damping

We develop, analyze and validate a new method for simulating fluid–structure interactions (FSIs), which is based on fictitious mass and fictitious damping in the structure equation. We employ a partitioned method for the fluid and structure motions in conjunction with sub-iteration and Aitken relaxation. In particular, the use of such fictitious parameters requires sub-iterations in order to reduce the induced error in addition to the local temporal truncation error. To this end, proper levels of tolerance for terminating the sub-iteration procedure have been obtained in order to recover the formal order of temporal accuracy. For the coupled FSI problem, these fictitious terms have a significant effect, leading to better convergence rate and hence substantially smaller number of sub-iterations. Through analysis we identify the proper range of these parameters, which we then verify by corresponding numerical tests. We implement the method in the context of spectral element discretization, which is more sensitive than low-order methods to numerical instabilities arising in the explicit FSI coupling. However, the method we present here is simple and general and hence applicable to FSI based on any other discretization. We demonstrate the effectiveness of the method in applications involving 2D vortex-induced vibrations (VIV) and in 3D flexible arteries with structural density close to blood density. We also present 3D results for a patient-specific aneurysmal flow under pulsatile flow conditions examining, in particular, the sensitivity of the results on different values of the fictitious parameters.     

Insights into the regulation mechanism of ring-shaped magnetoelectric energy harvesters via mechanical and magnetic conditions

This paper presents a theoretical model for predicting and tuning magnetoelectric(ME) effect of ring-shaped composites, in which stress boundary conditions are empoyed and the multi-field coupling property of giant magnetostrictive materials are taken into account. A linear analytical solutions for the closed-and open-circuit ME voltages are derived simultaneously using mechanical differential equations, interface and boundary conditions, and electrical equations. For nonlinear ME coupling effect, the nonlinear multi-field coupling constitutive equation is reduced to an equivalent form by expanding the strains as a Taylor series in the vicinity of bias magnetic field. Sequentially, the linear model is generalized to a nonlinear one involving the field-dependent material parameters. The results show that setting a stress-free condition is beneficial for reducing resonance frequency while applying clamped conditions on the inner and outer boundaries may improve the maximum output power density. In addition, performing stress conditions on one of the boundaries may enhance ME coupling significantly, without changing the corresponding resonance frequency and optimal resistance. When external stimuli like bias magnetic field and pre-stress are applied to the ring-shaped composites, a novel dual peak phenomenon in the ME voltage curve around resonance frequencies is revealed theoretically, indicating that strong ME coupling may be achieved within a wider bias field region. Eventually, the mutual coordination of the bias field and pre-stress may enhance ME coupling as well as tuning the resonance frequency, and thus is pivotal for tunable control of ME energy harvesters. The proposed model can be applied to design high-performance energy harvesters by manipulating the mechanical conditions and external stimuli.     

Elastic and anelastic relaxations in the relaxor ferroelectric Pb(Mg1/3Nb2/3)O3: I. Strain analysis and a static order parameter

The structural evolution of Pb(Mg(1/3)Nb(2/3))O(3) (PMN) has been reviewed in terms of characteristic temperatures, length scales and timescales, with a view to considering the overall relaxor behaviour from the perspectives of strain and elasticity. A conventional analysis of lattice parameter data in terms of spontaneous strain and strain/order parameter coupling shows that even though a normal phase transition does not occur the relaxor ordering process is accompanied by a significant volume strain which follows the pattern of a static order parameter evolving according to that expected for a tricritical phase transition with T(c) ≈ 350 K. This matches the evolution of the intensity of the elastic central peak in neutron scattering spectra, and reflects the development of static (or quasistatic) polar nanoregions (PNRs) as if by a mean-field phase transition. Use of a Landau free energy expansion, which includes Γ4(-) order parameter components to describe ferroelectric contributions and an R1(+) order parameter to describe cation ordering together with their formal coupling with strain, then allows the pattern of elastic softening expected for a cubic → rhombohedral phase transition to be anticipated. The extent to which observed softening differs from this static mean-field pattern serves to highlight the additional roles of local heterogeneity and relaxation dynamics in determining the relaxor properties of PMN.     

The Characterization of Absorptive Nonlinearities

The measurement of the order and magnitude of absorptive nonlinearities is essential to the development and improvement of a wide range of optical materials. This review presents and contrasts some of the most commonly used methods for making such measurements. As background, we introduce the history of, and the essential physics behind, multiphoton absorption (MPA) and phenomena that can interfere with its measurement. We review methods for the direct and indirect measurement of MPA. The direct methods covered include nonlinear transmission, absorption modulation, Z‐scan, and spectrally‐resolved, two‐beam coupling. The indirect methods reviewed, all of which rely on the measurement of multiphoton action cross sections, include nonlinear fluorescence excitation, thermal lensing, MPA‐induced photocurrent, and various techniques for determining the order of effective nonlinear absorption in multiphoton absorption polymerization.     

A theoretical analysis of mode-mixing in optical waveguides with nearest-neighbour mode coupling

Properties of the solutions of the coupled equations describing the flow of average mode power in optical waveguides are investigated. All of the important properties and their theoretical foundations are first reviewed for systems with arbitrary coupling and loss coefficients. It is then shown that, for the special class of systems with nearest-neighbour coupling, the solutions are expressible in terms of orthogonal polynomials. In particular, for systems obeying a quasi-uniform loss model and having coupling coefficients with simple dependence on mode number (uniform or linear), the solutions are associated with the classical polynomials. As illustrations of this analysis, the characteristics of optical power flow in three such systems are studied in detail. These include a slab-waveguide with uniform coupling, the corresponding cylindrical waveguide with uniform coupling, and a parabolic-index fibre with linear coupling dependence.     

基于原子光谱法的元素形态分析技术进展

元素的形态分析在环境和生物分析中极其重要,因为元素在生物体内的作用及其代谢过程在很大程度上取决于元素存在的化学形态,而不仅仅是元素的总量。形态分析是指测定样品中构成元素总量的单独物理化学形式的浓度。最先进的形态分析方法是色谱分离和光谱检测的联用技术,特别是色谱和电感耦合等离子-质谱(ICP-MS)的联用。但是,联用技术的设备投入大、运行成本高,难以在常规实验室中的推广应用。在许多情况下,采用非色谱分离方法对样品进行处理,也可以得到足够的元素形态信息。基于非色谱分离、原子光谱测定的元素形态分析方法的费用低、操作简单、易于推广应用。本文总结了元素形态分析的样品前处理方法,综述了基于原子光谱法的元素形态分析中比较常用的非色谱分离技术,对溶剂萃取、浊点萃取、单滴微萃取、分散液液萃取等分离技术的原理、应用和优缺点进行了评述,介绍了固相萃取中常用的吸附剂及其在元素形态分离中的应用,以及氢化物发生、共沉淀等分离方法。相比于色谱分离方法,非色谱方法是快速、灵敏、廉价的分离技术。     

Temporal summation of pulsate brain stimulation in normal and deafened cats.

Behaviorally measured, electrical-stimulation thresholds were obtained from 11 electrodes permanently positioned in the auditory system and other brain loci. Number of pulses and interpulse intervals were varied to determine how detection thresholds were affected by stimulation parameters. Detection thresholds generally decreased with increased number of pulses and with shorter interpulse intervals. A method is presented to describe the parametric threshold data for each electrode in terms of three constants: a single-pulse threshold which characterizes the sensitivity of the placement; a time constant of temporal summation; and a compression factor which describes the range of threshold variation. For three placements in the vicinity of cochlear nucleus, bilateral cochlear destruction permanently altered parametric thresholds. In particular, single-pulse threshold was lowered by 9.2 dB; time constant of temporal summation was reduced by a factor of 100; and the compression factor was increased. Classic strength-duration time constants were determined using behavioral methods and were shown to be equal in magnitude to the greatly reduced time constants for temporal summation in the deafened animals. This implies that capacity for temporal integration may be substantially reduced or lost in at least the lower level of the auditory system following deafness.