Comparing turbulence models for flow through a rigid glottal model

Flow through a rigid model of the human vocal tract featuring a divergent glottis was numerically modeled using the Reynolds-averaged Navier-Stokes approach. A number of different turbulence models, available in a widely used commercial computational fluid dynamics code, were tested to determine their ability to capture various flow features recently observed in laboratory experiments and large eddy simulation studies. The study reveals that results from unsteady simulations employing the k-omega shear stress transport model were in much better agreement with previous measurements and predictions with regard to the ability to predict glottal jet skewing due to the Coanda effect and the intraglottal pressure distribution or related skin friction coefficient, than either steady or unsteady simulations using the Spalart-Allmaras model or any other two-equation turbulence model investigated in this study.     

基于声带振动模型和声门波的嘶音研究

本文根据嘶音的主要病理表现为声带的病变、声门波反映了声带的运动状态,提出基于非对称四质量块声带振动模型及声门波分析一合成的嘶音研究方法。将嘶音信号波形与其声门波、声带振动模型联系起来,通过对正常语音和嘶哑病人术前术后语音信号的声门波和声带振动模型特征参数的研究,给出了常态语音和嘶音的声门波周期性、声带两侧参数对称性等参数的对比结果,分析了模型参数与嘶音生理与病理因素之间的关系。实验表明,基于声门波和声带振动模型的嘶音研究可以揭示嘶音的声学特征参数与病理因素的关系,为实现喉科疾病无接触诊断以及嘶音音质的改善提供理论和实验依据。     

A low-dimensional model of separation bubbles

In this work, motivated by the problem of model-based predictive control of separated flows, we identify the key variables and the requirements on a model-based observer and construct a prototype low-dimensional model to be embedded in control applications.Namely, using a phenomenological physics-based approach and dynamical systems and singularity theories, we uncover the low-dimensional nature of the complex dynamics of actuated separated flows and capture the crucial bifurcation and hysteresis inherent in separation phenomena. This new look at the problem naturally leads to several important implications, such as, firstly, uncovering the physical mechanisms for hysteresis, secondly, predicting a finite amplitude instability of the bubble, and, thirdly, to new issues to be studied theoretically and tested experimentally.     

Voice production model integrating boundary-layer analysis of glottal flow and source-filter coupling

A voice production model is created in this work by considering essential aerodynamic and acoustic phenomena in human voice production. A precise flow analysis is performed based on a boundary-layer approximation and the viscous-inviscid interaction between the boundary layer and the core flow. This flow analysis can supply information on the separation point of the glottal flow and the thickness of the boundary layer, both of which strongly depend on the glottal configuration and yield an effective prediction of the flow behavior. When the flow analysis is combined with the modified two-mass model of the vocal fold [Pelorson et al. (1994). J. Acoust. Soc. Am. 96, 3416-3431], the resulting acoustic wave travels through the vocal tract and a pressure change develops in the vicinity of the glottis. This change can affect the glottal flow and the motion of the vocal folds, causing source-filter coupling. The property of the acoustic feedback is explicitly expressed in the frequency domain by using an acoustic tube model, allowing a clear interpretation of the coupling. Numerical experiments show that the vocal-tract input impedance and frequency responses representing the source-filter coupling have dominant peaks corresponding to the fourth and fifth formants. Results of time-domain simulations also suggest the importance of these high-frequency peaks in voice production.     

Parameterization of the glottal area, glottal flow, and vocal fold contact area

A new set of parameters is described for analysis and synthesis of glottal area, vocal fold contact area, and glottal volume flow. Parameters are all nondimensionalized and consist of an abduction quotient, a shape quotient, a phase quotient, and a load quotient in addition to fundamental frequency and vibrational amplitude. The parameters show promise in interpretation of electroglottographic, photoglottographic, and inverse filtered volume velocity waveforms in terms of the glottal configuration. Some comparisons between modeled and measured glottographic waveforms are made.     

SIM--simultaneous inverse filtering and matching of a glottal flow model for acoustic speech signals

