Search NASA⌕ Search

SEARCH · Search NASA

Results for “boundary value problems”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

Analysis of longitudinal thickness-shear stiffness of a monolayer filamentary composite

The longitudinal thickness-shear behavior of a rectangular cross-section prism containing circular cross-section fibers embedded in a matrix is analyzed. The problem is formulated as a Saint-Venant flexure problem modeled by a tip-loaded cantilever composite beam and solution is obtained by the boundary-point least-squares method. It is shown that the longitudinal thickness-shear modulus is much greater than the in-plane shear modulus for a monofilament composite and therefore an independent analysis is necessary to predict the longitudinal thickness-shear behavior of the composite. Numerical results are presented in graphical form for various values of fiber volume fraction and material parameters typical of modern advanced composites.

Chang, S.↗
NASA NTRS1982

Finite-difference solution of the compressible stability eigenvalue problem

A compressible stability analysis computer code is developed. The code uses a matrix finite difference method for local eigenvalue solution when a good guess for the eigenvalue is available and is significantly more computationally efficient than the commonly used initial value approach. The local eigenvalue search procedure also results in eigenfunctions and, at little extra work, group velocities. A globally convergent eigenvalue procedure is also developed which may be used when no guess for the eigenvalue is available. The global problem is formulated in such a way that no unstable spurious modes appear so that the method is suitable for use in a black box stability code. Sample stability calculations are presented for the boundary layer profiles of a Laminar Flow Control (LFC) swept wing.

Malik, M. R.↗
NASA NTRS2000

Generation and Radiation of Acoustic Waves from a 2-D Shear Layer

A parallel numerical simulation of the radiation of sound from an acoustic source inside a 2-D jet is presented in this paper. This basic benchmark problem is used as a test case for scattering problems that are presently being solved by using the Impedance Mismatch Method (IMM). In this technique, a solid body in the domain is represented by setting the acoustic impedance of each medium, encountered by a wave, to a different value. This impedance discrepancy results in reflected and scattered waves with appropriate amplitudes. The great advantage of the use of this method is that no modifications to a simple Cartesian grid need to be made for complicated geometry bodies. Thus, high order finite difference schemes may be applied simply to all parts of the domain. In the IMM, the total perturbation field is split into incident and scattered fields. The incident pressure is assumed to be known and the equivalent sources for the scattered field are associated with the presence of the scattering body (through the impedance mismatch) and the propagation of the incident field through a non-uniform flow. An earlier version of the technique could only handle uniform flow in the vicinity of the source and at the outflow boundary. Scattering problems in non-uniform mean flow are of great practical importance (for example, scattering from a high lift device in a non-uniform mean flow or the effects of a fuselage boundary layer). The solution to this benchmark problem, which has an acoustic wave propagating through a non-uniform mean flow, serves as a test case for the extensions of the IMM technique.

Agarwal, Anurag↗
DOE OSTI2015

Mobility size and mass of nascent soot particles in a benchmark premixed ethylene flame

The burner stabilized stagnation flame technique coupled with micro-orifice probe sampling and mobility sizing has evolved into a useful tool for examining the evolution of the particle size distribution of nascent soot in laminar premixed flames. Several key aspects of this technique are examined through a multi-university collaborative study that involves both experimental measurement and computational modeling. Key issues examined include (a) data reproducibility and facility effects using four burners of different sizes and makers over three different facilities, (b) the mobility diameter and particle mass relationship, and (c) the degree to which the finite orifice flow rate affects the validity of the boundary condition in a pseudo one dimensional stagnation flow flame formulation. The results indicate that different burners across facilities yield nearly identical results after special attention is paid to a range of experimental details, including a proper selection of the sample dilution ratio and quantification of the experimental flame boundary conditions. The mobility size and mass relationship probed by tandem mass and mobility measurement shows that nascent soot with mobility diameter as small as 15 nm can deviate drastically from the spherical shape. Various non-spherical morphology models using a mass density value of 1.5 g/cm 3 can reconcile this discrepancy in nascent soot mass. Lastly, two-dimensional axisymmetric simulations of the experimental flame with and without the sample orifice flow reveal several problems of the pseudo one-dimensional stagnation flow flame approximation. The effect of the orifice flow on the flame and soot sampled, although small, is not negligible. Particular suggestions are provided as to how to treat the non-ideality of the experimental setup in experiment and model comparisons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1975

Steady subsonic flow around finite-thickness wings

