Weakly nonlinear acoustic wave propagation in a nonlinear orthotropic circular cylindrical waveguide

Nonlinear acoustic wave propagation is considered in an infinite orthotropic thin circular cylindrical waveguide. The modes are non-planar having small but finite amplitude. The fluid is assumed to be ideal and inviscid with no mean flow. The cylindrical waveguide is modeled using the Donnell's nonlinear theory for thin cylindrical shells. The approximate solutions for the acoustic velocity potential are found using the method of multiple scales (MMS) in space and time. The calculations are presented up to the third order of the small parameter. It is found that at some frequencies the amplitude modulation is governed by the Nonlinear Schrodinger Equation (NLSE). The first objective is to study the nonlinear term in the NLSE, as the sign of the nonlinear term determines the stability of the amplitude modulation. On the other hand, at other specific frequencies, interactions occur between the primary wave and its higher harmonics. Here, the objective is to identify the frequencies of the higher harmonic interactions. Lastly, the linear terms in the NLSE obtained using the MMS calculations are validated. All three objectives are met using an asymptotic analysis of the dispersion equation. (C) 2015 Acoustical Society of America.

On post-instability processes in adiabatic shear in hot rolled steel

EXPERIMENTS carried out using a split Hopkinson torsional bar have shown that only one shear band develops in specimens of hot rolled steel which break during testing. We observed, however, that in specimens which were not deformed to failure, several fine shear bands appeared. We believe that these formed during the loading cycle before the appearance of the final shear band and were not due to the effect of unloading. So we developed a numerical model to study the evolution of shear banding from several finite amplitude disturbances (FADs) in both temperature and strain rate. This numerical model reveals the detailed processes by which the FADs evolve into a fully developed shear band and suggests that beyond instability, the so-called shear banding process consists of two stages: inhomogeneous shearing and true shear-banding. The latter is characterized by the collapse of the stress and an abrupt increase of the local shear strain rate.

张力腿平台有限振幅运动的方程和数值解

论证了张力腿平台（TLP）在波浪作用下发生有限振幅运动时，所受惯性力、粘性力、浮力等载荷不仅与波浪场有关，还与瞬时响应有关，是响应的非线性函数；张力腿拉力也是各自由度位移的非线性函数．所以分析唧受力时必须考虑平台的瞬时加速度、速度和位移，在瞬时位置建立运动方程．据此推导出TLP发生有限振幅运动时的外力计算公式，建立了TLP 6自由度有限振幅运动非线性控制方程．其中考虑了由6自由度有限位移引起的多种非线性因素，如各自由度之间的耦合、瞬时湿表面、瞬时位置等；还包括自由表面效应、粘性力等因素引起的非线性．用数值方法求解所得到的非线性运动方程．对典型平台ISSC TLP进行了数值分析，求得该平台在规则波作用下的6自由度运动响应．用退化到线性范围的解与已有解进行了对比，吻合良好．数值结果表明，综合考虑非线性因素后响应有明显改变．

郭永怀文集

Steady capillary-gravity waves on deep water

The properties of capillary-gravity waves of permanent form on deep water are studied. Two different formulations to the problem are given. The theory of simple bifurcation is reviewed. For small amplitude waves a formal perturbation series is used. The Wilton ripple phenomenon is reexamined and shown to be associated with a bifurcation in which a wave of permanent form can double its period. It is shown further that Wilton's ripples are a special case of a more general phenomenon in which bifurcation into subharmonics and factorial higher harmonics can occur. Numerical procedures for the calculation of waves of finite amplitude are developed. Bifurcation and limit lines are calculated. Pure and combination waves are continued to maximum amplitude. It is found that the height is limited in all cases by the surface enclosing one or more bubbles. Results for the shape of gravity waves are obtained by solving an integra-differential equation. It is found that the family of solutions giving the waveheight or equivalent parameter has bifurcation points. Two bifurcation points and the branches emanating from them are found specifically, corresponding to a doubling and tripling of the wavelength. Solutions on the new branches are calculated.

