SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

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

At least 793 records · Page 44

All sources NASA NTRS DOE OSTI

NASA NTRS1982

Digital imaging techniques in experimental stress analysis

Digital imaging techniques are utilized as a measure of surface displacement components in laser speckle metrology. An image scanner which is interfaced to a computer records and stores in memory the laser speckle patterns of an object in a reference and deformed configuration. Subsets of the deformed images are numerically correlated with the references as a measure of surface displacements. Discrete values are determined around a closed contour for plane problems which then become input into a boundary integral equation method in order to calculate surface traction in the contour. Stresses are then calculated within this boundary. The solution procedure is illustrated by a numerical example of a case of uniform tension.

Peters, W. H.↗

NASA NTRS1991

Numerical recovery of material parameters in Euler-Bernoulli beam models

A fully Sinc-Galerkin method for recovering the spatially varying stiffness parameter in fourth-order time-dependence problems with fixed and cantilever boundary conditions is presented. The forward problems are discretized with a sinc basis in both the spatial and temporal domains. This yields an approximation solution which converges exponentially and is valid on the infinite time interval. When the forward methods are applied to parameter recovery problems, the resulting inverse problems are ill-posed. Tikhonov regularization is applied and the resulting minimization problems are solved via a quasi-Newton/trust region algorithm. The L-curve method is used to determine an appropriate value of the regularization parameter. Numerical results which highlight the method are given for problems with both fixed and cantilever boundary conditions.

Smith, R. C.↗

NASA NTRS2017

Wall-Resolved Large-Eddy Simulation of Flow Separation Over NASA Wall-Mounted Hump

This paper reports the findings from a study that applies wall-resolved large-eddy simulation to investigate flow separation over the NASA wall-mounted hump geometry. Despite its conceptually simple flow configuration, this benchmark problem has proven to be a challenging test case for various turbulence simulation methods that have attempted to predict flow separation arising from the adverse pressure gradient on the aft region of the hump. The momentum-thickness Reynolds number of the incoming boundary layer has a value that is near the upper limit achieved by recent direct numerical simulation and large-eddy simulation of incompressible turbulent boundary layers. The high Reynolds number of the problem necessitates a significant number of grid points for wall-resolved calculations. The present simulations show a significant improvement in the separation-bubble length prediction compared to Reynolds-Averaged Navier-Stokes calculations. The current simulations also provide good overall prediction of the skin-friction distribution, including the relaminarization observed over the front portion of the hump due to the strong favorable pressure gradient. We discuss a number of problems that were encountered during the course of this work and present possible solutions. A systematic study regarding the effect of domain span, subgrid-scale model, tunnel back pressure, upstream boundary layer conditions and grid refinement is performed. The predicted separation-bubble length is found to be sensitive to the span of the domain. Despite the large number of grid points used in the simulations, some differences between the predictions and experimental observations still exist (particularly for Reynolds stresses) in the case of the wide-span simulation, suggesting that additional grid resolution may be required.

Uzun, Ali↗

NASA NTRS1965

On the gravitational paradox

Solution of gravitational paradox problem using Poisson motion equations without boundary conditions

MOTION EQUATION↗

NASA NTRS1986

Acoustic-Liner Admittance in a Duct

Method calculates admittance from easily obtainable values. New method for calculating acoustic-liner admittance in rectangular duct with grazing flow based on finite-element discretization of acoustic field and reposing of unknown admittance value as linear eigenvalue problem on admittance value. Problem solved by Gaussian elimination. Unlike existing methods, present method extendable to mean flows with two-dimensional boundary layers as well. In presence of shear, results of method compared well with results of Runge-Kutta integration technique.

Watson, W. R.↗

NASA NTRS1993

Quantum noise and squeezing in optical parametric oscillator with arbitrary output coupling

The redistribution of intrinsic quantum noise in the quadratures of the field generated in a sub-threshold degenerate optical parametric oscillator exhibits interesting dependences on the individual output mirror transmittances, when they are included exactly. We present a physical picture of this problem, based on mirror boundary conditions, which is valid for arbitrary transmittances. Hence, our picture applies uniformly to all values of the cavity Q factor representing, in the opposite extremes, both perfect oscillator and amplifier configurations. Beginning with a classical second-harmonic pump, we shall generalize our analysis to the finite amplitude and phase fluctuations of the pump.

