Overlapping Domain Decomposition Algorithms for General Sparse Matrices

Domain decomposition methods for finite element problems using a partition based on the underlying finite element mesh have been extensively studied. In this paper, we discuss algebraic extensions of the class of overlapping domain decomposition algorithms for general sparse matrices. The subproblems are created with an overlapping partition of the graph corresponding to the sparsity structure of the matrix. These algebraic domain decomposition methods are especially useful for unstructured mesh problems. We also discuss some difficulties encountered in the algebraic extension, particularly the issues related to the coarse solver.     

The point source method for inverse scattering in the time domain

Many recent inverse scattering techniques have been designed for single frequency scattered fields in the frequency domain. In practice, however, the data is collected in the time domain. Frequency domain inverse scattering algorithms obviously apply to time‐harmonic scattering, or nearly time‐harmonic scattering, through application of the Fourier transform. Fourier transform techniques can also be applied to non‐time‐harmonic scattering from pulses. Our goal here is twofold: first, to establish conditions on the time‐dependent waves that provide a correspondence between time domain and frequency domain inverse scattering via Fourier transforms without recourse to the conventional limiting amplitude principle; secondly, we apply the analysis in the first part of this work toward the extension of a particular scattering technique, namely the point source method, to scattering from the requisite pulses. Numerical examples illustrate the method and suggest that reconstructions from admissible pulses deliver superior reconstructions compared to straight averaging of multi‐frequency data. Copyright © 2006 John Wiley & Sons, Ltd.     

COMPUTATION OF NONEQUILIBRIUM HYPERSONIC FLOW OVER CONCAVE CORNERS

1. IntroductionWe consider steady inviscid hyPersonic flow of air about a concave corner including chendcalreaction effects. Computation of nonequllibrium flow is duncult because of the steep gradielltsbehind the shock and an entroPy layer near the body For many cases the time scale of theChemical reactions is 1arger than, or close to, the time scale of the origina1 meChanical problemand the equllibrium chendstry mode1 is not realistic. Therefore, it is necessary to eva1uate thekinetics of c…     

A vector-parallel scheme for Navier-Stokes computations at multi-gigaflop performance rates

A class of vector-parallel schemes for solution of steady compressible or incompressible viscous flow is developed and performance studies carried out. The algorithms employ an artificial transient treatment that permits rapid integration to a steady state. In the present work a four-stage explicit Runge-Kutta scheme employing variable local step size is utilized for the ODE system integration. The RK-4 scheme is restructured to allow vectorization and enhance concurrency in the calculation for a streamfunction-vorticity formulation of the flow problem. The parameters of the resulting RK scheme can be selected to accelerate convergence of the RK recursion. Four main procedures are considered which permit vector-parallel solution: a Jacobi update, a hybrid of the Jacobi and Gauss-Seidel method, red-black ordering and domain decomposition. Numerical performance studies are conducted with a representative viscous incompressible flow calculation. Results indicate that a scheme involving domain decomposition with a Gauss-Seidel type of update for the RK four-stage scheme is most effective and provides performance in excess of 8 Gflops on the Cray C-90.     

CFD‐based optimization of aerofoils using radial basis functions for domain element parameterization and mesh deformation

A novel domain element shape parameterization method is presented for computational fluid dynamics‐based shape optimization. The method is to achieve two aims: (1) provide a generic ‘wrap‐around’ optimization tool that is independent of both flow solver and grid generation package and (2) provide a method that allows high‐fidelity aerodynamic optimization of two‐ and three‐dimensional bodies with a low number of design variables. The parameterization technique uses radial basis functions to transfer domain element movements into deformations of the design surface and corresponding aerodynamic mesh, thus allowing total independence from the grid generation package (structured or unstructured). Independence from the flow solver (either inviscid, viscous, aeroelastic) is achieved by obtaining sensitivity information for an advanced gradient‐based optimizer (feasible sequential quadratic programming) by finite‐differences. Results are presented for two‐dimensional aerofoil inverse design and drag optimization problems. Inverse design results demonstrate that a large proportion of the design space is feasible with a relatively low number of design variables using the domain element parameterization. Heavily constrained (in lift, volume, and moment) two‐dimensional aerofoil drag optimization has shown that significant improvements over existing designs can be achieved using this method, through the use of various objective functions. Copyright © 2008 John Wiley & Sons, Ltd.     

An X-band time domain electron paramagnetic resonance spectrometer

A computer-controlled X-band time domain electron paramagnetic resonance (EPR) spectrometer, with a time resolution of the order of 0.5μsec, has been constructed with many of the crucial microwave components designed and fabricated by the Microwave Engineering Group of TIFR. The spectrometer operates either in a microwave power pulsed mode for determination of spin-lattice relaxation times by the saturation recovery technique or in the kinetic mode for determination of the time dependence of EPR signal after laser excitation. It has an automatic frequency control, an automatic phase control and, most importantly, a field-frequency lock which ensures good stability of the EPR line positions enabling signal averaging for extended periods. The constructional details of the spectrometer and its performance in both the modes are described here by reporting results on certain typical systems.     

二元矩形软X射线衍射光栅的标量区衍射理论

利用标量区衍射理论分析了正交二元矩形软Ｘ射线衍射光栅的光学特性，给出了理论设计的初步结果。     

Propagation of elliptical flat-topped beam in spatial-frequency domain

Based on a tensor method, the propagation formulae of elliptical flat-topped beam through aligned and misaligned optical systems in spatial-frequency domain are derived analytically. The derived formulae provide a powerful tool for treating the propagation and transformation of elliptical flat-topped beam through complex optical systems in spatial-frequency domain. As a numerical example, the propagation properties of elliptical flat-topped Gaussian beam in free space are studied in spatial-frequency domain.     

Mie scattering calculation by FDTD employing a modified Debye model for Gold material

The parameters calculation for the modified Debye model in the electromagnetic finite difference time domain (FDTD) method is discussed. Validation for the selected parameters is done comparing the FDTD results with the Mie analytical solutions of a gold nanometer sphere under an optical plane wave incidence. Good agreement with an average relative error 6.7% has been demonstrated as to total electrical field amplitudes on the surface of the gold sphere. It indicates FDTD could be adopted in the optical region for a material even with a negative real relative permittivity particularly as metals.     

一类各向异性外问题的非重叠型区域分解算法

In this paper, based on the natural integral operator on elliptic boundary, a nonoverlapping domain decomposition method is presented for a kind of anisotropic elliptic problem with constant coefficients in an exterior domain, and the convergence of the method is analyzed. The choice of the relaxition factor is discussed.Some numerical examples are given. Theoretical analysis as well as numerical examples show that our method is performance.