Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CANONICAL PROBLEM”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Tetrahedral-Mesh Simulations of Shock-Turbulence Interaction

Despite decades of development of unstructured mesh methods, direct numerical simulations (DNS) of turbulent flows are still predominantly performed on structured or unstructured hexahedral meshes with high-order finite-difference methods, weighted essentially nonoscillatory (WENO) schemes, or hybrid schemes formed by their combinations. Tetrahedral meshes offer easy mesh generation and adaptation around complex geometries and the potential of an orientation-free grid that would benefit the isotropic nature of small-scale dissipation, as well as the solution accuracy of intermediate scales. To advance the state of the art of unstructured-mesh simulation capabilities for shock/turbulence interaction, DNS using pure tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method in this research. By its design, the CESE method is constructed based on a non-dissipative scheme and is a genuinely multidimensional numerical framework that is free from the use of an approximate Riemann-solver. The numerical framework also provides the ability to add numerical dissipation (the nondissipative scheme acts as the reference state like that of the reversible state in thermodynamics) when needed (with justification from mathematics/physics). The above-mentioned features along with the CESE method's consistent shock-capturing approach and strong enforcement of flux conservation in spacetime offers a novel method to accurately simulate turbulent flows and their interaction with shocks using tetrahedral meshes. Two canonical problems, namely, isotropic turbulence interaction with a normal shock and a Mach 2.9 turbulent boundary layer flow over a 24deg compression corner are investigated in this study. Computational results show reasonably good agreement with experimental data and results from structured-mesh, high-order simulations available in the literature. Successful validation of these canonical problems demonstrated here paves the way for future high-fidelity supersonic flow simulations involving complex-geometries.

Venkatachari, Balaji Shankar↗
NASA NTRS1989

Radar transponder antenna pattern analysis for the space shuttle

In order to improve tracking capability, radar transponder antennas will soon be mounted on the Shuttle solid rocket boosters (SRB). These four antennas, each being identical cavity-backed helices operating at 5.765 GHz, will be mounted near the top of the SRB's, adjacent to the intertank portion of the external tank. The purpose is to calculate the roll-plane pattern (the plane perpendicular to the SRB axes and containing the antennas) in the presence of this complex electromagnetic environment. The large electrical size of this problem mandates an optical (asymptotic) approach. Development of a specific code for this application is beyond the scope of a summer fellowship; thus a general purpose code, the Numerical Electromagnetics Code - Basic Scattering Code, was chosen as the computational tool. This code is based on the modern Geometrical Theory of Diffraction, and allows computation of scattering of bodies composed of canonical problems such as plates and elliptic cylinders. Apertures mounted on a curved surface (the SRB) cannot be accomplished by the code, so an antenna model consisting of wires excited by a method of moments current input was devised that approximated the actual performance of the antennas. The improvised antenna model matched well with measurements taken at the MSFC range. The SRB's, the external tank, and the shuttle nose were modeled as circular cylinders, and the code was able to produce what is thought to be a reasonable roll-plane pattern.

Radcliff, Roger↗
NASA NTRS1974

TM surface wave diffraction by a truncated dielectric slab recessed in a perfectly conducting surface

The diffraction of a TM sub o surface wave by a terminated dielectric slab which is flush mounted in a perfectly conducting surface is studied. The incident surface wave gives rise to waves reflected and diffracted by the termination; these reflected and diffracted fields may be expressed in terms of the geometrical theory of diffraction by introducing surface wave reflection and diffraction coefficients which are associated with the termination. In this investigation, the surface wave reflection and diffraction coefficients have been deduced from a formally exact solution to this canonical problem. The solution is obtained by a combination of the generalized scattering matrix technique and function theoretic methods.

Pathak, P. H.↗
NASA NTRS1974

The radiation from slots in truncated dielectric-covered surfaces

A theoretical approach based on the geometrical theory of diffraction is used to study the electromagnetic radiation from a narrow slot in a dielectric-covered perfectly-conducting surface terminated at an edge. The total far-zone field is composed of a geometrical optics field and a diffracted field. The geometrical optics field is the direct radiation from the slot to the field point. The slot also generates surface waves which are incident at the termination of the dielectric cover, where singly-diffracted rays and reflected surface waves are excited. The diffraction and reflection coefficients are obtained from the canonical problem of the diffraction of a surface wave by a right-angle wedge where the dielectric-covered surface is approximated by an impedance surface. This approximation is satisfactory for a very thin cover; however, the radiation from its vertical and faces cannot be neglected in treating the thicker dielectric cover. This is taken into account by using a Kirchhoff-type approximation, which contributes a second term to the diffraction coefficient previously obtained. The contributions from the geometrical optics field, the singly-diffracted rays and all significant multiply-diffracted rays are summed to give the total radiation. Calculated and measured patterns are found to be in good agreement.

