Search NASA⌕ Search

SEARCH · Search NASA

Results for “BOUNDARY CONDITIONS”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Computer-aided analysis of nonlinear problems in transport phenomena

The paper describes algorithms for equilibrium and steady-state problems with coefficients in the expansions derived by the Galerkin weighted residual method and calculated from the resulting sets of nonlinear algebraic equations by the Newton-Raphson method. Initial approximations are obtained from nearby solutions by continuation techniques as parameters are varied. The Newton-Raphson technique is preferred because the Jacobian of the solution is useful for continuation, for analyzing the stability of solutions, for detecting bifurcation of solution families, and for computing asymptotic estimates of the effects on any solution of small changes in parameters, boundary conditions, and boundary shape.

Brown, R. A.↗
NASA NTRS1980

Effect of grid system on finite element calculation

Detailed parametric studies of the effect of grid system on finite element calculation for potential flows were made. These studies led to the formulation of a design criteria for optimum mesh system and the development of two methods to generate the optimum mesh system. The guidelines for optimum mesh system are: (1) the mesh structure should be regular; (2) the element should be as regular and equilateral as possible; (3) the distribution of size of element should be consistent with that of flow variables to insure maximum uniformity in error distribution; (4) for non-Dirichlet boundary conditions, smaller boundary elements or higher order interpolation functions should be used; and (5) the mesh should accommodate the boundary geometry as accurately as possible. The results of the parametric studies are presented.

Lee, K. D.↗
NASA NTRS1985

Thermal stress analysis of a new turbine shroud seal concept

The thermal stress field of a two piece turbine shroud seal concept was analyzed and results compared to one piece designs by finite element analysis. The two piece seal has independently formed structure (substrate) and ceramic components that are assembled at ambient conditions. The boundary conditions used for analysis were hot gas surface temperatures of 1370 and 1650 C (2500 and 3000 F) and cooled surface temperature of 700 C (1285 F). The resulting thermal stress field, of the two piece seal when compared to the one piece seals in the region of all ceramic material, was reduced substantially.

Handschuh, R. F.↗
NASA NTRS1985

The role of isoprene oxidation in the tropospheric ozone budget in the tropics

A comprehensive chemical mechanism for the oxidation of isoprene (a hydrocarbon, C5H8 emitted primarily by vegetation) by OH and O3 in the troposphere was developed and incorporated into a one-dimensional steady-state photochemical model of the troposphere. Flux boundary conditions for NOx (NO + NO2), HNO3, O3, and CO were used to investigate the changes produced in the tropospheric concentrations and integrated column of ozone from including isoprene chemistry in the model. Two calculations were performed at 15 deg N latitude for annual conditions using identical flux boundary conditions for NOx, HNO3, O3, and CO; in one calculation, the chemistry describing isoprene oxidation was included while in the other it was not. Both sets of calculations included reactions describing the chemistry of anthropogenic nonmethane hydrocarbons. The calculations showed decreases in concentrations of ozone throughout the troposphere when isoprene chemistry was included. Concentrations of NOx and HNO3 increased in the lower troposphere and decreased in the upper troposphere while concentrations of CO and PAN increased throughout the troposphere when isoprene chemistry was included. Implications of this study to the budgets of these species in the tropics is discussed.

Brewer, D. A.↗
NASA NTRS1987

Electromagnetic backscattering by plates and disks

With the recent development of diffraction coefficients for imperfectly conducting half-planes, it has become possible to analyze a wide variety of problems for which the impedance surface boundary condition applies. This impedance boundary condition, while approximate, was utilized to extend the usefulness of the Uniform Geometrical Theory of Diffraction (UTD) beyond the perfectly conducting geometries. These half-plane diffraction coefficients are used to analyze patterns of an antenna in the presence of an imperfectly conducting flat polygonal plate. The Geometrical Theory of Diffraction (GTD) techniques were also used to investigate the backscattering from perfectly conducting plates. To further improve the soft polarization results for wide angles, a model for the creeping wave or circulating current on the edge of the disk was obtained and used to find an additional component of the backscattered field. The backscattering from a square plate was then analyzed using GTD. Backscattering in both the principal and off-principal planes was examined.

Balanis, Constantine A.↗
NASA NTRS1989

The morphology of strike-slip faults - Examples from the San Andreas Fault, California

The dilatational strains associated with vertical faults embedded in a horizontal plate are examined in the framework of fault kinematics and simple displacement boundary conditions. Using boundary element methods, a sequence of examples of dilatational strain fields associated with commonly occurring strike-slip fault zone features (bends, offsets, finite rupture lengths, and nonuniform slip distributions) is derived. The combinations of these strain fields are then used to examine the Parkfield region of the San Andreas fault system in central California.

Bilham, Roger↗
NASA NTRS1991

