Search NASA⌕ Search

SEARCH · Search NASA

Results for “boundary value 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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Linear Stability Analysis of Gravitational Effects on a Low-Density Gas Jet Injected into a High-Density Medium

The objective of this study was to determine the effects of buoyancy on the absolute instability of low-density gas jets injected into high-density gas mediums. Most of the existing analyses of low-density gas jets injected into a high-density ambient have been carried out neglecting effects of gravity. In order to investigate the influence of gravity on the near-injector development of the flow, a spatio-temporal stability analysis of a low-density round jet injected into a high-density ambient gas was performed. The flow was assumed to be isothermal and locally parallel; viscous and diffusive effects were ignored. The variables were represented as the sum of the mean value and a normal-mode small disturbance. An ordinary differential equation governing the amplitude of the pressure disturbance was derived. The velocity and density profiles in the shear layer, and the Froude number (signifying the effects of gravity) were the three important parameters in this equation. Together with the boundary conditions, an eigenvalue problem was formulated. Assuming that the velocity and density profiles in the shear layer to be represented by hyperbolic tangent functions, the eigenvalue problem was solved for various values of Froude number. The Briggs-Bers criterion was combined with the spatio-temporal stability analysis to determine the nature of the absolute instability of the jet whether absolutely or convectively unstable. The roles of the density ratio, Froude number, Schmidt number, and the lateral shift between the density and velocity profiles on the absolute instability of the jet were determined. Comparisons of the results with previous experimental studies show good agreement when the effects of these variables are combined together. Thus, the combination of these variables determines how absolutely unstable the jet will be.

Lawson, Anthony L.↗
NASA NTRS2001

Instability Analysis of a Low-Density Gas Jet Injected into a High-Density Gas

The objective of this study was to determine the effects of buoyancy on the absolute instability of low-density gas jets injected into high-density gas mediums. Most of the existing analyses of low-density gas jets injected into a high-density ambient have been carried out neglecting effects of gravity. In order to investigate the influence of gravity on the near-injector development of the flow, a linear temporal stability analysis and a spatio-temporal stability analysis of a low-density round jet injected into a high-density ambient gas were performed. The flow was assumed to be isothermal and locally parallel; viscous and diffusive effects were ignored. The variables were represented as the sum of the mean value and a normal-mode small disturbance. An ordinary differential equation governing the amplitude of the pressure disturbance was derived. The velocity and density profiles in the shear layer, and the Froude number (signifying the effects of gravity) were the three important parameters in this equation. Together with the boundary conditions, an eigenvalue problem was formulated. Assuming that the velocity and density profiles in the shear layer to be represented by hyperbolic tangent functions, the eigenvalue problem was solved for various values of Froude number. The temporal growth rates and the phase velocity of the disturbances were obtained. It was found that the presence of variable density within the shear layer resulted in an increase in the temporal amplification rate of the disturbances and an increase in the range of unstable frequencies, accompanied by a reduction in the phase velocities of the disturbances. Also, the temporal growth rates of the disturbances were increased as the Froude number was reduced (i.e. gravitational effects increased), indicating the destabilizing role played by gravity. The spatio-temporal stability analysis was performed to determine the nature of the absolute instability of the jet. The roles of the density ratio, Froude number, Schmidt number, and the lateral shift between the density and velocity profiles on the jet s absolute instability were determined. Comparisons of the results with previous experimental studies show good agreement when the effects of these variables are combined together. Thus, the combination of these variables determines how absolutely unstable the jet will be. Experiments were carried out to observe the qualitative differences between a round low-density gas jet injected into a high-density gas (helium jet injected into air) and a round constant density jet (air jet injected into air). Flow visualizations and velocity measurements in the near-injector region of the helium jet show more mixing and spreading of the helium jet than the air jet. The vortex structures develop and contribute to the jet spreading causing the helium jet to oscillate.

Lawson, Anthony Layiwola↗
NASA NTRS2002