A new method "simultaneous inverse filtering and model matching" (SIM) is proposed that allows one to calculate voice source measures without any user interaction. It is based on the discrete all-pole modeling (DAP) technique for inverse filtering (IF), which is modified to include a model of the glottal flow as integral part [LF model, Fant et al., STL-QPSR (Stockholm) 4/1985, 1-13 (1986)]. As the correct LF parameters are initially unknown, they are estimated in an iterative procedure using multi-dimensional optimization techniques that are initialized according to the results of an exhaustive search. The error criteria applied reflect how well the IF is performed after the spectral contribution of the glottal flow has been removed. The resulting optimal LF parameter constellation serves as the basis to calculate 11 voice source measures. The performance was evaluated using synthesized signals and recordings of natural utterances. For the synthesized signals, the accuracy to reproduce the original parameters was high (correlations exceeding 0.88) for measures where the starting point of the glottal cycle did not enter explicitly. Errors were smaller compared to conventional estimation methods where the measures were estimated from the IF signal. The analysis of natural utterances indicates that problems still exist with regard to robustness, but that under advantageous conditions the open quotient, the speed quotient, the closing quotient, the parabolic spectral parameter, and the negative peak amplitude of the glottal flow derivative can indeed be determined automatically by the SIM method.     

A new single-image fog removal algorithm based on physical model

Based on the physical model of atmospheric scattering and the optical reflectance imaging model, three major factors which affect the effect of fog removal are discussed in detail, dark channel phenomenon is explained via the optical model, and an approach for solving the parameter in the atmospheric scattering model is rigorously derived from a new perspective. Using gray-scale opening operation and fast joint bilateral filtering techniques, the proposed algorithm can effectively obtain the global atmospheric light and greatly improve the speed and accuracy of atmospheric scattering function solving. Finally, the scene albedo is recovered by inverting this model. Compared with existing algorithms, complexity of the proposed method is only a linear function of the number of input image pixels and this allows a very fast implementation. The simulation results show that the processing time of images with a resolution of 576*768 is only 1.7 s; Results on a variety of outdoor foggy images demonstrate that the proposed method achieves good restoration for contrast and color fidelity, resulting in a great improvement in image visibility.     

Numerical simulation of turbulence transition and sound radiation for flow through a rigid glottal model

Large eddy simulation (LES)-based computational aeroacoustics techniques were applied to a static model of the human glottis, idealized here as a planar channel with an orifice, to study flow-acoustic interactions related to speech. Rigid models of both converging and diverging glottal passages, each featuring a 20 deg included angle and a minimal glottal diameter of 0.04 cm, with an imposed transglottal pressure of 15 cm H2O, were studied. The Favre-filtered compressible Navier-Stokes equations were integrated for this low-Mach-number flow using an additive semi-implicit Runge-Kutta method and a high-order compact finite-difference scheme with characteristic-based nonreflecting boundary conditions and a multiblock approach. Flow asymmetries related to the Coanda effect and transition to turbulence, as well as the far-field sound, were captured. Acoustic-analogy-based far-field sound predictions were compared with direct simulations and showed that dipole sources, arising from unsteady flow forces exerted on the glottal walls, are primarily responsible for the tonal sound observed in the divergent glottis case.     

A low-dimensional model of binocular rivalry using winnerless competition

We discuss a novel minimal model for binocular rivalry (and more generally perceptual dominance) effects. The model has only three state variables, but nonetheless exhibits a wide range of input and noise-dependent switching. The model has two reciprocally inhibiting input variables that represent perceptual processes active during the recognition of one of the two possible states and a third variable that represents the perceived output. Sensory inputs only affect the input variables.We observe, for rivalry-inducing inputs, the appearance of winnerless competition in the perceptual system. This gives rise to a behaviour that conforms to well-known principles describing binocular rivalry (the Levelt propositions, in particular proposition IV: monotonic response of residence time as a function of image contrast) down to very low levels of stimulus intensity.     

A low-dimensional crystal growth model on an isotropic and quasi-free sustained substrate

A new crystal growth theoretical model is established for the low-dimensional nanocrystals on an isotropic and quasifree sustained substrate. The driven mechanism of the model is based on the competitive growth among the preferential growth directions of the crystals possessing anisotropic crystal structures, such as the hexagonal close-packed and wurtzite structures. The calculation results are in good agreement with the experimental findings in the growth process of the lowdimensional Zn nanocrystals on silicone oil surfaces. Our model shows a growth mechanism of various low-dimensional crystals on/in the isotropic substrates.     

