Iterative spectral methods and spectral solutions to compressible flows

A spectral multigrid scheme is described which can solve pseudospectral discretizations of self-adjoint elliptic problems in O(N log N) operations. An iterative technique for efficiently implementing semi-implicit time-stepping for pseudospectral discretizations of Navier-Stokes equations is discussed. This approach can handle variable coefficient terms in an effective manner. Pseudospectral solutions of compressible flow problems are presented. These include one dimensional problems and two dimensional Euler solutions. Results are given both for shock-capturing approaches and for shock-fitting ones

Suspended liquid particle disturbance on laser-induced blast wave and low density distribution

The impurity effect of suspended liquid particles on the laser-induced gas breakdown was experimentally investigated in quiescent gas. The focus of this study is the investigation of the influence of the impurities on the shock wave structure as well as the low density distribution. A 532 nm Nd:YAG laser beam with an 188 mJ/pulse was focused on the chamber filled with suspended liquid particles 0.9 ± 0.63 μm in diameter. Several shock waves are generated by multiple gas breakdowns along the beam path in the breakdown with particles. Four types of shock wave structures can be observed: (1) the dual blast waves with a similar shock radius, (2) the dual blast waves with a large shock radius at the lower breakdown, (3) the dual blast waves with a large shock radius at the upper breakdown, and (4) the triple blast waves. The independent blast waves interact with each other and enhance the shock strength behind the shock front in the lateral direction. The triple blast waves lead to the strongest shock wave in all cases. The shock wave front that propagates toward the opposite laser focal spot impinges on one another, and thereafter a transmitted shock wave (TSW) appears. The TSW interacts with the low density core called a kernel; the kernel then longitudinally expands quickly due to a Richtmyer-Meshkov-like instability. The laser-particle interaction causes an increase in the kernel volume which is approximately five times as large as that in the gas breakdown without particles. In addition, the laser-particle interaction can improve the laser energy efficiency

Spectral collocation methods

This review covers the theory and application of spectral collocation methods. Section 1 describes the fundamentals, and summarizes results pertaining to spectral approximations of functions. Some stability and convergence results are presented for simple elliptic, parabolic, and hyperbolic equations. Applications of these methods to fluid dynamics problems are discussed in Section 2

Hybrid Spectral Difference/Embedded Finite Volume Method for Conservation Laws

A novel hybrid spectral difference/embedded finite volume method is introduced in order to apply a discontinuous high-order method for large scale engineering applications involving discontinuities in the flows with complex geometries. In the proposed hybrid approach, the finite volume (FV) element, consisting of structured FV subcells, is embedded in the base hexahedral element containing discontinuity, and an FV based high-order shock-capturing scheme is employed to overcome the Gibbs phenomena. Thus, a discontinuity is captured at the resolution of FV subcells within an embedded FV element. In the smooth flow region, the SD element is used in the base hexahedral element. Then, the governing equations are solved by the SD method. The SD method is chosen for its low numerical dissipation and computational efficiency preserving high-order accurate solutions. The coupling between the SD element and the FV element is achieved by the globally conserved mortar method. In this paper, the 5th-order WENO scheme with the characteristic decomposition is employed as the shock-capturing scheme in the embedded FV element, and the 5th-order SD method is used in the smooth flow field. The order of accuracy study and various 1D and 2D test cases are carried out, which involve the discontinuities and vortex flows. Overall, it is shown that the proposed hybrid method results in comparable or better simulation results compared with the standalone WENO scheme when the same number of solution DOF is considered in both SD and FV elements.Comment: 27 pages, 17 figures, 2 tables, Accepted for publication in the Journal of Computational Physics, April 201

The Center of Excellence for Hypersonics Training and Research at the University of Texas at Austin

Over the period of this grant (1986-92), 23 graduate students were supported by the Center and received education and training in hypersonics through MS and Ph.D. programs. An additional 8 Ph.D. candidates and 2 MS candidates, with their own fellowship support, were attracted to The University of Texas and were recruited into the hypersonics program because of the Center. Their research, supervised by the 10 faculty involved in the Center, resulted in approximately 50 publications and presentations in journals and at national and international technical conferences. To provide broad-based training, a new hypersonics curriculum was created, enabling students to take 8 core classes in theoretical, computational, and experimental hypersonics, and other option classes over a two to four semester period. The Center also developed an active continuing education program. The Hypersonics Short Course was taught 3 times, twice in the USA and once in Europe. Approximately 300 persons were attracted to hear lectures by more than 25 of the leading experts in the field. In addition, a hypersonic aerodynamics short course was offered through AIAA, as well as short courses on computational fluid dynamics (CFD) and advanced CFD. The existence of the Center also enabled faculty to leverage a substantial volume of additional funds from other agencies, for research and graduate student training. Overall, this was a highly successful and highly visible program

Numerical solution of the Navier-Stokes equations for blunt nosed bodies in supersonic flows

A time dependent, two dimensional Navier-Stokes code employing the method of body fitted coordinate technique was developed for supersonic flows past blunt bodies of arbitrary shapes. The bow shock ahead of the body is obtained as part of the solution, viz., by shock capturing. A first attempt at mesh refinement in the shock region was made by using the forcing function in the coordinate generating equations as a linear function of the density gradients. The technique displaces a few lines from the neighboring region into the shock region. Numerical calculations for Mach numbers 2 and 4.6 and Reynolds numbers from 320 to 10,000 were performed for a circular cylinder with and without a fairing. Results of Mach number 4.6 and Reynolds number 10,000 for an isothermal wall temperature of 556 K are presented in detail

A CFD validation roadmap for hypersonic flows

A roadmap for computational fluid dynamics (CFD) code validation is developed. The elements of the roadmap are consistent with air-breathing vehicle design requirements and related to the important flow path components: forebody, inlet, combustor, and nozzle. Building block and benchmark validation experiments are identified along with their test conditions and measurements. Based on an evaluation criteria, recommendations for an initial CFD validation data base are given and gaps identified where future experiments would provide the needed validation data

Status of research at the Institute for Computer Applications in Science and Engineering (ICASE)

Research conducted at the Institute for Computer Applications in Science and Engineering in applied mathematics, numerical analysis and computer science is summarized