Linear Temporal Stability Analysis of a Low-Density Round Gas Jet Injected into a High-Density Gas

It has been observed in previous experimental studies that round helium jets injected into air display a repetitive structure for a long distance, somewhat similar to the buoyancy-induced flickering observed in diffusion flames. In order to investigate the influence of gravity on the near-injector development of the flow, a linear temporal stability analysis of a round helium jet injected into air was performed. The flow was assumed to be isothermal and locally parallel; viscous and diffusive effects were ignored. The variables were represented as the sum of the mean value and a normal-mode small disturbance. An ordinary differential equation governing the amplitude of the pressure disturbance was derived. The velocity and density profiles in the shear layer, and the Froude number (signifying the effects of gravity) were the three important parameters in this equation. Together with the boundary conditions, an eigenvalue problem was formulated. Assuming that the velocity and density profiles in the shear layer to be represented by hyperbolic tangent functions, the eigenvalue problem was solved for various values of Froude number. The temporal growth rates and the phase velocity of the disturbances were obtained. The temporal growth rates of the disturbances increased as the Froude number was reduced (i.e. gravitational effects increased), indicating the destabilizing role played by gravity.

Lawson, Anthony L.↗
NASA NTRS1984

Multigrid techniques for the numerical solution of the diffusion equation

An accurate numerical solution of diffusion problems containing large local gradients can be obtained with a significant reduction in computational time by using a multigrid computational scheme. The spatial domain is covered with sets of uniform square grids of different sizes. The finer grid patterns overlap the coarse grid patterns. The finite-difference expressions for each grid pattern are solved independently by iterative techniques. Two interpolation methods were used to establish the values of the potential function on the fine grid boundaries with information obtained from the coarse grid solution. The accuracy and computational requirements for solving a test problem by a simple multigrid and a multilevel-multigrid method were compared. The multilevel-multigrid method combined with a Taylor series interpolation scheme was found to be best.

Phillips, R. E.↗
NASA NTRS1990

Synoptic/planetary-scale interactions and blocking over the North Atlantic Ocean

Work was completed on the height tendency diagnoses of two extratropical cyclones that occurred upstream from the blocking event studied previously. One developed explosively over water 60 to 36 hours before the block first appeared, while the second developed explosively over the southeastern United States during the time of block formation. In both cases, both vorticity and temperature advection were consistently important forcing mechanisms. This is in contrast to the block itself, in which vorticity advection was easily the dominant forcing mechanism. Latent heat release was also significant, accounting for about 50 percent of the total height falls in the cyclone below 850 mb. Estimates of latent heat release were greatly enhanced by coupling parameterized estimates with values derived from GOES IR data using an algorithm developed by Marshall's F. R. Robertson. Among the difficulties encountered in this work was the identification of an appropriate lower boundary condition for the solution of the height tendency equation. The zero value currently used tends to yield underestimates of the lower troposphere height tendencies. To address this problem a new diagnostic technique was developed in cooperation with Dr. Peter Zwack of the University of Quebec at Montreal. Based on an equation Dr. Zwack had previously developed (the Zwack-Okossi development equation), researchers now have a relationship that is completely consistent with the height tendency equation and provides estimates of lower boundary geostrophic vorticity or height tendencies. Finally, comparison of the SAT (satellite data) and NOSAT (no satellite data) analyses is progressing well. The present focus is on both the new diagnostic technique and the SAT/NOSAT comparisons. The former is being tested on the southeastern United States cyclone case previously mentioned and compared with the height tendency diagnoses already completed. The latter are being examined for the blocking case described in the publications cited in this summary. In addition to obtaining statistics that will allow general comparison of the two analyses, it will be possible to determine whether conclusions about the dynamics of the block development are influenced by the analysis set used.

Smith, Phillip J.↗
NASA NTRS1977

Nonlinear singularly perturbed optimal control problems with singular arcs

