Methods for implementing the stream boundary condition in DSMC computations

In the direct simulation Monte‐Carlo (DSMC) method for simulating rarefied gas flows, the velocities of simulator particles that cross a simulation boundary and enter the simulation space are typically generated using the acceptance–rejection procedure that samples the velocities from a truncated theoretical velocity distribution that excludes low and high velocities. This paper analyses an alternative technique, where the velocities of entering particles are obtained by extending the simulation procedures to a region adjacent to the simulation space, and considering the movement of particles generated within that region during the simulation time step. The alternative method may be considered as a form of acceptance–rejection procedure, and permits the generation of all possible velocities, although the population of high velocities is depleted with respect to the theoretical distribution. Nevertheless, this is an improvement over the standard acceptance–rejection method. Previous implementations of the alternative method gave a number flux lower than the theoretical number required. Two methods for obtaining the correct number flux are presented. For upstream boundaries in high‐speed flows, the alternative method is more computationally efficient than the acceptance–rejection method. However, for downstream boundaries, the alternative method is extremely inefficient. The alternative method, with the correct theoretical number flux, should therefore be used in DSMC computations in favour of the acceptance–rejection method for upstream boundaries in high‐speed flows. Copyright © 2003 John Wiley & Sons, Ltd.     

Computing Zeta Functions of Kummer Curves via Multiplicative Characters

We present a practical polynomial-time algorithm for computing the zeta function of a Kummer curve over a finite field of small characteristic. Such algorithms have recently been obtained using a method of Kedlaya based upon Monsky–Washnitzer cohomology, and are of interest in cryptography. We take a different approach. The problem is reduced to that of computing the L-function of a multiplicative character sum. This latter task is achieved via a cohomological formula based upon the work of Dwork and Reich. We show, however, that our method and that of Kedlaya are very closely related.Dedicated to the memory of Gian-Carlo Rota     

LES of turbulent heat transfer: proper convection numerical schemes for temperature transport

Large eddy simulations of two basic configurations (decay of isotropic turbulence, and the academic plane channel flow) with heat transfer have been performed comparing several convection numerical schemes, in order to discuss their ability to evaluate temperature fluctuations properly. Results are compared with the available incompressible heat transfer direct numerical simulation data. It is shown that the use of regularizing schemes (such as high order upwind type schemes) for the temperature transport equation in combination with centered schemes for momentum transport equation gives better results than the use of centred schemes for both equations. Copyright © 2004 John Wiley & Sons, Ltd.     

Modelling of a single drop impact onto liquid film using particle method

The process of single liquid drop impact on thin liquid surface is numerically simulated with moving particle semi‐implicit method. The mathematical model involves gravity, viscosity and surface tension. The model is validated by the simulation of the experimental cases. It is found that the dynamic processes after impact are sensitive to the liquid pool depth and the initial drop velocity. In the cases that the initial drop velocity is low, the drop will be merged with the liquid pool and no big splash is seen. If the initial drop velocity is high enough, the dynamic process depends on the liquid depth. If the liquid film is very thin, a bowl‐shaped thin crown is formed immediately after the impact. The total crown subsequently expands outward and breaks into many tiny droplets. When the thickness of the liquid film increases, the direction of the liquid crown becomes normal to the surface and the crown propagates outward. It is also found that the radius of the crown is described by a square function of time: r_C = [c(t ? t₀)]^0.5. When the liquid film is thick enough, a crown and a deep cavity inside it are formed shortly after the impact. The bottom of the cavity is initially oblate and then the base grows downward to form a sharp corner and subsequently the corner moves downward. Copyright © 2004 John Wiley & Sons, Ltd.     

Photorefractive Four-Wave Mixing Switches Utilizing Pockels Effect -Longitudinal and Lateral Switches-

We carried out detailed calculations for photorefractive wave-mixing switches based on one of three crystals with high electro-optic coefficients, namely, BaTiO₃, Strontium Barium Niobate (SBN (0.75)), and Potasium Sodium Strontium Barium Niobate (KNSBN). A comparison of results for the three crystals shows that a 0_-cut BaTiO³ crystal is suitable for a longitudinal switch and requires a voltage of about 80 for a 2-mm-thick crystal to induce sufficient phase mismatch. The electrodes must be transparent for the incident and diffracted beams. A 45_-cut SBN (0.75) crystal, however, is suitable for a lateral switch and requires a voltage of about 150 for a 1-mm-wide crystal. The electrodes do not need to be transparent.     