Lithospheric extension near Lake Mead, Nevada - A model for ductile flow in the lower crust

Small variations in gravity anomalies and topographic elevation observed in areas that have undergone highly variable amounts of upper crustal thinning can be satisfactorily explained by ductile flow of lower crustal material under the proper conditions. The boundary between the unextended Colorado Plateau and a strongly extended domain in the Basin and Range Province in the Lake Mead (Nevada) region is examined. Finite element modeling of Newtonian flow and power law creep shows that flow over the length scale of the eastern Basin and Range (500) km or more) corresponding to upper crustal extension by a factor of 1.4-3 over 10 million years requires effective viscosities less than 10 to the 18th - 10 to the 20th Pa s for ductile channels 10-25 km thick. Modeling suggests that these effective viscosities may be sustained by lower crustal material deforming at laboratory-derived power law creep rates. The longer-scale flow may require elevated crustal temperatures (more than 700 C), depending on the composition and material properties assumed. Under the boundary conditions assumed in this study the linear viscous flow models yield a satisfactory approximation to deformation by power law creep. This work suggests that flow in the lower crust may be a viable mechanism for producing small variations in total crustal thickness between strongly extended and less extended regions, and thereby explaining the relative uniformity in gravity and topography between such regions.

Kruse, Sarah↗
NASA NTRS1993

Comparisons between DSMC and the Navier-Stokes equations for reentry flows

A detailed comparison is made between Navier-Stokes and DSMC calculations for flows near the continuum limit to assess the accuracy of the continuum equations in this regime. Meaningful comparisons require the use of similar physical models. This necessitates the inclusion of a separate rotational energy equation and use of slip boundary conditions. Inclusion of slip boundary conditions resulted in improved agreement between surface properties. Moreover, good agreement was obtained for the various temperatures in the nonequilibrium portion of the flow field that does not contain the shock region. Departures are noted in the shock region and in regions where thermal diffusion effects are important.

Olynick, David R.↗
NASA NTRS1994

Comparison of PDF and Moment Closure Methods in the Modeling of Turbulent Reacting Flows

In modeling turbulent reactive flows, Probability Density Function (PDF) methods have an advantage over the more traditional moment closure schemes in that the PDF formulation treats the chemical reaction source terms exactly, while moment closure methods are required to model the mean reaction rate. The common model used is the laminar chemistry approximation, where the effects of turbulence on the reaction are assumed negligible. For flows with low turbulence levels and fast chemistry, the difference between the two methods can be expected to be small. However for flows with finite rate chemistry and high turbulence levels, significant errors can be expected in the moment closure method. In this paper, the ability of the PDF method and the moment closure scheme to accurately model a turbulent reacting flow is tested. To accomplish this, both schemes were used to model a CO/H2/N2- air piloted diffusion flame near extinction. Identical thermochemistry, turbulence models, initial conditions and boundary conditions are employed to ensure a consistent comparison can be made. The results of the two methods are compared to experimental data as well as to each other. The comparison reveals that the PDF method provides good agreement with the experimental data, while the moment closure scheme incorrectly shows a broad, laminar-like flame structure.

Norris, Andrew T.↗
NASA NTRS1999

Preliminary Study of a Pull Plug Friction Weld

A pull plug friction weld, simply defined, comprises inserting a rotating cone-shaped cylinder into a somewhat cone-shaped hole in a plate. The rotating plug makes contact with the edge of the plate and the resulting friction generates heat. The temperature of the plate material eventually reaches a magnitude that will cause the plate material at the edge of the hole to flow. This can be termed a temperature dependent plastic flow. The rotation of the plug is terminated, additional pressure is applied and the metal at the interface of the two materials cools and welding occurs. This preliminary study addresses only three aspects of a complete analysis that is multi-faceted. The transient temperature distribution for different pull plug configurations has been studied in some detail even though the initial conditions and boundary conditions may still be deemed tentative. The stress distribution within the pull plug caused by the heating pressure was studied along with a preliminary analysis of the thermoelastic stress distribution caused by friction heating. There are no definitive results for the stress analysis. Additional study will be required.

Buchanan, George R.↗
NASA NTRS2002

New Boundary Constraints for Elliptic Systems used in Grid Generation Problems

This paper discusses new boundary constraints for elliptic partial differential equations as used in grid generation problems in generalized curvilinear coordinate systems. These constraints, based on the principle of local conservation of thermal energy in the vicinity of the boundaries, are derived using the Green's Theorem. They uniquely determine the so called decay parameters in the source terms of these elliptic systems. These constraints' are designed for boundary clustered grids where large gradients in physical quantities need to be resolved adequately. It is observed that the present formulation also works satisfactorily for mild clustering. Therefore, a closure for the decay parameter specification for elliptic grid generation problems has been provided resulting in a fully automated elliptic grid generation technique. Thus, there is no need for a parametric study of these decay parameters since the new constraints fix them uniquely. It is also shown that for Neumann type boundary conditions, these boundary constraints uniquely determine the solution to the internal elliptic problem thus eliminating the non-uniqueness of the solution of an internal Neumann boundary value grid generation problem.