A third order, nonlinear, singularly perturbed optimal control problem is considered under assumptions which assure that the full problem is singular and the reduced problem is nonsingular. The separation between the singular arc of the full problem and the optimal control law of the reduced one, both of which are hypersurfaces in state space, is of the same order as the small parameter of the problem. Boundary layer solutions are constructed which are stable and reach the outer solution in a finite time. A uniformly valid composite solution is then formed from the reduced and boundary layer solutions. The value of the approximate solution is that it is relatively easy to obtain and does not involve singular arcs. To illustrate the utility of the results, the technique is used to obtain an approximate solution of a simplified version of the aircraft minimum time-to-climb problem. A numerical example is included.

Ardema, M. D.↗
NASA NTRS1979

Nonlinear singularly perturbed optimal control problems with singular arcs

Singular perturbation techniques are studied for dealing with singular arc problems by analyzing a relatively low-order but otherwise general system. This system encompasses many flight mechanic problems including Goddard's problem and a version of the minimum time-to-climb problem. Boundary layer solutions are constructed which are stable and reach the outer solution in a finite time. A uniformly valid composite solution is then formed from the reduced and boundary layer solutions. The value of the approximate solution is that it is relatively easy to obtain and does not involve singular arcs. To illustrate the utility of the results, the technique is used to obtain an approximate solution of a simplified version of the aircraft minimum time-to-climb problem.

Ardema, M. D.↗
NASA NTRS1977

Stress-intensity factors for cracks emanating from the loaded fastener hole

Using a series approach and the Muskhelishvili formulation in the two-dimensional theory of elasticity, stress-intensity factors K are derived for problems in which cracks emanate radially from the boundary of an arbitrarily loaded internal circular hole in an infinite plate. Numerical values are obtained for K(I) and K(II) for radial cracks from a hole containing a loose-fitted pin or rivet that is pulled perpendicular to the crack direction in the plane of the plate. The method is a general one for determining K for a set of symmetrically emanating radial cracks for a variety of concentrated or distributed tractions on the circular hole.

Shivakumar, V.↗
NASA NTRS1947

Boundary-layer and stalling characteristics of two symmetrical NACA low-drag airfoil sections

Two symmetrical airfoils, an NACA 633-018 and an NACA 631-012, were investigated for the purpose of determining their stalling and boundary-layer characteristics with a view toward the eventual application of this information to the problem of boundary-layer control. Force measurements, pressure distributions, tuft studies, and boundary-layer-profile measurements were made at a value of 5,800,000 Reynolds number. It was found that the 18-percent-thick airfoil stalled progressively from the trailing edge because of separation of the turbulent boundary layer. In contrast, the12-percent-thick airfoil stalled abruptly from a separation of flow near the leading edge before the turbulent boundary layer became subject to separation. From this it was concluded that if high values of lift are to be obtained with thin, high-critical-speed sections by means of boundary-layer control, the work must be directed toward delaying the separation of flow near the leading edge. It was found that the presence of a nose flap on the 12-percent-thick section caused the airfoil to stall in a manner similar to that of the 18-percent-thick section.

Mccullough, George B↗
NASA NTRS1989

A finite difference treatment of Stokes-type flows: Preliminary report

The equations Laplacian operator omega = 0, (1.1a) and omega = Laplacian operator Chi, (1.1b) describe, in suitable units, 2-D Stokes flow of an incompressible fluid occupying a domain D in which omega is the vorticity and Chi is the stream function. The flow is uniquely determined by specifying the velocity on the boundary B of D, a condition which leads to specifying the stream function Chi and its normal derivative Chi sub n on B. A mathematically similar problem arises in describing the equilibrium of a flat plate in structural mechanics where a related 1-D problem by finite difference or finite element methods is to introduce effective methods for imposing the boundary conditions through which (1.1a) is coupled to (1.1b). These models thus provide a simple starting point for examining the general treatment of boundary conditions for more general time dependent Navier-Stokes incompressible flows. For the purpose of discussion it is assumed that D is a square domain. A standard finite difference method to solve (1.1) is to introduce a uniform grid and then use standard five point finite difference operators to express each equation in (1.1). At any point on the boundary B a value of Chi is specified by the boundary conditions but a value of omega at the same boundary mesh point will also be required to complete the computation. Methods are discussed which overcome the difficulty in solving these problems.

