SEARCH · Search NASA

Results for “corner singularities”

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.

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS1985

A critical assessment of viscous models of trench topography and corner flow

Stresses for Newtonian viscous flow in a simple geometry (e.g., corner flow, bending flow) are obtained in order to study the effect of imposed velocity boundary conditions. Stress for a delta function velocity boundary condition decays as 1/R(2); for a step function velocity, stress goes as 1/R; for a discontinuity in curvature, the stress singularity is logarithmic. For corner flow, which has a discontinuity of velocity at a certain point, the corresponding stress has a 1/R singularity. However, for a more realistic circular-slab model, the stress singularity becomes logarithmic. Thus the stress distribution is very sensitive to the boundary conditions, and in evaluating the applicability of viscous models of trench topography it is essential to use realistic geometries. Topography and seismicity data from northern Hoshu, Japan, were used to construct a finite element model, with flow assumed tangent to the top of the grid, for both Newtonian and non-Newtonian flow (power law 3 rheology). Normal stresses at the top of the grid are compared to the observed trench topography and gravity anomalies. There is poor agreement. Purely viscous models of subducting slables with specified velocity boundary conditions do not predict normal stress patterns compatible with observed topography and gravity. Elasticity and plasticity appear to be important for the subduction process.

Zhang, J.↗

NASA NTRS1995

Thermal barrier coating life modeling in aircraft gas turbine engines

Analytical models for predicting ceramic thermal barrier coating (TBC) spalling life in aircraft gas turbine engines are presented. Electron beam-physical vapor deposited (EB-PVD) and plasma sprayed TBC systems are discussed. An overview of the following TBC spalling mechanisms is presented: metal oxidation at the ceramic-metal interface, ceramic-metal interface stress singularities at edges and corners, ceramic-metal interface stresses caused by radius of curvature and interface roughness, material properties and mechanical behavior, temperature gradients, component design features and object impact damage. TBC spalling life analytical models are proposed based on observations of TBC spalling and plausible failure theories. TBC spalling was assumed to occur when the imposed stresses exceed the material strength (at or near the ceramic-metal interface). TBC failure knowledge gaps caused by lack of experimental evidence and analytical understanding are noted. The analytical models are considered initial engineering approaches that capture observed TBC failure trends.

Nissley, David M.↗

NASA NTRS1981

Solutions of the Helmholtz equation with boundary conditions for force-free magnetic fields

It is shown that the solution, with one ignorable coordinate, for the Taylor minimum energy state (resulting in a force-free magnetic field) in either a straight cylindrical or a toroidal geometry with arbitrary cross section can be reduced to the solution of either an inhomogeneous Helmholtz equation or a Grad-Shafranov equation with simple boundary conditions. Standard Green's function theory is, therefore, applicable. Detailed solutions are presented for the Taylor state in toroidal and cylindrical domains having a rectangular cross section. The focus is on solutions corresponding to the continuous eigenvalue spectra. Singular behavior at 90 deg corners is explored in detail.

Rasband, S. N.↗

NASA NTRS1981

A development of grid generation procedure for multicomponent aerodynamic configuration

Two approaches for solving the transonic flow in a multi-block grid were explored. The first approach examines a method involving "zonal decomposition" wherein block boundaries are treated as true boundary surfaces separating interfacing grids. The issues investigated involve techniques for matching solutions at a block boundary. A feasibility study was completed and the results are presented. The second approach involves overlapping grids for differencing across a block boundary near an artificially induced coordinate singularity occurring at a fictitious corner. This approach selects a set of neighboring nodes for the fictitious corner such that the resulting physical cells for a node are topologically the same as any other node on the airfoil surface.

Chen, H. C.↗

NASA NTRS2002

A Unified Overset Grid Generation Graphical Interface and New Concepts on Automatic Gridding Around Surface Discontinuities

