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

Convenient stability criteria for difference approximations of hyperbolic initial-boundary value problems

The purpose of this paper is to achieve more versatile, convenient stability criteria for a wide class of finite-difference approximations to initial boundary value problems associated with the hyperbolic system u sub t = au sub x + Bu + f in the quarter-plane x greater than or equal to 0, t greater than or equal to 0. With these criteria, stability is easily established for a large number of examples, thus incorporating and generalizing many of the cases studied in recent literature.

Goldberg, M.↗
NASA NTRS1973

A technique for increasing the accuracy of the numerical inversion of the Laplace transform with applications

A technique is introduced which extends the range of useful approximation of numerical inversion techniques to many cycles of an oscillatory function without requiring either the evaluation of the image function for many values of s or the computation of higher-order terms. The technique consists in reducing a given initial value problem defined over some interval into a sequence of initial value problems defined over a set of subintervals. Several numerical examples demonstrate the utility of the method.

Berger, B. S.↗
NASA NTRS1996

Quadrature-Free Implementation of the Discontinuous Galerkin Method for Hyperbolic Equations

A discontinuous Galerkin formulation that avoids the use of discrete quadrature formulas is described and applied to linear and nonlinear test problems in one and two space dimensions. This approach requires less computational time and storage than conventional implementations but preserves the compactness and robustness inherent to the discontinuous Galerkin method. Test problems include both linear and nonlinear one-dimensional scalar advection of botH smooth and discontinuous initial value problems, two-dimensional scalar advection of smooth initial value problems discretized by using unstructured grids with varying degrees of smoothness and regularity, and two-dimensional linear Euler solutions on unstructured grids.

Atkins, H. L.↗
NASA NTRS1986

Some insights into the stability of difference approximations for hyperbolic initial-boundary-value problems

This paper states a conjecture which relates Lax-Richmyer stability to the algebraic test of the stability theory of Gustafsson, Kreiss, and Sundstrom (1972), developed for difference approximations to initial boundary value problems where the matrix size J increases linearly with n as n goes to infinity. This corresponds to mesh refinement in both space and time for t = n x delta t = constant.

Warming, Robert F.↗
NASA NTRS1993

Adaptive wavelet collocation methods for initial value boundary problems of nonlinear PDE's

We have designed a cubic spline wavelet decomposition for the Sobolev space H(sup 2)(sub 0)(I) where I is a bounded interval. Based on a special 'point-wise orthogonality' of the wavelet basis functions, a fast Discrete Wavelet Transform (DWT) is constructed. This DWT transform will map discrete samples of a function to its wavelet expansion coefficients in O(N log N) operations. Using this transform, we propose a collocation method for the initial value boundary problem of nonlinear PDE's. Then, we test the efficiency of the DWT transform and apply the collocation method to solve linear and nonlinear PDE's.

Cai, Wei↗
NASA NTRS1973

Numerical integration of shell equations using the field method.

The 'field method' for the numerical solution of even-order linear boundary-value problems in ordinary differential equations is formulated. This method converts the boundary-value problem into two successive initial-value problems, which may be solved by standard forward integration techniques. The method has been implemented in a computer program to calculate the static response of ring-stiffened branched shells of revolution to asymmetric loads. For such problems the field method eliminates the well-known numerical problem of 'long subintervals' and also executes operations significantly faster than other numerical integration methods.

Cohen, G. A.↗
NASA NTRS1980

Evolution of solar magnetic fields - A new approach to MHD initial-boundary value problems by the method of nearcharacteristics

A method of analysis for the MHD initial-boundary problem is presented in which the model's formulation is based on the method of nearcharacteristics developed by Werner (1968) and modified by Shin and Kot (1978). With this method, the physical causality relationship can be traced from the perturbation to the response as in the method of characteristics, while achieving the advantage of a considerable reduction in mathematical procedures. The method offers the advantage of examining not only the evolution of nonforce free fields, but also the changes of physical conditions in the atmosphere accompanying the evolution of magnetic fields. The physical validity of the method is demonstrated with examples, and their significance in interpreting observations is discussed.

Nakagawa, Y.↗
NASA NTRS1984

Canonical solutions for unsteady flow fields

The initial value problem of one-dimensional gas-dynamics involving discontinuous, nonuniform initial data is discussed. Canonical solutions which are valid in a small x, t region aroung a discontinuity, and which include the first order effects of nonuniformities in the data, are derived explicitly. The theory is derived by considering a group of elementary piston problems. Solutions with a shock or with a centered expansion wave are worked out individually in order to relate initial flow properties and their gradients to the speed and acceleration of the discontinuity waves. They are then combined to represent the solution of a general initial value problem by regarding the piston path as a contact line. In addition, problems with chemical reaction are discussed in terms of elementary piston problems which involve strong detonation waves, Chapman-Jouguet detonation waves, and deflagration waves.