Hwang, Y. M.↗
NASA NTRS1974

The dyadic diffraction coefficient for a curved edge

A compact dyadic diffraction coefficient for electromagnetic waves obliquely incident on a curved edge formed by perfectly conducting curved or plane surfaces is obtained. This diffraction coefficent remains valid in the transition regions adjacent to shadow and reflection boundaries, where the diffraction coefficients of Keller's original theory fail. The method is on Keller's method of the canonical problem, which in this case is the perfectly conducting wedge illuminated by plane, cylindrical, conical, and spherical waves. When the proper ray fixed coordinate system is introduced, the dyadic diffraction coefficient for the wedge is found to be the sum of only two dyads, and it is shown that this is also true for the dyadic diffraction coefficients of higher order edges. One dyad contains the acoustic soft diffraction coefficient; the other dyad contains the acoustic hard diffraction coefficient. The expressions for the acoustic wedge diffraction coefficients contain Fresnel integrals, which ensure that the total field is continuous at shadow and reflection boundaries. The diffraction coefficients have the same form for the different types of edge illumination; only the arguments of the Fresnel integrals are different. Since diffraction is a local phenomenon, and locally the curved edge structure is wedge shaped, this result is readily extended to the curved edge.

Kouyoumjian, R. G.↗
NASA NTRS1974

A uniform geometrical theory of diffraction for an edge in a perfectly conducting surface

A compact dyadic diffraction coefficient for electromagnetic waves obliquely incident on a curved edge formed by perfectly conducting curved or plane surfaces is obtained. This diffraction coefficient remains valid in the transition regions adjacent to shadow and reflection boundaries, where the diffraction coefficients of Keller's original theory fail. Our method is based on Keller's method of the canonical problem, which in this case is the perfectly conducting wedge illuminated by plane, cylindrical, conical, and spherical waves. When the proper ray-fixed coordinate system is introduced, the dyadic diffraction coefficient for the wedge is found to be the sum of only two dyads, and it is shown that this is also true for the dyadic diffraction coefficients of higher order edges. One dyad contains the acoustic soft diffraction coefficient; the other dyad contains the acoustic hard diffraction coefficient. The expressions for the acoustic wedge diffraction coefficients contain Fresnel integrals, which ensure that the total field is continuous at shadow and reflection boundaries.

Kouyoumjian, R. G.↗
NASA NTRS1974

Computational procedures for mixed equations with shock waves

This paper discusses the procedures we have developed to treat a canonical problem involving a mixed nonlinear equation with boundary data that imply a discontinuous solution. This equation arises in various physical contexts and is basic to the description of the nonlinear acoustic behavior of a shock wave near a caustic. The numerical scheme developed is of second order, treats discontinuities as such by applying the appropriate jump conditions across them, and eliminates the numerical dissipation and dispersion associated with large gradients. Our results are compared with the results of a first-order scheme and with those of a second-order scheme we have developed. The algorithm used here can easily be generalized to more complicated problems, including transonic flows with imbedded shocks.

Yu, N. J.↗
NASA NTRS1979

Surface wave diffraction by a truncated dielectric slab recessed in a perfectly conducting surface

The diffraction of a TM sub 0 surface wave by a terminated dielectric slab which is flush mounted in a perfectly conducting surface is studied. The incident surface wave gives rise to waves reflected and diffracted by the termination; these reflected and diffracted fields may be expressed in terms of the geometrical theory of diffraction by introducing surface wave reflection and diffraction coefficients which are associated with the termination. In this investigation the surface wave reflection and diffraction coefficients have been deduced from a formally exact solution to this canonical problem. The solution is obtained via a combination of the generalized scattering matrix technique and the Wiener-Hopf procedure. Expressions for the reflection and diffraction coefficients contain integrals which can be evaluated numerically.

Pathak, P. H.↗
NASA NTRS1990

Scattering by an impedance sphere coated with a chiral layer

The scattering of a plane, linearly polarized electromagnetic wave by a sphere on whose surface an impedance boundary condition holds, and that is covered with a concentric layer of chiral material, is considered. Exact, explicit expressions are derived for the scattered field coefficients. The co-polarized and cross-polarized components of the far backscattered field are determined and discussed. The value of this canonical problem as a benchmark for computer codes is pointed out.

Uslenghi, Piergiorgio L. E.↗
NASA NTRS2005

Computation of Sound Propagation by Boundary Element Method