Normalized amplitude quotient for parametrization of the glottal flow

Normalized amplitude quotient (NAQ) is presented as a method to parametrize the glottal closing phase using two amplitude-domain measurements from waveforms estimated by inverse filtering. In this technique, the ratio between the amplitude of the ac flow and the negative peak amplitude of the flow derivative is first computed using the concept of equivalent rectangular pulse, a hypothetical signal located at the instant of the main excitation of the vocal tract. This ratio is then normalized with respect to the length of the fundamental period. Comparison between NAQ and its counterpart among the conventional time-domain parameters, the closing quotient, shows that the proposed parameter is more robust against distortion such as measurement noise that make the extraction of conventional time-based parameters of the glottal flow problematic. Experiments with breathy, normal, and pressed vowels indicate that NAQ is also able to separate the type of phonation effectively.     

Dynamic glottal pressures in an excised hemilarynx model

During phonation, air pressures act upon the vocal folds to help maintain their oscillation. The air pressures vary dynamically along the medial surface of the vocal folds, although no live human or excised studies have shown how those pressure profiles vary in time. The purpose of this study was to examine time-dependent glottal pressure profiles using a canine hemilarynx approach. The larynx tissue was cut in the midsaggital plane from the top to about 5 mm below the vocal folds. The right half was replaced with a Plexiglas pane with imbedded pressure taps. Simultaneous recordings were made of glottal pressure signals, subglottal pressure, particle velocity, and average airflow at various levels of adduction. The data indicate that the pressures in the glottis (on the Plexiglas) vary both vertically and longitudinally throughout the phonatory cycle. Pressures vary most widely near the location of maximum vibratory amplitude, and can include negative pressures during a portion of the cycle. Pressures anterior and posterior to the maximum amplitude location may have less variation and may remain positive throughout the cycle, giving rise to a new concept called dynamic bidirectional pressure gradients in the glottis. This is an important concept that may relate strongly to tissue health as well as basic oscillatory mechanics.     

A generalized critical-state model for low-dimensional superconductors with an edge barrier

A generalized model of the critical state in low-dimensional superconductors is proposed, which includes the effects of bulk and edge pinning of the magnetic flux. The conditions are determined for the penetration of Pearl-Abrikosov vortices into the superconductor, and metastable flux structures in these systems are described. A diagram of different vortex states is constructed for a complete cycle of variations of the external magnetic field.     

Measurement of the glottal impedance with a mechanical model

An approximate but realistic model of the human larynx was constructed to gain better knowledge of the complex glottal impedance and its dependence on glottal width, flow, and frequency. The glottal width was adjustable from 0 to 3 mm, the flow from 0 to 500 cm3/s. The model was fitted into a system of tubes, through which compressed air could be conducted. Supraglottally, a broadband signal was fed into the tube, and, with a two-microphone directional coupler, the complex glottal impedance at a given reference plane was directly determined as a function of frequency. Since the calculated impedance is sensitively dependent on the definition of the position of the reference plane, it is difficult to obtain quantitative statements about the frequency dependence. Nevertheless, in the presence of flow, it is possible to achieve reliable results by analysis of the relative position of the measured curves. On the one hand, the glottal inductance decreases linearly with increasing flow velocity; on the other hand, it diminishes nonlinearly with decreasing frequency. Finally, some difficulties in the definition of glottal impedance are pointed out.     

Critical exponents of the Ising model on low-dimensional fractal media

The critical behavior of the Ising model on fractal substrates with noninteger Hausdorff dimension d_H<2 and infinite ramification order is studied by means of the short-time critical dynamic scaling approach. Our determinations of the critical temperatures and critical exponents β, γ, and ν are compared to the predictions of the Wilson-Fisher expansion, the Wallace-Zia expansion, the transfer matrix method, and more recent Monte Carlo simulations using finite-size scaling analysis. We also determined the effective dimension (d_ef), which plays the role of the Euclidean dimension in the formulation of the dynamic scaling and in the hyperscaling relationship d_ef=2β/ν+γ/ν. Furthermore, we obtained the dynamic exponent z of the nonequilibrium correlation length and the exponent θ that governs the initial increase of the magnetization. Our results are consistent with the convergence of the lower-critical dimension towards d=1 for fractal substrates and suggest that the Hausdorff dimension may be different from the effective dimension.     

