Initial-value problem for relativistic plasma oscillations.
Initial value problem for collisionless relativistic Vlasov-Maxwell equations, noting new plasma oscillation modes and various phase velocities
Lerche, I.↗
SEARCH · Search NASA
Results for “Initial-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.
Chebyshev series solution of nonlinear ordinary differential equations - Initial-value problems
Chebyshev series solution of nonlinear ordinary differential equations - initial value problems
Byrd, P. F.↗
Formulation of plasma wave echoes as an initial-value problem.
High order plasma wave echoes formulation and relation to initial perturbations based on Landau solution of initial value problem for Vlasov equation
Nishikawa, K.↗
Solution of an initial-value problem in linear transport theory - Monoenergetic neutrons in a slab with infinite reflectors
Solution of initial value problem in linear transport theory obtained by monoenergetic neutrons migrating in thin slab with infinite reflector and isotropic scattering
Bowden, R. L.↗
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↗
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.↗
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.↗
THE INITIAL VALUE PROBLEM FOR LONGITUDINAL AND TRANSVERSAL PLASMA OSCILLATIONS
Solution of initial-value problem for linearized boltzmann equation for longitudinal and transversal plasma oscillations
INITIAL VALUE PROBLEM↗
Approximate solutions to the Navier-Stokes initial value problem
A Galerkin-Ritz procedure for any arbitrary system of field equations is shown to follow generically from 'two-functional' variational conditions. The initial-value problem for boundary-free incompressible Navier-Stokes flow is solved analytically in the two-parameter-function approximation.
Rosen, G.↗
Perfect spinfluid: A divergence-type approach
We present a new formulation of nondissipative relativistic spin hydrodynamics that incorporates spin degrees of freedom into the divergence-type theory framework. Due to the divergence-type structure, it is straightforward to enforce nonlinear causality and symmetric hyperbolicity of the equations of motion, ensuring local well-posedness of the initial-value problem and stability of the theory. Furthermore, in a specific realization based on spin kinetic theory, we prove that the equations of motion remain nonlinearly causal and symmetric-hyperbolic to all orders in the spin potential, provided a specific thermodynamic constraint is satisfied. Here, this framework can be applied for numerical simulations to study the dynamics of spin-polarized fluids, such as the quark-gluon plasma in heavy-ion collisions.
Chirality↗
Numerical Studies of the Non-linear Vlasov Equation
The subject of this dissertation is the numerical integration of the initial-value problem for the non-linear Vlasov equation. The Vlasov equation is used to describe the dynamics of a "collisionless", one-dimensional, classical electron gas confined between two perfectly reflecting boundaries. Only the long-range Coulomb interactions of the electrons are considered; effects associated with the discrete structure are neglected. The numerical results obtained for non-linear Landau damping compare well with similar results obtained by Knorr. A general statement of the results on stable initial conditions is: As the degree of non-linearity of the initial conditions is increased, the deviation from linear Landau damping appears sooner and is more severe. In some cases damping was observed to cease. Curves showing the time dependence of the damping decrement are derived and compared with predictions of non-linear theories. New results obtained in this study include the observation that for strongly non-linear cases, the damping of the electric field causes an initially Maxwellian fo (v, O) to develop a peak in the neighborhood of the phase velocity; strong growth of the second harmonic is seen after fo (v, t) develops such a peak. Also new in this study is the interpretation of the development of a certain class of strongly unstable initial conditions as approaching an inhomogeneous equilibrium.
NUMERICAL INTEGRATION↗
Numerical integration of ordinary differential equations of various orders
Report describes techniques for the numerical integration of differential equations of various orders. Modified multistep predictor-corrector methods for general initial-value problems are discussed and new methods are introduced.
Gear, C. W.↗
Sideband growth in nonlinear Landau wave-particle interaction.
The distortion of the electron velocity distribution caused by a large amplitude Landau wave is determined analytically for the initial-value problem. The resulting stability of electrostatic perturbations impressed on the evolving plasma is studied. Narrow sidebands of the applied frequency experience consecutive growths of large magnitude during the early stages of the nonlinear wave-particle interaction. The significance of the derived results to both wave propagation experiments and triggered VLF emissions in the magnetosphere is discussed.
Brinca, A. L.↗
On the Navier-Stokes equations with constant total temperature
For various applications in fluid dynamics, one can assume that the total temperature is constant. Therefore, the energy equations can be replaced by an algebraic relation. The resulting set of equations in the inviscid case is analyzed in this paper. It is shown that the system is strictly hyperbolic and well posed for the initial-value problem. Boundary conditions are described such that the linearized system is well posed. The hopscotch method is investigated and numerical results are presented.
Gottlieb, D.↗
Efficient and flexible multirate temporal adaptivity
In this work we present two new families of multirate time step adaptivity controllers, that are designed to work with embedded multirate infinitesimal (MRI) time integration methods for adapting time steps when solving problems with multiple time scales. We compare these controllers against competing approaches on two benchmark problems, showing that the proposed methods offer dramatically improved performance and flexibility. The combination of embedded MRI methods and the proposed controllers enable adaptive simulations of problems with a potentially arbitrary number of time scales, achieving high accuracy while maintaining low computational cost. Additionally, we introduce a new set of embeddings for the family of explicit multirate exponential Runge–Kutta (MERK) methods of orders 2 through 5, resulting in the first-ever fifth-order embedded MRI method. Finally, we compare the performance of a wide range of embedded MRI methods on our benchmark problems to provide guidance on how to select an appropriate MRI method and multirate controller.
97 MATHEMATICS AND COMPUTING↗
FLORA Equilibrium and Stability Code Archive
FLORA solves, in a 2-D domain (radial and axial dimensions) with a specified azimuthal Fourier mode, for the linearized stability of a long, thin, axisymmetric plasma equilibrium in an applied magnetic field. Before the stability equation is solved, FLORA solves a set of simple equations for pressure balance that specify the equilibrium magnetic field and plasma pressure in the long-thin limit given a simplified description of the magnetic coils. It uses an initial-value method for the linear stability problem in which an equilibrium is given an initial perturbation to its magnetic field, and the temporal behavior of the perturbation is followed. The perturbation has been Fourier expanded in the azimuthal coordinate; each azimuthal mode must be examined separately. The complex partial differential equation of motion for the perturbed radial displacement of the field lines is solved as a coupled system of two real p.d.e.'s and the solution consists of two parts, the real part and the imaginary part. The system is solved by bringing the coupling terms in each equation to the right side and using an iterative technique.
Cohen, Bruce [Lawrence Livermore National Laborato↗
A high order accurate finite element algorithm for high Reynolds number flow prediction
A Galerkin-weighted residuals formulation is employed to establish an implicit finite element solution algorithm for generally nonlinear initial-boundary value problems. Solution accuracy, and convergence rate with discretization refinement, are quantized in several error norms, by a systematic study of numerical solutions to several nonlinear parabolic and a hyperbolic partial differential equation characteristic of the equations governing fluid flows. Solutions are generated using selective linear, quadratic and cubic basis functions. Richardson extrapolation is employed to generate a higher-order accurate solution to facilitate isolation of truncation error in all norms. Extension of the mathematical theory underlying accuracy and convergence concepts for linear elliptic equations is predicted for equations characteristic of laminar and turbulent fluid flows at nonmodest Reynolds number. The nondiagonal initial-value matrix structure introduced by the finite element theory is determined intrinsic to improved solution accuracy and convergence. A factored Jacobian iteration algorithm is derived and evaluated to yield a consequential reduction in both computer storage and execution CPU requirements while retaining solution accuracy.
Baker, A. J.↗
Reduction of the Dirichlet problem to an initial value problem.
Although the derivation is concerned with solutions for plane regions with prescribed boundary values, the approach presented could by easily generalized to higher dimensions. The initial-value method is derived by a combination of invariant imbedding techniques and the Fredholm integral equation method of representation of the potential as a function of a dilayer distribution on the boundary of the region in question.
Kalaba, R.↗