The general method for analyzing steady subsonic potential aerodynamic flow around a lifting body having arbitrary shape is presented. By using the Green function method, an integral representation for the potential is obtained. Under small perturbation assumption, the potential at any point, P, in the field depends only upon the values of the potential and its normal derivative on the surface of the body. Hence if the point P approaches the surface of the body, the representation reduces to an integral equation relating the potential and its normal derivative (which is known from the boundary conditions) on the surface. The question of uniqueness is examined and it is shown that, for thin wings, the operator becomes singular as the thickness approaches zero. This fact may yield numerical problems for very thin wings. However, numerical results obtained for a rectangular wing in subsonic flow show that these problems do not appear even for thickness ratio tau = .001. Comparison with existing results shows that the proposed method is at least as fast and accurate as the lifting surface theories.

Kuo, C. C.↗
DOE OSTI2014

Steady fall of isothermal, resistive-viscous, compressible fluid across magnetic field

This is a basic MHD study of the steady fall of an infinite, vertical slab of isothermal, resistive-viscous, compressible fluid across a dipped magnetic field in uniform gravity. This double-diffusion steady flow in unbounded space poses a nonlinear but numerically tractable, one-dimensional (1D) free-boundary problem, assuming constant coefficients of resistivity and viscosity. The steady flow is determined by a dimensionless number μ{sub 1} proportional to the triple product of the two diffusion coefficients and the square of the linear total mass. For a sufficiently large μ{sub 1}, the Lorentz, viscous, fluid-pressure, and gravitational forces pack and collimate the fluid into a steady flow of a finite width defined by the two zero-pressure free-boundaries of the slab with vacuum. The viscous force is essential in this collimation effect. The study conjectures that in the regime μ{sub 1}→0, the 1D steady state exists only for μ{sub 1}∈Ω, a spectrum of an infinite number of discrete values, including μ{sub 1} = 0 that corresponds to two steady states, the classical zero-resistivity static slab of Kippenhahn and Schlüter [R. Kippenhahn and A. Schlüter, Z. Astrophys. 43, 36 (1957)] and its recent generalization [B. C. Low et al., Astrophys. J. 755, 34 (2012)] to admit an inviscid resistive flow. The pair of zero-pressure boundaries of each of the μ{sub 1}→0 steady-state slabs are located at infinity. Computational evidence suggests that the Ω steady-states are densely distributed around μ{sub 1} = 0, as an accumulation point, but are sparsely separated by open intervals of μ{sub 1}-values for which the slab must be either time-dependent or spatially multi-dimensional. The widths of these intervals are vanishingly small as μ{sub 1}→0. This topological structure of physical states is similar to that described by Landau and Liftshitz [L. D. Landau and E. M. Lifshitz, Fluid Mechanics (Addison-Wesley, Reading, MA, 1959)] to explain the onset of hydrodynamic turbulence. The implications of this MHD study are discussed, with an interest in the prominences in the solar atmosphere and the interstellar clouds in the Galaxy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2000

Supercomputing Aspects for Simulating Incompressible Flow

The primary objective of this research is to support the design of liquid rocket systems for the Advanced Space Transportation System. Since the space launch systems in the near future are likely to rely on liquid rocket engines, increasing the efficiency and reliability of the engine components is an important task. One of the major problems in the liquid rocket engine is to understand fluid dynamics of fuel and oxidizer flows from the fuel tank to plume. Understanding the flow through the entire turbo-pump geometry through numerical simulation will be of significant value toward design. One of the milestones of this effort is to develop, apply and demonstrate the capability and accuracy of 3D CFD methods as efficient design analysis tools on high performance computer platforms. The development of the Message Passage Interface (MPI) and Multi Level Parallel (MLP) versions of the INS3D code is currently underway. The serial version of INS3D code is a multidimensional incompressible Navier-Stokes solver based on overset grid technology, INS3D-MPI is based on the explicit massage-passing interface across processors and is primarily suited for distributed memory systems. INS3D-MLP is based on multi-level parallel method and is suitable for distributed-shared memory systems. For the entire turbo-pump simulations, moving boundary capability and efficient time-accurate integration methods are built in the flow solver, To handle the geometric complexity and moving boundary problems, an overset grid scheme is incorporated with the solver so that new connectivity data will be obtained at each time step. The Chimera overlapped grid scheme allows subdomains move relative to each other, and provides a great flexibility when the boundary movement creates large displacements. Two numerical procedures, one based on artificial compressibility method and the other pressure projection method, are outlined for obtaining time-accurate solutions of the incompressible Navier-Stokes equations. The performance of the two methods is compared by obtaining unsteady solutions for the evolution of twin vortices behind a flat plate. Calculated results are compared with experimental and other numerical results. For an unsteady flow, which requires small physical time step, the pressure projection method was found to be computationally efficient since it does not require any subiteration procedure. It was observed that the artificial compressibility method requires a fast convergence scheme at each physical time step in order to satisfy the incompressibility condition. This was obtained by using a GMRES-ILU(0) solver in present computations. When a line-relaxation scheme was used, the time accuracy was degraded and time-accurate computations became very expensive.