A non-rigid cardiac image registration method based on an optical flow model

According to non-rigid medical image registration, new method of classification registration is proposed. First, Feature points are extracted based on SIFT (Scale Invariant Feature Transform) from reference images and floating images to match feature points. And the coarse registration is performed using the least square method. Then the precise registration is achieved using the optical flow model algorithm. SIFT algorithm is based on local image features that are with good scale, rotation and illumination invariance. Optical flow algorithm does not extract features and use the image gray information directly, and its registration speed is faster. The both algorithms are complementary. SIFT algorithm is used for improving the convergence speed of optical flow algorithm, and optical flow algorithm makes the registration result more accurate. The experimental results prove that the algorithm can improve the accuracy of the non-rigid medical image registration and enhance the convergence speed. Therefore, the algorithm has some advantages in the image registration.     

Streaminglike diffusion in the low-dimensional stochastic pump model

In this paper we analyze a diffusion phenomenon in a few-dimensional Hamiltonian system of coupled mappings in which the principal component of diffusion occurs along resonances. The result is that the diffusion can have power-law dependence in coupling parameter mu and be independent of the stochastic parameter K. For the same range of parameters, the usual analytical Arnold diffusion across resonances is dependent on K and can be much smaller than resonance streaming diffusion. The results are used to qualitatively explain recent results in multidimensional coupled standard maps.     

Experimental investigation of the influence of a posterior gap on glottal flow and sound

The influence of a posterior gap on the airflow through the human glottis was investigated using a driven synthetic model. Instantaneous orifice discharge coefficient of a glottal shaped orifice was obtained from the time-varying orifice area and the velocity distribution of the pulsated jet measured on the axial plane using a single hot-wire probe. Instantaneous orifice discharge coefficient values were found to undergo a cyclic hysteresis loop when plotted versus Reynolds number and time, indicating a pressure head increase and a net energy transfer from the air flow to the orifice wall. The net energy transferred was estimated to be around 10% of the value presumably required to achieve self-sustained oscillation. The radiated sound pressure was measured to characterize the influence of the minimum flow through the posterior gap on the broadband component of the radiated sound. The presence of a posterior gap was found to significantly increase the broadband sound level produced over the frequency range in which human hearing is most sensitive.     

一个用简单物理模型构建的加密系统

一些简单的物理模型可产生复杂的行为,并有很强的逻辑表达能力.基于旋转元器件和格子气两个物理模型设计了一个加密系统,具有可逆、并行、简单、高效的特性.实验结果证实该系统具有良好的随机性和敏感性.研究结果表明,采用简单物理模型构建加密系统是一种很有潜力的方法. 关键词： 物理模型 加密系统 并行 效率     

低维原子/分子晶体材料的可控生长、物性调控和原理性应用

本文介绍了高鸿钧课题组在物理所20年来的部分代表性工作.研究的主要方向为低维纳米功能材料的分子束外延可控制备、生长机制、物性调控及其在未来信息技术中的原理性应用.从材料的可控制备入手,结合第一性原理的理论计算,阐明材料生长机制和结构与物性的关系,进而实现物性调控和原理性应用.主要内容有:1)纳米尺度"海马"分形结构的形成及其生长机制;2)STM分辨率的提高及最高分辨Si(111)-7×7原子图像的获得;3)固体表面上功能分子的吸附、组装及其机制;4)稳定、重复、可逆的纳米尺度电导转变与超高密度信息存储;5)固体表面上单分子自旋态的量子调控及其原理性应用;6)原子尺度上朗德g因子的空间分辨及其空间分布不均匀性的发现;7)晶圆尺寸、高质量、单晶石墨烯的制备及原位硅插层绝缘化;8)几种新型二维原子晶体材料的可控构筑及其物性调控;9)"自然图案化"的新型二维原子晶体材料及其功能化.这些工作为低维量子结构的构造、物性调控及其原理性应用奠定了基础.