Search NASA⌕ Search

SEARCH · Search NASA

Results for “BOUNDARY VALUE”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

Asymptotic solution of the problem for a thin axisymmetric cavern

The boundary value problem which describes the axisymmetric separation of the flow around a body by a stationary infinite stream is considered. It is understood that the cavitation number varies over the length of the cavern. Using the asymptotic expansions for the potential of a thin body, the orders of magnitude of terms in the equations of the problem are estimated. Neglecting small quantities, a simplified boundary value problem is obtained.

Serebriakov, V. V.↗
NASA NTRS1984

Shape design sensitivity analysis of built-up structures

Selection of the best shape of a fillet in a tension bar such that no yielding occurs has long attracted the attention of engineers. Dimensions and notations for the bar and fillet are shown. With symmetry, only the upper half of the bar is considered. The optimal design problem is to find a boundary shape to minimize the total area of the fillet such that no yielding occurs. The classical boundary value problem is reduced to a variational or energy related problem which not only has excellent properties of existence and uniqueness but also provides the mathematical foundation for finite element analysis. The variational formulation may be viewed as the principle of virtual work and the finite element method as a application of the Galerkin method to the variational equation for approximate solution of the boundary value problem.

Choi, K. K.↗
NASA NTRS2010

Determining Shear Stress Distribution in a Laminate

A "simplified shear solution" method approximates the through-thickness shear stress distribution within a composite laminate based on an extension of laminated beam theory. The method does not consider the solution of a particular boundary value problem; rather, it requires only knowledge of the global shear loading, geometry, and material properties of the laminate or panel. It is thus analogous to lamination theory in that ply-level stresses can be efficiently determined from global load resultants at a given location in a structure and used to evaluate the margin of safety on a ply-by-ply basis. The simplified shear solution stress distribution is zero at free surfaces, continuous at ply boundaries, and integrates to the applied shear load. The method has been incorporated within the HyperSizer commercial structural sizing software to improve its predictive capability for designing composite structures. The HyperSizer structural sizing software is used extensively by NASA to design composite structures. In the case of through-thickness shear loading on panels, HyperSizer previously included a basic, industry-standard, method for approximating the resulting shear stress distribution in sandwich panels. However, no such method was employed for solid laminate panels. The purpose of the innovation is to provide an approximation of the through-thickness shear stresses in a solid laminate given the through-thickness shear loads (Qx and Qy) on the panel. The method was needed for implementation within the HyperSizer structural sizing software so that the approximated ply-level shear stresses could be utilized in a failure theory to assess the adequacy of a panel design. The simplified shear solution method was developed based on extending and generalizing bi-material beam theory to plate-like structures. It is assumed that the through-thickness shear stresses arise due to local bending of the laminate induced by the through-thickness shear load, and by imposing equilibrium both vertically and horizontally, the through-thickness shear stress distribution can be calculated. The resulting shear stresses integrate to the applied shear load, are continuous at the ply interfaces, and are zero at the laminate-free surfaces. If both Qx and Qy shear loads are present, it is assumed that they act independently and that their effects can be superposed. The calculated shear stresses can be rotated within each ply to the principal material coordinates for use in a ply-level failure criterion. The novelty of the simplified shear solution method is its simplicity and the fact that it does not require solution of a particular boundary value problem. The advantages of the innovation are that an approximation of the though-thickness shear stress distribution can be quickly determined for any solid laminate or solid laminate region within a stiffened panel.

Bednarcyk, Brett A.↗
NASA NTRS1992

Computational aspects of helicopter trim analysis and damping levels from Floquet theory

