轴向斜缝机匣处理与转子通道之间耦合流动分析

轴向斜缝机匣处理能够有效地扩大压气机的稳定工作范围,本文对轴向斜缝处理机匣与亚音速压气机转子通道之间的耦合流动进行了时间精确的全三维数值模拟,数值结果所获得的总性能与试验结果符合较好。详细分析了处理机匣与转子通道之间非定常干涉现象以及处理机匣结构引入后转子顶部区域不同的流场结构,数值计算结果表明,对具有高损失高阻塞效应的间隙泄漏流及其间隙泄漏涡的控制是轴向斜缝机匣处理扩大压气机稳定工作裕度的主要原因。     

Dual-Band Bi-Directional Pattern Reconfigurable Fractal Patch Antenna for Millimeter Wave Application

A novel dual-band bi-directional reconfigurable antenna based on Koch patch is presented in this paper. By controlling the switches in the slots etched on the Koch patch, different far field bi-directional patterns at the dual-band around 60 GHz/80 GHz could be achieved. And full azimuthal coverage at this dual-band is possible by electronically controlling the switches in the slots of the patch. This work was supported by the National Natural Science Foundation of China (No. 60371008, No. 90505001) and the CRT Program of UESTC.     

A stochastic model for solidification

A 3-dimensional (2-space, 1-time) model relating the diffusion of heat and mass to the kinetic processes at the solid-liquid interface, using a stochastic approach is presented in this paper. This paper is divided in two parts. In the first part the basic set of equations describing solidification alongwith their analysis and solution are given. The process of solidification has a stochastic character and depends on the net probability of transfer of atoms from liquid to the solid phase. This has been modeled by a Markov process in which knowledge of the parameters at the initial time only is needed to evaluate the time evolution of the system. Solidification process is expressed in terms of four coupled equations, namely, the diffusion equations for heat and mass, the equations for concentration of the solid phase and for rate of growth of the solid-liquid interface. The position of the solid-liquid interface is represented with the help of a delta function and it is defined as the surface at which latent heat is evolved. A numerical method is used to solve the equations appearing in the model. In the second part the results i.e. the time evolution of the solid-liquid interface shape and its concentration, rate of growth and temperature are given.     

Space–time discontinuous Galerkin finite element method for two-fluid flows

A novel numerical method for two-fluid flow computations is presented, which combines the space–time discontinuous Galerkin finite element discretization with the level set method and cut-cell based interface tracking. The space–time discontinuous Galerkin (STDG) finite element method offers high accuracy, an inherent ability to handle discontinuities and a very local stencil, making it relatively easy to combine with local hp-refinement. The front tracking is incorporated via cut-cell mesh refinement to ensure a sharp interface between the fluids. To compute the interface dynamics the level set method (LSM) is used because of its ability to deal with merging and breakup. Also, the LSM is easy to extend to higher dimensions. Small cells arising from the cut-cell refinement are merged to improve the stability and performance. The interface conditions are incorporated in the numerical flux at the interface and the STDG discretization ensures that the scheme is conservative as long as the numerical fluxes are conservative. The numerical method is applied to one and two dimensional two-fluid test problems using the Euler equations.     

Time-Domain Numerical Solutions of Maxwell Interface Problems with Discontinuous Electromagnetic Waves

This paper is devoted to time domain numerical solutions of two-dimensional (2D) material interface problems governed by the transverse magnetic (TM) and transverse electric (TE) Maxwell's equations with discontinuous electromagnetic solutions. Due to the discontinuity in wave solutions across the interface, the usual numerical methods will converge slowly or even fail to converge. This calls for the development of advanced interface treatments for popular Maxwell solvers. We will investigate such interface treatments by considering two typical Maxwell solvers – one based on collocation formulation and the other based on Galerkin formulation. To restore the accuracy reduction of the collocation finite-difference time-domain (FDTD) algorithm near an interface, the physical jump conditions relating discontinuous wave solutions on both sides of the interface must be rigorously enforced. For this purpose, a novel matched interface and boundary (MIB) scheme is proposed in this work, in which new jump conditions are derived so that the discontinuous and staggered features of electric and magnetic field components can be accommodated. The resulting MIB time-domain (MIBTD) scheme satisfies the jump conditions locally and suppresses the staircase approximation errors completely over the Yee lattices. In the discontinuous Galerkin time-domain (DGTD) algorithm – a popular Galerkin Maxwell solver, a proper numerical flux can be designed to accurately capture the jumps in the electromagnetic waves across the interface and automatically preserves the discontinuity in the explicit time integration. The DGTD solution to Maxwell interface problems is explored in this work, by considering a nodal based high order discontinuous Galerkin method. In benchmark TM and TE tests with analytical solutions, both MIBTD and DGTD schemes achieve the second order of accuracy in solving circular interfaces. In comparison, the numerical convergence of the MIBTD method is slightly more uniform, while the DGTD method is more flexible and robust.     

