Nonlinear wave propagation in fluids
Waves nonlinear propagation in fluids due to convection, reviewing various solution methods
Lick, W.↗
SEARCH · Search NASA
Results for “nonlinear propagation”
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.
High-order Two-way Artificial Boundary Conditions for Nonlinear Wave Propagation with Backscattering
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.
Fibich, Gadi↗
Nonlinear Acoustic Propagation in Low-Density Media
Nonlinear acoustic propagation in the atmosphere is usually modeled using an augmented Burgers equation accounting for atmospheric absorption and weak nonlinearity. However, the weak-nonlinearity assumption may not apply to acoustic signals propagating from the lower atmosphere that are subsequently refracted downward from the upper atmosphere (e.g., stratosphere and thermosphere) due to the decreasing air density with increasing altitude [Lonzaga, et al., Geophysical Journal International, 200(3), pp.1347-1361]. Consequently, this paper discusses the effects of a strong nonlinearity that lead to an amplitude-dependent increase in signal propagation speed. This propagation speed is obtained using a perturbation expansion where the leading-order term is simply the small-amplitude sound speed while the first-order term is the existing expression that gives rise to waveform steepening and stretching. Furthermore, the second-order term is proportional to the spectral density of the acoustic signal and is inversely proportional to the local density of the medium. Consequently, for an impulsive signal such as a sonic boom or an infrasound, the second-order effect causes a dispersion of the signal similar to the observed dispersion of acoustic signals from supersonic Concorde as well as from large explosions. This paper also discusses the extension of these results in general low-density propagation media.
sonic booms↗
Calculating Sonic-Boom Propagation
Nonlinear effects included, enabling more-realistic modeling. Modified Method of Characteristics Sonic Boom Extrapolation Program (MMOC) is computer program for sonic-boom propagation that includes shock coalescence and incorporates effects of asymmetries due to volume and lift. Numerically integrates nonlinear governing equations using data on initial data line approximately one body length from aircraft and yields sonic-boom pressure at ground as function of time or of position at given time. MMOC written in FORTRAN IV.
Darden, C. M.↗
Similarity rules for nonlinear acoustic propagation through a caustic
Similarity rules for ray analysis in nonlinear acoustic wave propagation
Hayes, W. D.↗
Nonlinear acoustic propagation in rectangular ducts
The method of multiple scales is used to obtain a second-order uniformly valid expansion for nonlinear acoustic wave propagation in a rectangular duct whose walls are treated with a nonlinear acoustic material. The wave propagation in the duct is characterized by the unsteady nonlinear Euler equations. The results show that nonlinear materials attenuate sound more than linear materials except at high acoustic frequencies. The nonlinear materials produce higher and combination tones which have higher attenuation rates than the fundamentals. Moreover, the attenuation rates of the fundamentals increase with increasing amplitude.
Nayfeh, A. H.↗
Attenuation characteristics of nonlinear pressure waves propagating in pipes
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.
Shih, C. C.↗
Nonlinear acoustic propagation in two-dimensional ducts
The method of multiple scales is used to obtain a second-order uniformly valid expansion for the nonlinear acoustic wave propagation in a two-dimensional duct whose walls are treated with a nonlinear acoustic material. The wave propagation in the duct is characterized by the unsteady nonlinear Euler equations. The results show that nonlinear effects tend to flatten and broaden the absorption versus frequency curve, in qualitative agreement with the experimental observations. Moreover, the effect of the gas nonlinearity increases with increasing sound frequency, whereas the effect of the material nonlinearity decreases with increasing sound frequency.
Nayfeh, A. H.↗
Nonlinear wave propagation using three different finite difference schemes (category 2 application)
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.
Pope, D. Stuart↗
Rao-Blackwellization for Adaptive Gaussian Sum Nonlinear Model Propagation
When dealing with imperfect data and general models of dynamic systems, the best estimate is always sought in the presence of uncertainty or unknown parameters. In many cases, as the first attempt, the Extended Kalman filter (EKF) provides sufficient solutions to handling issues arising from nonlinear and non-Gaussian estimation problems. But these issues may lead unacceptable performance and even divergence. In order to accurately capture the nonlinearities of most real-world dynamic systems, advanced filtering methods have been created to reduce filter divergence while enhancing performance. Approaches, such as Gaussian sum filtering, grid based Bayesian methods and particle filters are well-known examples of advanced methods used to represent and recursively reproduce an approximation to the state probability density function (pdf). Some of these filtering methods were conceptually developed years before their widespread uses were realized. Advanced nonlinear filtering methods currently benefit from the computing advancements in computational speeds, memory, and parallel processing. Grid based methods, multiple-model approaches and Gaussian sum filtering are numerical solutions that take advantage of different state coordinates or multiple-model methods that reduced the amount of approximations used. Choosing an efficient grid is very difficult for multi-dimensional state spaces, and oftentimes expensive computations must be done at each point. For the original Gaussian sum filter, a weighted sum of Gaussian density functions approximates the pdf but suffers at the update step for the individual component weight selections. In order to improve upon the original Gaussian sum filter, Ref. [2] introduces a weight update approach at the filter propagation stage instead of the measurement update stage. This weight update is performed by minimizing the integral square difference between the true forecast pdf and its Gaussian sum approximation. By adaptively updating each component weight during the nonlinear propagation stage an approximation of the true pdf can be successfully reconstructed. Particle filtering (PF) methods have gained popularity recently for solving nonlinear estimation problems due to their straightforward approach and the processing capabilities mentioned above. The basic concept behind PF is to represent any pdf as a set of random samples. As the number of samples increases, they will theoretically converge to the exact, equivalent representation of the desired pdf. When the estimated qth moment is needed, the samples are used for its construction allowing further analysis of the pdf characteristics. However, filter performance deteriorates as the dimension of the state vector increases. To overcome this problem Ref. [5] applies a marginalization technique for PF methods, decreasing complexity of the system to one linear and another nonlinear state estimation problem. The marginalization theory was originally developed by Rao and Blackwell independently. According to Ref. [6] it improves any given estimator under every convex loss function. The improvement comes from calculating a conditional expected value, often involving integrating out a supportive statistic. In other words, Rao-Blackwellization allows for smaller but separate computations to be carried out while reaching the main objective of the estimator. In the case of improving an estimator's variance, any supporting statistic can be removed and its variance determined. Next, any other information that dependents on the supporting statistic is found along with its respective variance. A new approach is developed here by utilizing the strengths of the adaptive Gaussian sum propagation in Ref. [2] and a marginalization approach used for PF methods found in Ref. [7]. In the following sections a modified filtering approach is presented based on a special state-space model within nonlinear systems to reduce the dimensionality of the optimization problem in Ref. [2]. First, the adaptive Gaussian sum propagation is explained and then the new marginalized adaptive Gaussian sum propagation is derived. Finally, an example simulation is presented.
state estimation↗
Strongly nonlinear wave propagation in elasto-plastic metamaterials: Low-order dynamic modeling
Nonlinear elastic metamaterials are known to support a variety of dynamic phenomena that enhance our capacity to manipulate elastic waves. Since these properties stem from complex, subwavelength geometry, full-scale dynamic simulations are often prohibitively expensive at scales of interest. Prior studies have therefore utilized low-order effective medium models, such as discrete mass-spring lattices, to capture essential properties in the long-wavelength limit. While models of this type have been successfully implemented for a wide variety of nonlinear elastic systems, they have predominantly considered dynamics depending only on the instantaneous kinematics of the lattice, neglecting history-dependent effects, such as wear and plasticity. Here, to address this limitation, the present study develops a lattice-based modeling framework for nonlinear elastic metamaterials undergoing plastic deformation. Due to the history- and rate-dependent nature of plasticity, the framework generally yields a system of differential-algebraic equations whose computational cost is significantly greater than an elastic system of comparable size. We demonstrate the method using several models inspired by classical lattice dynamics and continuum plasticity theory and explore means to obtain empirical plasticity models for general geometries, thereby gaining insight into the influence of microstructural plasticity on effective material performance, which can be used to improve the design of nonlinear mechanical metamaterials.
Dynamic simulation↗
An explicit algorithm for a fluid approach to nonlinear optics propagation using splitting and rezoning techniques
The paper presents a three-dimensional analysis of the nonlinear light matter interaction in a hydrodynamic context. It is reported that the resulting equations are a generalization of the Navier-Stokes equations subjected to an internal potential which depends solely upon the fluid density. In addition, three numerical approaches are presented to solve the governing equations using an extension of McCormack predict-corrector scheme. These are a uniform grid, a dynamic rezoned grid, and a splitting technique. It is concluded that the use of adaptive mapping and splitting techniques with McCormack two-level predictor-corrector scheme results in an efficient and reliable code whose storage requirements are modest compared with other second order methods of equal accuracy.
Mattar, F. P.↗
Microscopic Lagrangian description of warm plasmas. I - Linear wave propagation. II - Nonlinear wave interactions
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.
Kim, H.↗
Computational Modeling of Ultrafast Pulse Propagation in Nonlinear Optical Materials
There is an emerging technology of photonic (or optoelectronic) integrated circuits (PICs or OEICs). In PICs, optical and electronic components are grown together on the same chip. rib build such devices and subsystems, one needs to model the entire chip. Accurate computer modeling of electromagnetic wave propagation in semiconductors is necessary for the successful development of PICs. More specifically, these computer codes would enable the modeling of such devices, including their subsystems, such as semiconductor lasers and semiconductor amplifiers in which there is femtosecond pulse propagation. Here, the computer simulations are made by solving the full vector, nonlinear, Maxwell's equations, coupled with the semiconductor Bloch equations, without any approximations. The carrier is retained in the description of the optical pulse, (i.e. the envelope approximation is not made in the Maxwell's equations), and the rotating wave approximation is not made in the Bloch equations. These coupled equations are solved to simulate the propagation of femtosecond optical pulses in semiconductor materials. The simulations describe the dynamics of the optical pulses, as well as the interband and intraband.
Goorjian, Peter M.↗
Robust direct laser acceleration of electrons with flying-focus laser pulses
Direct laser acceleration (DLA) offers a compact source of high-charge, energetic electrons for generating secondary radiation or neutrons. While DLA in high-density plasma optimizes the energy transfer from a laser pulse to electrons, it exacerbates nonlinear propagation effects, such as filamentation, that can disrupt the acceleration process. Here, we show that superluminal flying-focus pulses (FFPs) mitigate nonlinear propagation, thereby enhancing the number of high-energy electrons and resulting x-ray yield. Three-dimensional particle-in-cell simulations show that, compared to a Gaussian pulse of equal energy (1 J) and intensity (2 × 10 20 W/cm 2 ), an FFP produces 80 × more electrons above 100 MeV, increases the electron cutoff energy by 20%, triples the high-energy x-ray yield, and improves x-ray collimation. These results illustrate the ability of spatiotemporally structured laser pulses to provide additional control in the highly nonlinear, relativistic regime of laser-plasma interactions.
Laser-produced plasmas↗
Propagation of nonlinear, radiatively damped longitudinal waves along magnetic flux tubes in the solar atmosphere
For solar magnetic flux tubes three types of waves are compared: longitudinal MHD tube waves, acoustic tube waves propagating in the same tube geometry but with rigid walls and ordinary acoustic waves in plane geometry. It is found that the effect of the distensibility of the tube is small and that longitudinal waves are essentially acoustic tube waves. Due to the tube geometry there is considerable difference between longitudinal waves or acoustic tube waves and ordinary acoustic waves. Longitudinal waves as well as acoustic tube waves show a smaller amplitude growth, larger shock formation heights, smaller mean chromospheric temperature but a steeper dependence of the temperature gradient on wave period.
Herbold, G.↗
Farfield viscous effects in nonlinear noise propagation
Discussion of the method of parametric differentiation in application to predictions of farfield noise propagation in both lossless and dissipative media. It is shown that, in the lossless medium, the governing equation, transformed to parameter space, reduces to a wave equation in the farfield. In the dissipative medium, the system of nonlinear partial differential equations, transformed to parameter space, reduces to a linear partial differential equation of the propagating type which contains a third-order derivative as well as the wave operator.
Harris, W. L., Sr.↗
The effect of non-linear propagation in jet noise
An experimental investigation of the nonlinear propagation effects which occur in the noise radiated from low and moderate Reynolds number supersonic jets has been performed. An array of three condenser microphones was used to measure the waveforms propagated by axisymmetric, cold model jets of Mach numbers 2.1 and 2.5. Relatively low Reynolds numbers were obtained by exhausting the jets into a low pressure anechoic test chamber. The results show that phenomena normally associated with nonlinear acoustic propagation, such as wave steepening, harmonic generation and wave merging, are measurable in the high speed model jets.
Gallagher, J. A.↗