Engineering Papers⌕ Search

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.

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

At least 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Modal and non-modal baroclinic waves

Solution of the initial-value problem for the Eady model is presented. In the presence of boundaries, normal mode waves as well as non-modal waves exist. Energy extracted from the mean flow during the initial development of a perturbation is found to excite the persistent normal modes. It is suggested that this process may be important to cyclogenesis and in providing energy to neutral or near-neutral normal modes. In particular, the Petterssen criterion for cyclogenesis is clarified.

Farrell, B.↗
NASA NTRS1994

Pseudospectral collocation methods for fourth order differential equations

Collocation schemes are presented for solving linear fourth order differential equations in one and two dimensions. The variational formulation of the model fourth order problem is discretized by approximating the integrals by a Gaussian quadrature rule generalized to include the values of the derivative of the integrand at the boundary points. Collocation schemes are derived which are equivalent to this discrete variational problem. An efficient preconditioner based on a low-order finite difference approximation to the same differential operator is presented. The corresponding multidomain problem is also considered and interface conditions are derived. Pseudospectral approximations which are C1 continuous at the interfaces are used in each subdomain to approximate the solution. The approximations are also shown to be C3 continuous at the interfaces asymptotically. A complete analysis of the collocation scheme for the multidomain problem is provided. The extension of the method to the biharmonic equation in two dimensions is discussed and results are presented for a problem defined in a nonrectangular domain.

Malek, Alaeddin↗
NASA NTRS1973

Application of matched asymptotic expansions to lunar and interplanetary trajectories. Volume 2: Derivations of second-order asymptotic boundary value solutions

Previously published asymptotic solutions for lunar and interplanetery trajectories have been modified and combined to formulate a general analytical solution to the problem of N-bodies. The earlier first-order solutions, derived by the method of matched asymptotic expansions, have been extended to second order for the purpose of obtaining increased accuracy. The complete derivation of the second-order solution, including the application of a regorous matching principle, is given. It is shown that the outer and inner expansions can be matched in a region of order mu to the alpha power, where 2/5 alpha 1/2, and mu (the moon/earth or planet/sun mass ratio) is much less than one. The second-order asymptotic solution has been used as a basis for formulating a number of analytical two-point boundary value solutions. These include earth-to-moon, one- and two-impulse moon-to-Earth, and interplanetary solutions. Each is presented as an explicit analytical solution which does not require iterative steps to satisfy the boundary conditions. The complete derivation of each solution is shown, as well as instructions for numerical evaluation. For Vol. 1, see N73-27738.

Lancaster, J. E.↗
NASA NTRS1987

Goertler instability on an airfoil

An effective computational scheme was developed to study the growth/damping of Goertler vortices along walls of variable curvature. Computational experiments indicate that when the amplification rates for the u-, v-, and w-perturbations are the same, the finite difference approach to solve the initial value problem and the normal mode approach give identical results for the Blasius boundary layer on constant curvature concave walls. The growth of Goertler vortices was rapid in the concave regions and was followed by sharp damping in the convex region. However, multiple sets of counter-rotating vortices were formed and remained far downstream in the convex region. The current computational scheme can be easily extended to more realistic problems including variable pressure gradients and suction effects.

Kalburgi, Vijay↗
NASA NTRS1951

Analysis of plane-plastic stress problems with axial symmetry in strain-hardening range

A simple method is developed for solving plane-plastic-stress problems with axial symmetry in the strain-hardening range which is based on the deformation theory of plasticity employing the finite-strain concept. The equations defining the problems are first reduced to two simultaneous nonlinear differential equations involving two dependent variables: (a) the octahedral shear strain, and (b) a parameter indicating the ratio of principal stresses. By multiplying the load and dividing the radius by an arbitrary constant, it is possible to solve these problems without iteration for any value of the modified load. The constant is determined by the boundary condition. This method is applied to a circular membrane under pressure, a rotating disk without and with a central hole, and an infinite plate with a circular hole. Two materials, inconel x and 16-25-6, the octahedral shear stress-strain relations of which do not follow the power law, are used. Distributions of octahedral shear strain, as well as of principal stresses and strains, are obtained. These results are compared with the results of the same problems in the elastic range.

Wu, M H Lee↗
NASA NTRS1958

Drag Minimization for Wings and Bodies in Supersonic Flow

The minimization of inviscid fluid drag is studied for aerodynamic shapes satisfying the conditions of linearized theory, and subject to imposed constraints on lift, pitching moment, base area, or volume. The problem is transformed to one of determining two-dimensional potential flows satisfying either Laplace's or Poisson's equations with boundary values fixed by the imposed conditions. A general method for determining integral relations between perturbation velocity components is developed. This analysis is not restricted in application to optimum cases; it may be used for any supersonic wing problem.

Heaslet, Max A↗
NASA NTRS2000

Direct Numerical Simulation of Automobile Cavity Tones

