AN ADAPTIVE MULTI-SCALE CONJUGATE GRADIENT METHOD FOR DISTRIBUTED PARAMETER ESTIMATION OF 2-D WAVE EQUATION

An adaptive multi-scale conjugate gradient method for distributed parameter estimations (or inverse problems) of wave equation is presented. The identification of the coefficients of wave equations in two dimensions is considered. First, the conjugate gradient method for optimization is adopted to solve the inverse problems. Second, the idea of multi-scale inversion and the necessary conditions that the optimal solution should be the fixed point of multi-scale inversion method is considered. An adaptive multi-scale inversion method for the inoerse problem is developed in conjunction with the conjugate gradient method. Finally, some numerical results are shown to indicate the robustness and effectiveness of our method.     

ON INTERACTION OF SHOCK AND SOUND WAVE （I）

This paper studies the interaction of shock and gradient wave (sound wave) of solutions to the system of inviscid isentropic gas dynamics as a model for the corresponding problems for nonlinear hyperbolic systems. The problem can be reduced to a boundary value problem in a wedged dormain, By using the method of constructing asymptotic solutions and Newton‘siteration process it is proved that if a weak shock hits a gradient wave, then the grandient wave will split into two gradient waves, while the shock continuses propagating. In this paper the author reduces the problem to a standard form and constructs asymptotic solution of the problem. The existence of the genuine solution will he given in the following paper.     

Continuous compression waves in the two-dimensional Riemann problem

The interaction between a plane shock wave in a plate and a wedge is considered within the framework of the nondissipative compressible fluid dynamic equations. The wedge is filled with a material that may differ from that of the plate. Based on the numerical solution of the original equations, self-similar solutions are obtained for several versions of the problem with an iron plate and a wedge filled with aluminum and for the interaction of a shock wave in air with a rigid wedge. The behavior of the solids at high pressures is approximately described by a two-term equation of state. In all the problems, a two-dimensional continuous compression wave develops as a wave reflected from the wedge or as a wave adjacent to the reflected shock. In contrast to a gradient catastrophe typical of one-dimensional continuous compression waves, the spatial gradient of a two-dimensional compression wave decreases over time due to the self-similarity of the solution. It is conjectured that a phenomenon opposite to the gradient catastrophe can occur in an actual flow with dissipative processes like viscosity and heat conduction. Specifically, an initial shock wave is transformed over time into a continuous compression wave of the same amplitude.     

ON INTERACTION OF SHOCK AND SOUND WAVE (Ⅰ)

