General-domain compressible Navier-Stokes solvers exhibiting quasi-unconditional stability and high-order accuracy in space and time

This thesis presents a new class of solvers for the subsonic compressible Navier-Stokes equations in general two- and three-dimensional spatial domains. The proposed methodology incorporates: 1) A novel linear-cost implicit solver based on use of higher-order backward differentiation formulae (BDF) and the alternating direction implicit approach (ADI); 2) A fast explicit solver; 3) Dispersionless spectral spatial discretizations; and 4) A domain decomposition strategy that negotiates the interactions between the implicit and explicit domains. In particular, the implicit methodology is quasi-unconditionally stable (it does not suffer from CFL constraints for adequately resolved flows), and it can deliver orders of time accuracy between two and six in the presence of general boundary conditions. In fact this thesis presents, for the first time in the literature, high-order time-convergence curves for Navier-Stokes solvers based on the ADI strategy---previous ADI solvers for the Navier-Stokes equations have not demonstrated orders of temporal accuracy higher than one. An extended discussion is presented in this thesis which places on a solid theoretical basis the observed quasi-unconditional stability of the methods of orders two through six. The performance of the proposed solvers is favorable. For example, a two-dimensional rough-surface configuration including boundary layer effects at Reynolds number equal to one million and Mach number 0.85 (with a well-resolved boundary layer, run up to a sufficiently long time that single vortices travel the entire spatial extent of the domain, and with spatial mesh sizes near the wall of the order of one hundred-thousandth the length of the domain) was successfully tackled in a relatively short (approximately thirty-hour) single-core run; for such discretizations an explicit solver would require truly prohibitive computing times. As demonstrated via a variety of numerical experiments in two- and three-dimensions, further, the proposed multi-domain parallel implicit-explicit implementations exhibit high-order convergence in space and time, useful stability properties, limited dispersion, and high parallel efficiency.

Observing the fractional Fourier transform by free-space Fresnel diffraction

The fractional Fourier transform of an object can be observed in the free-space Fresnel diffraction pattern of the object. (C) 1997 Optical Society of America

Scaled fractional Fourier transform and its optical implementation

The scaled fractional Fourier transform is suggested and is implemented optically by one lens for different values of phi and output scale. In addition, physically it relates the FRT with the general lens transform-the optical diffraction between two asymmetrically positioned planes before and after a lens. (C) 1997 Optical Society of America.

Extended fractional Fourier transforms

The concept of an extended fractional Fourier transform (FRT) is suggested. Previous PBT's and complex FRT's are only its subclasses. Then, through this concept and its method, we explain the physical meaning of any optical Fresnel diffraction through a lens: It is just an extended FRT; a lens-cascaded system can equivalently be simplified to a simple analyzer of the FRT; the two-independent-parameter FRT of an object illuminated with a plane wave can be readily implemented by a lens of arbitrary focal length; when cascading, the Function of each lens unit and the relationship between the adjacent ones are clear and simple; and more parameters and fewer restrictions on cascading make the optical design easy. (C) 1997 Optical Society of America.

Highly sensitive wave-front sensor with a non-zero-order phase plate

We propose a novel highly sensitive wave front detection method for a quick check of a flat wave front by taking advantage of a non-zero-order pi phase plate that yields a non-zero-order diffraction pattern. When a light beam with a flat wave front illuminates a phase plate, the zero-order intensity is zero. When there is a slight distortion of the wave front, the zero-order intensity increases. The ratio of first-order intensity to that of zero-order intensity is used as the criterion with which to judge whether the wave front under test is flat, eliminating the influence of background light. Experimental results demonstrate that this method is efficient, robust, and cost-effective and should be highly interesting for a quick check of a flat wave front of a large-aperture laser beam and adaptive optical systems. (c) 2005 Optical Society of America.

A search for fractionally charged particles in cosmic rays and a theoretical interpretation of the results

An array of two spark chambers and six trays of plastic scintillation counters was used to search for unaccompanied fractionally charged particles in cosmic rays near sea level. No acceptable events were found with energy losses by ionization between 0.04 and 0.7 that of unit-charged minimum-ionizing particles. New 90%-confidence upper limits were thereby established for the fluxes of fractionally charged particles in cosmic rays, namely, (1.04 ± 0.07)x10^-10 and (2.03 ± 0.16)x10^-10 cm^-2sr^-1sec^-1 for minimum-ionizing particles with charges 1/3 and 2/3, respectively.

