Nonlinear wave propagation in fluids
Waves nonlinear propagation in fluids due to convection, reviewing various solution methods
SEARCH · Search NASA
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.
Waves nonlinear propagation in fluids due to convection, reviewing various solution methods
Three common finite difference schemes are used to examine the computation of one-dimensional nonlinear wave propagation. The schemes are studied for their responses to numerical parameters such as time step selection, boundary condition implementation, and discretization of governing equations. The performance of the schemes is compared and various numerical phenomena peculiar to each is discussed.
When solving linear scattering problems, one typically first solves for the impinging wave in the absence of obstacles. Then, by linear superposition, the original problem is reduced to one that involves only the scattered waves driven by the values of the impinging field at the surface of the obstacles. In addition, when the original domain is unbounded, special artificial boundary conditions (ABCs) that would guarantee the reflectionless propagation of waves have to be set at the outer boundary of the finite computational domain. The situation becomes conceptually different when the propagation equation is nonlinear. In this case the impinging and scattered waves can no longer be separated, and the problem has to be solved in its entirety. In particular, the boundary on which the incoming field values are prescribed, should transmit the given incoming waves in one direction and simultaneously be transparent to all the outgoing waves that travel in the opposite direction. We call this type of boundary conditions two-way ABCs. In the paper, we construct the two-way ABCs for the nonlinear Helmholtz equation that models the laser beam propagation in a medium with nonlinear index of refraction. In this case, the forward propagation is accompanied by backscattering, i.e., generation of waves in the direction opposite to that of the incoming signal. Our two-way ABCs generate no reflection of the backscattered waves and at the same time impose the correct values of the incoming wave. The ABCs are obtained for a fourth-order accurate discretization to the Helmholtz operator; the fourth-order grid convergence is corroborated experimentally by solving linear model problems. We also present solutions in the nonlinear case using the two-way ABC which, unlike the traditional Dirichlet boundary condition, allows for direct calculation of the magnitude of backscattering.
It is pointed out that the conventional iterative analysis of nonlinear plasma wave phenomena, which involves a direct use of Maxwell's equations and the equations describing the particle dynamics, leads to formidable theoretical and algebraic complexities, especially for warm plasmas. As an effective alternative, the Lagrangian method may be applied. It is shown how this method may be used in the microscopic description of small-signal wave propagation and in the study of nonlinear wave interactions. The linear theory is developed for an infinite, homogeneous, collisionless, warm magnetoplasma. A summary is presented of a perturbation expansion scheme described by Galloway and Kim (1971), and Lagrangians to third order in perturbation are considered. Attention is given to the averaged-Lagrangian density, the action-transfer and coupled-mode equations, and the general solution of the coupled-mode equations.
A series of experiments was conducted to investigate temporal and spatial velocity distributions of fluid flow in 3-in. open-end pipes of various lengths up to 210 ft, produced by the propagation of nonlinear pressure waves of various intensities. Velocity profiles across each of five sections along the pipes were measured as a function of time with the use of hot-film and hot-wire anemometers for two pressure waves produced by a piston. Peculiar configurations of the velocity profiles across the pipe section were noted, which are uncommon for steady pipe flow. Theoretical consideration was given to this phenomenon of higher velocity near the pipe wall for qualitative confirmation. Experimentally time-dependent velocity distributions along the pipe axis were compared with one-dimensional theoretical results obtained by the method of characteristics with or without diffusion term for the purpose of determining the attenuation characteristics of the nonlinear wave propagation in the pipes.
Nonlinear electromagnetic wave propagation in magnetoactive finite temperature plasma
The investigation of the acoustic wave propagation theory and numerical implementation for the situation of an isothermal atmosphere is described. A one-dimensional model to validate an asymptotic theory and a 3-D situation to relate to a realistic situation are considered. In addition, nonlinear wave propagation and the numerical treatment are included. It is known that the gravitational effects play a crucial role in the low frequency acoustic wave propagation. They propagate large distances and, as such, the numerical treatment of those problems become difficult in terms of posing boundary conditions which are valid for all frequencies.
Consideration is given to properties of obliquely propagating (with respect to the large-scale magnetic field) nonlinear Alfven and fast magnetosonic waves as described by the derivative nonlinear Schroedinger equation. The focus is on those waves which in the parallel-propagating limit are circularly polarized, exact solutions to the MHD equations. Results are presented of a study of the changes in the structure of these waves as the propagation angle is increased. Attention is also given to the evolution of initial conditions which do not correspond to stationary wave propagations. It is found that, as the propagation angle with respect to the mean field increases, the functional forms of the wave magnetic field components in the two directions orthogonal to the propagation direction become increasingly dissimilar. The angle about the mean field within which stationary waves resemble parallel-propagating, circularly polarized waves becomes smaller as the wavelength increases. If the 'initial condition' wave train is close to parallel-propagating, then there is little subsequent evolution in its form.
A model problem that simulates an atmospheric acoustic wave propagation situation that is nonlinear is considered. The model is derived from the basic Euler equations for the atmospheric flow and from the regular perturbations for the acoustic part. The nonlinear effects are studied by obtaining two successive linear problems in which the second one involves the solution of the first problem. Well posedness of these problems is discussed and approximations of the radiation boundary conditions that can be used in numerical simulations are presented.
In this paper a model problem is considered that simulates an atmospheric acoustic wave propagation situation that is nonlinear. The model is derived from the basic Euler equations for the atmospheric flow and from the regular perturbations for the acoustic part. The nonlinear effects are studied by obtaining two successive linear problems in which the second one involves the solution of the first problem. Well-posedness of these problems is discussed and approximations of the radiation boundary conditions that can be used in numerical simulations are presented.
In this paper a model problem is considered that simulates an atmospheric acoustic wave propagation situation that is nonlinear. The model is derived from the basic Euler equations for the atmospheric flow and from the regular perturbations for the acoustic part. The nonlinear effects are studied by obtaining two successive linear problems in which the second one involves the solution of the first problem. Well-posedness of these problems is discussed and approximations of the radiation boundary conditions that can be used in numerical simulations are presented.
The development of analytical techniques and numerical methods for the prediction of the stability behavior of liquid propellant rocket combustors is discussed. The studies involved nonlinear wave propagation effects and nonlinear combustion zone, liner, and nozzle responses. A second primary extension of the basic theory was directed at the inclusion of the effects of injector face baffles. Attention was restricted to the linear (small amplitude) problem, finite Mach numbers, three dimensional oscillations, and distributed combustion effects.
By the use of a Hamiltonian formulation, a basic group velocity is defined as the derivative of frequency with respect to wavenumber keeping action density constant, and is shown to represent an incremental action velocity in the general nonlinear case. The stability treatment of Whitham and Lighthill is extended to several dimensions. The water-wave analysis of Whitham (1967) is extended to two space dimensions, and is shown to predict oblique-mode instabilities for kh smaller than 1.36. A treatment of Lighthill's (1965) solution in the one-dimensional elliptic case resolves the problem of the energy distribution in the solution past the critical time.
The observations made during the encounter with comet Giacobini Zinner show that the character of MHD turbulence is governed by the magnetosonic (MS) waves generated by the pickup ions via a resonant cyclotron instability. The interaction of cometary ions in the mass loaded solar wind with MS waves propagating away from the comet and oblique to the interplanetary magnetic field (IMF) is investigated using the test particle approach. Ion trajectories, distribution functions, widths of pitch angle scattering and energy diffusion are obtained. We have also investigated the relationship between pitch angle scattering and the three parameters, namely, alpha, the angle between the solar wind flow direction and the ambient magnetic field, theta(sub BK), the angle between the ambient magnetic field and the wave propagation vector, and the ion injection velocity.
Elliptic and hyperbolic problems in unbounded regions are considered. These problems, when one wants to solve them numerically, have the difficulty of prescribing boundary conditions at infinity. Computationally, one needs a finite region in which to solve these problems. The corresponding conditions at infinity imposed on the finite distance boundaries should dictate the boundary conditions at infinity and be accurate with respect to the interior numerical scheme. The treatment of these boundary conditions for wave-like equations is discussed.
The observations made during the encounter with comet Giacobini Zinner show that the character of MHD turbulence is governed by the magnetosonic (MS) waves generated by the pickup ions via a resonant cyclotron instability. We have performed the test particle calculations to study wave particle interactions via cyclotron resonance of water group cometary ions with (1) linearly polarized monochromatic MS waves, (2) circularly polarized monochromatic MS waves, and (3) magnetic fluctuations consisting of MS waves with a power spectrum which varies as 1 / k 2 .
Coronal holes are well know sources of high speed solar wind, however, the exact acceleration mechanism of the wind is still unknown. Interplanetary scintillation (IPS) observations indicate that the fast solar wind reaches an average velocity of 800 km s(exp -1) within several solar radii with large velocity fluctuations. However, the origin of the IPS velocity spread below 10 solar radii is unclear. A previously developed coronal home model with a more realistic initial state is applied, and time-dependent, nonlinear, resistive 2.5-DMHD equations are numerically solved. It is found that nonlinear solitary-like waves with a supersonic phase speed are generated in coronal holes by torisonal Alfven waves in the radial flow velocity. The outward propagating nonlinear waves are similar in properties to sound solitons. When these waves are present, the solar wind speed and density fluctuate considerably on a time scale of an hour and on spatial scales of several solar radii in addition to the Alfvenic fluctuations. This is in qualitative agreement with the IPS velocity observations beyond 10 solar radii.
The selection of a satisfactory numerical method for calculating the propagation of steep fronted shock life waveforms in a solid rocket motor combustion chamber is discussed. A number of different numerical schemes were evaluated by comparing the results obtained for three problems: the shock tube problems; the linear wave equation, and nonlinear wave propagation in a closed tube. The most promising method--a combination of the Lax-Wendroff, Hybrid and Artificial Compression techniques, was incorporated into an existing nonlinear instability program. The capability of the modified program to treat steep fronted wave instabilities in low smoke tactical motors was verified by solving a number of motor test cases with disturbance amplitudes as high as 80% of the mean pressure.