Passive systems theory with narrow-band and linear constraints: Part II -- Temporal diversity

This paper is part of a series of three papers studying passive tracking problems arising in navigation and positioning applications. The basic question here lies with the determination of the position and dynamics of a point source being tracked by an omnidirectional observer, through demodulation of the Doppler effect induced on the radiated signals by the relative motions. A simple model, fitting a finite parameter nonlinear estimation context, is developed, the receiver designed, and its mean-square error performance studied. It is shown that, besides the speed and angle estimation, simultaneous global range passive tracking is possible. The signal model precludes range acquisition from synchronous measurement of the absolute phase reference: the global range estimation is attained by processing the higher order temporal modulations (varying Doppler). Quantifying the statistical and geometric performance tradeoffs, the work presents simple expressions and graphical displays that can be used as design tools in practical passive tracking problems. A subsequent paper considers the space/ time coupling issues, generalizing the study to the context where a moving source is tracked by a directional array.     

Passive systems theory with narrow-band and linear constraints: Part I--Spatial diversity

This paper studies passive problems where the receiver extracts from the source radiated signature information concerning the parameters defining the relative source/receiver geometry. A model encompassing the fundamental global and local characteristics for passive positioning and navigation is presented. It considers narrow-band signals, imposes linear constraints on the geometry, and exhibits explicitly the symmetry between the space and time aspects. The analysis concentrates on questions of global geometry identifiability, emphasizing the passive global range acquisition. The maximum-likelihood processor is analyzed by studying the ambiguity structure associated with inhomogeneous passive narrowband tracking. Bounds on the global and local mean-square error performance are studied and tested via Monte Carlo simulations. By considering two limiting geometries, a distant and a close observer, simple approximate expressions for the mean-square errors are presented and compared to the exact bounds. Herein the study is restricted to stationary geometries where the source is located by an extended array (spatial diversity). Subsequent papers generalize the study to moving sources (temporal diversity) and to coupled geometries.     

基于相空间重构与模糊神经网络的海温垂直分布预测模型

在海温预报中引入混沌理论,将相空间重构与模糊神经网络相结合,提出了海温垂直建模预测模型。通过相空间重构,把海温时间序列拓展为多维序列,而多维序列包含着各态历经的信息,从而挖掘出了丰富的海温变化空间的信息,有利于模糊神经网络的训练。利用建立好的模糊神经网络模型,较好地对海温的垂直结构进行了建模、训练和预测。实际的预测结果表明,该模型预报精度较高,超前1~5个月的预测值的相对误差均控制在10%以内,预测结果可以为业务工作提供一定的参考与借鉴。     

3D空间航片纹理源的纹理映射问题

随着数字三维城市建设的快速推进,传统的纹理获取方法已不能满足快速建模的要求。提出了一种基于航片纹理源的3D空间纹理映射方法,基于摄影测量学相关理论,通过设计算法,实现了建筑物侧面纹理的自动匹配。实验证明,这种方法简单、高效,适用于大面积建筑物3D建模。     

On the Development of a Second-Order Bistatic Radar Cross Section of the Ocean Surface: A High-Frequency Result for a Finite Scattering Patch

The development of a model for the second-order bistatic high-frequency (HF) radar cross section on an ocean surface patch remote from the transmitter and receiver is addressed. A new approach is taken that allows a direct comparison with existing monostatic cross sections for finite regions of the ocean surface. The derivation starts with a general expression for the bistatically received second-order electric field in which the scattering surface is assumed to be of small height and slope. The source field is taken to be that of a vertically polarized dipole, and it is assumed that the ocean surface can be described, as is usually done, by a Fourier series in which the coefficients are zero-mean Gaussian random variables. Subsequently, a bistatic cross section of the surface, normalized to patch area, is derived. The result is verified by the following two means: 1) the complete form of the bistatic HF radar cross section in backscattering case is shown to contain an earlier monostatic result that has, itself, been used extensively in radio oceanography applications; and 2) the bistatic electromagnetic coupling coefficient is shown to reduce exactly to the monostatic result when backscattering geometry is imposed. The model is also depicted and discussed based on simulated data     

A modeling study of acoustic propagation through movingshallow-water solitary wave packets