Some problems of edge waves and standing waves on beaches

Some problems of edge waves and standing waves on beaches are examined.

The nonlinear interaction of a wave normally incident on a sloping beach with a subharmonic edge wave is studied. A two-timing expansion is used in the full nonlinear theory to obtain the modulation equations which describe the evolution of the waves. It is shown how large amplitude edge waves are produced; and the results of the theory are compared with some recent laboratory experiments.

Traveling edge waves are considered in two situations. First, the full linear theory is examined to find the finite depth effect on the edge waves produced by a moving pressure disturbance. In the second situation, a Stokes' expansion is used to discuss the nonlinear effects in shallow water edge waves traveling over a bottom of arbitrary shape. The results are compared with the ones of the full theory for a uniformly sloping bottom.

The finite amplitude effects for waves incident on a sloping beach, with perfect reflection, are considered. A Stokes' expansion is used in the full nonlinear theory to find the corrections to the dispersion relation for the cases of normal and oblique incidence.

Finally, an abstract formulation of the linear water waves problem is given in terms of a self adjoint but nonlocal operator. The appropriate spectral representations are developed for two particular cases.

Theoretical investigation of turbulent boundary layer over a wavy surface

The important features of the two-dimensional incompressible turbulent flow over a wavy surface of wavelength comparable with the boundary layer thickness are analyzed.

A turbulent field method using model equation for turbulent shear stress similar to the scheme of Bradshaw, Ferriss and Atwell (1967) is employed with suitable modification to cover the viscous sublayer. The governing differential equations are linearized based on the small but finite amplitude to wavelength ratio. An orthogonal wavy coordinate system, accurate to the second order in the amplitude ratio, is adopted to avoid the severe restriction to the validity of linearization due to the large mean velocity gradient near the wall. Analytic solution up to the second order is obtained by using the method of matched-asymptotic-expansion based on the large Reynolds number and hence the small skin friction coefficient.

In the outer part of the layer, the perturbed flow is practically "inviscid." Solutions for the velocity, Reynolds stress and also the wall pressure distributions agree well with the experimental measurement. In the wall region where the perturbed Reynolds stress plays an important role in the process of momentum transport, only a qualitative agreement is obtained. The results also show that the nonlinear second-order effect is negligible for amplitude ratio of 0.03. The discrepancies in the detailed structure of the velocity, shear stress, and skin friction distributions near the wall suggest modifications to the model are required to describe the present problem.

Strongly nonlinear long interfacial waves in uniform basic flows and their dispersion relationship

In this paper, long interfacial waves of finite amplitude in uniform basic flows are considered with the assumption that the aspect ratio between wavelength and water depth is small. A new model is derived using the velocities at arbitrary distances from the still water level as the velocity variables instead of the commonly used depth-averaged velocities. This significantly improves the dispersion properties and makes them applicable to a wider range of water depths. Since its derivation requires no assumption on wave amplitude, the model thus can be used to describe waves with arbitrary amplitude.

Algorithm Studies on How to Obtain a Conditional Nonlinear Optimal Perturbation (CNOP)

The conditional nonlinear optimal perturbation (CNOP), which is a nonlinear generalization of the linear singular vector (LSV), is applied in important problems of atmospheric and oceanic sciences, including ENSO predictability, targeted observations, and ensemble forecast. In this study, we investigate the computational cost of obtaining the CNOP by several methods. Differences and similarities, in terms of the computational error and cost in obtaining the CNOP, are compared among the sequential quadratic programming (SQP) algorithm, the limited memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) algorithm, and the spectral projected gradients (SPG2) algorithm. A theoretical grassland ecosystem model and the classical Lorenz model are used as examples. Numerical results demonstrate that the computational error is acceptable with all three algorithms. The computational cost to obtain the CNOP is reduced by using the SQP algorithm. The experimental results also reveal that the L-BFGS algorithm is the most effective algorithm among the three optimization algorithms for obtaining the CNOP. The numerical results suggest a new approach and algorithm for obtaining the CNOP for a large-scale optimization problem.

