Search NASA⌕ Search

SEARCH · Search NASA

Results for “weak solutions”

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 451 records · Page 25

Progress on a Taylor weak statement finite element algorithm for high-speed aerodynamic flows

A new finite element numerical Computational Fluid Dynamics (CFD) algorithm has matured to the point of efficiently solving two-dimensional high speed real-gas compressible flow problems in generalized coordinates on modern vector computer systems. The algorithm employs a Taylor Weak Statement classical Galerkin formulation, a variably implicit Newton iteration, and a tensor matrix product factorization of the linear algebra Jacobian under a generalized coordinate transformation. Allowing for a general two-dimensional conservation law system, the algorithm has been exercised on the Euler and laminar forms of the Navier-Stokes equations. Real-gas fluid properties are admitted, and numerical results verify solution accuracy, efficiency, and stability over a range of test problem parameters.

Baker, A. J.↗

The Nonlinear Evolution of the Inviscid Secondary Instability of Streamwise Vortex Structures

The weakly nonlinear evolution of an inviscid marginally unstable wave growing on a boundary layer supporting a streamwise vortex structure is investigated. The nonlinear growth of the wave is found to be controlled by the diffusion layer located at the edge of the critical layer associated with the wave. The evolution equation is found to depend on the upstream history of the wave and the solution of the equation suggests that the wave either restructures the mean state so as to make it stable or develops a singularity at a finite distance downstream of the point of neutral stability.

Blackaby, Nicholas D.↗

Spectroscopy of Pluto, 380-930 Nm at Six Longitudes

We have obtained spectra of the Pluto-Charon pair (unresolved) in the wavelength range 380-930 nm with resolution approx..450 at six roughly equally spaced longitudes. The data were taken in May and June, 2014, with the 4.2-m Isaac Newton Telescope at Roque de Los Muchachos Observatory in the Canary Islands, using the ACAM (auxiliary-port camera) in spectrometer mode, and using two solar analog stars. The new spectra clearly show absorption bands of solid CH4 at 620, 728, and 850-910 nm, which were known from earlier work. The 620-nm CH4 band is intrinsically very weak, and its appearance indicates a long optical path-length through the ice. This is especially true if it arises from CH4 dissolved in N2 ice. Earlier work (Owen et al. Science 261, 745, 1993) on the near-infrared spectrum of Pluto (1-2.5 microns) has shown that the CH4 bands are shifted to shorter wavelengths because the CH4 occurs as a solute in beta-phase crystalline N2. The optical path-length through the N2 crystals must be on the order of several cm to produce the N2 band observed at 2.15 microns. The new spectra exhibit a pronounced red slope across the entire wavelength range; the slope is variable with longitude, and differs in a small but significant way from that measured at comparable longitudes by Grundy & Fink (Icarus 124, 329, 1996) in their 15-year study of Pluto's spectrum (500-1000 nm). The new spectra will provide an independent means for calibrating the color filter bands on the Multispectral Visible Imaging Camera (MVIC) (Reuter et al. Space Sci. Rev. 140, 129, 2008) on the New Horizons spacecraft, which will encounter the Pluto-Charon system in mid-2015. They will also form the basis of modeling the spectrum of Pluto at different longitudes to help establish the nature of the non-ice component(s) of Pluto's surface. It is presumed that the non-ice component is the source of the yellow-red coloration of Pluto, which is known to be variable across the surface.

Spectroscopy↗

Effect of metallurgical structure and properties on adhesion and friction behavior of cobalt alloys

The metallurgical structure and some of the mechanical properties of two cobalt alloys, cobalt-50% iron and cobalt-25% molybdenum-10% chromium, were determined under various heat treated conditions. The mechanical properties of the bcc disordered Co-50Fe alloy, which was found to be very brittle, indicated an exceedingly low fracture strength, low hardness, and very weak grain boundary strength. Ordering by suitable heat treatment only produced a more brittle material with a lower fracture strength and a slightly higher hardness value. Work hardening was found to produce a finer grain structure and a greater grain boundary strength. Tensile properties were examined. It was found that the Co-25Mo-10Cr alloy was difficult to place in the alpha Co solid solution condition, which limited the ability to use precipitation as a hardening reaction. Over two hundred adhesion cycles from zero contact load, to maximum load, to fracture were conducted between couples for each of the above alloys in an ultrahigh vacuum system which would permit the sample surfaces to be cleaned of all contaminant layers. In the Co-50Fe case, the calculated fracture stress from the adhesion tests showed values in the range of 80 to 150 k.s.i., which is about ten times greater than the values from tension tests.

Keller, D. V., Jr.↗

Hybrid approach to steady transonic normal-shock compressible laminar boundary-layer interactions over airfoils