Helicopter trim settings of periodic initial state and control inputs are investigated for convergence of Newton iteration in computing the settings sequentially and in parallel. The trim analysis uses a shooting method and a weak version of two temporal finite element methods with displacement formulation and with mixed formulation of displacements and momenta. These three methods broadly represent two main approaches of trim analysis: adaptation of initial-value and finite element boundary-value codes to periodic boundary conditions, particularly for unstable and marginally stable systems. In each method, both the sequential and in-parallel schemes are used and the resulting nonlinear algebraic equations are solved by damped Newton iteration with an optimally selected damping parameter. The impact of damped Newton iteration, including earlier-observed divergence problems in trim analysis, is demonstrated by the maximum condition number of the Jacobian matrices of the iterative scheme and by virtual elimination of divergence. The advantages of the in-parallel scheme over the conventional sequential scheme are also demonstrated.

Gaonkar, Gopal H.↗
NASA NTRS1993

Helicopter trim analysis by shooting and finite element methods with optimally damped Newton iterations

Helicopter trim settings of periodic initial state and control inputs are investigated for convergence of Newton iteration in computing the settings sequentially and in parallel. The trim analysis uses a shooting method and a weak version of two temporal finite element methods with displacement formulation and with mixed formulation of displacements and momenta. These three methods broadly represent two main approaches of trim analysis: adaptation of initial-value and finite element boundary-value codes to periodic boundary conditions, particularly for unstable and marginally stable systems. In each method, both the sequential and in-parallel schemes are used, and the resulting nonlinear algebraic equations are solved by damped Newton iteration with an optimally selected damping parameter. The impact of damped Newton iteration, including earlier-observed divergence problems in trim analysis, is demonstrated by the maximum condition number of the Jacobian matrices of the iterative scheme and by virtual elimination of divergence. The advantages of the in-parallel scheme over the conventional sequential scheme are also demonstrated.

Achar, N. S.↗
NASA NTRS2003

Time as an Observable in Nonrelativistic Quantum Mechanics

The argument follows from the viewpoint that quantum mechanics is taken not in the usual form involving vectors and linear operators in Hilbert spaces, but as a boundary value problem for a special class of partial differential equations-in the present work, the nonrelativistic Schrodinger equation for motion of a structureless particle in four- dimensional space-time in the presence of a potential energy distribution that can be time-as well as space-dependent. The domain of interest is taken to be one of two semi-infinite boxes, one bounded by two t=constant planes and the other by two t=constant planes. Each gives rise to a characteristic boundary value problem: one in which the initial, input values on one t=constant wall are given, with zero asymptotic wavefunction values in all spatial directions, the output being the values on the second t=constant wall; the second with certain input values given on both z=constant walls, with zero asymptotic values in all directions involving time and the other spatial coordinates, the output being the complementary values on the z=constant walls. The first problem corresponds to ordinary quantum mechanics; the second, to a fully time-dependent version of a problem normally considered only for the steady state (time-independent Schrodinger equation). The second problem is formulated in detail. A conserved indefinite metric is associated with space-like propagation, where the sign of the norm of a unidirectional state corresponds to its spatial direction of travel.

Hahne, G. E.↗
NASA NTRS1994

Linear and nonlinear dynamic analysis of redundant load path bearingless rotor systems

The goal of this research is to develop the transfer matrix method to treat nonlinear autonomous boundary value problems with multiple branches. The application is the complete nonlinear aeroelastic analysis of multiple-branched rotor blades. Once the development is complete, it can be incorporated into the existing transfer matrix analyses. There are several difficulties to be overcome in reaching this objective. The conventional transfer matrix method is limited in that it is applicable only to linear branch chain-like structures, but consideration of multiple branch modeling is important for bearingless rotors. Also, hingeless and bearingless rotor blade dynamic characteristics (particularly their aeroelasticity problems) are inherently nonlinear. The nonlinear equations of motion and the multiple-branched boundary value problem are treated together using a direct transfer matrix method. First, the formulation is applied to a nonlinear single-branch blade to validate the nonlinear portion of the formulation. The nonlinear system of equations is iteratively solved using a form of Newton-Raphson iteration scheme developed for differential equations of continuous systems. The formulation is then applied to determine the nonlinear steady state trim and aeroelastic stability of a rotor blade in hover with two branches at the root. A comprehensive computer program is developed and is used to obtain numerical results for the (1) free vibration, (2) nonlinearly deformed steady state, (3) free vibration about the nonlinearly deformed steady state, and (4) aeroelastic stability tasks. The numerical results obtained by the present method agree with results from other methods.