This report documents the development of a Boundary Element Method (BEM) code for the computation of sound propagation in uniform mean flows. The basic formulation and implementation follow the standard BEM methodology; the convective wave equation and the boundary conditions on the surfaces of the bodies in the flow are formulated into an integral equation and the method of collocation is used to discretize this equation into a matrix equation to be solved numerically. New features discussed here include the formulation of the additional terms due to the effects of the mean flow and the treatment of the numerical singularities in the implementation by the method of collocation. The effects of mean flows introduce terms in the integral equation that contain the gradients of the unknown, which is undesirable if the gradients are treated as additional unknowns, greatly increasing the sizes of the matrix equation, or if numerical differentiation is used to approximate the gradients, introducing numerical error in the computation. It is shown that these terms can be reformulated in terms of the unknown itself, making the integral equation very similar to the case without mean flows and simple for numerical implementation. To avoid asymptotic analysis in the treatment of numerical singularities in the method of collocation, as is conventionally done, we perform the surface integrations in the integral equation by using sub-triangles so that the field point never coincide with the evaluation points on the surfaces. This simplifies the formulation and greatly facilitates the implementation. To validate the method and the code, three canonic problems are studied. They are respectively the sound scattering by a sphere, the sound reflection by a plate in uniform mean flows and the sound propagation over a hump of irregular shape in uniform flows. The first two have analytical solutions and the third is solved by the method of Computational Aeroacoustics (CAA), all of which are used to compare the BEM solutions. The comparisons show very good agreements and validate the accuracy of the BEM approach implemented here.

Guo, Yueping↗
NASA NTRS2016

Numerical Simulation of Bolide Entry with Ground Footprint Prediction

As they decelerate through the atmosphere, meteors deposit mass, momentum and energy into the surrounding air at tremendous rates. Trauma from the entry of such bolides produces strong blast waves that can propagate hundreds of kilometers and cause substantial terrestrial damage even when no ground impact occurs. We present a new simulation technique for airburst blast prediction using a fully-conservative, Cartesian mesh, finite-volume solver and investigate the ability of this method to model far- field propagation over hundreds of kilometers. The work develops mathematical models for the deposition of mass, momentum and energy into the atmosphere and presents verification and validation through canonical problems and the comparison of surface overpressures, and blast arrival times with actual results in the literature for known bolides. The discussion also examines the effects of various approximations to the physics of bolide entry that can substantially decrease the computational expense of these simulations. We present parametric studies to quantify the influence of entry-angle, burst-height and other parameters on the ground footprint of the airburst, and these values are related to predictions from analytic and handbook-methods.

Simulation↗
NASA NTRS2016

Canonical Statistical Model for Maximum Expected Immission of Wire Conductor in an Aperture Enclosure

Prediction of the maximum expected electromagnetic pick-up of conductors inside a realistic shielding enclosure is an important canonical problem for system-level EMC design of space craft, launch vehicles, aircraft and automobiles. This paper introduces a simple statistical power balance model for prediction of the maximum expected current in a wire conductor inside an aperture enclosure. It calculates both the statistical mean and variance of the immission from the physical design parameters of the problem. Familiar probability density functions can then be used to predict the maximum expected immission for deign purposes. The statistical power balance model requires minimal EMC design information and solves orders of magnitude faster than existing numerical models, making it ultimately viable for scaled-up, full system-level modeling. Both experimental test results and full wave simulation results are used to validate the foundational model.

Computational electromagnetics↗
NASA NTRS2019

Weighted Least-Squares Cell-Average Gradient Construction Methods for the VULCAN-CFD Second-Order Accurate Unstructured-Grid Cell-Centered Finite-Volume Solver

The ability to solve the equations governing the hypersonic turbulent flow of a real gas on unstructured grids using a spatially-elliptic, 2nd-order accurate, cell-centered, finite-volume method has been recently implemented in the VULCAN-CFD code. The construction of cell-average gradients using a weighted linear least-squares method and the use of these gradients in the construction of the inviscid fluxes is the focus of this paper. A comparison of least-squares stencil construction methodologies is presented and approaches to augment the number of cells participating in the stencil while preserving accuracy are explored. Due to our interest in hypersonic flow, a robust multidimensional cell-average gradient limiter procedure that is consistent with the stencil used to construct the cell-average gradients is described and investigated. Canonical problems are computed to illustrate the challenges and investigate the accuracy, robustness and convergence behavior of the cell-average gradient methods on unstructured cell-centered finite-volume grids. Finally, thermally perfect, chemically frozen, Mach 8 turbulent flow of air around a blunt wedge is computed to demonstrate the robustness and convergence behavior of the new method for constructing stencils of use in a weighted linear least-squares gradient method for a hypersonic flow.