A simulation of steady transonic viscous-inviscid interactions through a combination of zonal solution methods (involving different equation sets and numerical regions for various flow regions) is used in an investigation of a fluid mechanics problem in which the pressure distribution is determined by the interaction between transonic inviscid flow and an inner laminar viscous layer. The results obtained are in good agreement with laminar experimental data. Both attached and separated boundary layer flows are considered for either weak or strong interactions.

Ram, R. B.↗

Ignition and combustion characteristics of metallized propellants

During this reporting period, theoretical work on the secondary atomization process was continued and the experimental apparatus was improved. A one-dimensional model of a rocket combustor, incorporating multiple droplet size classes, slurry combustion, secondary atomization, radiation heat transfer, and two-phase slip between slurry droplets and the gas flow was derived and a computer code was written to implement this model. The STANJAN chemical equilibrium solver was coupled with this code to yield gas temperature, density, and composition as functions of axial location. Preliminary results indicate that the model is performing correctly, given current model assumptions. Radiation heat transfer in the combustion chamber is treated as an optically-thick participating media problem requiring a solution of the radiative transfer equation. A cylindrical P sub 1 approximation was employed to yield an analytical expression for chamber-wall heat flux at each axial location. The code exercised to determine the effects of secondary atomization intensity, defined as the number of secondary drops produced per initial drop, on chamber burnout distance and final Al2O3 agglomerate diameter. These results indicate that only weak secondary atomization is required to significantly reduce these two parameters. Stronger atomization intensities were found to yield decreasing marginal benefits. The experimental apparatus was improved to reduce building vibration effects on the optical system alignment. This was accomplished by mounting the burner and the transmitting/receiving optics on a single frame supported by vibration-isolation legs. Calibration and shakedown tests indicate that vibration problems were eliminated and that the system is performing correctly.

Mueller, D. C.↗

Analysis of the interaction of a weak normal shock wave with a turbulent boundary layer

The method of matched asymptotic expansions is used to analyze the interaction of a normal shock wave with an unseparated turbulent boundary layer on a flat surface at transonic speeds. The theory leads to a three-layer description of the interaction in the double limit of Reynolds number approaching infinity and Mach number approaching unity. The interaction involves an outer, inviscid rotational layer, a constant shear-stress wall layer, and a blending region between them. The pressure distribution is obtained from a numerical solution of the outer-layer equations by a mixed-flow relaxation procedure. An analytic solution for the skin friction is determined from the inner-layer equations. The significance of the mathematical model is discussed with reference to existing experimental data.

Melnik, R. E.↗

Deconvolution of a pre-outburst picture of SN 1987A

The block iterative method (Young 1971) of algebraic image restoration is applied to a photographic plate of the field around SN 1987A obtained 4 yr before outburst by Chu at the CTIO 4-m telescope. By setting appropriate constraints for the solution on a pixel-by-pixel basis, the image of a starfield centered on Sk -69 202, the precursor to the supernova, is restored. It is found that this star has two companions, which others (e.g., Walborn et al., 1987 and West et al., 1987) have also noted and designated as Star 2 and Star 3. Also found is marginal evidence for Star 4, a weak source at theta = 222 deg, r = 1.9 arcsec. Correction for the contributions of these companions yields an estimate of the magnitude for Sk -69 202 alone of V = 12.37.

Heap, Sara R.↗

Transmission of singularities through a shock wave and the sound generation

The interaction of a plane shock wave of finite strength with a vortex line, point vortex, doublet or quadrupole of weak strength is studied. Based upon the physical condition that a free vortex line cannot support a pressure difference, rules are established which define the change of the linear intensity of the segment of the vortex line after its passage through the shock. The rules for point vortex, doublet, and quadrupole are then established as limiting cases. These rules can be useful for the construction of the solution of the entire flow field and for its physical interpretation. However, the solution can be obtained directly by the technique developed for shock diffraction problems. Explicit solutions and the associated sound generation are obtained for the passage of a point vortex through the shock wave.

Ting, L.↗

Monte Carlo simulation of reentry flows with ionization

The Direct Simulation Monte Carlo method is applied to a rarefied, weakly ionized, hypersonic flow over a blunt axisymmetric body. An ionization model based on the concept of ambipolar diffusion is used and a model for the sheath is presented. The effects of the new modeling techniques are investigated for flow over the Project Fire II configuration at 11.37 km/s at an altitude of 84.6 km. The calculated results are presented and compared with both experimental data and solutions where ionization effects were not included. In general, the calculated results overpredict the experimental values by about 15-20 percent.

Taylor, Jeff C.↗

Symmetries of the TDNLS equations for weakly nonlinear dispersive MHD waves

