Search NASA⌕ Search

SEARCH · Search NASA

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

An Application of the Difference Potentials Method to Solving External Problems in CFD

Numerical solution of infinite-domain boundary-value problems requires some special techniques that would make the problem available for treatment on the computer. Indeed, the problem must be discretized in a way that the computer operates with only finite amount of information. Therefore, the original infinite-domain formulation must be altered and/or augmented so that on one hand the solution is not changed (or changed slightly) and on the other hand the finite discrete formulation becomes available. One widely used approach to constructing such discretizations consists of truncating the unbounded original domain and then setting the artificial boundary conditions (ABC's) at the newly formed external boundary. The role of the ABC's is to close the truncated problem and at the same time to ensure that the solution found inside the finite computational domain would be maximally close to (in the ideal case, exactly the same as) the corresponding fragment of the original infinite-domain solution. Let us emphasize that the proper treatment of artificial boundaries may have a profound impact on the overall quality and performance of numerical algorithms. The latter statement is corroborated by the numerous computational experiments and especially concerns the area of CFD, in which external problems present a wide class of practically important formulations. In this paper, we review some work that has been done over the recent years on constructing highly accurate nonlocal ABC's for calculation of compressible external flows. The approach is based on implementation of the generalized potentials and pseudodifferential boundary projection operators analogous to those proposed first by Calderon. The difference potentials method (DPM) by Ryaben'kii is used for the effective computation of the generalized potentials and projections. The resulting ABC's clearly outperform the existing methods from the standpoints of accuracy and robustness, in many cases noticeably speed up the multigrid convergence, and at the same time are quite comparable to other methods from the standpoints of geometric universality and simplicity of implementation.

Ryaben 'Kii, Victor S.↗
NASA NTRS1974

Control optimization of a lifting body entry problem by an improved and a modified method of perturbation function

A study of the solution problem of a complex entry optimization was studied. The problem was transformed into a two-point boundary value problem by using classical calculus of variation methods. Two perturbation methods were devised. These methods attempted to desensitize the contingency of the solution of this type of problem on the required initial co-state estimates. Also numerical results are presented for the optimal solution resulting from a number of different initial co-states estimates. The perturbation methods were compared. It is found that they are an improvement over existing methods.

Garcia, F., Jr.↗
NASA NTRS1978

Excitation of instability waves in a two-dimensional shear layer by sound

The excitation of instability waves in a plane compressible shear layer by sound waves is studied. The problem is formulated mathematically as an inhomogeneous boundary-value problem. A general solution for abitrary incident sound wave is found by first constructing the Green's function of the problem. Numerical values of the coupling constants between incident sound waves and excited instability waves for a range of flow Mach number are calculated. The effect of the angle of incidence in the case of a beam of acoustic waves is analyzed. It is found that for moderate subsonic Mach numbers a narrow beam aiming at an angle between 50 to 80 deg to the flow direction is most effective in exciting instability waves.

Tam, C. K. W.↗
NASA NTRS1980

Vertical slender jets

The shape of a vertical slender jet of fluid falling steadily under the force of gravity is studied. The problem is formulated as a nonlinear free boundary-value problem for the potential. Surface tension effects are neglected. The use of perturbation expansions results in a system of equations that can be solved by an efficient numerical procedure. Computations were made for jets issuing from orifices in various shapes including an ellipse, a rectangle, and an equilateral triangle. Computational results are presented illustrating the propagation of discontinuities and the formation of thin sheets of fluid.

Geer, J. F.↗
NASA NTRS1994

Hybrid near-optimal aeroassisted orbit transfer plane change trajectories

In this paper, a hybrid methodology is used to determine optimal open loop controls for the atmospheric portion of the aeroassisted plane change problem. The method is hybrid in the sense that it combines the features of numerical collocation with the analytically tractable portions of the problem which result when the two-point boundary value problem is cast in the form of a regular perturbation problem. Various levels of approximation are introduced by eliminating particular collocation parameters and their effect upon problem complexity and required number of nodes is discussed. The results include plane changes of 10, 20, and 30 degrees for a given vehicle.

Calise, Anthony J.↗
DOE OSTI2025

Analysis of the weighted shifted boundary method for the Poisson and Stokes problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. Accuracy is maintained by properly shifting the location and values of the boundary conditions. This avoids integration over cut cells and the associated implementation issues. Recently, the Weighted SBM (WSBM) was proposed for the Navier-Stokes equations with free surfaces and the Stokes flow with moving boundaries. The attribute “weighted” in the name WSBM stems from the fact that its variational form is weighted with the elemental volume fraction of active fluid. The motivation for the development of the WSBM was the preservation of the volume of active fluid to a higher degree of accuracy, which in turn resulted in improved stability and robustness characteristics in moving-boundary, time-dependent simulations. In this article, we present the numerical analysis of the WSBM formulations for the Poisson and Stokes problems. We give mathematical conditions under which the bilinear forms defining the discrete variational formulations are uniformly coercive (Poisson problem) or inf-sup stable (Stokes problem). By these results, stability and optimal convergence is proven in the natural norm; L2-error estimates can also be derived.

Approximate domain boundaries↗
NASA NTRS1973

Unsteady, combined radiation and conduction in an absorbing, scattering, and emitting medium.

The transient cooldown of a gray, absorbing, isotropic scattering, emitting, and conducting medium bounded by gray, diffusely emitting and reflecting parallel plates is considered. Numerical solutions are obtained for the initial boundary-value problem with a discontinuous decrease in temperature at one boundary. The quasi-steady equation of radiative transfer is solved using Gaussian quadrature and a matrix eigenvector technique together with explicit numerical solution of the unsteady energy equation. Temperature and energy flux distributions are presented for variations of optical thickness, boundary emissivity, albedo, and conduction-radiation parameter.

Weston, K. C.↗
NASA NTRS1976

Hybrid finite element methods

The purpose of this paper is to show how Lagrange multipliers can be used with finite elements to achieve a number of desirable properties in the underlying approximation. For elliptic boundary value problems, variational principles can be developed in which all boundary conditions are natural. In fluid flow problems, one can endow the approximations with physically essential conservation laws.

Fix, G. M.↗
NASA NTRS1983

Parallel architectures for iterative methods on adaptive, block structured grids

A parallel computer architecture well suited to the solution of partial differential equations in complicated geometries is proposed. Algorithms for partial differential equations contain a great deal of parallelism. But this parallelism can be difficult to exploit, particularly on complex problems. One approach to extraction of this parallelism is the use of special purpose architectures tuned to a given problem class. The architecture proposed here is tuned to boundary value problems on complex domains. An adaptive elliptic algorithm which maps effectively onto the proposed architecture is considered in detail. Two levels of parallelism are exploited by the proposed architecture. First, by making use of the freedom one has in grid generation, one can construct grids which are locally regular, permitting a one to one mapping of grids to systolic style processor arrays, at least over small regions. All local parallelism can be extracted by this approach. Second, though there may be a regular global structure to the grids constructed, there will be parallelism at this level. One approach to finding and exploiting this parallelism is to use an architecture having a number of processor clusters connected by a switching network. The use of such a network creates a highly flexible architecture which automatically configures to the problem being solved.

Gannon, D.↗
NASA NTRS1993

Robust neighboring extremal guidance for the advanced launch system

With the availability of modern flight computers, realtime neighboring extremal guidance seems feasible. To overcome sensitivity to unknown system parameters and environmental uncertainties, a robust neighboring extremal guidance scheme is proposed. About the optimal trajectory, the accessory problem in the calculus of variations is formed, generating a quadratic cost criterion in the perturbed states and controls. By formulating a disturbance attenuation problem based upon the second variation cost criterion, a differential game is formulated. The game theoretic cost criterion is minimized with respect to the perturbed control but maximized with respect to the unknown parameters in the linearized dynamics. The resulting differential game problem gives rise to a two-point boundary-value problem solved using the sweep method. The sweep method solution provides a linear robust neighboring extremal guidance scheme that is applied to the Advanced Launch System.

Bain, John↗
NASA NTRS2007

Low-Thrust Control of a Lunar Mapping Orbit

A method is presented for generating and maintaining a lunar mapping orbit using continuous low-thrust hardware. Optimal control theory is used to maintain a lunar orbit that is low-altitude, near-polar, and Sun-synchronous; three typical requirements for a successful lunar mapping mission. The analysis of the optimal control problem leads to the commonly seen two-point boundary value problem, which is solved using a simple indirect shooting algorithm. Simulations are presented for a 50-day mapping duration, in which it is shown that a very tight control is achieved with thrust levels below 1 N for a 1000 kg spacecraft. A straightforward approach for using the method presented to compute missions of any duration is also discussed.

Harl, Nathan↗
NASA NTRS1972

The computation of optimal control programmes using a modified successive sweep method.

A second-order method for numerically solving control optimization problems has been developed. The method, referred to as the modified sweep method (MSM), differs from the successive sweep method (SSM) proposed by McReynolds and Bryson (1965) in that the conditions for local control optimality are used to determine the control as an explicit function of the state variables and time. The control is eliminated from the problem and the solution to the resulting two-point boundary value problem can be obtained by linear perturbation methods. The modified sweep method proposed here uncouples the perturbation equations for the state variables and the Lagrange multipliers by using a generalized matrix-Riccati transformation of variables. The resulting algorithm for the numerical iteration process is concerned with determining the initial values of a set of Lagrange multipliers rather than correcting a numerical control programme over the entire time interval of interest.

Colunga, D.↗
NASA NTRS1979

Some new results on spherical dynamos

Small-scale motions in the earth's liquid core are likely to be highly anisotropic because of the effects of rotation. Guided by physical considerations, the 'alpha-effect' described by an anisotropic tensor alpha sub ik is formulated and the corresponding boundary value problem for a sphere is solved for a variety of boundary conditions. A converged solution has been obtained only in the case of the Fermi condition of an infinitely conducting exterior of the sphere. Some remarks are made on the hypothetical upper bound on magnetic field strengths in planetary cores originally proposed by Busse (1976).

Busse, F. H.↗