Light Diffraction of Aligned Polymer Fibers Periodically Dispersed by Phase Separation of Liquid Crystal and Polymer

We have confirmed light diffraction of aligned polymer fibers obtained by a phase separation of an anisotropic-phase solution of liquid crystal and polymer. He—Ne laser light passing through the polymer fibers was scattered in the axis vertical to the fibers, and had two peaks of light intensity symmetrical to the center of the transmitting laser spot. The two peaks were found to be caused by light diffraction due to the periodic polymer-fiber dispersion because the peaks corresponded to values calculated by intervals between the fibers. The periodical fiber networks are considered to be formed by anisotropic spinodal decomposition. This effect can be used to measure the dispersion order of the polymer fibers. © 2004 The Optical Society of Japan     

Measurements of Time-Varying Characteristics of Optical Signals and Inhomogeneities of Tested Substances with Additional Phase Modulation

The problem of investigation of the amplitude and phase structure of a time-varying probing optical signal and the structure of time-varying inhomogeneities of a substance tested by this signal is considered. The analysis is concerned, in particular, with determination of the structure of signals and processes with resolution in the pico- and femtosecond range. The scheme used for the analysis is based on registration of four spatially separated spectra of the studied radiation. The spectra are formed in a four-channel scheme with a twin-wave Michelson interferometer and a spectral device. Modulators based on electrooptical crystals (perovskites) are placed in the channels. The sum spectra are formed: without modulators, with the effect of either of the modulators, and with both of them affecting the radiation. The effect of the studied substance implies either modulating the radiation (in this case it is described by multiplication) or redistributing the radiation (then it is described by convolution).     

Physical Interpretation of the Motion of a Laser Pulse within the Framework of the Four-Oscillator Model

A physical multidimensional local model of the propagation process of a laser pulse in air having no analogy to famous models in laser physics is constructed. It is based on the representations of classical wave mechanics characteristic of the theory of scalar wave fields.     

New liquid‐crystalline poly(ester amide)s: The role of nitro groups in the phase behavior

New hydrogen‐bonded liquid‐crystalline poly(ester amide)s (PEA)s were obtained from 1,4‐terephthaloyl[bis‐(3‐nitro‐N‐anthranilic acid)] (5) or 1,4‐terephthaloyl[bis‐(N‐anthranilic acid)] (6), with or without nitro groups, respectively, through the separate condensation of each with hydroquinone or dihydroxynaphthalene. The dicarboxylic monomers were synthesized from 2‐aminobenzoic acid. The phase behavior of the monomers and polymers were studied with differential scanning calorimetry, polarized light microscopy, and wide‐angle X‐ray diffraction methods. Monomer 5, containing nitro groups, exhibited a smectic liquid‐crystalline phase, whereas the texture of monomer 6 without nitro groups appeared to be nematic. The PEAs containing nitro groups exhibited polymorphism (smectic and nematic), whereas those without nitro groups exhibited only one phase transition (a nematic threaded texture). The changes occurring in the phase behavior of the polymers were explained by the introduction of nitro groups. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 1289–1298, 2004     

Shear banding in strained semicrystalline polyamide 6 films as revealed by atomic force microscopy: Role of the amorphous phase

Atomic force microscopy (AFM) has been used to visualize the plastic deformation mechanisms that are responsible for the yielding of semicrystalline polymers of low degree of crystallinity (<50%). Indeed, AFM, if operated in suitable conditions, is able to image both the amorphous and the crystalline phases. Polyamide 6 films have been drawn at temperatures T < 160 °C. Postmortem AFM observations show that, at yield, shear bands nucleate and propagate in the amorphous phase. They cross the crystalline lamellae and run over the whole surface of the sample. By crossing the lamellae, they form nanoblocks of uniform size. Neither the size of the nanoblocks nor the angle between the tensile axis and the shear bands can be explained in terms of crystal plasticity. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 687–701, 2004