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 91 records · Page 5

Thermal microwave emission from a random inhomogeneous layer over a homogeneous medium using the method of invariant imbedding

The paper studies thermal microwave emission from an inhomogeneous slab of a random medium, with possible nonuniform absorption, scattering, and temperature profiles, bounded by different dielectrics on both sides. The invariant imbedding method is used to cast the boundary value problem of the radiative transfer equations into an initial value problem at zero slab thickness. As a numerical example, the angular and polarization variations of brightness temperatures for ice over water are considered.

Tsang, L.↗

The theoretical accuracy of Runge-Kutta time discretizations for the initial boundary value problem: A careful study of the boundary error

The conventional method of imposing time dependent boundary conditions for Runge-Kutta (RK) time advancement reduces the formal accuracy of the space-time method to first order locally, and second order globally, independently of the spatial operator. This counter intuitive result is analyzed in this paper. Two methods of eliminating this problem are proposed for the linear constant coefficient case: (1) impose the exact boundary condition only at the end of the complete RK cycle, (2) impose consistent intermediate boundary conditions derived from the physical boundary condition and its derivatives. The first method, while retaining the RK accuracy in all cases, results in a scheme with much reduced CFL condition, rendering the RK scheme less attractive. The second method retains the same allowable time step as the periodic problem. However it is a general remedy only for the linear case. For non-linear hyperbolic equations the second method is effective only for for RK schemes of third order accuracy or less. Numerical studies are presented to verify the efficacy of each approach.

Carpenter, Mark H.↗

Numerical Simulation of Rotation-Driven Plasma Transport In the Jovian Magnetosphere

A Jupiter version of the Rice Convection Model (RCM-J) was developed with support of an earlier NASA SR&T grant. The conversion from Earth to Jupiter included adding currents driven by centrifugal force, reversing the planetary magnetic field, and rescaling various parameters. A series of informative runs was carried out, all of them solving initial value problems. The simulations followed an initial plasma torus configuration as it fell apart by interchange instability. Some conclusions from the simulations were the following: 1. We confirmed that, for conventional values of the torus density and ionospheric conductance, the torus disintegrates by interchange instability on a time scale of approx. one day, which is 1-2 orders of magnitude shorter than the best estimates of the average residence time of plasma in the torus. 2. In the model, the instability could be slowed to an arbitrary degree by the addition of sufficient impounding energetic particles, as suggested earlier by Siscoe et al (1981). However, the observed energetic particles do not seem sufficient to guarantee impoundment (e.g., Mauk et al., 1996). 3. Whether inhibited by impoundment or not, the interchange was found to proceed by the formation of long fingers, which get thinner as they get longer. This picture differed dramatically from the conventional radial-diffusion picture (e.g., Siscoe and Summers (1981)), more superficially with the outward-moving-blob picture (Pontius and Hill, 1989). The obvious limitation of the original RCM-J was that it could not represent a plasma source. We could represent the decay of a pre-existing torus, but we could not represent the way ionization of material from Io continually replenishes the plasma. We consequently were precluded from studying a whole set of fundamental issues of torus theory, including whether the system can come to a steady state.

Wolf, Richard A.↗

Numerical integration of kinetic equations.

Relaxation to Maxwellian of Balescu-Lenard equation solved numerically as initial value problem for isotropic distributions of electrons in background of positive charge

MAXWELL DISTRIBUTION↗

Applications of estimation theory to numerical weather prediction

Numerical weather prediction (NWP) is an initial value problem for a system of nonlinear partial differential equations in which the initial values are known only incompletely and inaccurately. Data at initial time can be supplemented, however, by observations of the system distributed over a time interval preceding it. Estimation theory was successful in approaching such problems for models governed by systems of ordinary differential equations and of linear PDEs. Estimation-theoretic methods for NWP are developed. A model exhibiting many features of large scale atmospheric flow important in NWP is the one governed by the shallow fluid equations. The estimation problem for a linearized formulation of these equations is studied. A finite difference version of the equations is used as a forecast model to simulate the numerical models used in NWP.

Cohn, S.↗

Crack propagation and arrest in pressurized containers

The problem of crack propagation and arrest in a finite volume cylindrical container filled with pressurized gas is considered. It is assumed that the cylinder contains a symmetrically located longitudinal part-through crack with a relatively small net ligament. The net ligament suddenly ruptures initiating the process of fracture propagation and depressurization in the cylinder. Thus the problem is a coupled gas dynamics and solid mechanics problem the exact formulation of which does not seem to be possible. The problem is reduced to a proper initial value problem by introducing a dynamic fracture criterion which relates the crack acceleration to the difference between a load factor and the corresponding strength parameter. The results indicate that generally in gas filled cylinders fracture arrest is not possible unless the material behaves in a ductile manner and the container is relatively long.

Erdogan, F.↗

Crack propagation and arrest in pressurized containers

The problem of crack propagation and arrest in a finite volume cylindrical container filled with pressurized gas is considered. It is assumed that the cylinder contains a symmetrically located longitudinal part-through crack with a relatively small net ligament. The net ligament suddenly ruptures initiating the process of fracture propagation and depressurization in the cylinder. The problem is formulated by making two major assumptions, namely, that the shell problem is quasi-static and the gas dynamics problem is one-dimensional. The problem is reduced to a proper initial value problem by introducing a dynamic fracture criterion which relates the crack acceleration to the difference between a load factor and the corresponding strength parameter. The main results are demonstrated by considering two examples, an aluminum cylinder which may behave in a quasi-brittle manner, and a steel cylinder which undergoes ductile fracture. The results indicate that generally in gas-filled cylinders fracture arrest is not possible unless the material behaves in a ductile manner and the container is relatively long.

Erdogan, F.↗