Nonlinear Acoustics in Underwater and Biomedical Applications: Array Performance Degradation and Time Reversal Invariance

This dissertation describes a model for acoustic propagation in inhomogeneous flu- ids, and explores the focusing by arrays onto targets under various conditions. The work explores the use of arrays, in particular the time reversal array, for underwater and biomedical applications. Aspects of propagation and phasing which can lead to reduced focusing effectiveness are described. An acoustic wave equation was derived for the propagation of finite-amplitude waves in lossy time-varying inhomogeneous fluid media. The equation was solved numerically in both Cartesian and cylindrical geometries using the finite-difference time-domain (FDTD) method. It was found that time reversal arrays are sensitive to several debilitating factors. Focusing ability was determined to be adequate in the presence of temporal jitter in the time reversed signal only up to about one-sixth of a period. Thermoviscous absorption also had a debilitating effect on focal pressure for both linear and nonlinear propagation. It was also found that nonlinearity leads to degradation of focal pressure through amplification of the received signal at the array, and enhanced absorption in the shocked waveforms. This dissertation also examined the heating effects of focused ultrasound in a tissue-like medium. The application considered is therapeutic heating for hyperther- mia. The acoustic model and a thermal model for tissue were coupled to solve for transient and steady temperature profiles in tissue-like media. The Pennes bioheat equation was solved using the FDTD method to calculate the temperature fields in tissue-like media from focused acoustic sources. It was found that the temperature-dependence of the medium's background prop- erties can play an important role in the temperature predictions. Finite-amplitude effects contributed excess heat when source conditions were provided for nonlinear ef- fects to manifest themselves. The effect of medium heterogeneity was also found to be important in redistributing the acoustic and temperature fields, creating regions with hotter and colder temperatures than the mean by local scattering and lensing action. These temperature excursions from the mean were found to increase monotonically with increasing contrast in the medium's properties.

PROPAGATION OF SONIC BOOMS THROUGH A REAL, STRATIFIED ATMOSPHERE

Sonic boom propagation in a quiet) stratified) lossy atmosphere is the subject of this dissertation. Two questions are considered in detail: (1) Does waveform freezing occur? (2) Are sonic booms shocks in steady state? Both assumptions have been invoked in the past to predict sonic boom waveforms at the ground. A very general form of the Burgers equation is derived and used as the model for the problem. The derivation begins with the basic conservation equations. The effects of nonlinearity) attenuation and dispersion due to multiple relaxations) viscosity) and heat conduction) geometrical spreading) and stratification of the medium are included. When the absorption and dispersion terms are neglected) an analytical solution is available. The analytical solution is used to answer the first question. Geometrical spreading and stratification of the medium are found to slow down the nonlinear distortion of finite-amplitude waves. In certain cases the distortion reaches an absolute limit) a phenomenon called waveform freezing. Judging by the maturity of the distortion mechanism, sonic booms generated by aircraft at 18 km altitude are not frozen when they reach the ground. On the other hand, judging by the approach of the waveform to its asymptotic shape, N waves generated by aircraft at 18 km altitude are frozen when they reach the ground. To answer the second question we solve the full Burgers equation and for this purpose develop a new computer code, THOR. The code is based on an algorithm by Lee and Hamilton (J. Acoust. Soc. Am. 97, 906-917, 1995) and has the novel feature that all its calculations are done in the time domain, including absorption and dispersion. Results from the code compare very well with analytical solutions. In a NASA exercise to compare sonic boom computer programs, THOR gave results that agree well with those of other participants and ran faster. We show that sonic booms are not steady state waves because they travel through a varying medium, suffer spreading, and fail to approximate step shocks closely enough. Although developed to predict sonic boom propagation, THOR can solve other problems for which the extended Burgers equation is a good propagation model.