Murthy, V. R.↗
NASA NTRS1990

Global pressure relaxation for laminar two-dimensional internal flow

This study extends the reduced Navier-Stokes (RNS) global pressure relaxation procedure developed by Rubin and co-workers for external flow to internal flow applications. The streamwise pressure gradient is split into a backward-differenced or initial value component, as in boundary layer marching, and a forward-differenced or boundary value component that represents the elliptic downstream effects. The streamwise convection terms are upwind-differenced and all other streamwise derivatives are backward-differenced. A standard boundary layer marching technique imbedded in a conventional line relaxation technique is obtained. For compressible flow the pressure iteration determines the interior flow interaction as well as the inlet mass flux that is consistent with the outflow pressure boundary condition. Results have been computed for incompressible flow in both rectangular and curved channels, and for subsonic compressible flow in the simulation of an aerofoil in a wind tunnel. Converged solutions were obtained over a range of Reynolds numbers generating small to moderately large separation bubbles.

Rosenbaum, D.↗
NASA NTRS1974

User's guide to four-body and three-body trajectory optimization programs

A collection of computer programs and subroutines written in FORTRAN to calculate 4-body (sun-earth-moon-space) and 3-body (earth-moon-space) optimal trajectories is presented. The programs incorporate a variable step integration technique and a quadrature formula to correct single step errors. The programs provide capability to solve initial value problem, two point boundary value problem of a transfer from a given initial position to a given final position in fixed time, optimal 2-impulse transfer from an earth parking orbit of given inclination to a given final position and velocity in fixed time and optimal 3-impulse transfer from a given position to a given final position and velocity in fixed time.

Pu, C. L.↗
NASA NTRS1988

Solving time-dependent two-dimensional eddy current problems

Results of transient eddy current calculations are reported. For simplicity, a two-dimensional transverse magnetic field which is incident on an infinitely long conductor is considered. The conductor is assumed to be a good but not perfect conductor. The resulting problem is an interface initial boundary value problem with the boundary of the conductor being the interface. A finite difference method is used to march the solution explicitly in time. The method is shown. Treatment of appropriate radiation conditions is given special consideration. Results are validated with approximate analytic solutions. Two stringent test cases of high and low frequency incident waves are considered to validate the results.

Lee, Min Eig↗
NASA NTRS1975

The determination of orbits using Picard iteration

The determination of orbits by using Picard iteration is reported. This is a direct extension of the classical method of Picard that has been used in finding approximate solutions of nonlinear differential equations for a variety of problems. The application of the Picard method of successive approximations to the initial value and the two point boundary value problems is given.

Mikkilineni, R. P.↗
NASA NTRS1984

Navier-Stokes calculations for unsteady three-dimensional vortical flows in unbounded domains

Finite-difference Navier-Stokes calculations for unsteady, three-dimensional, incompressible, viscous flows induced by initial vorticity distributions are presented and discussed in this paper. The initial vorticity distributions are assumed to be embedded in a flow field of infinite extent that is quiescent at infinity. These vorticity distributions are typical of vortex rings and other closed vortical tubes or structures. Such structures are important elements in fluid flows such as jets, atmospheric convection and the far-field wakes of aircraft; studies of their interaction may aid in an understanding of complex fluid flows. The calculations employ a method recently proposed by Ting to approximate the infinite-domain boundary value problem with a finite boundary computational domain, and this method is shown to yield accurate three-dimensional results for reasonable expenditures of computer time. Because of the efficiency of the boundary condition technique and the resulting Navier-Stokes code, a 16-bit minicomputer with virtual memory was capable of performing the calculations for the unsteady motion of two obliquely colliding vortex rings. The results of these calculations are presented in the paper.

Chamberlain, J. P.↗