Kaul, Upender K.↗
NASA NTRS2022

Test data analysis of the thermodynamic vent system-augmented top spray injector liquid nitrogen transfer experiments

Traditionally, a cryogenic tank must be pre-chilled to some “target” temperature before the main vent valve can be closed to attempt a non-vented fill (NVF) of cryogenic liquid propellant. This methodology is particularly attractive for performing in-space transfer of cryogens due to the unknown location of the liquid/vapor interface in microgravity and the high likelihood of venting liquid if the vent valve is opened during transfer. This paper presents in-depth test data analysis of a Thermodynamic Vent System (TVS) augmented injector used for cryogenic tank chilldown and fill experiments of a thin-walled Titanium tank. Eight tests were conducted using liquid nitrogen across a range of inlet conditions and boundary conditions, and three different chilldown/fill methods. For four of the tests, the injector sprays liquid into the tank as normal, but also uses a TVS heat exchanger to cool the metallic injector itself as well as the main incoming liquid stream. Results show that using the TVS augmented injector simplifies transfer operation via enhanced condensation at the injector surface at the cost of sacrificing only a small amount of propellant.

No-vent Fill↗
NASA NTRS1991

The asymptotic spectra of banded Toeplitz and quasi-Toeplitz matrices

Toeplitz matrices occur in many mathematical, as well as, scientific and engineering investigations. This paper considers the spectra of banded Toeplitz and quasi-Toeplitz matrices with emphasis on non-normal matrices of arbitrarily large order and relatively small bandwidth. These are the type of matrices that appear in the investigation of stability and convergence of difference approximations to partial differential equations. Quasi-Toeplitz matrices are the result of non-Dirichlet boundary conditions for the difference approximations. The eigenvalue problem for a banded Toeplitz or quasi-Toeplitz matrix of large order is, in general, analytically intractable and (for non-normal matrices) numerically unreliable. An asymptotic (matrix order approaches infinity) approach partitions the eigenvalue analysis of a quasi-Toeplitz matrix into two parts, namely the analysis for the boundary condition independent spectrum and the analysis for the boundary condition dependent spectrum. The boundary condition independent spectrum is the same as the pure Toeplitz matrix spectrum. Algorithms for computing both parts of the spectrum are presented. Examples are used to demonstrate the utility of the algorithms, to present some interesting spectra, and to point out some of the numerical difficulties encountered when conventional matrix eigenvalue routines are employed for non-normal matrices of large order. The analysis for the Toeplitz spectrum also leads to a diagonal similarity transformation that improves conventional numerical eigenvalue computations. Finally, the algorithm for the asymptotic spectrum is extended to the Toeplitz generalized eigenvalue problem which occurs, for example, in the stability of Pade type difference approximations to differential equations.

Beam, Richard M.↗
NASA NTRS1994

Three-dimensional fully spectral numerical method for mantle convection with depth-dependent properties

A semi-implicit fully spectral collocation method for the simulation of three-dimensional mantle convection with depth-dependent thermo-dynamic and transport properties is presented. The variable property Navier-Stokes equation expressed in terms of the primitive variable velocity and pressure is solved with the mass continuity and temperature equations. The periodic horizontal boundary conditions allow a Fourier expansion for the two horizontal directions. The stress-free, impermeable isothermal boundary conditions along with the depth dependent coefficients are handled with a Chebyshev expansion in the vertical direction. In the limit of an infinite Prandtl number appropriate to mantle convection, the inertial terms in the momentum equation are unimportant. In this case an explicit solution of a Poisson equation for pressure can be avoided; instead a fourth-order equation for vertical velocity can be solved. Simultaneous imposition of both impermeable and continuity boundary conditions during the vertical velocity evaluation is discussed. The pressure distributions on the top and bottom bounding planes were determined by means of an influence matrix technique. The numerical method employed here avoids time-splitting errors and enforces velocity boundary conditions and continuity over the entire domain, including the boundaries, to machine accuracy. Strongly time-dependent three-dimensional solutions up to a surface Rayleigh number of 1 x 10(exp 7) have been obtained. Strong upwellings, pulsating chaotically, are formed by the collective merging of cylindrical plumes.

Balachandar, S.↗
NASA NTRS1998

Inlet-Compressor Analysis Using Coupled CFD Codes