Propagation of 400-Hz sound through continental-shelf internal solitary wave packets is shown by numerical simulation to be strongly influenced by coupling of normal modes. Coupling in a packet is controlled by the mode coefficients at the point where sound enters the packet, the dimensions of the waves and packet, and the ambient depth structures of temperature and salinity. In the case of a moving packet, changes of phases of the incident modes with respect to each other dominate over the other factors, altering the coupling over time and thus inducing signal fluctuations. The phasing within a moving packet varies with time scales of minutes, causing coupling and signal fluctuations with comparable time scales. The directionality of energy flux between high-order acoustic modes and (less attenuated) low-order modes determines a gain factor for long-range propagation. A significant finding is that energy flux toward low-order modes through the effect of a packet near a source favoring high-order modes will give net amplification at distant ranges. Conversely, a packet far from a source sends energy into otherwise quiet higher modes. The intermittency of the coupling and of high-mode attenuation via bottom interaction means that signal energy fluctuations and modal diversity fluctuations at a distant receiver are complementary, with energy fluctuations suggesting a source-region packet and mode fluctuations suggesting a receiver-region packet. Simulations entailing 33-km propagation are used in the analyses, imitating the SWARM experiment geometry, allowing comparison with observations     

潜艇操纵控制方法的现状与发展

在大量相关文献的基础上,对潜艇空间机动中出现的非线性、两平面运动间的强耦合、参数的时变特性以及近水面时的定深控制等潜艇操纵控制的难点问题及其处理方法进行了分析、探讨和综述。     

Multistatic Sonar Localization

This paper develops bistatic contact-localization statistics as a function of source-target-receiver geometry, source and receiver location errors, sound-speed errors, time and bearing errors, and array-heading errors. We verify the analysis through Monte Carlo simulation and illustrate localization accuracy for some representative choices of source-target-receiver geometry and measurement error assumptions. We examine the improved localization accuracy of multistatic contacts that result from the fusion of multiple sonar contacts     

Design of a 3D Chirp Sub-bottom Imaging System

Chirp sub-bottom profilers are marine acoustic devices that use a known and repeatable source signature (1–24 kHz) to produce decimetre vertical resolution cross-sections of the sub-seabed. Here the design and development of the first true 3D Chirp system is described. When developing the design, critical factors that had to be considered included spatial aliasing, and precise positioning of sources and receivers. Full 3D numerical modelling of the combined source and receiver directivity was completed to determine optimal source and receiver geometries. The design incorporates four source transducers (1.5–13 kHz) that can be arranged into different configurations, including Maltese Cross, a square and two separated pairs. The receive array comprises 240 hydrophones in 60 groups whose group-centres are separated by 25 cm in both horizontal directions, with each hydrophone group containing four individual elements and a pre-amplifier. After careful consideration, it was concluded that the only way to determine with sufficient accuracy the source–receiver geometry, was to fix the sources and receivers within a rigid array. Positional information for the array is given by a Real Time Kinematic GPS and attitude system incorporating four antennas to give position, heading, pitch and roll. It is shown that this system offers vertical positioning accuracy with a root-mean-square (rms) error less than 2.6 cm, while the horizontal positioning rms error was less than 2.0 cm. The system is configured so that the Chirp source signature can be chosen by software aboard the acquisition vessel. The complete system is described and initial navigational and seismic data results are presented. These data demonstrate that the approach of using fixed source-receiver geometry combined with RTK navigation can provide complete 3D imaging of the sub-surface.     

Parametric Geometric Model and Hydrodynamic Shape Optimization of A Flying-Wing Structure Underwater Glider

Combining high precision numerical analysis methods with optimization algorithms to make a systematic exploration of a design space has become an important topic in the modern design methods. During the design process of an underwater glider''s flying-wing structure, a surrogate model is introduced to decrease the computation time for a high precision analysis. By these means, the contradiction between precision and efficiency is solved effectively. Based on the parametric geometry modeling, mesh generation and computational fluid dynamics analysis, a surrogate model is constructed by adopting the design of experiment (DOE) theory to solve the multi-objects design optimization problem of the underwater glider. The procedure of a surrogate model construction is presented, and the Gaussian kernel function is specifically discussed. The Particle Swarm Optimization (PSO) algorithm is applied to hydrodynamic design optimization. The hydrodynamic performance of the optimized flying-wing structure underwater glider increases by 9.1%.     