Liu, G. C.↗
NASA NTRS1979

Possible inverse cascade behavior for drift-wave turbulence

The turbulent spectral properties of the dynamical equation of Hasegawa and Mima (1978) governing the evolution of the electrostatic potential in drift-wave turbulence is investigated for two formulations of the problem: (1) as a nondissipative initial value problem, with the potential represented by a truncated Fourier series with large number of terms, and (2) as a dissipative problem with a small viscous dissipation at very short spatial scales, and a long wavelength forcing term at longer wavelengths. It is found that Hasegawa and Mima's prediction for the nondissipative, truncated initial value modal problem is accurate, but substantial differences exist for the forced dissipative case between computer results and analytical predictions based on a wave kinetic equation of Kadomtsev. Much better agreement is found with a simple dual-cascade model based on Kraichnan's generalization of Kolmogorov's cascade arguments.

Fyfe, D.↗
NASA NTRS1999

Towards Understanding the Mechanism of Receptivity and Bypass Dynamics in Laminar Boundary Layers

Three problems concerning laminar-turbulent transition are addressed by solving a series of initial value problems. The first problem is the calculation of resonance within the continuous spectrum of the Blasius boundary layer. The second is calculation of the growth of Tollmien-Schlichting waves that are a direct result of disturbances that only lie outside of the boundary layer. And, the third problem is the calculation of non-parallel effects. Together, these problems represent a unified approach to the study of freestream disturbance effects that could lead to transition. Solutions to the temporal, initial-value problem with an inhomogeneous forcing term imposed upon the flow is sought. By solving a series of problems, it is shown that: A transient disturbance lying completely outside of the boundary layer can lead to the growth of an unstable Tollmien-Schlichting wave. A resonance with the continuous spectrum leads to strong amplification that may provide a mechanism for bypass transition once nonlinear effects are considered. A disturbance with a very weak unstable Tollmien-Schlichting wave can lead to a much stronger Tollmien-Schlichting wave downstream, if the original disturbance has a significant portion of its energy in the continuum modes.

Lasseigne, D. G.↗
NASA NTRS2008

Numerical Determination of Critical Conditions for Thermal Ignition

The determination of ignition or thermal explosion in an oxidizing porous body of material, as described by a dimensionless reaction-diffusion equation of the form .tu = .2u + .e-1/u over the bounded region O, is critically reexamined from a modern perspective using numerical methodologies. First, the classic stationary model is revisited to establish the proper reference frame for the steady-state solution space, and it is demonstrated how the resulting nonlinear two-point boundary value problem can be reexpressed as an initial value problem for a system of first-order differential equations, which may be readily solved using standard algorithms. Then, the numerical procedure is implemented and thoroughly validated against previous computational results based on sophisticated path-following techniques. Next, the transient nonstationary model is attacked, and the full nonlinear form of the reaction-diffusion equation, including a generalized convective boundary condition, is discretized and expressed as a system of linear algebraic equations. The numerical methodology is implemented as a computer algorithm, and validation computations are carried out as a prelude to a broad-ranging evaluation of the assembly problem and identification of the watershed critical initial temperature conditions for thermal ignition. This numerical methodology is then used as the basis for studying the relationship between the shape of the critical initial temperature distribution and the corresponding spatial moments of its energy content integral and an attempt to forge a fundamental conjecture governing this relation. Finally, the effects of dynamic boundary conditions on the classic storage problem are investigated and the groundwork is laid for the development of an approximate solution methodology based on adaptation of the standard stationary model.

Luo, W.↗
NASA NTRS1975

Computer analysis of multicircuit shells of revolution by the field method

The field method, presented previously for the solution of even-order linear boundary value problems defined on one-dimensional open branch domains, is extended to boundary value problems defined on one-dimensional domains containing circuits. This method converts the boundary value problem into two successive numerically stable initial value problems, which may be solved by standard forward integration techniques. In addition, a new method for the treatment of singular boundary conditions is presented. This method, which amounts to a partial interchange of the roles of force and displacement variables, is problem independent with respect to both accuracy and speed of execution. This method was implemented in a computer program to calculate the static response of ring stiffened orthotropic multicircuit shells of revolution to asymmetric loads. Solutions are presented for sample problems which illustrate the accuracy and efficiency of the method.

Cohen, G. A.↗
NASA NTRS2003

Aeroelastic System Development Using Proper Orthogonal Decomposition and Volterra Theory