For many years, generation of overset grids for complex configurations has required the use of a number of different independently developed software utilities. Results created by each step were then visualized using a separate visualization tool before moving on to the next. A new software tool called OVERGRID was developed which allows the user to perform all the grid generation steps and visualization under one environment. OVERGRID provides grid diagnostic functions such as surface tangent and normal checks as well as grid manipulation functions such as extraction, extrapolation, concatenation, redistribution, smoothing, and projection. Moreover, it also contains hyperbolic surface and volume grid generation modules that are specifically suited for overset grid generation. It is the first time that such a unified interface existed for the creation of overset grids for complex geometries. New concepts on automatic overset surface grid generation around surface discontinuities will also be briefly presented. Special control curves on the surface such as intersection curves, sharp edges, open boundaries, are called seam curves. The seam curves are first automatically extracted from a multiple panel network description of the surface. Points where three or more seam curves meet are automatically identified and are called seam corners. Seam corner surface grids are automatically generated using a singular axis topology. Hyperbolic surface grids are then grown from the seam curves that are automatically trimmed away from the seam corners.

Chan, William M.↗

NASA NTRS1986

On the inception lengths of swept shock-wave/turbulent boundary-layer interactions

Experimental data are shown for the inception lengths of swept compression corner-generated and fin-generated shock/boundary layer interactions at Mach 2.95. These results are found to correlate on the basis of three different flow regimes. The inception lengths of these flows are dominated by a singularity at the cylindrical/conical boundary for swept corners and by an elongation due to shock wave sweepback for fin interactions. Similarity rules for both Re(delta)and shock generator geometry effects on inception lengths are demonstrated.

Settles, G. S.↗

NASA NTRS1979

A careful numerical study of flowfields about external conical corners. I - Symmetric configurations

A numerical study of the flowfield about symmetrical external axial corners formed by the juncture of swept compressive wedges is presented. The geometrical configuration under investigation allows a unified treatment of external corners typical of delta wings and of rectangular inlets. Comparisons are made with other numerical results. For the delta wing configuration, the occurrence of an anomalous shock behavior predicted by Gonor is discussed. A parametric study shows the singular behavior of the pressure at the corner as a function of the external corner angle, and a second parametric study shows the effect of finite corner radius on the location of the crossflow stagnation point. Previous theoretical predictions of certain flow features typical of corner flows agree well with present numerical results.

Salas, M. D.↗

NASA NTRS1979

A three-dimensional viscous/potential flow interaction analysis method for multi-element wings: Modifications to the potential flow code to allow part-span, high-lift devices and close-interference calculations

The description of the modified code includes details of a doublet subpanel technique in which panels that are close to a velocity calculation point are replaced by a subpanel set. This treatment gives the effect of a higher panel density without increasing the number of unknowns. In particular, the technique removes the close approach problem of the earlier singularity model in which distortions occur in the detailed pressure calculation near panel corners. Removal of this problem allowed a complete wake relaxation and roll-up iterative procedure to be installed in the code. The geometry package developed for the new technique and also for the more general configurations is based on a multiple patch scheme. Each patch has a regular array of panels, but arbitrary relationships are allowed between neighboring panels at the edges of adjacent patches. This provides great versatility for treating general configurations.

Maskew, B.↗

NASA NTRS1985

A clamped rectangular plate containing a crack

The general problem of a rectangular plate clamped along two parallel sides and containing a crack parallel to the clamps is considered. The problem is formulated in terms of a system of singular integral equations and the asymptotic behavior of the stress state near the corners is investigated. Numerical examples are considered for a clamped plate without a crack and with a centrally located crack, and the stress intensity factors and the stresses along the clamps are calculated.

Tang, R.↗

NASA NTRS1985

A clamped rectangular plate containing a crack

Tang, R.↗

NASA NTRS1977

Numerical determination of the fundamental eigenvalue for the Laplace operator on a spherical domain

Methods for obtaining approximate solutions for the fundamental eigenvalue of the Laplace-Beltrami operator (i.e., the membrane eignevalue problem for the vibration equation) on the unit spherical surface are developed. Two types of spherical surface domains are considered: the interior of a spherical triangle, and the exterior of a great circle arc extending for less than pi radians (a spherical surface with a slit). In both cases, zero boundary conditions are imposed. In order to solve the resulting second-order elliptic partial differential equations in two independent variables, a finite difference approximation is employed. The fundamental eigenvalue is approximated by iteration utilizing the power method and point successive overrelaxation. Some numerical results are given and compared, in certain special cases, with analytical solutions to the eigenvalue problem. The significance of the numerical eigenvalue results is discussed in terms of the singularities in the solution of three-dimensional boundary-value problems near a polyhedral corner of the domain.