Analysis of the radiative transfer equation with highly assymetric phase function

This paper considers a scalar radiative transfer problem with high scattering anisotropy. Two computational methods are presented based on decomposition of the diffuse light field into a regular and anisotropic part. The first algorithm (DOMAS) singles out the anisotropic radiance in the forward scattering peak using the Small-Angle Modification of RTE. The second algorithm (DOM2+) separates the single scattering radiance as an anisotropic part, which largely defines the fine detail of the total radiance in the backscattering directions. In both cases, the anisotropic part is represented analytically. With anisotropy subtraction, the regular part of the signal, which requires a numerical solution, is essentially smoothed as a function of angles. Further, the transport equation is obtained for the regular part that contains an additional source function from the anisotropic part of the signal. This equation is solved with the discrete ordinates method. A conducted numerical analysis of this work showed that algorithm DOMAS has a strong advantage as compared to the standard discrete ordinates method for simulation of the radiance transmission, and DOM2+ is the best of the three for the reflection computations. Both algorithms offer at least a factor of three acceleration of convergence of the azimuthal series for highly anisotropic phase functions.     

Boundary integral solutions of coupled Stokes and Darcy flows

An accurate computational method based on the boundary integral formulation is presented for solving boundary value problems for Stokes and Darcy flows. The method also applies to problems where the equations are coupled across an interface through appropriate boundary conditions. The adopted technique consists of first reformulating the singular integrals for the fluid quantities as single and double layer potentials. Then the layer potentials are regularized and discretized using standard quadratures. As a final step, the leading term in the regularization error is eliminated in order to gain one more order of accuracy. The numerical examples demonstrate the increase of the convergence rate from first to second order and show a decrease in magnitude of the error. The coupled problems require the computation of the gradient of the Stokes velocity at the common interface. This boundary condition is also written as a combination of single and double layer potentials so that the same approach can be used to compute it accurately. Extensive numerical examples show the increased accuracy gained by the correction terms.     

Laser transmission welding of composites – Part B: Experimental validation of numerical model

In this paper, experimental measurements are performed to confirm the global analytical model (refraction and absorption phenomena) presented in the previous work. Initially, an experimental approach to the estimation of the laser beam spread in a semi-transparent composite and at interface (width of the output beam) subjected to an incident heat flux, was presented. This parameter represents a fundamental input data for the global analytical model (refraction and absorption) during the numerical simulation of transmission infrared welding. Then, an experimental setup for the temperature measurement is performed using infrared camera, during infrared transmission welding of materials joints to validate the welding simulation results (a transient numerical model, based both on conduction and radiation mode heat transfer) with the developed analytical model. The commercial FEM software COMSOL Multiphysics® is used to compute temperature distribution by implementing a radiative source term. Numerical simulations are compared with experimental data. The agreement between simulations and experiments is fair, which gives confidence to use the developed model with acceptable accuracy.     

Analysis of mechanically induced processes in the Langmuir film

The understanding of processes in the monolayer at the air-water interface induced by mechanical compression is important as a part of basic research of the system with reduced dimensionality as well as for the investigation of processes during the Langmuir-Blodgett deposition. The Maxwell displacement current technique provides a substantial contribution for the study of structural and electrical properties. Analysis based on imperfect gas approximation with semi-empirical intermolecular potentials is used. Detail theoretical study of molecular tilt in a continual lateral compression and dielectric relaxation phenomena (step-compression) is presented. Obtained results are confronted with standard surface pressure analysis and surface potential measurement.     

平面上方二维介质目标对高斯波束的电磁散射研究

基于矩量法、互易性定理及镜像理论,提出了一种新的混合方法用于研究水平分界面上方二维介质目标对垂直入射高斯波束的差值散射场.应用镜像理论,介质水平分界面可被原始目标相对于该分界面的镜像目标所替代,从而给出散射问题的等效模型.在等效模型中,应用矩量法求解了原始目标及镜像目标对高斯波束的散射场,同时结合互易性定理得到了原目标与其镜像目标之间的耦合散射场.数值计算结果与相关文献方法及MoM所得结果进行了比较,验证了该混合方法的有效性. 关键词： 互易性定理 电磁散射 高斯波束 二维目标     

Pulse laser acoustics for the characterization of inhomogeneities at interfaces of microstructures