Propulsion performance and operability are key factors in the development of a successful aircraft. For high-speed supersonic aircraft, mixed-compression inlets offer high performance but are susceptible to an instability referred to as unstart. An unstart occurs when a disturbance originating in the atmosphere or the engine causes the shock system to be expelled from the inlet. This event can have adverse effects on control of the aircraft, which is unacceptable for a passenger plane such as the high speed civil transport (HSCT). The ability to predict the transient response of such inlets to flow perturbations is, therefore, important to the proper design of the inlet and the control measures used to prevent unstart. Computational fluid dynamics (CFD) is having an increasing role in the analysis of individual propulsion components. isolated inlet studies are relatively easy to perform, but a major uncertainty is the boundary condition used at the inlet exit to represent the engine - the so-called compressor face boundary condition. A one-dimensional (I-D) Euler inlet simulation (ref. 1) showed that the predicted inlet unstart tolerance to free-stream pressure perturbations can vary by as much as a factor of about six, depending on the boundary condition used. Obviously, a thorough understanding of dynamic interactions between inlets and compressors/fans is required to provide the proper boundary condition.

Cole, Gary↗
NASA NTRS1995

Review of Combustion-acoustic Instabilities

Combustion-acoustic instabilities occur when the acoustic energy increase due to the unsteady heat release of the flame is greater than the losses of acoustic energy from the system. The problem of combustion-acoustic instability is a concern in many devices for various reasons, as each device may have a unique mechanism causing unsteady heat release rates and many have unique boundary conditions. To accurately predict and quantify combustion-acoustic stabilities, the unsteady heat release rate and boundary conditions need to be accurately determined. The present review brings together work performed on a variety of practical combustion devices. Many theoretical and experimental investigations of the unsteady heat release rate have been performed, some based on perturbations in the fuel delivery system particularly for rocket instabilities, while others are based on hydrodynamic processes as in ramjet dump combustors. The boundary conditions for rocket engines have been analyzed and measured extensively. However, less work has been done to measure acoustic boundary conditions in many other combustion systems.

Oyediran, Ayo↗
NASA NTRS1993

A method for the modelling of porous and solid wind tunnel walls in computational fluid dynamics codes

Porous wall wind tunnels have been used for several decades and have proven effective in reducing wall interference effects in both low speed and transonic testing. They allow for testing through Mach 1, reduce blockage effects and reduce shock wave reflections in the test section. Their usefulness in developing computational fluid dynamics (CFD) codes has been limited, however, by the difficulties associated with modelling the effect of a porous wall in CFD codes. Previous approaches to modelling porous wall effects have depended either upon a simplified linear boundary condition, which has proven inadequate, or upon detailed measurements of the normal velocity near the wall, which require extensive wind tunnel time. The current work was initiated in an effort to find a simple, accurate method of modelling a porous wall boundary condition in CFD codes. The development of such a method would allow data from porous wall wind tunnels to be used more readily in validating CFD codes. This would be beneficial when transonic validations are desired, or when large models are used to achieve high Reynolds numbers in testing. A computational and experimental study was undertaken to investigate a new method of modelling solid and porous wall boundary conditions in CFD codes. The method utilized experimental measurements at the walls to develop a flow field solution based on the method of singularities. This flow field solution was then imposed as a pressure boundary condition in a CFD simulation of the internal flow field. The effectiveness of this method in describing the effect of porosity changes on the wall was investigated. Also, the effectiveness of this method when only sparse experimental measurements were available has been investigated. The current work demonstrated this approach for low speed flows and compared the results with experimental data obtained from a heavily instrumented variable porosity test section. The approach developed was simple, computationally inexpensive, and did not require extensive or intrusive measurements of the boundary conditions during the wind tunnel test. It may be applied to both solid and porous wall wind tunnel tests.

Beutner, Thomas John↗
NASA NTRS1986

Application of a low order panel method to complex three-dimensional internal flow problems

An evaluation of the ability of a low order panel method to predict complex three-dimensional internal flow fields was made. The computer code VSAERO was used as a basis for the evaluation. Guidelines for modeling internal flow geometries were determined and the effects of varying the boundary conditions and the use of numerical approximations on the solutions accuracy were studied. Several test cases were run and the results were compared with theoretical or experimental results. Modeling an internal flow geometry as a closed box with normal velocities specified on an inlet and exit face provided accurate results and gave the user control over the boundary conditions. The values of the boundary conditions greatly influenced the amount of leakage an internal flow geometry suffered and could be adjusted to eliminate leakage. The use of the far-field approximation to reduce computation time influenced the accuracy of a solution and was coupled with the values of the boundary conditions needed to eliminate leakage. The error induced in the influence coefficients by using the far-field approximation was found to be dependent on the type of influence coefficient, the far-field radius, and the aspect ratio of the panels.

Ashby, D. L.↗