Rose, M. E.↗
NASA NTRS1981

Analyses of shuttle orbiter approach and landing conditions

A study of one shuttle orbiter approach and landing conditions are summarized. Causes of observed PIO like flight deficiencies are identified and potential cures are examined. Closed loop pilot/vehicle analyses are described and path/attitude stability boundaries defined. The latter novel technique proved of great value in delineating and illustrating the basic causes of this multiloop pilot control problem. The analytical results are shown to be consistent with flight test and fixed base simulation. Conclusions are drawn relating to possible improvements of the shuttle orbiter/digital flight control system.

Teper, G. L.↗
DOE OSTI2017

An Approach to Quad Meshing Based On Cross Valued Maps and the Ginzburg-Landau Theory

A generalization of vector fields, referred to as N-direction fields or cross fields when N=4, has been recently introduced and studied for geometry processing, with applications in quadrilateral (quad) meshing, texture mapping, and parameterization. We make the observation that cross field design for two-dimensional quad meshing is related to the well-known Ginzburg-Landau problem from mathematical physics. This identification yields a variety of theoretical tools for efficiently computing boundary-aligned quad meshes, with provable guarantees on the resulting mesh, for example, the number of mesh defects and bounds on the defect locations. The procedure for generating the quad mesh is to (i) find a complex-valued "representation" field that minimizes the Dirichlet energy subject to a boundary constraint, (ii) convert the representation field into a boundary-aligned, smooth cross field, (iii) use separatrices of the cross field to partition the domain into four sided regions, and (iv) mesh each of these four-sided regions using standard techniques. Under certain assumptions on the geometry of the domain, we prove that this procedure can be used to produce a cross field whose separatrices partition the domain into four sided regions. To solve the energy minimization problem for the representation field, we use an extension of the Merriman-Bence-Osher (MBO) threshold dynamics method, originally conceived as an algorithm to simulate motion by mean curvature, to minimize the Ginzburg-Landau energy for the optimal representation field. Lastly, we demonstrate the method on a variety of test domains.

Computational Geometry↗
NASA NTRS1993

Flutter analysis using transversality theory

A new method of calculating flutter boundaries of undamped aeronautical structures is presented. The method is an application of the weak transversality theorem used in catastrophe theory. In the first instance, the flutter problem is cast in matrix form using a frequency domain method, leading to an eigenvalue matrix. The characteristic polynomial resulting from this matrix usually has a smooth dependence on the system's parameters. As these parameters change with operating conditions, certain critical values are reached at which flutter sets in. Our approach is to use the transversality theorem in locating such flutter boundaries using this criterion: at a flutter boundary, the characteristic polynomial does not intersect the axis of the abscissa transversally. Formulas for computing the flutter boundaries and flutter frequencies of structures with two degrees of freedom are presented, and extension to multi-degree of freedom systems is indicated. The formulas have obvious applications in, for instance, problems of panel flutter at supersonic Mach numbers.

Afolabi, D.↗
NASA NTRS1994

Evaluation of water vapor distribution in general circulation models using satellite observations