Walden, H.↗

NASA NTRS1987

A new formulation of electromagnetic wave scattering using an on-surface radiation boundary condition approach

A new formulation of electromagnetic wave scattering by convex, two-dimensional conducting bodies is reported. This formulation, called the on-surface radiation condition (OSRC) approach, is based upon an expansion of the radiation condition applied directly on the surface of a scatterer. It is now shown that application of a suitable radiation condition directly on the surface of a convex conducting scatterer can lead to substantial simplification of the frequency-domain integral equation for the scattered field, which is reduced to just a line integral. For the transverse magnetic case, the integrand is known explicitly. For the transverse electric case, the integrand can be easily constructed by solving an ordinary differential equation around the scatterer surface contour. Examples are provided which show that OSRC yields computed near and far fields which approach the exact results for canonical shapes such as the circular cylinder, square cylinder, and strip. Electrical sizes for the examples are ka = 5 and ka = 10. The new OSRC formulation of scattering may present a useful alternative to present integral equation and uniform high-frequency approaches for convex cylinders larger than ka = 1. Structures with edges or corners can also be analyzed, although more work is needed to incorporate the physics of singular currents at these discontinuities. Convex dielectric structures can also be treated using OSRC.

Kriegsmann, Gregory A.↗

NASA NTRS1982

Altitude transitions in energy climbs

The aircraft energy-climb trajectory for configurations with a sharp transonic drag rise is well known to possess two branches in the altitude/Mach-number plane. Transition in altitude between the two branches occurs instantaneously, a 'corner' in the minimum-time solution obtained with the energy-state model. If the initial and final values of altitude do not lie on the energy-climb trajectory, then additional jumps (crude approximations to dives and zooms) are required at the initial and terminal points. With a singular-perturbation approach, a 'boundary-layer' correction is obtained for each altitude jump, the transonic jump being a so-called 'internal' boundary layer, different in character from the initial and terminal layers. The determination of this internal boundary layer is examined and some computational results for an example presented.

Weston, A. R.↗

NASA NTRS1994

Equilibrium fluid interface behavior under low- and zero-gravity conditions

We describe here some of our recent mathematical work, which forms a basis for the Interface Configuration Experiment scheduled for USML-2. The work relates to the design of apparatus that exploits microgravity conditions for accurate determination of contact angle. The underlying motivation for the procedures rests on a discontinuous dependence of the capillary free surface interface S on the contact angle gamma, in a cylindrical capillary tube whose section (base) omega contains a protruding corner with opening angle 2 alpha. Specifically, in a gravity-free environment, omega can be chosen so that, for all sufficiently large fluid volume, the height of S is uniquely determined as a (single-valued) function mu(x,y) entirely covering the base; the height mu is bounded over omega uniformly in gamma throughout the range absolute value of (gamma -(pion/2)) less than or equal to alpha, while for absolute value of (gamma - (pion/2)) greater than alpha fluid will necessarily move to the corner and uncover the base, rising to infinity (or falling to negative infinity) at the vertex, regardless of volume. We mention here only that procedures based on the phenomenon promise excellent accuracy when gamma is close pion/2 but may be subject to experimental error when gamma is close to zero (orpion), as the 'singular' part of the domain over which the fluid accumulates (or disappears) when a critical angle gamma theta is crossed then becomes very small and may be difficult to observe. We ignore the trivial case gamma is equal to pion/2 (planar free surface), to simplify the discussion.

Concus, Paul↗

NASA NTRS1977

Mathematical and computational studies of equilibrium capillary free surfaces