Interfaces between neighbouring materials are often subjected to diffusion processes which cause layers having gradually varying mechanical properties--like densities, Young's moduli or shear moduli--perpendicular to the surface or interface. In this investigation particular interest is drawn on the question how the propagation characteristics of bulk acoustic waves are affected by diffusion layers. The reflection and transmission behavior of bulk acoustic waves encountering a continuum having a spatially dependent sound velocity is discussed based on numerical simulations as well as on experimental verifications. The simulated results are part of an on-going project in which material properties of MEMS devices are investigated by short pulse laser acoustic methods. Mechanical waves are excited and detected thermoelastically using laser pulses of 70 fs duration. For metals this leads to wavelengths of 10-20 nm and the corresponding frequencies amount to 0.3-0.6 THz. In contrast to previous work done in this field in which diffusion effects are generally considered as undesirable phenomena, the deliberate realization of microstructures having well defined gradually varying material properties in one or more dimensions represents a goal of this investigation. For metallic thin film multilayers thermally induced diffusion processes have shown to be an easy and reliable technique for the realization of layered structures having continuously varying mechanical properties within several 10 nm. Among the experimental methods suitable for the in-depth profiling of submicron metallic thin films providing resolutions of several nanometers, are short pulse laser acoustic methods, Rutherford backscattering spectroscopy (RBS), and glow discharge optical emission spectroscopy (GDOES). Short pulse laser acoustic methods and RBS have the advantage to be nondestructive. The short pulse laser acoustic method is described in detail and RBS measurements are presented for verification purposes. Finally potential engineering applications like micro-machined spectrum analyzers, acoustic isolation layers, and band pass filters, operating at very high frequencies are presented.     

Development of Lattice Boltzmann Flux Solver for Simulation of Incompressible Flows

A lattice Boltzmann flux solver (LBFS) is presented in this work for simulation of incompressible viscous and inviscid flows. The new solver is based on Chapman-Enskog expansion analysis, which is the bridge to link Navier-Stokes (N-S) equations and lattice Boltzmann equation (LBE). The macroscopic differential equations are discretized by the finite volume method, where the flux at the cell interface is evaluated by local reconstruction of lattice Boltzmann solution from macroscopic flow variables at cell centers. The new solver removes the drawbacks of conventional lattice Boltzmann method such as limitation to uniform mesh, tie-up of mesh spacing and time interval, limitation to viscous flows. LBFS is validated by its application to simulate the viscous decaying vortex flow, the driven cavity flow, the viscous flow past a circular cylinder, and the inviscid flow past a circular cylinder. The obtained numerical results compare very well with available data in the literature, which show that LBFS has the second order of accuracy in space, and can be well applied to viscous and inviscid flow problems with non-uniform mesh and curved boundary.     

A new modeling approach for mixture fraction statistics based on dissipation elements

The probabilty density function (PDF) of the mixture fraction is of integral importance to a large number of combustion models. Here, a novel modelling approach for the PDF of the mixture fraction is proposed which employs dissipation elements. While being restricted to the commonly used mean and variance of the mixture fraction, this model approach individually considers contributions of the laminar regions as well as the turbulent core and the turbulent/non-turbulent interface region. The later region poses a highly intermittent part of the flow which is of high relevance to the non-premixed combustion of pure hydrocarbon fuels. The model assumptions are justified by means of the gradient trajectory based analysis of high fidelity direct numerical simulation (DNS) datasets of two turbulent inert configurations and a turbulent non-premixed jet flame. The new dissipation element based model is validated against the DNS datasets and a comparison with the beta PDF is presented.     

Gas bubble pulsation in a semiconfined space subjected to ultrasound

In the case of ultrasound application in biological tissues, gas bubbles might form and collapse within cells, in the intercellular spaces and on tissue surfaces. In this work the effect of confined space on the behavior of the gas bubble in the presence of ultrasonic field is studied. A numerical model for bubble pulsation in a planar liquid layer, bounded by two rigid walls, is developed. Surface tension at the interface between the host liquid and the gas in the bubble is considered as well. A mathematical statement and solution technique based on the boundary integral method are presented. In some cases, the bubble divides into two symmetrical parts and high-velocity jets are generated, aimed at the walls. The final velocity of the jets strongly depends on the surface tension of the host liquid. Two new parameters that predict the occurrence of jet formation are developed.     

An RKDG finite element method for the one-dimensional inviscid compressible gas dynamics equations in a Lagrangian coordinate

In this paper,Runge-Kutta Discontinuous Galerkin(RKDG) finite element method is presented to solve the onedimensional inviscid compressible gas dynamic equations in a Lagrangian coordinate.The equations are discretized by the DG method in space and the temporal discretization is accomplished by the total variation diminishing Runge-Kutta method.A limiter based on the characteristic field decomposition is applied to maintain stability and non-oscillatory property of the RKDG method.For multi-medium fluid simulation,the two cells adjacent to the interface are treated differently from other cells.At first,a linear Riemann solver is applied to calculate the numerical ?ux at the interface.Numerical examples show that there is some oscillation in the vicinity of the interface.Then a nonlinear Riemann solver based on the characteristic formulation of the equation and the discontinuity relations is adopted to calculate the numerical ?ux at the interface,which suppresses the oscillation successfully.Several single-medium and multi-medium fluid examples are given to demonstrate the reliability and efficiency of the algorithm.     