Prasad, Sudhakar↗

NASA NTRS2010

Does the Hall Effect Solve the Flux Pileup Saturation Problem?

It is well known that magnetic flux pileup can significantly speed up the rate of magnetic reconnection in high Lundquist number resistive MHD,allowing reconnection to proceed at a rate which is insensitive to the plasma resistivity over a wide range of Lundquist number. Hence, pileup is a possible solution to the Sweet-Parker time scale problem. Unfortunately, pileup tends to saturate above a critical value of the Lundquist number, S_c, where the value ofS_c depends on initial and boundary conditions, with Sweet-Parker scaling returning above S_c. It has been argued (see Dorelli and Bim [2003] and Dorelli [2003]) that the Hall effect can allow flux pileup to saturate (when the scale of the current sheet approaches ion inertial scale, di) before the reconnection rate begins to stall. However, the resulting saturated reconnection rate, while insensitive to the plasma resistivity, was found to depend strongly on the di. In this presentation, we revisit the problem of magnetic island coalescence (which is a well known example of flux pileup reconnection), addressing the dependence of the maximum coalescence rate on the ratio of di in the "large island" limit in which the following inequality is always satisfied: l_eta di lambda, where I_eta is the resistive diffusion length and lambda is the island wavelength.

Dorelli, John C.↗

NASA NTRS1998

Exact Solution to Stationary Onset of Convection Due to Surface Tension Variation in a Multicomponent Fluid Layer With Interfacial Deformation

Stationary onset of convection due to surface tension variation in an unbounded multicomponent fluid layer is considered. Surface deformation is included and general flux boundary conditions are imposed on the stratifying agencies (temperature/composition) disturbance equations. Exact solutions are obtained to the general N-component problem for both finite and infinitesimal wavenumbers. Long wavelength instability may coexist with a finite wavelength instability for certain sets of parameter values, often referred to as frontier points. For an impermeable/insulated upper boundary and a permeable/conductive lower boundary, frontier boundaries are computed in the space of Bond number, Bo, versus Crispation number, Cr, over the range 5 x 10(exp -7) less than or equal to Bo less than or equal to 1. The loci of frontier points in (Bo, Cr) space for different values of N, diffusivity ratios, and, Marangoni numbers, collapsed to a single curve in (Bo, D(dimensional variable)Cr) space, where D(dimensional variable) is a Marangoni number weighted diffusivity ratio.

Skarda, J. Raymond Lee↗

NASA NTRS1972

Numerical method for the solution of large systems of differential equations of the boundary layer type

A numerical method for the solution of large systems of nonlinear differential equations of the boundary-layer type is described. The method is a modification of the technique for satisfying asymptotic boundary conditions. The present method employs inverse interpolation instead of the Newton method to adjust the initial conditions of the related initial-value problem. This eliminates the so-called perturbation equations. The elimination of the perturbation equations not only reduces the user's preliminary work in the application of the method, but also reduces the number of time-consuming initial-value problems to be numerically solved at each iteration. For further ease of application, the solution of the overdetermined system for the unknown initial conditions is obtained automatically by applying Golub's linear least-squares algorithm. The relative ease of application of the proposed numerical method increases directly as the order of the differential-equation system increases. Hence, the method is especially attractive for the solution of large-order systems. After the method is described, it is applied to a fifth-order problem from boundary-layer theory.

Green, M. J.↗

NASA NTRS1995

Wave propagation and scattering in computational aeroacoustics