In order to be certain that the spark chambers could have functioned for the low levels of ionization expected from particles with small fractional charges, tests were conducted to estimate the efficiency of the chambers as they had been used in this experiment. These tests showed that the spark-chamber system with the track-selection criteria used might have been over 99% efficient for the entire range of energy losses considered.

Lower limits were then obtained for the mass of a quark by considering the above flux limits and a particular model for the production of quarks in cosmic rays. In this model, which is one involving the multi-peripheral Regge hypothesis, the production cross section and a corresponding mass limit are critically dependent on the Regge trajectory assigned to a quark. If quarks are "elementary'' with a flat trajectory, the mass of a quark can be expected to be at least 6 ± 2 BeV/c². If quarks have a trajectory with unit slope, just as the existing hadrons do, the mass of a quark might be as small as 1.3 ± 0.2 BeV/c². For a trajectory with unit slope and a mass larger than a couple of BeV/c², the production cross section may be so low that quarks might never be observed in nature.

The effect of the group delay ripple of chirped fiber grating on composite second-order in optical fiber CATV system

The effect of group delay ripple of chirped fiber gratings on composite second-order (CSO) performance in optical fiber CATV system is investigated. We analyze the system CSO performances for different ripple amplitudes, periods and residual dispersion amounts in detail. It is found that the large ripple amplitude and small ripple period will deteriorate the system CSO performance seriously. Additionally, the residual dispersion amount has considerable effect on CSO performance in the case of small ripple amplitude and large ripple period. (c) 2004 Elsevier B.V. All rights reserved.

Composite second-order performance improvement in optical fibre CATV transmission system using chirped fibre grating

Theoretically, we analyse the dispersion compensation characteristics of the chirped fibre grating (CFG) in an optical fibre cable television (CATV) system and obtain the analytic expression of the composite second-order (CSO) distortion using the time-domain form of the field envelope wave equation. The obtained result is in good agreement with the numerical simulation result. Experimentally, we verify the result by making use of the tunable characteristics of CFG to change the dispersion compensation amount and obtain an optimal CSO performance in a 125km fibre transmission link. Both the theoretical and experimental results show that the CSO performance can be improved by properly choosing the dispersion compensation amount for a certain fibre transmission link.

Analysis of a second-order polarization-independent filter made of a single ring resonator and a Sagnac interferometer

A novel second-order polarization-independent filter made of a single ring resonator and a Sagnac interferometer (SRRSI) is proposed, and its filtering characteristics are investigated. By using birefringence in waveguide, a single ring resonator can be used to synthesize a filter with second-order response. Analytical formulas are derived for characteristics of the SRRSI varied with waveguide parameters.. such as the coupling coefficient; and the critical condition of a second-order Butterworth filter is given. The influence of loss in the ring resonator is also analyzed. (c) 2005 Elsevier B.V. All rights reserved.

A generalized treatment of the order-disorder transformation in alloys and its effects on their magnetic properties

A theory of the order-disorder transformation is developed in complete generality. The general theory is used to calculate long range order parameters, short range order parameters, energy, and phase diagrams for a face centered cubic binary alloy. The theoretical results are compared to the experimental determination of the copper-gold system, Values for the two adjustable parameters are obtained.

An explanation for the behavior of magnetic alloys is developed, Curie temperatures and magnetic moments of the first transition series elements and their alloys in both the ordered and disordered states are predicted. Experimental agreement is excellent in most cases. It is predicted that the state of order can effect the magnetic properties of an alloy to a considerable extent in alloys such as Ni₃Mn. The values of the adjustable parameter used to fix the level of the Curie temperature, and the adjustable parameter that expresses the effect of ordering on the Curie temperature are obtained.

I. Numerical solution of exact pair equations. II. Numerical solution of first-order pair equations

Part I

Numerical solutions to the S-limit equations for the helium ground state and excited triplet state and the hydride ion ground state are obtained with the second and fourth difference approximations. The results for the ground states are superior to previously reported values. The coupled equations resulting from the partial wave expansion of the exact helium atom wavefunction were solved giving accurate S-, P-, D-, F-, and G-limits. The G-limit is -2.90351 a.u. compared to the exact value of the energy of -2.90372 a.u.

Part II

The pair functions which determine the exact first-order wavefunction for the ground state of the three-electron atom are found with the matrix finite difference method. The second- and third-order energies for the (1s1s)¹S, (1s2s)³S, and (1s2s)¹S states of the two-electron atom are presented along with contour and perspective plots of the pair functions. The total energy for the three-electron atom with a nuclear charge Z is found to be E(Z) = -1.125•Z² +1.022805•Z-0.408138-0.025515•(1/Z)+O(1/Z²)a.u.