The results of several independent studies are presented. The general question is considered of whether a wetting liquid always rises higher in a small capillary tube than in a larger one, when both are dipped vertically into an infinite reservoir. An analytical investigation is initiated to determine the qualitative behavior of the family of solutions of the equilibrium capillary free-surface equation that correspond to rotationally symmetric pendent liquid drops and the relationship of these solutions to the singular solution, which corresponds to an infinite spike of liquid extending downward to infinity. The block successive overrelaxation-Newton method and the generalized conjugate gradient method are investigated for solving the capillary equation on a uniform square mesh in a square domain, including the case for which the solution is unbounded at the corners. Capillary surfaces are calculated on the ellipse, on a circle with reentrant notches, and on other irregularly shaped domains using JASON, a general purpose program for solving nonlinear elliptic equations on a nonuniform quadrilaterial mesh. Analytical estimates for the nonexistence of solutions of the equilibrium capillary free-surface equation on the ellipse in zero gravity are evaluated.

Albright, N.↗

NASA NTRS1976

Computation of the inviscid supersonic flow over an external axial corner

A second-order finite-difference procedure is used to evaluate the inviscid supersonic flowfield surrounding an external axial corner composed of swept planar compression surfaces and representing the inlets on existing high-speed aircraft. The governing partial differential equations in conservation-law form are hyperbolic with respect to the axial coordinate and are solved iteratively by means of MacCormack's algorithm. The procedure treats both the peripheral shock wave and vortical singularities as discontinuities. Numerical results are presented for two parametric studies regarding the effects on the flowfield of varying the free-stream Mach number and the leading edge sweep of the horizontal wedge. Results of parametric Mach number study agree with the Mach number independence principle in that as the Mach number increases, such characteristics as shock shape, cross-flow sonic line location, and vortical singularity position approach an asymptote.

Kutler, P.↗

NASA NTRS1976

Diffraction of a shock wave by a compression corner. I - Regular reflection

The unsteady, two-dimensional flowfield resulting from the interaction of a moving planar shock wave with a compression corner is determined using a second-order, discontinuity-fitting, finite-difference approach. The time-dependent Euler equations are transformed to normalize the distance between the body and peripheral shock and to include the existing self-similar property of the flow. The resulting set of partial differential equations in conservation-law form is then solved in a time-dependent fashion using MacCormack's scheme. The vortical singularity, which lies on the body surface, and the single reflected shock are both treated as discontinuities in the numerical procedure. The results of the numerical simulation compare quite favorably with existing experimental interferograms and yield better flowfield resolution than previous first-order, shock-capturing, numerical solutions.

Kutler, P.↗

NASA NTRS2003

Gibbs-Thomson Law for Singular Step Segments: Thermodynamics Versus Kinetics

Classical Burton-Cabrera-Frank theory presumes that thermal fluctuations are so fast that at any time density of kinks on a step is comparable with the reciprocal intermolecular distance, so that the step rate is about isotropic within the crystal plane. Such azimuthal isotropy is, however, often not the case: Kink density may be much lower. In particular, it was recently found on the (010) face of orthorhombic lysozyme that interkink distance may exceed 500-600 intermolecular distances. Under such conditions, Gibbs-Thomson law (GTL) may not be applicable: On a straight step segment between two corners, communication between the comers occurs exclusively by kink exchange. Annihilation between kinks of opposite sign generated at the comers results in the grain in step energy entering GTL. If the step segment length l much greater than D/v, where D and v are the kink diffusivity and propagation rate, respectively, the opposite kinks have practically no chance to annihilate and GTL is not applicable. The opposite condition of the GTL applicability, l much less than D/v, is equivalent to the requirement that relative supersaturation Delta(sub mu)/kT much less than alpha/l, where alpha is molecular size. Thus, GTL may be applied to a segment of 10(exp 3)alpha approx. 3 x 10(exp -5)cm approx 0.3 micron only if supersaturation is less than 0.1%, while practically used driving forces for crystallization are much larger. Relationships alternative to the GTL for different, but low, kink density have been discussed. They confirm experimental evidences that the Burton-Cabrera-Frank theory of spiral growth is growth rates twice as low as compared to the observed figures. Also, application of GTL results in unrealistic step energy while suggested kinetic law give reasonable figures.

Chernov, A. A.↗