This paper presents a comparison of the water vapor distribution obtained from two general circulation models, the European Centre for Medium-Range Weather Forecasts (ECMWF) model and the National Center for Atmospheric Research (NCAR) Community Climate Model (CCM), with satellite observations of total precipitable water (TPW) from Special Sensor Microwave/Imager (SSM/I) and upper tropospheric relative humidity (UTH) from GOES. Overall, both models are successful in capturing the primary features of the observed water vapor distribution and its seasonal variation. For the ECMWF model, however, a systematic moist bias in TPW is noted over well-known stratocumulus regions in the eastern subtropical oceans. Comparison with radiosonde profiles suggests that this problem is attributable to difficulties in modeling the shallowness of the boundary layer and large vertical water vapor gradients which characterize these regions. In comparison, the CCM is more successful in capturing the low values of TPW in the stratocumulus regions, although it tends to exhibit a dry bias over the eastern half of the subtropical oceans and a corresponding moist bias in the western half. The CCM also significantly overestimates the daily variability of the moisture fields in convective regions, suggesting a problem in simulating the temporal nature of moisture transport by deep convection. Comparison of the monthly mean UTH distribution indicates generally larger discrepancies than were noted for TPW owing to the greater influence of large-scale dynamical processes in determining the distribution of UTH. In particular, the ECMWF model exhibits a distinct dry bias along the Intertropical Convergence Zone (ITCZ) and a moist bias over the subtropical descending branches of the Hadley cell, suggesting an underprediction in the strength of the Hadley circulation. The CCM, on the other hand, demonstrates greater discrepancies in UTH than are observed for the ECMWF model, but none that are as clearly correlated with well-known features of the large-scale circulation.

Soden, Brian J.↗
NASA NTRS1986

Kelvin-Helmholtz instability at the magnetopause boundary

An MHD model is defined for the convective Kelvin-Helmholtz instability (CKHI) which satellite data indicates produces vortices in the near-earth region of the plasma sheet through turbulence caused by shear flow. The MHD equations are solved as an initial value problem, which reveals that the CKHI growth rate is the same as that of the periodic KHI (PKHI). However, the CKHI is not as limited in the range of amplitude as is the PKHI, which experiences stability in nonlinear growth because of the periodic characteristics of the growth. A perpendicularly moving CKHI can produce a large vortex flow and shocks off the vortex boundary. The satellite data indicate that the vortices are in the size range of 10-20 earth radii. Additional effects can include the expulsion of magnetic flux and the formation of the boundary layer. The success of the model up to the point of formation of the boundary layer is taken as an incentive for further model development using more realistic magnetospheric data.

Wu, C. C.↗
DOE OSTI2018

Summer term 2018 work with IMC radiative transport in binary stochastic media

In the quest for new and improved methods to solve thermal radiation transport in stochastic media, this summer project was to characterize radation transport in a deterministic lattice problem that represents a real, binary stochastic media. A representative geometry for a binary stochastic media consisting of a cube of a high volume fraction media containing constant-diameter non-overlappying spheres randomly distributed within was run previously. Previous thermal radiation wave problems in binary stochastic media displayed three characteristic thermal energy deposition wavefronts in the following order, front of the spheres, bulk material, and the backs of the spheres. Several physical aspects of the representative geometry for a stochastic binary media using an optically thin background with optically thick spheres were analyzed. This analysis includes wave speed propagation using a one dimensional diffusion approximation, energy fluence transmission, and material energy deposition. The approaches to analysis included representing one sphere of the high opacity material in a diffusive background of the low opacity material, as well as repeating spheres in a column con guration with reflective boundary conditions to approximate a lattice. Additionally, both atomic mix and Levermore-Pomraning homogenization schemes were compared against the representative cube geometry as well as the sphere column lattice geometry. The wave speed analysis with a one dimensional Marshak diffusion wave was largely ineffective at determining an effective opacity value to which a problem could be homogenized, likely due to low mean free path resolution within the problem. However, the sphere column lattice geometry showed promising results for approximating the energy fluence transmitted through the problem domain. One concern with the sphere column lattice geometry indicated the high forward-bias effect of fluence reaching the problem boundary due to streaming paths around the spheres in the case of small mean free path resolution of the background material.

36 MATERIALS SCIENCE↗