Correlation-based decision-feedback equalizer for underwater acoustic communications

The purpose of this paper is to develop a decision-feedback equalizer (DFE) using a fixed set of parameters applicable to most shallow oceans with minimal user supervision (i.e., a turn key system). This work is motivated by the superior performance [bit error rate (BER)] of the multichannel DFE compared with other methods, such as passive-phase conjugation (PPC), at the same time noting its sensitivity to different acoustic environments. The approach is to couple PPC, utilizing its adaptability to different environments, with a single-channel DFE. This coupling forms an optimal processor for acoustic communications in theory, but it has never been implemented in practice. By coupling with DFE, the method achieves the same spatial diversity as conventional multichannel DFE, without requiring a large number of receivers as does PPC. The correlation-based DFE in terms of the autocorrelation functions of the channel impulse responses summed over the receiver channels (the Q function) is derived. This paper shows in terms of waveguide physics, further supported by real data, the many desirable features of the Q function that suggest, given adequate sampling of the water column, a general applicability of the correlation-based equalizer to different environments, irrespective of the sound speed profiles, bottom properties, and source-receiver ranges/depths. This property can be expected to hold approximately for a small number of receivers with spatial diversity. This paper demonstrates the robustness of the new equalizer with moving source data despite the range change (which modifies the impulse response) and symbol phase change due to time-varying Doppler.     

A spectral barotropic model of the wind-driven world ocean

A global spectral barotropic ocean model is introduced to describe the depth-averaged flow. The equations are based on vorticity and divergence (instead of horizontal momentum); continents exert a nearly infinite drag on the fluid. The coding follows that of spectral atmospheric general circulation models using triangular truncation and implicit time integration to provide a first step for seamless coupling to spectral atmospheric global circulation models and an efficient method for filtering of ocean wave dynamics. Five experiments demonstrate the model performance: (i) Bounded by an idealized basin geometry and driven by a zonally uniform wind stress, the ocean circulation shows close similarity with Munk’s analytical solution. (ii) With a real land–sea mask the model is capable of reproducing the spin-up, location and magnitudes of depth-averaged barotropic ocean currents. (iii) The ocean wave-dynamics of equatorial waves, excited by a height perturbation at the equator, shows wave dispersion and reflection at eastern and western coastal boundaries. (iv) The model reproduces propagation times of observed surface gravity waves in the Pacific with real bathymetry. (v) Advection of tracers can be simulated reasonably by the spectral method or a semi-Langrangian transport scheme. This spectral barotropic model may serve as a first step towards an intermediate complexity spectral atmosphere–ocean model for studying atmosphere–ocean interactions in idealized setups and long term climate variability beyond millennia.     

A WKB analysis of the surface signature and vertical structure of long extratropical baroclinic Rossby waves over topography

WKB theory has been used recently to construct global approximations of purely periodic long extratropical baroclinic Rossby waves in a continuously stratified ocean, with the goal of building a global theoretical framework that can serve to interpret observed features of the waves, such as sea surface height wave activity and anomalous propagation. This study adopts the same approach in the idealized context of an ocean with constant buoyancy frequency N and longitudinally varying Gaussian topography, to gain insight into issues that have received little or no attention so far, namely: (a) the nature of the links between the vertical structure of the wave field and its surface signature, and the extent to which constraints on the interior dynamics can be derived solely from observing the surface; (b) which of the phase/amplitude variations determine the visual impression of westward propagation, addressed by constructing longitude/time plots of the signal in the attempt to mimic more closely the way the satellite SSH data of TOPEX/Poseidon are analyzed; (c) the validity and accuracy of the classical leading order WKB theory, addressed by estimating the residual a posteriori by computing the next-order term of the formal WKB series expansion. To a lesser extent, this also serves to assess the validity of the planetary geostrophic equations used to describe the dynamics of long Rossby waves.The main results are that: (a) faster propagation is unambiguously related to the surface intensification of the waves, while slower propagation is associated with a vertical structure intermediate between that of the first and second standard baroclinic modes; (b) westward propagation is dominated by the phase variations; (c) the residual is inversely proportional to the frequency (or equivalently the wavelength by dividing by the phase speed), and is found to vary strongly with position. The hilltop is where the residual is the highest, and hence where WKB is the most likely to breakdown, in agreement with recent published predictions of mode conversion theory in a two-layer model setting.     