White, Jeffery A.↗
NASA NTRS2019

Efficient and Robust Weighted Least-Squares Cell-Average Gradient Construction Methods for the Simulation of Scramjet Flows

The ability to solve the equations governing the hypersonic turbulent flow of a real gas on unstructured grids using a spatially-elliptic, 2nd-order accurate, cell-centered, finite-volume method has been recently implemented in the VULCAN-CFD code. The construction of cell-average gradients using a weighted linear least-squares method and the use of these gradients in the construction of the inviscid fluxes is the focus of this paper. A comparison of least-squares stencil construction methodologies is presented and approaches designed to minimize the number of cells used to augment/stabilize the least-squares stencil while preserving accuracy are explored. Due to our interest in hypersonic flow, a robust multidimensional cell-average gradient limiter procedure that is consistent with the stencil used to construct the cellaverage gradients is described. Canonical problems are computed to illustrate the challenges and investigate the accuracy, robustness and convergence behavior of the cell-average gradient methods on unstructured cell-centered finite-volume grids. Finally, thermally perfect, chemically frozen, Mach 7.8 turbulent flow of air through a scramjet engine flowpath is computed and compared with experimental data to demonstrate the robustness, accuracy and convergence behavior of the preferred gradient method for a realistic 3-D geometry on a non-hex-dominant grid.

White, Jeffery A.↗
NASA NTRSDate not supplied

The Effect of Reduced Pressure on the Characteristics of Spreading Flames

Flame spread over solid fuels is a canonical problem in fire science, due to its direct implications on material flammability and importance in fire development. In a microgravity environment, such as onboard a spacecraft, flames can behave very differently than on Earth. This is concerning for spaceflight life safety, especially in higher-oxygen environments. Due to the difficulties associated with microgravity testing, low-pressure environments have been proposed as an alternative to approximately replicate the burning behavior of solid fuels observed in reduced gravity conditions because of similar diffusion and heat transfer mechanisms. However, the roles played by gravity and pressure on flame length, standoff distance, and flame spread rate vary with the burning configuration. In concurrent flame spread, the buoyant flow enhances the spread rate by bringing the flame closer to the fuel surface and increasing the heating of the solid fuel. In opposed flame spread, the sample is preheated by the flame ahead of the flame leading edge, which is strongly affected by the surrounding flow field. In this work, we consider flames spreading over thin cotton samples in both downward (opposed) and upward (concurrent) configurations to investigate the effect of pressure (30-100 kPa) on flame characteristics, such as spread rate and standoff distance. A small forced flow is induced upward so that the flames are exposed to a mixed (forced and free) flow. By reducing pressure, flames become less bright, their standoff distance increases, and their spread rates decrease similar to what is observed in low-gravity environments. These results could in help understanding the differences between flames spreading at low pressure and low gravity environments for these similar, yet very different, spreading configurations. They could also provide more information about potential Earth-based flammability testing of materials for spacecraft applications.

fire safety↗
NASA NTRSDate not supplied

The effect of reduced pressure on the characteristics of spreading flames

Flame spread over solid fuels is a canonical problem in fire science, due to its direct implications on material flammability and importance in fire development. In a microgravity environment, such as onboard a spacecraft, flames can behave very differently than on Earth. This is concerning for spaceflight life safety, especially in higher-oxygen environments. Due to the difficulties associated with microgravity testing, low-pressure environments have been proposed as an alternative to approximately replicate the burning behavior of solid fuels observed in reduced gravity conditions because of similar diffusion and heat transfer mechanisms. However, the roles played by gravity and pressure on flame length, standoff distance, and flame spread rate vary with the burning configuration. In concurrent flame spread, the buoyant flow enhances the spread rate by bringing the flame closer to the fuel surface and increasing the heating of the solid fuel. In opposed flame spread, the sample is preheated by the flame ahead of the flame leading edge, which is strongly affected by the surrounding flow field. In this work, we consider flames spreading over thin cotton samples in both downward (opposed) and upward (concurrent) configurations to investigate the effect of pressure (30-100 kPa) on flame characteristics, such as spread rate and standoff distance. A small forced flow is induced upward so that the flames are exposed to a mixed (forced and free) flow. By reducing pressure, flames become less bright, their standoff distance increases, and their spread rates decrease similar to what is observed in low-gravity environments. These results could in help understanding the differences between flames spreading at low pressure and low gravity environments for these similar, yet very different, spreading configurations. They could also provide more information about potential Earth-based flammability testing of materials for spacecraft applications.

fire safety↗