Non-reflective boundary conditions for non-uniform mean velocity inlet and outlet flows

The numerical solution of problems in unbounded physical space requires a truncation of the computational domain to a reasonable size. As a result, the conditions on the artificial boundaries are generally unknown. Assumptions like constant pressure or velocities are only valid in the far field and lead to spurious reflections if applied on the boundaries of the truncated domain. A number of attempts have been made over the past decades to design conditions that prevent such reflections. One approach is based on characteristics. The standard analysis assumes a spatially uniform mean flow field but this is often impractical. In the present paper we show how to extend the formulation to the more general case of a non-uniform mean velocity field. A number of test cases are provided and our results compare favourably with other boundary conditions. In principle the present approach can be extended to include non-uniformities in all variables.

Shock propagation and MPT noise from a transonic rotor in non-uniform flow

One of the major challenges in hig4h-speed fan stages used in compact, embedded propulsion systems is inlet distortion noise. A body-force-based approach for the prediction of multiple-pure-tone (MPT) noise was previously introduced and validated. In this paper, it is employed with the objective of quantifying the effects of non-uniform flow on the generation and propagation of MPT noise. First-of-their-kind back-to-back coupled aero-acoustic computations were carried out using the new approach for conventional and serpentine inlets. Both inlets delivered flow to the same NASA/GE R4 fan rotor at equal corrected mass flow rates. Although the source strength at the fan is increased by 45 dB in sound power level due to the non-uniform inflow, farfield noise for the serpentine inlet duct is increased on average by only 3.1 dBA overall sound pressure level in the forward arc. This is due to the redistribution of acoustic energy to frequencies below 11 times the shaft frequency and the apparent cut-off of tones at higher frequencies including blade-passing tones. The circumferential extent of the inlet swirl distortion at the fan was found to be 2 blade pitches, or 1/11th of the circumference, suggesting a relationship between the circumferential extent of the inlet distortion and the apparent cut-off frequency perceived in the far field. A first-principles-based model of the generation of shock waves from a transonic rotor in non-uniform flow showed that the effects of non-uniform flow on acoustic wave propagation, which cannot be captured by the simplified model, are more dominant than those of inlet flow distortion on source noise. It demonstrated that non-linear, coupled aerodynamic and aeroacoustic computations, such as those presented in this paper, are necessary to assess the propagation through non-uniform mean flow. A parametric study of serpentine inlet designs is underway to quantify these propagation effects. Copyright © 2011 by ASME.

Analysis of modes in a freestanding microsquare resonator by 3-D finite-difference time-domain

Modes in a microsquare resonator slab with strong vertical waveguide consisting of air/semiconductor/air are analyzed by three-dimensional (3-D) finite-difference time-domain simulation, and compared with that of two-dimensional (2-D) simulation under effective index approximation. Mode frequencies and field distributions inside the resonator obtained by the 3-D simulation are in good agreement with those of the 2-D approximation. However, field distributions at the boundary of the resonator obtained by 3-D simulation are different from that of the 2-D simulation, especially the vertical field distribution near the boundary is great different from that of the slab waveguide, which is used in the effective index approximation. Furthermore the quality factors obtained by 3-D simulation are much larger. than that by 2-D simulation for the square resonator slab with the strong vertical waveguide.

Finite-difference time-domain analysis of deformed square cavity filters with a traveling-wave-like filtering response by mode coupling

An add-drop filter based on a perfect square resonator can realize a maximum of only 25% power dropping because the confined modes are standing-wave modes. By means of mode coupling between two modes with inverse symmetry properties, a traveling-wave-like filtering response is obtained in a two-dimensional single square cavity filter with cut or circular corners by finite-difference time-domain simulation. The optimized deformation parameters for an add-drop filter can be accurately predicted as the overlapping point of the two coupling modes in an isolated deformed square cavity. More than 80% power dropping can be obtained in a deformed square cavity filter with a side length of 3.01 mu m. The free spectral region is decided by the mode spacing between modes, with the sum of the mode indices differing by 1. (c) 2007 Optical Society of America.

Investigation of mode characteristics for a square cavity with a pedestal by a three-dimensional finite-difference time-domain technique