This research combines Volterra theory and proper orthogonal decomposition (POD) into a hybrid methodology for reduced-order modeling of aeroelastic systems. The out-come of the method is a set of linear ordinary differential equations (ODEs) describing the modal amplitudes associated with both the structural modes and the POD basis functions for the uid. For this research, the structural modes are sine waves of varying frequency, and the Volterra-POD approach is applied to the fluid dynamics equations. The structural modes are treated as forcing terms which are impulsed as part of the uid model realization. Using this approach, structural and uid operators are coupled into a single aeroelastic operator. This coupling converts a free boundary uid problem into an initial value problem, while preserving the parameter (or parameters) of interest for sensitivity analysis. The approach is applied to an elastic panel in supersonic cross ow. The hybrid Volterra-POD approach provides a low-order uid model in state-space form. The linear uid model is tightly coupled with a nonlinear panel model using an implicit integration scheme. The resulting aeroelastic model provides correct limit-cycle oscillation prediction over a wide range of panel dynamic pressure values. Time integration of the reduced-order aeroelastic model is four orders of magnitude faster than the high-order solution procedure developed for this research using traditional uid and structural solvers.

Lucia, David J.↗
NASA NTRS1972

Numerical method for the solution of large systems of differential equations of the boundary layer type

A numerical method for the solution of large systems of nonlinear differential equations of the boundary-layer type is described. The method is a modification of the technique for satisfying asymptotic boundary conditions. The present method employs inverse interpolation instead of the Newton method to adjust the initial conditions of the related initial-value problem. This eliminates the so-called perturbation equations. The elimination of the perturbation equations not only reduces the user's preliminary work in the application of the method, but also reduces the number of time-consuming initial-value problems to be numerically solved at each iteration. For further ease of application, the solution of the overdetermined system for the unknown initial conditions is obtained automatically by applying Golub's linear least-squares algorithm. The relative ease of application of the proposed numerical method increases directly as the order of the differential-equation system increases. Hence, the method is especially attractive for the solution of large-order systems. After the method is described, it is applied to a fifth-order problem from boundary-layer theory.

Green, M. J.↗
NASA NTRS1992

Convective instability in differentially rotating disks

A normal mode analysis for nonaxisymmetric perturbations in a thin, differentially rotating disk with a vertical structure that is isothermal and convectively unstable is performed. The vertical gravity is assumed to be external and constant. The perturbation scale is assumed to be much shorter than the radius of the disk but comparable to or less than the thickness. The initial value problem is formulated in shearing coordinates. Dispersion relations are obtained for the three limiting cases of zero shear, axisymmetric perturbations, and small radial wavelengths. The full effects of shear are studied by integrating numerically the initial value problem. Nonaxisymmetric local Fourier modes are found to have a radial wavenumber that increases linearly with time in proportion to the shear times the azimuthal wavenumber. While Coriolis forces exert stabilizing effects on the convective modes, reducing their growth rate and the range of unstable wavelengths, shear has destablizing effects inasmuch as it reduces the epicyclic frequency at a given angular velocity. In a Keplerian disk, perturbations with azimuthal wavelengths about 2 times smaller than vertical wavelengths grow exponentially.

Ryu, Dongsu↗
DOE OSTI2026

An efficient explicit implementation of a near-optimal quantum algorithm for simulating linear dissipative differential equations

We propose an efficient block-encoding technique for the implementation of the Linear Combination of Hamiltonian Simulations (LCHS) for simulating dissipative initial-value problems. This algorithm approximates a target nonunitary operator as a weighted sum of Hamiltonian evolutions, thereby emulating a dissipative problem by mixing various time scales. We introduce an efficient encoding of the LCHS into a quantum circuit based on a simple coordinate transformation that turns the dependence on the summation index into a trigonometric function. Classically, this method is equivalent to the use of a highly accurate Fejér-Clenshaw-Curtis quadrature formula. Quantumly, this significantly simplifies block-encoding of a dissipative problem and allows one to perform an exponential number of Hamiltonian simulations by a single Quantum Signal Processing (QSP) circuit. The resulting LCHS circuit has high success probability and the selector scales logarithmically with the number of terms in the LCHS sum and linearly with time. Careful analysis of error convergence proves that this method is more efficient than other LCHS circuits that have recently appeared in the literature. We verify the quantum circuit and its scaling by simulating it on a digital emulator of fault-tolerant quantum computers and, as a test problem, solve the advection-diffusion equation. The proposed algorithm can be used for simulating a wide class of nonunitary initial-value problems including the Liouville equation with added dissipation and linear embeddings of nonlinear systems, such as the Koopman-von Neumann and Carleman embeddings.

Novikau, I [Lawrence Livermore National Laboratory↗