Search NASA⌕ Search

SEARCH · Search NASA

Results for “INITIAL 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
DOE OSTI2025

A provably stable numerical method for the anisotropic diffusion equation in confined magnetic fields

We present a novel numerical method for solving the anisotropic diffusion equation in magnetic fields confined to a periodic box which is accurate and provably stable. We derive energy estimates of the solution of the continuous initial boundary value problem. A discrete formulation is presented using operator splitting in time with the summation by parts finite difference approximation of spatial derivatives for the perpendicular diffusion operator. Weak penalty procedures are derived for implementing both boundary conditions and parallel diffusion operator obtained by field line tracing. We prove that the fully-discrete approximation is unconditionally stable. Discrete energy estimates are shown to match the continuous energy estimate given the correct choice of penalty parameters. A nonlinear penalty parameter is shown to provide an effective method for tuning the parallel diffusion penalty and significantly minimises rounding errors. Several numerical experiments, using manufactured solutions, the “NIMROD benchmark” problem and a single island problem, are presented to verify numerical accuracy, convergence, and asymptotic preserving properties of the method. Finally, we present a magnetic field with chaotic regions and islands and show the contours of the anisotropic diffusion equation reproduce key features in the field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Explicit Runge–Kutta Methods that Alleviate Order Reduction

Explicit Runge–Kutta (RK) methods are susceptible to a reduction in the observed order of convergence when applied to an initial boundary value problem with time-dependent boundary conditions. We study conditions on explicit RK methods that guarantee high order convergence for linear problems; we refer to these conditions as weak stage order conditions. We prove a general relationship between the method’s order, weak stage order, and number of stages. Furthermore, we derive explicit RK methods with high weak stage order and demonstrate, through numerical tests, that they avoid the order reduction phenomenon up to any order for linear problems and up to order three for nonlinear problems.

explicit Runge–Kutta↗
NASA NTRS1962

The Restricted Three-Body-Problem as a Perturbation of Euler's Problem of Two Fixed Centers and Its Application to Lunar Trajectories

The restricted Three-Body-Problem considers the motion of an infinitesimal mass under the gravitational attraction of two finite masses, which revolve about their common center of gravity in coplanar circles. It is well known that Euler's problem of two fixed centers, consisting of the motion of an infinitesimal mass under the gravitational attraction of two finite masses fixed in space, can be solved by elliptic functions. The idea presented here is to take the solution of Euler's problem as the solution of the restricted Three-Body-Problem by allowing the initial values to be functions of time now. Differential equations for the perturbed initial values are established. These equations can be given in closed form by using the fact that the transformation to the perturbed initial values of Euler's problem is canonical. Thus, an approximation can be obtained for the solution of the restricted Three-Body-Problem. The method can also be used to represent classes of neighboring trajectories for guidance purposes.

Euler equation↗
NASA NTRS1965

A users' guide for SEARCH and ITER

Computer program for calculating initial values of trajectory and iterative program for running multiple cases with varied parameters

INITIAL VALUE PROBLEM↗
NASA NTRS1966

Dynamic behavior of soap films.

Dynamic behavior of soap film stretched between two coaxial rings, assuming potential energy of this surface to be proportional to surface area

FILM↗
NASA NTRS1973

Transient ion neutralization by electrons.

The nonlinear initial-boundary-value problems describing the lateral neutralization of ion beams for the cases that (1) an auxiliary electric field accelerates the electrons into the ion space, and (2) the electrons are injected into the ion space at a prescribed current density are treated. Analytical solutions are derived which give the position and speed of the neutralization front as a function of time, and the temporal development of the electron density, velocity, and electric fields during the neutralization process.

Wilhelm, H. E.↗
NASA NTRS1972

A reciprocal theorem for a mixture theory

A dynamic reciprocal theorem for a linearized theory of interacting media is developed. The constituents of the mixture are a linear elastic solid and a linearly viscous fluid. In addition to Steel's field equations, boundary conditions and inequalities on the material constants that have been shown by Atkin, Chadwick and Steel to be sufficient to guarantee uniqueness of solution to initial-boundary value problems are used. The elements of the theory are given and two different boundary value problems are considered. The reciprocal theorem is derived with the aid of the Laplace transform and the divergence theorem and this section is concluded with a discussion of the special cases which arise when one of the constituents of the mixture is absent.

Martin, C. J.↗
NASA NTRS1973

A highly stable explicit integration technique for computational continuum mechanics

A user-oriented subroutine package is built around a highly stable explicit integration algorithm for solution of large order systems of ordinary differential equations, as result from discretization of initial-boundary value problems in continuum mechanics. Fast and accurate solutions, for problems in laminar and turbulent, two and three-dimensional viscous flow fields and multi-dimensional transient heat transfer, are presented using this algorithm, as embodied within a general purpose finite element computer program.

Baker, A. J.↗
NASA NTRS1975

A reciprocal theorem for a mixture of interacting continuous media

Using a linearized theory of interacting continuous media composed of linear elastic solid and linear viscous fluid, a reciprocal theorem is derived which relates the solution of one initial-boundary value problem to another. As an illustration of the reciprocal theorem, we obtain the early time displacement and velocity field of a mixture occupying an infinite region and subjected to an impulsively applied moving point load acting on the solid constituent.

Lee, Y. M.↗
NASA NTRS1979

Comptonization of X-rays by low-temperature electrons

A method is described for calculating the spectrum that results from the Compton scattering of a monochromatic source of X-rays by low-temperature electrons, both for initial-value relaxation problems and for steady-state spatial diffusion problems. The method gives an exact solution of the inital-value problem for evolution of the spectrum in an infinite homogeneous medium if Klein-Nishina corrections to the Thomson cross section are neglected. This, together with approximate solutions for problems in which Klein-Nishina corrections are significant and/or spatial diffusion occurs, shows spectral structure near the original photon wavelength that may be used to infer physical conditions in cosmic X-ray sources. Explicit results, shown for examples of time relaxation in an infinite medium and spatial diffusion through a uniform sphere, are compared with results obtained by Monte Carlo calculations and by solving the appropriate Fokker-Planck equation.

Illarionov, A.↗
NASA NTRS1980

Computation of three-dimensional viscous supersonic flow in inlets

A new approach has been developed for the computation of the three-dimensional viscous supersonic flow with embedded subsonic regions adjacent to solid boundaries and is applied to a mixed-compression supersonic inlet typical of current designs. The approach uses a reduced form of the three-dimensional Navier-Stokes equations so that the resultant equations can be treated as an initial boundary value problem and thus be solved by non-iterative forward marching in space. The numerical procedure utilizes an efficient consistently-split linearized block implicit technique to solve the finite difference analogues to the set of governing partial differential equations.

Buggeln, R. C.↗
NASA NTRS1981

Numerical approximation of boundary conditions with applications to inviscid equations of gas dynamics

A comprehensive overview of the state of the art of well-posedness and stability analysis of difference approximations for initial boundary value problems of the hyperbolic type is presented. The applicability of recent theoretical development to practical calculations for nonlinear gas dynamics is examined. The one dimensional inviscid gas dynamics equations in conservation law form are selected for numerical experiments. The class of implicit schemes developed from linear multistep methods in ordinary differential equations is chosen and the use of linear extrapolation as an explicit or implicit boundary scheme is emphasized. Specification of boundary data in the primitive variables and computation in terms of the conservative variables in the interior is discussed. Some numerical examples for the quasi-one-dimensional nozzle are given.

Yee, H. C.↗