GRID3D-v2: An updated version of the GRID2D/3D computer program for generating grid systems in complex-shaped three-dimensional spatial domains

In order to generate good quality systems for complicated three-dimensional spatial domains, the grid-generation method used must be able to exert rather precise controls over grid-point distributions. Several techniques are presented that enhance control of grid-point distribution for a class of algebraic grid-generation methods known as the two-, four-, and six-boundary methods. These techniques include variable stretching functions from bilinear interpolation, interpolating functions based on tension splines, and normalized K-factors. The techniques developed in this study were incorporated into a new version of GRID3D called GRID3D-v2. The usefulness of GRID3D-v2 was demonstrated by using it to generate a three-dimensional grid system in the coolent passage of a radial turbine blade with serpentine channels and pin fins

A two-dimensional numerical study of the flow inside the combustion chambers of a motored rotary engine

A numerical study was performed to investigate the unsteady, multidimensional flow inside the combustion chambers of an idealized, two-dimensional, rotary engine under motored conditions. The numerical study was based on the time-dependent, two-dimensional, density-weighted, ensemble-averaged conservation equations of mass, species, momentum, and total energy valid for two-component ideal gas mixtures. The ensemble-averaged conservation equations were closed by a K-epsilon model of turbulence. This K-epsilon model of turbulence was modified to account for some of the effects of compressibility, streamline curvature, low-Reynolds number, and preferential stress dissipation. Numerical solutions to the conservation equations were obtained by the highly efficient implicit-factored method of Beam and Warming. The grid system needed to obtain solutions were generated by an algebraic grid generation technique based on transfinite interpolation. Results of the numerical study are presented in graphical form illustrating the flow patterns during intake, compression, gaseous fuel injection, expansion, and exhaust

Numerical simulation of the flow and fuel-air mixing in an axisymmetric piston-cylinder arrangement

The implicit factored method of Beam and Warming was employed to describe the flow and the fuel-air mixing in an axisymmetric piston-cylinder configuration during the intake and compression strokes. The governing equations were established on the basis of laminar flow. The increased mixing due to turbulence was simulated by appropriately chosen effective transport properties. Calculations were performed for single-component gases and for two-component gases and for two-component gas mixtures. The flow field was calculated as functions of time and position for different geometries, piston speeds, intake-charge-to-residual-gas-pressure ratios, and species mass fractions of the intake charge. Results are presented in graphical form which show the formation, growth, and break-up of those vortices which form during the intake stroke and the mixing of fuel and air throughout the intake and compression strokes. It is shown that at bore-to-stroke ratio of less than unity, the vortices may break-up during the intake stroke. It is also shown that vortices which do not break-up during the intake stroke coalesce during the compression stroke. The results generated were compared to existing numerical solutions and to available experimental data

Vortex motion in axisymmetric piston-cylinder configurations

By using the Beam and Warming implicit-factored method of solution of the Navier-Stokes equations, velocities were calculated inside axisymmetric piston cylinder configurations during the intake and compression strokes. Results are presented in graphical form which show the formation, growth and breakup of those vortices which form during the intake stroke by the jet issuing from the valve. It is shown that at bore-to-stroke ratio of less than unity, the vortices may breakup during the intake stroke. It is also shown that vortices which do not breakup during the intake stroke coalesce during the compression stroke

Numerical Study of Boundary-Layer in Aerodynamics

The accomplishments made in the following three tasks are described: (1) The first task was to study shock-wave boundary-layer interactions with bleed - this study is relevant to boundary-layer control in external and mixed-compression inlets of supersonic aircraft; (2) The second task was to test RAAKE, a code developed for computing turbulence quantities; and (3) The third task was to compute flow around the Ames ER-2 aircraft that has been retrofitted with containers over its wings and fuselage. The appendices include two reports submitted to AIAA for publication

Intense Field Multiphoton Ionization via Complex Dressed States: Application to the H Atom

This is the publisher's version, also available electronically from http://journals.aps.org/prl/abstract/10.1103/PhysRevLett.39.1195.Extension of Floquet theory to include continuum as well as bound atomic states yields a practical technique for computation of multiphoton ionization rates in the region where rms field strengths approach the strength of the internal atomic fields

GRID2D/3D: A computer program for generating grid systems in complex-shaped two- and three-dimensional spatial domains. Part 2: User's manual and program listing