Methods for ultimate limit state assessment of ships and ship-shaped offshore structures: Part II stiffened panels

The present paper is Part II of a series of three papers on methods useful for the ultimate limit state assessment of ships and ship-shaped offshore structures. In contrast to Part I [Paik et al., 2007a. Methods for ultimate limit state assessment of ships and ship-shaped offshore structures: Part I unstiffened plates, Ocean Engineering, doi:10.1016/j.oceaneng.2007.08.004] that deals with unstiffened plates, the present paper (Part II) is focused on methods for the ultimate limit state assessment of stiffened plate structures under combined biaxial compression and lateral pressure actions. The object structure is the bottom part of an AFRAMAX-class hypothetical double-hull oil tanker structure designed by IACS common structural rules (CSR) method, that is the same ship studied in Part I. Three candidate methods, namely ANSYS nonlinear finite element method, DNV PULS method, and ALPS/ULSAP method, are employed for the present study. The results and insights developed from the present study are summarized in terms of ultimate strength characteristics of bottom-stiffened plate structures.     

A spectral barotropic model of the wind-driven world ocean

A global spectral barotropic ocean model is introduced to describe the depth-averaged flow. The equations are based on vorticity and divergence (instead of horizontal momentum); continents exert a nearly infinite drag on the fluid. The coding follows that of spectral atmospheric general circulation models using triangular truncation and implicit time integration to provide a first step for seamless coupling to spectral atmospheric global circulation models and an efficient method for filtering of ocean wave dynamics. Five experiments demonstrate the model performance: (i) Bounded by an idealized basin geometry and driven by a zonally uniform wind stress, the ocean circulation shows close similarity with Munk’s analytical solution. (ii) With a real land–sea mask the model is capable of reproducing the spin-up, location and magnitudes of depth-averaged barotropic ocean currents. (iii) The ocean wave-dynamics of equatorial waves, excited by a height perturbation at the equator, shows wave dispersion and reflection at eastern and western coastal boundaries. (iv) The model reproduces propagation times of observed surface gravity waves in the Pacific with real bathymetry. (v) Advection of tracers can be simulated reasonably by the spectral method or a semi-Langrangian transport scheme. This spectral barotropic model may serve as a first step towards an intermediate complexity spectral atmosphere–ocean model for studying atmosphere–ocean interactions in idealized setups and long term climate variability beyond millennia.     

Study of the mutual influence of the El Niño-Southern Oscillation processes and the North Atlantic and Arctic Oscillations

On the basis of the nonlinear techniques for the estimation of coupling between oscillatory systems from time series, we investigate the dynamics of climatic modes characterizing global and Northern Hemisphere (NH) processes. The North Atlantic Oscillation (NAO) and Arctic Oscillation indices and the El Niño-Southern Oscillation (ENSO) indices are analyzed in terms of the most reliable data from 1950 through 2004 and earlier data since the 19th century. These indices characterize changes in NH atmospheric pressure (specifically, sea-level pressure in the North Atlantic and NH extratropical latitudes as a whole) and in equatorial Pacific sea-surface temperature and sea-level pressure to which the strongest variations of global surface temperature and global climate on interannual time scales and of regional climatic anomalies in the NH are linked. The methods used are based on phase-dynamics modeling and nonlinear prediction models (a nonlinear version of Granger causality). From both methods and various ENSO indices, the inference about the ENSO effect on the NAO during the latter half of the 20th century and in the early 21st century is made with confidence probability of at least 0.95. The influence is characterized by a time delay of about two years. No inverse influence is found with a similar degree of reliability. Results of estimating the coupling between the ENSO and the NAO depend on the type of index that is used to describe the NAO. The ENSO effect on the NAO is detected with sufficient confidence when the NAO index is chosen to be a larger scale characteristic. However, when a more local index of the NAO is used, no statistically significant coupling to the ENSO is found. Increasing the length of the analyzed ENSO and NAO series (over more than 100 yr) does not lead to any more reliable detection of coupling. Analysis of the data for different time intervals during 1950–2004 has revealed a strengthening of the ENSO effect on the NAO, although this inference is not reliable.     