In this paper we consider the symmetries and conservation laws for the TDNLS equations derived by Hada (1993) and Brio, Hunter and Johnson, to describe the propagation of weakly nonlinear dispersive MHD waves in beta approximately 1 plasmas. The equations describe the interaction of the Alfven and magnetoacoustic modes near the triple umbilic, where the fast magnetosonic, slow magnetosonic and Alfven speeds coincide and a(g)(exp 2) = V(A)(exp 2) where a(g) is the gas sound speed and V(A) is the Alfven speed. We discuss Lagrangian and Hamiltonian formulations, and similarity solutions for the equations.

Webb, G. M.↗

2-D/Axisymmetric Formulation of Multi-dimensional Upwind Scheme

A multi-dimensional upwind discretization of the two-dimensional/axisymmetric Navier-Stokes equations is detailed for unstructured meshes. The algorithm is an extension of the fluctuation splitting scheme of Sidilkover. Boundary conditions are implemented weakly so that all nodes are updated using the base scheme, and eigen-value limiting is incorporated to suppress expansion shocks. Test cases for Mach numbers ranging from 0.1-17 are considered, with results compared against an unstructured upwind finite volume scheme. The fluctuation splitting inviscid distribution requires fewer operations than the finite volume routine, and is seen to produce less artificial dissipation, leading to generally improved solution accuracy.

Wood, William A.↗

A multiple-scales model of the shock-cell structure of imperfectly expanded supersonic jets

The present investigation is concerned with the development of an analytical model of the quasi-periodic shock-cell structure of an imperfectly expanded supersonic jet. The investigation represents a part of a program to develop a mathematical theory of broadband shock-associated noise of supersonic jets. Tam and Tanna (1982) have suggested that this type of noise is generated by the weak interaction between the quasi-periodic shock cells and the downstream-propagating large turbulence structures in the mixing layer of the jet. In the model developed in this paper, the effect of turbulence in the mixing layer of the jet is simulated by the addition of turbulent eddy-viscosity terms to the momentum equation. Attention is given to the mean-flow profile and the numerical solution, and a comparison of the numerical results with experimental data.

Tam, C. K. W.↗

Diskoseismology: Probing accretion disks. I - Trapped adiabatic oscillations

The normal modes of acoustic oscillations within thin accretion disks which are terminated by an innermost stable orbit around a slowly rotating black hole or weakly magnetized compact neutron star are analyzed. The dominant relativistic effects which allow modes to be trapped within the inner region of the disk are approximated via a modified Newtonian potential. A general formalism is developed for investigating the adiabatic oscillations of arbitrary unperturbed disk models. The generic behavior is explored by way of an expansion of the Lagrangian displacement about the plane of symmetry and by assuming separable solutions with the same radial wavelength for the horizontal and vertical perturbations. The lowest eigenfrequencies and eigenfunctions of a particular set of radial and quadrupole modes which have minimum motion normal for the plane are obtained. These modes correspond to the standard dispersion relation of disk theory.

Nowak, Michael A.↗

Tollmien-Schlichting/vortex interactions in compressible boundary layer flows

The weakly nonlinear interaction of oblique Tollmien-Schlichting waves and longitudinal vortices in compressible, high Reynolds number, boundary-layer flow over a flat plate is considered for all ranges of the Mach number. The interaction equations comprise of equations for the vortex which is indirectly forced by the waves via a boundary condition, whereas a vortex term appears in the amplitude equation for the wave pressure. The downstream solution properties of interaction equations are found to depend on the sign of an interaction coefficient. Compressibility is found to have a significant effect on the interaction properties; principally through its impact on the waves and their governing mechanism, the triple-deck structure. It is found that, in general, the flow quantities will grow slowly with increasing downstream co-ordinate; i.e. in general, solutions do not terminate in abrupt, finite-distance 'break-ups'.

Blackaby, Nicholas D.↗

Nucleation and Growth According to Lysozyme