Kwak, Dochan↗
NASA NTRS1994

Coupled Marangoni-Benard/Rayleigh-Benard Instability with Temperature Dependent Viscosity

The onset of convection induced by coupled surface tension gradient and buoyancy forces is investigated with temperature dependent viscosity. Both surface tension and viscosity are assumed to vary linearly with temperature. The limiting case, Ma=0, is the buoyancy driven convection problem typically referred to as the Rayleigh-Benard stability problem. The other limiting case, Ra=0, is the surface tension (gradient) driven flow problem referred to as the Marangoni-Benard problem. The equations and boundary conditions obtained from the linear analysis are solved numerically as a generalized eigenvalue problem. Neutral stability curves for different viscosity slopes have been generated for the Marangoni-Benard and Rayleigh-Benard problems. It is shown that the curves can be collapsed to a single curve by appropriately scaling the results for each of the limiting cases. The critical Marangoni number is determined as a function of the slope of the viscosity temperature variation, epsilon, for different values of the Rayleigh number. When the viscosity decreases linearly with temperature, the coupled buoyancy-surface tension problem, including the limiting cases of Ra=0 and Ma=0, is found to be more stable than the constant viscosity case.

Skarda, J. Raymond Lee↗
NASA NTRS2000

Towards Flange-to-Flange Turbopump Simulations for Liquid Rocket Engines

The primary objective of this research is to support the design of liquid rocket systems for the Advanced Space Transportation System. Since the space launch systems in the near future are likely to rely on liquid rocket engines, increasing the efficiency and reliability of the engine components is an important task. One of the major problems in the liquid rocket engine is to understand fluid dynamics of fuel and oxidizer flows from the fuel tank to plume. Understanding the flow through the entire turbopump geometry through numerical simulation will be of significant value toward design. This will help to improve safety of future space missions. One of the milestones of this effort is to develop, apply and demonstrate the capability and accuracy of 3D CFD methods as efficient design analysis tools on high performance computer platforms. The development of the MPI and MLP versions of the INS3D code is currently underway. The serial version of INS3D code is a multidimensional incompressible Navier-Stokes solver based on overset grid technology. INS3D-MPI is based on the explicit massage-passing interface across processors and is primarily suited for distributed memory systems. INS3D-MLP is based on multi-level parallel method and is suitable for distributed-shared memory systems. For the entire turbopump simulations, moving boundary capability and an efficient time-accurate integration methods are build in the flow solver. To handle the geometric complexity and moving boundary problems, overset grid scheme is incorporated with the solver that new connectivity data will be obtained at each time step. The Chimera overlapped grid scheme allows subdomains move relative to each other, and provides a great flexibility when the boundary movement creates large displacements. The performance of the two time integration schemes for time-accurate computations is investigated. For an unsteady flow which requires small physical time step, the pressure projection method was found to be computationally efficient since it does not require any subiterations procedure. It was observed that the artificial compressibility method requires a fast convergence scheme at each physical time step in order to satisfy incompressibility condition. This was obtained by using a GMRES-ILU(0) solver in our computations. When a line-relaxation scheme was used, the time accuracy was degraded and time-accurate computations became very expensive. The current geometry for the LOX boost turbopump has various rotating and stationary components, such as inducer, stators, kicker, hydrolic turbine, where the flow is extremely unsteady. Figure 1 shows the geometry and computed surface pressure of the inducer. The inducer and the hydrolic turbine rotate in different rotational speed.

Kiris, Cetin↗
NASA NTRS1979

Liquid oxygen/liquid hydrogen auxiliary power system thruster investigation