1.IntroductionRecentlythestudyofdiscontinuoussolutionforthesystemofconservationlawsinhigherdimensionalspacehasbeenconsiderablydeveloped.In[1,W12]thelocalekistenceof8olutionforsuclisystemwithdiscontinuityinvolvingsingleshock,rarefactionwaveorsoundwav(gradientwave)hasbeenestablished.In[2Jand[14]theproblemsoninteractionoftwoshocksorinteractionofweaksingularitiesarealsoconsidered.Itisnaturaltoaskwhatabouttheresultwhenashockisinteractedbyawavewitliweakersillgularities,particularly,forthenbynsystem…     

On Identification of Memory Kernels in Linear Viscoelasticity

The inverse problem of identification of the memory kernel in the linear constitutive stress-strain-relation of Boltzmann type is reduced to an optimal control problem for an initial-boundary-value problem of the related wave equation for the displacement. For the control problem the existence of an optimal control is proved, where both classical and generalized solutions of the equation are dealt with. Further the existence and an expression for the gradient of the cost functional are derived.     

Recovering a space-dependent source term in a time-fractional diffusion wave equation

This work is concerned with identifying a space-dependent source function from noisy final time measured data in a time-fractional diffusion wave equation by a variational regularization approach. We provide a regularity of direct problem as well as the existence and uniqueness of adjoint problem. The uniqueness of the inverse source problem is discussed. Using the Tikhonov regularization method, the inverse source problem is formulated into a variational problem and a conjugate gradient algorithm is proposed to solve it. The efficiency and robust of the proposed method are supported by some numerical experiments.     

Non-linear peristaltic transport in an inclined asymmetric channel

The problem of peristaltic transport of a hydromagnetic (electrically conducting) viscous incompressible fluid in an inclined planar channel having electrically insulated walls has been investigated under long-wavelength and low-Reynolds number assumptions. The channel asymmetry is produced by choosing the peristaltic wave train on the walls to have different amplitudes and phase. The flow is investigated in a wave frame of reference moving with the velocity of the wave. Expressions for velocity field, shear stress and pressure gradient on the wall are obtained. The effects of different parameters entering into the problem are discussed numerically and explained graphically.     

Controllability method for acoustic scattering with spectral elements

We formulate the Helmholtz equation as an exact controllability problem for the time-dependent wave equation. The problem is then discretized in time domain with central finite difference scheme and in space domain with spectral elements. This approach leads to high accuracy in spatial discretization. Moreover, the spectral element method results in diagonal mass matrices, which makes the time integration of the wave equation highly efficient. After discretization, the exact controllability problem is reformulated as a least-squares problem, which is solved by the conjugate gradient method. We illustrate the method with some numerical experiments, which demonstrate the significant improvements in efficiency due to the higher order spectral elements. For a given accuracy, the controllability technique with spectral element method requires fewer computational operations than with conventional finite element method. In addition, by using higher order polynomial basis the influence of the pollution effect is reduced.     

Peristaltic transport of a Newtonian fluid in a vertical asymmetric channel with heat transfer and porous medium

The problem of peristaltic flow of a Newtonian fluid with heat transfer in a vertical asymmetric channel through porous medium is studied under long-wavelength and low-Reynolds number assumptions. The flow is examined in a wave frame of reference moving with the velocity of the wave. The channel asymmetry is produced by choosing the peristaltic wave train on the walls to have different amplitudes and phase. The analytical solution has been obtained in the form of temperature from which an axial velocity, stream function and pressure gradient have been derived. The effects of permeability parameter, Grashof number, heat source/sink parameter, phase difference, varying channel width and wave amplitudes on the pressure gradient, velocity, pressure drop, the phenomenon of trapping and shear stress are discussed numerically and explained graphically.     

A double smoothing technique for solving unconstrained nondifferentiable convex optimization problems

The aim of this paper is to develop an efficient algorithm for solving a class of unconstrained nondifferentiable convex optimization problems in finite dimensional spaces. To this end we formulate first its Fenchel dual problem and regularize it in two steps into a differentiable strongly convex one with Lipschitz continuous gradient. The doubly regularized dual problem is then solved via a fast gradient method with the aim of accelerating the resulting convergence scheme. The theoretical results are finally applied to an l ₁ regularization problem arising in image processing.     

The efficient solution of electromagnetic scattering for inhomogeneous media

We consider an electromagnetic scattering problem for inhomogeneous media. In particular, we focus on the numerical computation of the electromagnetic scattered wave generated by the interaction of an electromagnetic plane wave and an inhomogeneity in the corresponding propagation medium. This problem is studied in the VV polarization case, where some special symmetry requirements for the incident wave and for the inhomogeneity are assumed. This problem is reformulated as a Fredholm integral equation of second kind, which is discretized by a linear system having a special form. This allows to compute efficiently an approximate solution of the scattering problem by using iterative techniques for linear systems. Some numerical examples are reported.     

Optimal design of the damping set for the stabilization of the wave equation

We consider the problem of optimizing the shape and position of the damping set for the internal stabilization of the linear wave equation in RN, N=1,2. In a first theoretical part, we reformulate the problem into an equivalent non-convex vector variational one using a characterization of divergence-free vector fields. Then, by means of gradient Young measures, we obtain a relaxed formulation of the problem in which the original cost density is replaced by its constrained quasi-convexification. This implies that the new relaxed problem is well-posed in the sense that there exists a minimizer and, in addition, the infimum of the original problem coincides with the minimum of the relaxed one. In a second numerical part, we address the resolution of the relaxed problem using a first-order gradient descent method. We present some numerical experiments which highlight the influence of the over-damping phenomena and show that for large values of the damping potential the original problem has no minimizer. We then propose a penalization technique to recover the minimizing sequences of the original problem from the optimal solution of the relaxed one.     

Solving the Inverse Generalized Problem of Vertical Seismic Profiling

The dynamic inverse seismics problem is considered in a generalized setting. We investigate whether the wave propagation problem in a vertically nonhomogeneous medium is well-posed. We show that the regular part of the solution is an L ₂ function and the inverse problem, i.e., the determination of the reflection coefficient, is thus reducible to minimizing the error functional. The gradient of the functional is obtained in explicit form from the conjugate problem, and approximate formulas for its evaluation are derived. A regularization algorithm for the solution of the inverse problem is considered; simulation results using various excitation sources are reported.     

Martingale solutions of a stochastic wave equation with reflection

We discuss an initial boundary value problem for a one-dimensional stochastic wave equation with reflection. For stochastic parabolic equations with reflection, there are some well-known results. However, there seems to be no existence result for a stochastic wave equation with reflection. Even for a deterministic wave equation, the problem has not been completely resolved. Our goal is to establish the existence of a martingale solution for this problem.     

Asynchronous gradient algorithms for a class of convex separable network flow problems

We consider the single commodity strictly convex network flow problem. The dual of this problem is unconstrained, differentiable, and well suited for solution via distributed or parallel iterative methods. We present and prove convergence of gradient and asynchronous gradient algorithms for solving the dual problem. Computational results are given and analysed.     

On bounded solutions of a nonlinear Schrödinger equation on the half-line

We consider the problem on nonzero solutions of the Schrödinger equation on the half-line with potential that implicitly depends on the wave function via a nonlinear ordinary differential equation of the second order under zero boundary conditions for the wave function and the condition that the potential is zero at the beginning of the interval and its derivative is zero at infinity. The problem is reduced to the analysis and investigation of solutions of the Cauchy problem for a system of two nonlinear second-order ordinary differential equations with initial conditions depending on two parameters. We show that if the solution of the Cauchy problem for some parameter values can be extended to the entire half-line, then there exists a nonzero solution of the original problem with finitely many zeros.     

Stability of Rarefaction Wave to the 1-D Piston Problem for the Pressure-Gradient Equations

The 1-D piston problem for the pressure gradient equations arising from the flux-splitting of the compressible Euler equations is considered. When the total variations of the initial data and the velocity of the piston are both sufficiently small, the author establishes the global existence of entropy solutions including a strong rarefaction wave without restriction on the strength by employing a modified wave front tracking method.     

The impedance boundary-value problem of diffraction by a strip

We consider the impedance boundary-value problem for the Helmholtz equation originated by the problem of wave diffraction by an infinite strip with imperfect conductivity. The two possible different situations of real and complex wave numbers are considered. Bessel potential spaces are used to deal with the problem, and the identification of corresponding operators of single and double layer potentials allow a reformulation of the problem into a system of integral equations. The well-posedness of the problem is obtained for a set of impedance parameters (and wave numbers), after the incorporation of some compatibility conditions on the data. At the end, an improvement of the regularity of the solution is derived for the same set of parameters previously considered.     

Dynamic modeling and simulations of the wave glider

Wave glider is a new wave-powered autonomous marine vehicle, which is composed of a submerged glider connected to a surface floater via a tether. Such an advanced vessel is designed to harvest energy in ocean surface waves to generate forward thrust. Equipped with solar panel and battery as well as some dedicated sensors, the wave glider is able to achieve long duration missions via sea-side control. In this paper, a 4-DOF (degree-of-freedom) mathematical model of the wave glider is established using Newton–Euler approach. The second-order wave drift force on the horizontal plane and the first order wave force on the vertical direction are considered. The hydrodynamic parameters were calculated using the potential flow theory and empirical formula. Furthermore, motion simulation of the wave glider with respect to the sensitivity analysis to some key environmental factors and the heading control ability is conducted. The simulation results are presented and discussed in detail, which provides theoretical guidance and reference for wave glider design.     

Effect of Hall currents on interaction of pulsatile and peristaltic transport induced flows of a particle-fluid suspension

The interaction of purely periodic mean flow with a peristaltic induced flow is investigated within the framework of a two-dimensional analogue. The mathematical model considers a viscous incompressible fluid under the effect of transverse magnetic field, taking into account the effect of Hall currents for a magneto-fluid with suspended particles between infinite parallel walls on which a sinusoidal traveling wave is imposed. A perturbation solution to the complete set of Navier-Stokes equations is found for the case in which the frequency of the traveling wave and that of the imposed pressure gradient are equal. The ratio of the traveling wave amplitude to channel width is assumed to be small. For this case a first order steady flow is found to exist, as contrasted to a second order effect in the absence of the imposed periodic pressure gradient. The effect of Hall parameter, Hartmann number and the various parameters included in the problem are discussed numerically.