Mode characteristics of a strongly confined square cavity suspended in air via a pedestal on the substrate are investigated by a three-dimensional finite-difference time-domain technique. The mode wavelengths and mode quality factors (Q factors) are calculated as the functions of the size of the pedestal and the slope angle 0 of the sidewalls of the square slab, respectively For the square slab with side length of 2 mu m, thickness of 0.2 mu m, and refractive index of 3.4, on a square pedestal with refractive index of 3.17, the Q factor of the whispering-gallery (WG)-like mode transverse-electric TE(3.5)o first increases with the side length b of the square pedestal and then quickly decreases as b > 0.4 mu m, but the Q factor of the WG-like mode TE(4.6)o drops down quickly as b > 0.2 mu m, owing to their different symmetries. The results indicate that the pedestal can also result in mode selection in the WG-like modes. In addition, the numerical results show that the Q factors decrease 50% as the slope angle of the sidewalls varies from 90 degrees to 80 degrees. The mode characteristics of WG-like modes in the square cavity with a rectangular pedestal are also discussed. The results show that the nonsquare pedestal largely degrades the WG-like modes. (c) 2006 Optical Society of America

Investigation of mode characteristics for microdisk resonators by S-matrix and three-dimensional finite-difference time-domain technique

The mode characteristics of a three-dimensional (3D) microdisk with a vertical refractive index distribution of n(2)/3.4/n(2) are investigated by the S-matrix method and 3D finite-difference time-domain (FDTD) technique. For the microdisk with a thickness of 0.2 mu m. and a radius of 1 mu m, the mode wavelengths and quality factors for the HE7,1 mode obtained by 3D FDTD simulation and the S-matrix method are in good agreement as n(2) increases from 1.0 to 2.6. But the Q factor obtained by the 3D FDTD rapidly decreases from 1.12 X 10(4) to 379 as n2 increases from 2.65 to 2.8 owing to the vertical radiation losses, which cannot be predicted by the proposed S-matrix method. The comparisons also show that quality factors obtained from the analytical solution of two-dimensional microdisks under the effective index approximation are five to seven times smaller than those of the 3D FDTD as n(2) = 1 and R = 1 mu m. (c) 2006 Optical Society of America.

Mode-coupling analysis of three-dimensional microdisk resonators by the finite-difference time-domain technique

Quality factor enhancement due to mode coupling is observed in a three-dimensional microdisk resonator. The microdisk, which is vertically sandwiched between air and a substrate, with a radius of 1 mu m, a thickness of 0.2 mu m, and a refractive index of 3.4, is considered in a finite-difference time-domain (FDTD) numerical simulation. The mode quality factor of the fundamental mode HE71 decreases with an increase of the refractive index of the substrate, n(sub), from 2.0 to 3.17. However, the mode quality factor of the first-order mode HE72 reaches a peak value at n(sub) = 2.7 because of the mode coupling between the fundamental and the first-order modes. The variation of mode field distributions due to the mode coupling is also observed. This mechanism may be used to realize high-quality-factor modes in microdisks with high-refractive-index substrates. (c) 2006 Optical Society of America.

Digital simulation in a thin layer spectroelectrochemical cell with minigrid platinum electrode by microregion approximation explicit finite difference method

The microregion approximation explicit finite difference method is used to simulate cyclic voltammetry of an electrochemical reversible system in a three-dimensional thin layer cell with minigrid platinum electrode. The simulated CV curve and potential scan-absorbance curve were in very good accordance with the experimental results, which differed from those at a plate electrode. The influences of sweep rate, thickness of the thin layer, and mesh size on the peak current and peak separation were also studied by numerical analysis, which give some instruction for choosing experimental conditions or designing a thin layer cell. The critical ratio (1.33) of the diffusion path inside the mesh hole and across the thin layer was also obtained. If the ratio is greater than 1.33 by means of reducing the thickness of a thin layer, the electrochemical property will be far away from the thin layer property.

地震波模拟：非规则网格离散的数值解法、吸收边界及其实际应用