The design, fabrication, and demonstration of a 111 newton (25 lb) thrust, integrated auxiliary propulsion system (IAPS) thruster for use with LH2/LO2 propellants is described. Hydrogen was supplied at a temperature range of 22 to 33 K (40 to 60 R), and oxygen from 89 to 122 K (160 to 220 R). The thruster was designed to operate in both pulse mode and steady-state modes for vehicle attitude control, space maneuvering, and as an abort backup in the event of failure of the main propulsion system. A dual-sleeve, tri-axial injection system was designed that utilizes a primary injector/combustor where 100 percent of the oxygen and 8 percent of the hydrogen is introduced; a secondary injector/combustor where 45 percent of the hydrogen is introduced to mix with the primary combustor gases; and a boundary layer injector that uses the remaining 45 percent of the hydrogen to cool the thrust throat/nozzle design. Hot-fire evaluation of this thruster with a BLC injection distance of 2.79 cm (1.10 in.) indicated that a specific impulse value of 390 sec can be attained using a coated molybdenum thrust chamber. Pulse mode tests indicated that a chamber pressure buildup to 90 percent thrust can be achieved in a time on the order of 48 msec. Some problems were encountered in achieving ignition of each pulse during pulse trains. This was interpreted to indicate that a higher delivered spark energy level ( 100 mJ) would be required to maintain ignition reliability of the plasma torch ignition system under the extra 'cold' conditions resulting during pulsing.

Eberle, E. E.↗
NASA NTRS1985

A method for estimating soil moisture availability

A method for estimating values of soil moisture based on measurements of infrared surface temperature is discussed. A central element in the method is a boundary layer model. Although it has been shown that soil moistures determined by this method using satellite measurements do correspond in a coarse fashion to the antecedent precipitation, the accuracy and exact physical interpretation (with respect to ground water amounts) are not well known. This area of ignorance, which currently impedes the practical application of the method to problems in hydrology, meteorology and agriculture, is largely due to the absence of corresponding surface measurements. Preliminary field measurements made over France have led to the development of a promising vegetation formulation (Taconet et al., 1985), which has been incorporated in the model. It is necessary, however, to test the vegetation component, and the entire method, over a wide variety of surface conditions and crop canopies.

Carlson, T. N.↗
NASA NTRS1991

Transformation of two and three-dimensional regions by elliptic systems

A reliable linear system is presented for grid generation in 2-D and 3-D. The method is robust in the sense that convergence is guaranteed but is not as reliable as other nonlinear elliptic methods in generating nonfolding grids. The construction of nonfolding grids depends on having reasonable approximations of cell aspect ratios and an appropriate distribution of grid points on the boundary of the region. Some guidelines are included on approximating the aspect ratios, but little help is offered on setting up the boundary grid other than to say that in 2-D the boundary correspondence should be close to that generated by a conformal mapping. It is assumed that the functions which control the grid distribution depend only on the computational variables and not on the physical variables. Whether this is actually the case depends on how the grid is constructed. In a dynamic adaptive procedure where the grid is constructed in the process of solving a fluid flow problem, the grid is usually updated at fixed iteration counts using the current value of the control function. Since the control function is not being updated during the iteration of the grid equations, the grid construction is a linear procedure. However, in the case of a static adaptive procedure where a trial solution is computed and used to construct an adaptive grid, the control functions may be recomputed at every step of the grid iteration.

Mastin, C. Wayne↗
NASA NTRS1994

Multilayer Thin Film Polarizer Design for Far Ultraviolet using Induced Transmission and Absorption Technique

Good theoretical designs of far ultraviolet polarizers have been reported using a MgF2/Al/MgF2 three layer structure on a thick Al layer as a substrate. The thicknesses were determined to induce transmission and absorption of p-polarized light. In these designs Al optical constants were used from films produced in ultrahigh vacuum (UHV: 10(exp -10) torr). Reflectance values for polarizers fabricated in a conventional high vacuum (p approx. 10(exp -6 torr)) using the UHV design parameters differed dramatically from the design predictions. Al is a highly reactive material and is oxidized even in a high vacuum chamber. In order to solve the problem other metals have been studied. It is found that a larger reflectance difference is closely related to higher amplitude and larger phase difference of Fresnel reflection coefficients between two polarizations at the boundary of MgF2/metal. It is also found that for one material a larger angle of incidence from the surface normal brings larger amplitude and phase difference. Be and Mo are found good materials to replace Al. Polarizers designed for 121.6 nm with Be at 60 deg and with Mo at 70 deg are shown as examples.

Kim, Jongmin↗
NASA NTRS2015

Analysis of Advanced Thermoelectric Materials and Their Functional Limits