How does one take a molecule, strongly asymmetric in both shape and charge distribution, and assemble it into a crystal? We propose a model for the nucleation and crystal growth process for tetragonal lysozyme that may be very germane to other monomeric proteins. The first species formed is postulated to be a dimer. Through repeating associations involving the same intermolecular interactions this becomes the 4(sub 3) helix, that in turn serves as the basic unit for nucleation and crystal growth. High salt attenuates surface charges while promoting hydrophobic interactions. Symmetry facilitates helix self-association. Assembly stability is enhanced when a four helix structure is obtained, with each bound to two neighbors. Only two unique interactions are required. The first are those for helix formation, where the dominant interaction is the intermolecular bridging anion. The second is the anti-parallel side-by-side helix-helix interaction, guided by alternating pairs of symmetry related salt bridges along each side. At this stage all eight unique positions of the P4(sub 3)2(sub 1)2(sub 1) unit cell are filled. From the above, the process is one of a) attenuating the most strongly interacting groups, such that b) the molecules begin to self-associate in defined patterns, so that c) symmetry is obtained, which d) propagates as a growing crystal. Simple and conceptually obvious in hindsight, this tells much about what we are empirically doing when we crystallize macromolecules. By adjusting the solution parameters we are empirically balancing the intermolecular interactions, preferentially attenuating the dominant strong (for lysozyme the charged groups) while strengthening the lesser strong (hydrophobic) interactions. Lysozyme is atypical in the breadth of its crystallization conditions; many proteins only crystallize under narrowly defined conditions, pointing to the criticality of the empirical balancing process. Lack of a singularly defined association pathway leads to formation of multiple species, i.e., amorphous precipitation. Weak interactions, such as hydrogen bonds, are promiscuous, serving to strengthen rather than define specific interactions. Participation in an interaction sequesters that surface from subsequent interactions, and we expect the strongest bonds to form first. When two molecules self associate the resulting species will have an axis of symmetry. Subsequent interactions between two associated species having equivalent interactions will also have symmetry. Only a few unique sets of interactions are required to give any of the commonly found space groups for monomeric proteins. This model and what it suggests will be discussed.

Pusey, Marc L.↗

Fundamental Mixing and Combustion Experiments for Propelled Hypersonic Flight

The first experiment is a study of a coaxial jet discharging into stagnant laboratory air, with center jet of a mixture of 5% oxygen and 95% helium by volume and coflow jet of air. The exit flow pressure of both center-jet and coflow nozzles is 1 atmosphere. The presence of oxygen in the center jet is to allow the use of an oxygen flow-tagging technique (RELIEF4) to obtain non-intrusive velocity measurements. Both jets are nominally Mach 1.8, but, because of the greater speed of sound, the center jet velocity is more than twice that of the coflow. The mixing layer which forms between the center jet and the coflow near the nozzle exit is compressible, with a calculated convective Mach number of approximately 0.7. This geometry has several advantages: The streamwise development of the flow is generally dominated by turbulent stresses (rather than pressure forces), and thus calculations are sensitive to turbulence modeling. It includes features present in supersonic combustors, including a compressible mixing layer near the nozzle exit and a light-gas/air plume downstream. Since it is a free jet, it provides easy access for both optical instrumentation and probes. Since it is axisymmetric, it requires fewer experimental measurements to fully characterize, and calculations can be performed with more modest computer resources. However, weak shock waves formed at the nozzle exit strengthen and turn normal as they approach the axis, complicating the flow. Care is thus taken in the design of the facility to provide as near as possible to 1-D flow at the exit of both center and coflow nozzles, and to minimize the strength of waves generated at the nozzle exit. Results from this experiment are compared to CFD solutions obtained by VULCAN, a previously developed code used in engine analysis. The second experiment is a study of a supersonic combustor consisting of a diverging duct with single downstream-angled wall injector. Thus, the geometry is relatively simple and large regions of subsonic recirculating flow are avoided. The nominal entrance Mach number is 2 and the enthalpy of the test gas (hot air "simulant") is nominally that of Mach 7 flight. It was believed, on the basis of calculations performed that this would produce mixing-limited flow, that is to say, one for which chemical reaction to equilibrium proceeds at a much greater rate than mixing. It later proved that this was not the case. The primary experimental technique employed is coherent anti-Stokes Raman spectroscopy, known by its acronym CARS. The species probed is molecular nitrogen and the quantity measured is temperature. Intrusive probes, such as Pitot, total temperature, hot-wire, etc., are not used due to access difficulty and high heat flux in the combustor, and because they may alter the flow. CARS has several advantages over other optical methods. It is a relatively mature and well-understood technique. Signal levels are relatively high and the signal is in the form of a coherent (laser) beam that can be collected through small windows. Incoherent (non-CARS) interferences are rejected by spatial filtering.

Diskin, G. S.↗

A globally well-posed finite element algorithm for aerodynamics applications

A finite element CFD algorithm is developed for Euler and Navier-Stokes aerodynamic applications. For the linear basis, the resultant approximation is at least second-order-accurate in time and space for synergistic use of three procedures: (1) a Taylor weak statement, which provides for derivation of companion conservation law systems with embedded dispersion-error control mechanisms; (2) a stiffly stable second-order-accurate implicit Rosenbrock-Runge-Kutta temporal algorithm; and (3) a matrix tensor product factorization that permits efficient numerical linear algebra handling of the terminal large-matrix statement. Thorough analyses are presented regarding well-posed boundary conditions for inviscid and viscous flow specifications. Numerical solutions are generated and compared for critical evaluation of quasi-one- and two-dimensional Euler and Navier-Stokes benchmark test problems.

Iannelli, G. S.↗