This paper describes predictions of model problems in computational aeroacoustics. Two problem classes are considered: The first (Workshop Category 3, Problem 1 and 2) considers two-dimensional wave propagation and non-reflecting boundary conditions in the presence of a mean flow. The second (Workshop Category 4, Problem 1) examines wall boundary conditions. For the last problem we introduce the Impedance Mismatch Method (IMM) to treat the solid wall boundaries. In this method the solid wall is simulated using a wall region in which the characteristic impedance is set to a different value to that in the fluid region. This method has advantages over traditional solid wall boundary conditions including simplicity of coding, speed of computations, and the ability to treat curved boundaries efficiently. Several numerical examples are given in addition to the Workshop Problems. The discretization of the Euler equations is performed in all cases with a Dispersion-Relation-Preserving (DRP) algorithm. The numerical results are compared with either analytical solutions or solutions obtained using traditional solid wall boundary conditions.

Chung, Cathy↗

NASA NTRS1993

A direct approach to finding unknown boundary conditions in steady heat conduction

The capability of the boundary element method (BEM) in determining thermal boundary conditions on surfaces of a conducting solid where such quantities are unknown was demonstrated. The method uses a non-iterative direct approach in solving what is usually called the inverse heat conduction problem (IHCP). Given any over-specified thermal boundary conditions such as a combination of temperature and heat flux on a surface where such data is readily available, the algorithm computes the temperature field within the object and any unknown thermal boundary conditions on surfaces where thermal boundary values are unavailable. A two-dimensional, steady-state BEM program was developed and was tested on several simple geometries where the analytic solution was known. Results obtained with the BEM were in excellent agreement with the analytic values. The algorithm is highly flexible in treating complex geometries, mixed thermal boundary conditions, and temperature-dependent material properties and is presently being extended to three-dimensional and unsteady heat conduction problems. The accuracy and reliability of this technique was very good but tended to deteriorate when the known surface conditions were only slightly over-specified and far from the inaccessible surface.

Martin, Thomas J.↗

NASA NTRS1980

A nonreflecting outflow boundary condition for subsonic Navier-Stokes calculations

A nonreflecting boundary condition is presented for the numerical solution of the time-dependent compressible Navier-Stokes equations when these equations are used to obtain a steady state. This boundary condition is shown to be effective in reducing reflections at a subsonic outflow boundary. Numerical calculations using a model problem were made to compare this boundary condition with other outflow boundary conditions. The nonreflecting boundary condition contains a parameter whose optimal value is estimated using the analysis of a simplified set of equations.

Rudy, D. H.↗

NASA NTRS1990

On the instability of hypersonic flow past a flat plate

The instability of hypersonic boundary-layer flows over flat plates is considered. The viscosity of the fluid is taken to be governed by Sutherland's law, which gives a much more accurate representation of the temperature dependence of fluid viscosity at hypersonic speeds than Chapman's approximate linear law; although at lower speeds the temperature variation of the mean state is less pronounced so that the Chapman law can be used with some confidence. Attention is focussed on the so-called (vorticity) mode of instability of the viscous hypersonic boundary layer. This is thought to be the fastest growing inviscid disturbance at hypersonic speeds; it is also believed to have an asymptotically larger growth rate than any viscous or centrifugal instability. As a starting point the instability of the hypersonic boundary layer which exists far downstream from the leading edge of the plate is investigated. In this regime the shock that is attached to the leading edge of the plate plays no role, so that the basic boundary layer is non-interactive. It is shown that the vorticity mode of instability of this flow operates on a significantly different lengthscale than that obtained if a Chapman viscosity law is assumed. In particular, it is found that the growth rate predicted by a linear viscosity law overestimates the size of the growth rate by O(M(exp 2). Next, the development of the vorticity mode as the wavenumber decreases is described, and it is shown that acoustic modes emerge when the wavenumber has decreased from it's O(1) initial value to O(M (exp -3/2). Finally, the inviscid instability of the boundary layer near the leading edge in the interaction zone is discussed and particular attention is focussed on the strong interaction region which occurs sufficiently close to the leading edge. It is found that the vorticity mode in this regime is again unstable, and that it is concentrated in the transition layer at the edge of the boundary layer where the temperature adjusts from its large, O(M(exp 2), value in the viscous boundary layer, to its O(1) free stream value. The existence of the shock indirectly, but significantly, influences the instability problem by modifying the basic flow structure in this layer.

Blackaby, Nicholas↗

NASA NTRS2005

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

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.↗