A RKDG finite element method for the one-dimensional inviscid compressible gas dynamics equations in Lagrangian coordinate

In this paper, Runge-Kutta Discontinuous Galerkin (RKDG) finite element method is presented to solve the one-dimensional inviscid compressible gas dynamic equations in Lagrangian coordinate. The equations are discretized by the DG method in space and the temporal discretization is accomplished by the total variation diminishing Runge-Kutta method. A limiter based on the characteristic field decomposition is applied to maintain stability and non-oscillatory property of the RKDG method. For multi-medium fluid simulation, the two cells adjacent to the interface are treated differently from other cells. At first, a linear Riemann solver is applied to calculate the numerical flux at the interface. Numerical examples show that there is some oscillation in the vicinity of the interface. Then a nonlinear Riemann solver based on the characteristic formulation of the equation and the discontinuity relations is adopted to calculate the numerical flux at the interface, which suppress the oscillation successfully. Several single-medium and multi-medium fluid examples are given to demonstrate the reliability and efficiency of the algorithm.     

Accelerating flames in tubes—an analysis

Flame acceleration in tubes is studied. A tube filled with flammable mixture is closed at one end and open to the atmosphere at its second end. When ignition takes place near the closed end, it is well-known from experiments that the flame may accelerate, oscillate and eventually reach considerable speeds. A one-dimensional analysis is presented, based upon the assumption that the flame front propagates at a speed that is small compared to the speed of sound. The analysis leads to a construction of the complete unsteady solution. Results from the analysis and from a numerical simulation are compared. They are similar enough to validate the analysis. The tube acoustics are set in motion by the expansion of the fluid due to ignition at the closed end. Subsequently, both spectrum and amplitude evolve because of the motion of the temperature interface, and because of forcing by the flame front, which the analysis precisely quantifies. Oscillations in the front position are strong enough to result in flow reversal. In addition, the induced periodic acoustic acceleration of the temperature and density interface will periodically make the flame front Rayleigh–Taylor unstable, which should result in the dramatic increase in the propagation speed seen in experiments.     

Elastoacoustic model with uncertain mechanical properties for ultrasonic wave velocity prediction: application to cortical bone evaluation

The axial transmission technique can measure the longitudinal wave velocity of an immersed solid. An elementary model of the technique is developed with a set of source and receivers placed in a semi-infinite fluid coupled at a plane interface with a semi-infinite solid. The acoustic fluid is homogeneous. The solid is homogeneous, isotropic, and linearly elastic. The work is focused on the prediction of the measured velocity (apparent velocity) when the solid is considered to have random material properties. The probability density functions of the random variables modeling each mechanical parameter of the solid are derived following the maximum entropy principle. Specific attention is paid to the modeling of Poisson's ratio so that the second-order moments of the velocities remain finite. The stochastic solver is based on a Monte Carlo numerical simulation and uses an exact semianalytic expression of the acoustic response derived with the Cagniard-de Hoop method. Results are presented for a solid with the material properties of cortical bone. The estimated mean values and confidence regions of the apparent velocity are presented for various dispersion levels of the random parameters. A sensibility analysis with respect to the source and receivers locations is presented.     

Multiscale molecular dynamics using the matched interface and boundary method

The Poisson-Boltzmann (PB) equation is an established multiscale model for electrostatic analysis of biomolecules and other dielectric systems. PB based molecular dynamics (MD) approach has a potential to tackle large biological systems. Obstacles that hinder the current development of PB based MD methods are concerns in accuracy, stability, efficiency and reliability. The presence of complex solvent-solute interface, geometric singularities and charge singularities leads to challenges in the numerical solution of the PB equation and electrostatic force evaluation in PB based MD methods. Recently, the matched interface and boundary (MIB) method has been utilized to develop the first second order accurate PB solver that is numerically stable in dealing with discontinuous dielectric coefficients, complex geometric singularities and singular source charges. The present work develops the PB based MD approach using the MIB method. New formulation of electrostatic forces is derived to allow the use of sharp molecular surfaces. Accurate reaction field forces are obtained by directly differentiating the electrostatic potential. Dielectric boundary forces are evaluated at the solvent-solute interface using an accurate Cartesian-grid surface integration method. The electrostatic forces located at reentrant surfaces are appropriately assigned to related atoms. Extensive numerical tests are carried out to validate the accuracy and stability of the present electrostatic force calculation. The new PB based MD method is implemented in conjunction with the AMBER package. MIB based MD simulations of biomolecules are demonstrated via a few example systems.