The world's demand for energy is increasing dramatically, but the best energy conversion systems operate at approximately 30% efficiency. One way to decrease energy loss is in the recovery of waste heat using thermoelectric (TE) generators. A TE generator is device that generates electricity by exploiting heat flow across a thermal gradient. The efficiency of a TE material for power generation and cooling is determined by the dimensionless Figure of Merit (ZT): ZT = S(exp. 2)sigmaT/κ: where S is the Seebeck coefficient, sigma is the electrical conductivity, T is the absolute temperature, and κ is the thermal conductivity. The parameters are not physically independent, but intrinsically coupled since they are a function of the transport properties of electrons. Traditional research on TE materials has focused on synthesizing bulk semiconductor-type materials that have low thermal conductivity and high electrical conductivity affording ZT values of 1. The optimization of the σ/κ ratio is difficult to achieve using current material formats, as these material constants are complementary. Recent areas of research are focusing on using nanostructural artifacts that introduce specific dislocations and boundary conditions that scatter the phonons. This disrupts the physical link between thermal (phonon) and electrical (electron) transport. The result is that κ is decreased without decreasing σ. These material formats give ZT values of up to 2 which represent approximately 18% energy gain from waste heat recovery. The next challenge in developing the next generation of TE materials with superior performance is to tailor the interconnected thermoelectric physical parameters of the material system. In order to approach this problem, the fundamental physics of each parameter S, sigma, and κ need to be physically understood in their context of electron/phonon interaction for the construction of new high ZT thermoelectric devices. Is it possible to overcome the physical limit imposed by of the effect of phonon lattice oscillation and energetic electrons towards thermal conductivity? Is the Seebeck coefficient, based on the difference in voltage over temperature gradient ( deltaV/deltaT), an intrinsic parameter of each material? All these parameters were manipulated using nano-bridge and twin-lattice structural concepts at the NASA Langley Research Center. This talk will review the current trend of TE research to optimize the ZT and discuss about new approaches on increasing ZT within functional limits of each parameter.

Kim, Hyun Jung↗
NASA NTRS2021

Evaluation of Stereology for Snow Microstructure Measurement and Microwave Emission Modeling: A Case Study

Reliable microstructure measurement of snow is a requirement for microwave radiative transfer model validation. Snow specific surface area (SSA) can be measured using stereological methods, in which snow samples are cast in the field and photographed in the laboratory. Processing stereology photographs manually by counting intersections of test cycloids with air–ice boundaries reduces the problems in binary segmentation. This paper is a case study to evaluate the repeatability of the manually stereology interpretation by two independent research groups. We further assessed how uncertainty in snow SSA influences simulated brightness temperature(TB) driven by the Microwave Emission Model of Layered Snowpacks (MEMLS), and how stereology compares to Near Infrared (NIR) camera and hand lens. Data was obtained from two alpine snow profiles from Steamboat Springs, Colorado. Results showed that stereological SSA values measured by two groups are highly consistent, and the ground radiometer measured TBat 19 and 37 GHz was successfully predicted (RMSE<3.8 K);simulations using NIR SSA and hand-lens geometric grain size (Dg)measurements have larger errors. This conclusion was not sensitive to uncertainty in the free parameters of TB modeling.

Microwave radiometry↗
DOE OSTI2024

An Efficient Numerical Algorithm for Solving Coupled Time-Dependent Ginzburg-Landau Equation for Superconductivity and Elasticity

A decoupled finite element algorithm is developed for simulating the vortex dynamics on an elastic superconductor which couples the time-dependent Ginzburg- Landau equation with the complex-valued superconducting order parameter and the vector-valued magnetic potential, and the elasticity equation. We present an iterative algorithm for the decoupled system arising from the time and spatial discretization using a combination of preconditioner, algebraic multigrid method (AMG) and preconditioned conjugate gradient method (PCG). The iterative algorithm allows us to perform large-scale three-dimensional simulations of mesoscale pattern formation during superconducting phase transitions with arbitrary elastic boundary conditions. Here, the performance and efficiency of the algorithm are numerically verified by several benchmark problems, exhibiting up to two orders of magnitude improvement depending on the scale of discrete system compared to the exact solver.

Efficiency↗
DOE OSTI2017

Influence of Sedimentation Length on the Convective Stability of a Colloidal Suspension

The convective stability of a colloidal suspension is studied in the case when the vertical dimension of the cavity is less than or comparable with the sedimentation length of nanoparticles. The analysis is carried out within the Boussinesq approximation on the basis of a modified model that takes into account the dependence of a thermodiffusion flow on the local value of impurity concentration. A new parameter of the problem is the ratio of the sedimentation length to the vertical dimension of the cavity. For a quiescent colloidal suspension, exact and approximate (in the case of small concentrations) solutions are obtained that describe the distributions of nanoparticles. A transformation is obtained that allows one to investigate the convective stability of a colloidal suspension stratified in the gravitational field by the Galerkin method with a set of simple trial functions. Instability boundaries and the characteristics of critical perturbations are determined. It is shown that, in the case of negative thermodiffusion, a decrease in the sedimentation length leads to a decrease in the convection threshold and the frequency of neutral oscillations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