The Navier Stokes equation is solved computationally by the Dispersion-Relation-Preserving (DRP) scheme for the flow and acoustic fields associated with a laminar boundary layer flow over an automobile door cavity. In this work, the flow Reynolds number is restricted to R(sub delta*) < 3400; the range of Reynolds number for which laminar flow may be maintained. This investigation focuses on two aspects of the problem, namely, the effect of boundary layer thickness on the cavity tone frequency and intensity and the effect of the size of the computation domain on the accuracy of the numerical simulation. It is found that the tone frequency decreases with an increase in boundary layer thickness. When the boundary layer is thicker than a certain critical value, depending on the flow speed, no tone is emitted by the cavity. Computationally, solutions of aeroacoustics problems are known to be sensitive to the size of the computation domain. Numerical experiments indicate that the use of a small domain could result in normal mode type acoustic oscillations in the entire computation domain leading to an increase in tone frequency and intensity. When the computation domain is expanded so that the boundaries are at least one wavelength away from the noise source, the computed tone frequency and intensity are found to be computation domain size independent.

Kurbatskii, Konstantin↗
NASA NTRS1992

Multiple scattering technique lidar

The Bernouilli-Ricatti equation is based on the single scattering description of the lidar backscatter return. In practice, especially in low visibility conditions, the effects of multiple scattering can be significant. Instead of considering these multiple scattering effects as a nuisance, we propose here to use them to help resolve the problems of having to assume a backscatter-to-extinction relation and specifying a boundary value for a position far remote from the lidar station. To this end, we have built a four-field-of-view lidar receiver to measure the multiple scattering contributions. The system has been described in a number of publications that also discuss preliminary results illustrating the multiple scattering effects for various environmental conditions. Reported here are recent advances made in the development of a method of inverting the multiple scattering data for the determination of the aerosol scattering coefficient.

Bissonnette, Luc R.↗
NASA NTRS1974

Interaction between the atmospheric and oceanic boundary layers

The two-layer system of an atmosphere over water bodies is reduced to a single-layer problem. Values of the interfacial quantities, such as the friction velocity, the surface velocity, the angles, alpha and beta, between the surface shear stress and the geostrophic wind velocity and the surface wind velocity, respectively, and the surface roughness, all of which depend upon external parameters, such as the geostrophic wind and stratifications, are obtained. The geostrophic drag coefficient, the geostrophic wind coefficient, and the angles alpha, and beta, of the turbulent flow at the sea-air interface are functions of a dimensionless number, mfG/kg, with S sub 1 and S sub 2 as two free stratification parameters. The surface roughness is uniquely determined from the geostrophic wind rather than from the wind profile in the boundary layer.

Yeh, G.-T.↗
NASA NTRS1984

Multiple steady states for characteristic initial value problems

The time dependent, isentropic, quasi-one-dimensional equations of gas dynamics and other model equations are considered under the constraint of characteristic boundary conditions. Analysis of the time evolution shows how different initial data may lead to different steady states and how seemingly anamolous behavior of the solution may be resolved. Numerical experimentation using time consistent explicit algorithms verifies the conclusions of the analysis. The use of implicit schemes with very large time steps leads to erroneous results.

Salas, M. D.↗
NASA NTRS1991

Measurement and analysis of grain boundary grooving by volume diffusion

Experimental measurements of isothermal grain boundary grooving by volume diffusion are carried out for Sn bicrystals in the Sn-Pb system near the eutectic temperature. The dimensions of the groove increase with a temporal exponent of 1/3, and measurement of the associated rate constant allows the determination of the product of the liquid diffusion coefficient D and the capillarity length Gamma associated with the interfacial free energy of the crystal-melt interface. The small-slope theory of Mullins is generalized to the entire range of dihedral angles by using a boundary integral formulation of the associated free boundary problem, and excellent agreement with experimental groove shapes is obtained. By using the diffusivity measured by Jordon and Hunt, the present measured values of Gamma are found to agree to within 5 percent with the values obtained from experiments by Gunduz and Hunt on grain boundary grooving in a temperature gradient.

Hardy, S. C.↗
NASA NTRS1997

Low-Dispersion Scheme for Nonlinear Acoustic Waves in Nonuniform Flow

The linear dispersion-relation-preserving scheme and its boundary conditions have been extended to the nonlinear Euler equations. This allowed computing, a nonuniform flowfield and a nonlinear acoustic wave propagation in such a medium, by the same scheme. By casting all the equations, boundary conditions, and the solution scheme in generalized curvilinear coordinates, the solutions were made possible for non-Cartesian domains and, for the better deployment of the grid points, nonuniform grid step sizes could be used. It has been tested for a number of simple initial-value and periodic-source problems. A simple demonstration of the difference between a linear and nonlinear propagation was conducted. The wall boundary condition, derived from the momentum equations and implemented through a pressure at a ghost point, and the radiation boundary condition, derived from the asymptotic solution to the Euler equations, have proven to be effective for the nonlinear equations and nonuniform flows. The nonreflective characteristic boundary conditions also have shown success but limited to the nonlinear waves in no mean flow, and failed for nonlinear waves in nonuniform flow.

Baysal, Oktay↗
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↗