Methods for ultimate limit state assessment of ships and ship-shaped offshore structures: Part I—Unstiffened plates

The present paper is Part I of a series of three papers prepared by the authors on the methods useful for ultimate limit state assessment of marine structures, that have been developed in the literature during the last few decades. It is considered that such methods are now mature enough to enter day-by-day design and strength assessment practice. The aims of the three papers are to conduct some benchmark studies of such methods on ultimate limit state assessment of (unstiffened) plates, stiffened panels, and hull girders of ships and ship-shaped offshore structures, using some candidate methods such as ANSYS nonlinear finite element analysis (FEA), DNV PULS, ALPS/ULSAP, ALPS/HULL, and IACS common structural rules (CSR) methods. As an illustrative example, an AFRAMAX-class hypothetical double hull oil tanker structure designed by CSR method is studied. In the present paper (Part I), the ultimate limit state assessment of unstiffened plates under combined biaxial compression and lateral pressure loads is emphasized using ANSYS, DNV PULS, and ALPS/ULSAP methods, and their resulting computations are compared. Part II will deal with methods for the ultimate limit state assessment of stiffened panels under combined biaxial compression and lateral pressure using ANSYS, DNV PULS, and ALPS/ULSAP methods, and Part III will treat methods for the progressive collapse analysis of the hull structure using ANSYS, ALPS/HULL, and IACS CSR methods.     

Height variability from the MIROC – IPCC model for the 20th century compared to that of the TOPEX/POSEIDON altimeter

Planetary waves are key to large-scale dynamical adjustment in the global ocean as they transfer energy from the east to the west side of oceanic basins; they connect the forcing in the ocean interior with the variability at its boundaries; and they change the local heat content, thus coupling oceanic, atmospheric, and biological processes. Planetary waves, mostly of the first baroclinic mode, are observed as distinctive patterns in global time series of sea surface height anomaly (SSHA) and heat storage. The goal of this study is to compare and validate large-scale SSHA signals from coupled ocean-atmosphere general circulation Model for Interdisciplinary Research on Climate (MIROC) with TOPEX/POSEIDON satellite altimeter observations. The last decade of the models’ time series is selected for comparison with the altimeter data. The wave patterns are separated from the meso- and large-scale SSHA signals by digital filters calibrated to select the same spectral bands in both model and altimeter data. The band-wise comparison allows for an assessment of the model skill to simulate the dynamical components of the observed wave field. Comparisons regarding both the seasonal cycle and the Rossby wave field differ significantly among basins. When carried within the same basin, differences can occur between equal latitudes in opposite hemispheres. Furthermore, at some latitudes the MIROC reproduces biannual, annual and semiannual planetary waves with phase speeds and average amplitudes similar to those observed by the altimeter, but with significant differences in phase.     

Second-order coupling of numerical and physical wave tanks for 2D irregular waves. Part II: Experimental validation in two-dimensions

This paper provides an experimental validation of the second-order coupling theory outlined by Yang et al. (Z. Yang, S. Liu, H.B. Bingham and J. Li., 2013. Second-order coupling of numerical and physical wave tanks for 2D irregular waves. Part I: Formulation, implementation and numerical properties, submitted for publication) using 2D irregular waves. This work provides a second-order dispersive correction for the physical wavemaker signal which improves the nonlinear transfer of information between the numerical and physical models compared to the first-order method of Zhang et al. (2007). The important nonlinear parameters and numerical performance were theoretically investigated in Part I. In the present Part II, careful experimental validation is carried out using a sequence of progressively more complex analytical and numerical target waves. The results demonstrate clearly that improved performance is achieved by using the second-order correction. When controlling with a second-order coupling signal, two key points are notable: (i) The higher harmonics underlying the numerical waves are accurately captured and transferred into the physical model. (ii) The second-order behavior leads to an unwanted spurious freely propagating second harmonic that is substantially reduced when compared to an identical wave paddle operating with a first-order coupling signal. Using nonlinear regular (monochromatic), bi-chromatic and irregular wave cases as well as varying coupled wave tank bathymetries, both these aspects are verified over a broad range of wave frequencies and shown to be extensively applicable to physical wave tanks.