This dissertation presents a series of irregular-grid based numerical technique for modeling seismic wave propagation in heterogeneous media. The study involves the generation of the irregular numerical mesh corresponding to the irregular grid scheme, the discretized version of motion equations under the unstructured mesh, and irregular-grid absorbing boundary conditions. The resulting numerical technique has been used in generating the synthetic data sets on the realistic complex geologic models that can examine the migration schemes. The motion equation discretization and modeling are based on Grid Method. The key idea is to use the integral equilibrium principle to replace the operator at each grid in Finite Difference scheme and variational formulation in Finite Element Method. The irregular grids of complex geologic model is generated by the Paving Method, which allow varying grid spacing according to meshing constraints. The grids have great quality at domain boundaries and contain equal quantities of nodes at interfaces, which avoids the interpolation of parameters and variables. The irregular grid absorbing boundary conditions is developed by extending the Perfectly Matched Layer method to the rotated local coordinates. The splitted PML equations of the first-order system is derived by using integral equilibrium principle. The proposed scheme can build PML boundary of arbitrary geometry in the computational domain, avoiding the special treatment at corners in a standard PML method and saving considerable memory and computation cost. The numerical implementation demonstrates the desired qualities of irregular grid based modeling technique. In particular, (1) smaller memory requirements and computational time are needed by changing the grid spacing according to local velocity; (2) Arbitrary surfaces and interface topographies are described accurately, thus removing the artificial reflection resulting from the stair approximation of the curved or dipping interfaces; (3) computational domain is significantly reduced by flexibly building the curved artificial boundaries using the irregular-grid absorbing boundary conditions. The proposed irregular grid approach is apply to reverse time migration as the extrapolation algorithm. It can discretize the smoothed velocity model by irregular grid of variable scale, which contributes to reduce the computation cost. The topography. It can also handle data set of arbitrary topography and no field correction is needed.

Fully developed laminar flow and heat transfer in smooth trapezoidal microchannel

Numerical analysis of fully developed laminar slip flow and heat transfer in trapezoidal micro-channels has been studied with uniform wall heat flux boundary conditions. Through coordinate transformation, the governing equations are transformed from physical plane to computational domain, and the resulting equations are solved by a finite-difference scheme. The influences of velocity slip and temperature jump on friction coefficient and Nusselt number are investigated in detail. The calculation also shows that the aspect ratio and base angle have significant effect on flow and heat transfer in trapezoidal micro-channel. (c) 2005 Elsevier Ltd. All rights reserved.

Non-Uniform Reductions

We study properties of non-uniform reductions and related completeness notions. We strengthen several results of Hitchcock and Pavan and give a trade-off between the amount of advice needed for a reduction and its honesty on NEXP. We construct an oracle relative to which this trade-off is optimal. We show, in a more systematic study of non-uniform reductions, that among other things non-uniformity can be removed at the cost of more queries. In line with Post's program for complexity theory we connect such 'uniformization' properties to the separation of complexity classes.

Coupled 3-D finite difference time domain and finite volume methods for solving microwave heating in porous media

Computational results for the microwave heating of a porous material are presented in this paper. Combined finite difference time domain and finite volume methods were used to solve equations that describe the electromagnetic field and heat and mass transfer in porous media. The coupling between the two schemes is through a change in dielectric properties which were assumed to be dependent both on temperature and moisture content. The model was able to reflect the evolution of temperature and moisture fields as the moisture in the porous medium evaporates. Moisture movement results from internal pressure gradients produced by the internal heating and phase change.

Numerical simulation of flow-induced cavity noise in self-sustained oscillations

The generation and near-field radiation of aerodynamic sound from a low-speed unsteady flow over a two-dimensional automobile door cavity is simulated by using a source-extraction-based coupling method. In the coupling procedure, the unsteady cavity flow field is first computed solving the Reynolds averaged Navier–Stokes (RANS) equations. The radiated sound is then calculated by using a set of acoustic perturbation equations with acoustic source terms which are extracted from the time-dependent solutions of the unsteady flow. The aerodynamic and its resulting acoustic field are computed for the Reynolds number of 53,266 based on the base length of the cavity. The free stream flow velocity is taken to be 50.9m/s. As first stage of the numerical investigation of flow-induced cavity noise, laminar flow is assumed. The CFD solver is based on a cell-centered finite volume method. A dispersion-relation-preserving (DRP), optimized, fourth-order finite difference scheme with fully staggered-grid implementation is used in the acoustic solver