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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Multigrid-sinc methods

A Galerkin method using Whittaker cardinal or 'sinc' functions as basis functions is described for the solution of boundary-value problems. When the solution is analytic in the interior of the domain, the error of approximation using 2N + 1 points is O(e exp /-gamma sq rt N/) even if derivatives of the solution are singular at the boundaries. A multigrid method with overall complexity O(N log N) is used to solve the discrete equations. This paper contains a description of the multigrid-sinc algorithm along with some preliminary numerical results for two-point boundary-value problems.

Schaffer, Steve↗
NASA NTRSDate not supplied

Extending Explicit Guidance Methods to Higher Dimensions, Additional Conditions, and Higher Order Integration

Guidance functions play critical roles in autonomy to steer vehicles and aircraft to the intended target or destination. Explicit guidance (E Guidance) solves the two-point boundary value problem with initial and final conditions for position and velocity. The original formulation of E Guidance involves translational acceleration commands with a direct relationship to time, and it is possible to modify E Guidance for rotational acceleration. Other extensions for E Guidance include higher dimensions, additional conditions, and higher-order integration of the linearly independent E Guidance functions. The most promising extension involves higher-order integration of the E Guidance functions, but it may be physically impractical by initially moving away from the target. This paper provides a brief overview of some methods that extend E Guidance to higher dimensions, utilize additional conditions, or perform higher-order integration, and if they satisfy the two-point boundary value problem.

explicit guidance↗
NASA NTRSDate not supplied

Extending Explicit Guidance Methods to Higher Dimensions, Additional Conditions, and Higher Order Integration

Guidance functions play critical roles in autonomy to steer vehicles and aircraft to the intended target or destination. Explicit guidance (E Guidance) solves the two-point boundary value problem with initial and final conditions for position and velocity. The original formulation of explicit guidance involves translational acceleration commands with a direct relationship to time, and it is possible to modify E Guidance for rotational acceleration. Other extensions for E Guidance include higher dimensions, additional conditions, and higher-order integration of the linearly independent E Guidance functions. The most promising extension involves higher-order integration of the E Guidance functions, but it may be physically impractical by initially moving away from the target. This paper provides a brief overview of some methods that extend E Guidance to higher dimensions, utilize additional conditions, or perform higher-order integration, and if they satisfy the two-point boundary value problem.

explicit guidance↗
DOE OSTI2021

Source term method for binary neutron stars initial data

The initial condition problem for a binary neutron star system requires a Poisson equation solver for the velocity potential with a Neumann-like boundary condition on the surface of the star. Difficulties that arise in this boundary value problem are: (a) the boundary is not known a priori, but constitutes part of the solution of the problem; (b) various terms become singular at the boundary. In this work, we present a new method to solve the fluid Poisson equation for irrotational/spinning binary neutron stars. The advantage of the new method is that it does not require complex fluid surface fitted coordinates and it can be implemented in a Cartesian grid, which is a standard choice in numerical relativity calculations. This is accomplished by employing the source term method proposed by Towers, where the boundary condition is treated as a jump condition and is incorporated as additional source terms in the Poisson equation, which is then solved iteratively. The issue of singular terms caused by vanishing density on the surface is resolved with an additional separation that shifts the computation boundary to the interior of the star. We present two-dimensional tests to show the convergence of the source term method, and we further apply this solver to a realistic three-dimensional binary neutron star problem. By comparing our solution with the one coming from the initial data solver cocal, we demonstrate agreement to approximately 1%. We report our method can be used in other problems with non-smooth solutions like in magnetized neutron stars.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1981

The choice of numerical boundary conditions for hyperbolic systems

Two fundamental problems for mixed initial boundary value problems with applications in fluid mechanics are discussed. First, different stability properties are discussed which are of importance for long time integrations and steady state calculations. Secondly, a new numerical technique for problems with an artificial boundary is introduced.

Gustafson, B.↗
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.↗
DOE OSTI2018

Spectrum and Stabilization in Hyperbolic Problems

We study the connection between the stabilization of solutions of a mixed hyperbolic problem and spectral properties of the corresponding elliptic boundary value problem. We consider the first mixed problem for the wave equation in bounded and unbounded domains in ℝ{sup n}, determine the class of its energy solutions, and represent the solutions in terms of the Bochner–Stieltjes integral. We study how the spectrum of the elliptic operator affects the behavior of local energy of a solution and describe a method which allows us to study the stabilization of solutions with the help of estimates in the spectral parameter for solutions of the stationary problem on the upper half-plane.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1983

Finite difference numerical methods for boundary control problems governed by hyperbolic partial differential equations

This paper briefly reviews convergent finite difference schemes for hyperbolic initial boundary value problems and their applications to boundary control systems of hyperbolic type which arise in the modelling of vibrations. These difference schemes are combined with the primal and the dual approaches to compute the optimal control in the unconstrained case, as well as the case when the control is subject to inequality constraints. Some of the preliminary numerical results are also presented.

Chen, G.↗
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 NTRS1973

Resonant oscillations of intermediate frequency in a stratified atmosphere.

A class of solutions to a model of forced oscillations in a rotating stratified atmospheric layer is derived and analyzed. The basic model is found to reduce to a boundary value problem with a second-order linear partial differential equation of the hyperbolic type for this range of forcing frequencies. The forced solutions are shown to exhibit resonances with the normal modes of oscillation of the layer. The characteristics of the resonant modes are analyzed in terms of mean tropospheric values of temperature, temperature lapse, wind speed, horizontal and vertical wind shears, latitude, and the frequency and horizontal wavelength of the forcing mechanism. These solutions are compared with solutions to the model for a different (subinertial) range of forcing frequencies. This comparison leads to an elliptic boundary value problem. The solutions in that case do not exhibit the same type of resonance and generally decay away from the region of forcing.

Thompson, O. E.↗
NASA NTRS1984

Stability of the laminar boundary layer in a streamwise corner

The stability of viscous, incompressible flow along a streamwise corner, often called the corner boundary layer problem is examined. The semi-infinite boundary value problem satisfied by small amplitude disturbances in the "bending boundary layer' region is obtained. The mean secondary flow induced by the corner exhibits a flow reversal in this region. Uniformly valid "first approximations' to solutions of the governing differential equations are derived. Uniformity at infinity is achieved by a suitable choice of the large parameter and use of an approximate Langer variable. Approximations to solutions of balanced type have a phase shift across the critical layer which is associated with instabilities in the case of two dimensional boundary layer profiles.

Lakin, W. D.↗
NASA NTRS1984

Stability of the laminar boundary layer in a streamwise corner

The stability of viscous, incompressible flow along a streamwise corner, often called the corner boundary layer problem is examined. The semi-infinite boundary value problem satisfied by small amplitude disturbances in the 'bending boundary layer' region is obtained. The mean secondary flow induced by the corner exhibits a flow reversal in this region. Uniformly valid 'first approximations' to solutions of the governing differential equations are derived. Uniformity at infinity is achieved by a suitable choice of the large parameter and use of an approximate Langer variable. Approximations to solutions of balanced type have a phase shift across the critical layer which is associated with instabilities in the case of two dimensional boundary layer profiles. Previously announced in STAR as N84-17532

Lakin, W. D.↗
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.↗