Broadband chaos generated by an optoelectronic oscillator.

We study an optoelectronic time-delay oscillator that displays high-speed chaotic behavior with a flat, broad power spectrum. The chaotic state coexists with a linearly stable fixed point, which, when subjected to a finite-amplitude perturbation, loses stability initially via a periodic train of ultrafast pulses. We derive approximate mappings that do an excellent job of capturing the observed instability. The oscillator provides a simple device for fundamental studies of time-delay dynamical systems and can be used as a building block for ultrawide-band sensor networks.

Ion-acoustic waves in a plasma consisting of adiabatic warm ions, non-isothermal electrons and a weakly relativistic electron beam: linear and higher-order nonlinear effects

The nonlinear propagation of finite amplitude ion acoustic solitary waves in a plasma consisting of adiabatic warm ions, nonisothermal electrons, and a weakly relativistic electron beam is studied via a two-fluid model. A multiple scales technique is employed to investigate the nonlinear regime. The existence of the electron beam gives rise to four linear ion acoustic modes, which propagate at different phase speeds. The numerical analysis shows that the propagation speed of two of these modes may become complex-valued (i.e., waves cannot occur) under conditions which depend on values of the beam-to-background-electron density ratio , the ion-to-free-electron temperature ratio , and the electron beam velocity v0; the remaining two modes remain real in all cases. The basic set of fluid equations are reduced to a Schamel-type equation and a linear inhomogeneous equation for the first and second-order potential perturbations, respectively. Stationary solutions of the coupled equations are derived using a renormalization method. Higher-order nonlinearity is thus shown to modify the solitary wave amplitude and may also deform its shape, even possibly transforming a simple pulse into a W-type curve for one of the modes. The dependence of the excitation amplitude and of the higher-order nonlinearity potential correction on the parameters , , and v0 is numerically investigated.

Dust ion acoustic solitons in a plasma with kappa-distributed electrons

Dust ion acoustic solitons in an unmagnetized dusty plasma comprising cold dust particles, adiabatic fluid ions, and electrons satisfying a kappa distribution are investigated using both small amplitude and arbitrary amplitude techniques. Their existence domain is discussed in the parameter space of Mach number M and electron density fraction f over a wide range of values of kappa. For all kappa > 3/2, including the Maxwellian distribution, negative dust supports solitons of both polarities over a range in f. In that region of parameter space solitary structures of finite amplitude can be obtained even at the lowest Mach number, the acoustic speed, for all kappa. These cannot be found from small amplitude theories. This surprising behavior is investigated, and it is shown that f(c), the value of f at which the KdV coefficient A vanishes, plays a critical role. In the presence of positive dust, only positive potential solitons are found. (C) 2010 American Institute of Physics. [doi: 10.1063/1.3400229]

Modeling relativistic soliton interactions in overdense plasmas: A perturbed nonlinear Schrodinger equation framework

We investigate the dynamics of localized solutions of the relativistic cold-fluid plasma model in the small but finite amplitude limit, for slightly overcritical plasma density. Adopting a multiple scale analysis, we derive a perturbed nonlinear Schrodinger equation that describes the evolution of the envelope of circularly polarized electromagnetic field. Retaining terms up to fifth order in the small perturbation parameter, we derive a self-consistent framework for the description of the plasma response in the presence of localized electromagnetic field. The formalism is applied to standing electromagnetic soliton interactions and the results are validated by simulations of the full cold-fluid model. To lowest order, a cubic nonlinear Schrodinger equation with a focusing nonlinearity is recovered. Classical quasiparticle theory is used to obtain analytical estimates for the collision time and minimum distance of approach between solitons. For larger soliton amplitudes the inclusion of the fifth-order terms is essential for a qualitatively correct description of soliton interactions. The defocusing quintic nonlinearity leads to inelastic soliton collisions, while bound states of solitons do not persist under perturbations in the initial phase or amplitude