An efficient computer program, called GRID2D/3D, was developed to generate single and composite grid systems within geometrically complex two- and three-dimensional (2- and 3-D) spatial domains that can deform with time. GRID2D/3D generates single grid systems by using algebraic grid generation methods based on transfinite interpolation in which the distribution of grid points within the spatial domain is controlled by stretching functions. All single grid systems generated by GRID2D/3D can have grid lines that are continuous and differentiable everywhere up to the second-order. Also, grid lines can intersect boundaries of the spatial domain orthogonally. GRID2D/3D generates composite grid systems by patching together two or more single grid systems. The patching can be discontinuous or continuous. For continuous composite grid systems, the grid lines are continuous and differentiable everywhere up to the second-order except at interfaces where different single grid systems meet. At interfaces where different single grid systems meet, the grid lines are only differentiable up to the first-order. For 2-D spatial domains, the boundary curves are described by using either cubic or tension spline interpolation. For 3-D spatial domains, the boundary surfaces are described by using either linear Coon's interpolation, bi-hyperbolic spline interpolation, or a new technique referred to as 3-D bi-directional Hermite interpolation. Since grid systems generated by algebraic methods can have grid lines that overlap one another, GRID2D/3D contains a graphics package for evaluating the grid systems generated. With the graphics package, the user can generate grid systems in an interactive manner with the grid generation part of GRID2D/3D. GRID2D/3D is written in FORTRAN 77 and can be run on any IBM PC, XT, or AT compatible computer. In order to use GRID2D/3D on workstations or mainframe computers, some minor modifications must be made in the graphics part of the program; no modifications are needed in the grid generation part of the program. The theory and method used in GRID2D/3D is described

GRID2D/3D: A computer program for generating grid systems in complex-shaped two- and three-dimensional spatial domains. Part 1: Theory and method

An efficient computer program, called GRID2D/3D was developed to generate single and composite grid systems within geometrically complex two- and three-dimensional (2- and 3-D) spatial domains that can deform with time. GRID2D/3D generates single grid systems by using algebraic grid generation methods based on transfinite interpolation in which the distribution of grid points within the spatial domain is controlled by stretching functions. All single grid systems generated by GRID2D/3D can have grid lines that are continuous and differentiable everywhere up to the second-order. Also, grid lines can intersect boundaries of the spatial domain orthogonally. GRID2D/3D generates composite grid systems by patching together two or more single grid systems. The patching can be discontinuous or continuous. For continuous composite grid systems, the grid lines are continuous and differentiable everywhere up to the second-order except at interfaces where different single grid systems meet. At interfaces where different single grid systems meet, the grid lines are only differentiable up to the first-order. For 2-D spatial domains, the boundary curves are described by using either cubic or tension spline interpolation. For 3-D spatial domains, the boundary surfaces are described by using either linear Coon's interpolation, bi-hyperbolic spline interpolation, or a new technique referred to as 3-D bi-directional Hermite interpolation. Since grid systems generated by algebraic methods can have grid lines that overlap one another, GRID2D/3D contains a graphics package for evaluating the grid systems generated. With the graphics package, the user can generate grid systems in an interactive manner with the grid generation part of GRID2D/3D. GRID2D/3D is written in FORTRAN 77 and can be run on any IBM PC, XT, or AT compatible computer. In order to use GRID2D/3D on workstations or mainframe computers, some minor modifications must be made in the graphics part of the program; no modifications are needed in the grid generation part of the program. This technical memorandum describes the theory and method used in GRID2D/3D

Stark ionization in dc and ac fields: An L2 complex-coordinate approach

This is the published version, also available here: http://dx.doi.org/10.1103/PhysRevA.27.2946.A finite-dimensional-matrix technique valid for computation of complex eigenvalues and eigenfunctions useful for discussing time evolution in both dc and ac Stark fields is presented. The complex eigenvalue parameters are those of appropriately analytically continued, time-independent Stark Hamiltonians as obtained via the complex scale transformation r→reiθ. Such a transformation distorts the continuous spectrum away from the real axis, exposing the Stark resonances, and also allowing use of finite variational expansions employing L2 basis functions chosen from a complete discrete basis. The structure of the dc and ac Stark Hamiltonians is discussed and extensive convergence studies performed in both the dc and ac cases to fully document the utility of the method. Sudden and adiabatic dc Stark time evolution is used to illustrate the power of finite-dimensional-matrix methods in describing complex, multiple-time-scale time evolution. The relationship between the ac Stark Hamiltonian used (a time-independent truncated Floquet Hamiltonian) and continued-fraction perturbation theory follows easily via use of matrix partitioning, and provides a particularly straightforward derivation of these results. Finally, some illustrative calculations of off-resonant generalized cross sections are given at low and high intensities, indicating that the method works satisfactorily at intensities the order of internal atomic field strengths. A more detailed discussion of time evolution in two-, three-, and four-photon ionization processes appears in the following paper by Holt, Raymer, and Reinhardt

Vortex motion in axisymmetric piston-cylinder configurations

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/76950/1/AIAA-8430-329.pd