Shallow water tidal constituents in the Bohai Sea and the Yellow Sea from a numerical adjoint model with TOPEX/POSEIDON altimeter data

A numerical adjoint model with TOPEX/POSEIDON (T/P) altimeter data was set up to investigate the shallow water tidal constituents in the Bohai Sea and the Yellow Sea. Shallow water tidal constituents W-4, MS4 and M-6) in the Bohai Sea and the Yellow Sea were first extracted from nearly 10 years of T/P data and then assimilated into a nonlinear barotropic tidal model by using adjoint method in order to fully describe the tides in this area. The general patterns of M-4 and MS4 solutions were in good agreement with those of Kang et al. (Cont. Shelf. Res. IS (1998) 739.) and Lefevre et al., (J. Geophys. Res. 105 (2000) 8707.). The RMS values for the principal constituents and coastal constituents were obviously less than those calculated by Kang et al. (1998) and Lefevre et al. (2000). It was shown that the calculated tidal constituents charts obtained in the present study were more accurate than those in other models. In the future the model will be applied to other coastal areas and some semi-enclosed seas. (C) 2004 Elsevier Ltd. All rights reserved.

Numerical modeling of water wave interaction with a soft mud bed

A vertical 2-D numerical model is presented for simulating the interaction between water waves and a soft mud bed. Taking into account nonlinear rheology, a semi-empirical rheological model is applied to this water-mud model, reflecting the combined visco-elasto-plastic properties of soft mud under such oscillatory external forces as water waves. In order to increase the resolution of the flow in the neighborhood of both sides of the inter-surface, a logarithmic grid in the vertical direction is employed for numerical treatment. Model verifications are given through comparisons between the calculated and the measured mud mass transport velocities as well as wave height changes.

Numerical wave flume study on wave motion around submerged plates

Nonlinear interaction between surface waves and a submerged horizontal plate is investigated in the absorbed numerical wave flume developed based on the volume of fluid (VOF) method. The governing equations of the numerical model are the continuity equation and the Reynolds-Averaged Navier-Stokes (RANS) equations with the k-epsilon turbulence equations. Incident waves are generated by an absorbing wave-maker that eliminates the waves reflected from structures. Results are obtained for a range of parameters, with consideration of the condition under which the reflection coefficient becomes maximal and the transmission coefficient minimal. Wave breaking over the plate, vortex shedding downwave, and pulsating flow below the plate are observed. Time-averaged hydrodynamic force reveals a negative drift force. All these characteristics provide a reference for construction of submerged plate breakwaters.

Numerical simulation of principal tidal constituents in the South China Sea, Gulf of Tonkin and Gulf of Thailand

The principal tidal constituents M-2, S-2, K-1 and O-1 in the South China Sea, Gulf of Tonkin and Gulf of Thailand are simulated simultaneously using the numerical scheme of Kwok et al. (1995 Proceedings of the 1st Asian Computational Fluid Dynamics Conference, pp. 16-19). The average differences between the computed and observed harmonic constants are mostly within 5 cm and 10 degrees for amplitudes and phase-lags, respectively. The simulated tidal regimes in the present model are believed to be more accurate than the previous numerical results. Our studies confirm that a clockwise rotating M-2 amphidromic system lies in the southeast of the Gulf of Thailand and an S-2 amphidromic system at the near-shore area of the northeast South China Sea. The linear tidal energy equation developed by Garrett (1975 Deep-Sea Research 22, 23-35) is generalized to the nonlinear case. Based on the numerical results, the energy budgets in the South China Sea and its subareas, namely the Taiwan Strait, the Gulf of Tonkin, the Gulf of Thailand and the remaining area are investigated. The tidal motion in the Taiwan Strait is maintained mainly by the energy fluxes from the East China Sea for both semidiurnal and diurnal species and partially from the Luzon Strait for semidiurnal species. For the other parts of the South China Sea, the tidal motion is mainly maintained by the energy fluxes through the Luzon Strait. The energy inputs from the tide-generating force are negative for semidiurnal species and positive for diurnal species. (C) 1999 Elsevier Science Ltd. All rights reserved.

A HYBRID METHOD OF FRACTIONAL STEPS WITH PREDICTOR-CORRECTOR DIFFERENCE-PSEUDOSPECTRUM FOR NUMERICAL-SOLUTION OF THE CONVECTION-DOMINATED FLOW PROBLEMS

A fractional-step method of predictor-corrector difference-pseudospectrum with unconditional L(2)-stability and exponential convergence is presented. The stability and convergence of this method is strictly proved mathematically for a nonlinear convection-dominated flow. The error estimation is given and the superiority of this method is verified by numerical test.

Applications of improved 2-D numerical wave tanks in wave breaking

Song and Banner (2002, henceforth referred to as SB02) used a numerical wave tank (developed by Drimer and Agnon, and further refined by Segre, henceforth referred to as DAS) to study the wave breaking in the deep water, and proposed a dimensionless breaking threshold that based on the behaviour of the wave energy modulation and focusing during the evolution of the wave group. In this paper, two modified DAS models are used to further test the SB02's results, the first one (referred to MDAS1) corrected many integral calculation errors appeared in the DAS code, and the second one (referred to MDAS2) replaced the linear boundary element approximation of DAS into the cubic element on the free surface. Researches show that the results of MDAS1 are the same with those of DAS for the simulations of deep water wave breaking, but, the different values of the wavemaker amplitude, the breaking time and the maximum local average energy growth rate delta(max) for the marginal breaking cases are founded by MDAS2 and MDAS1. However, MDAS2 still satisfies the SB02' s breaking threshold. Furthermore, MDAS1 is utilized to study the marginal breaking case in the intermediate water depth when wave passes over a submerged slope, where the slope is given by 1 : 500, 1 : 300, 1 : 150 or 1 : 100. It is found that the maximum local energy density U increases significantly if the slope becomes steeper, and the delta(max) decreases weakly and increases intensively for the marginal recurrence case and marginal breaking case respectively. SB02's breaking threshold is still valid for the wave passing over a submerged slope gentler than 1 : 100 in the intermediate water depth.

Numerical Simulation of Waves Generated by Seafloor Movements

Waves generated by vertical seafloor movements are simulated by use of a fully nonlinear two-dimensional numerical wave tank. In the source region, the seafloor lifts to a designated height by a generation function. The numerical tests show that file linear theory is only valid for estimating the wave behaviors induced by the seafloor movements with a small amplitude, and the fully nonlinear numerical model should be adopted in the simulation of the wave generation by the large amplitude seafloor movements. Without the background surface waves, many numerical tests on the stable maximum elevations eta(max)(0) are carried out by both the linear theory and the fully nonlinear model. The results of two models are compared and analyzed. For the fully nonlinear model, the influences of the amplitudes and the horizontal lengths on eta(max)(0) are stronger than that of the characteristic duration times. Furthermore, results reveal that there are significant differences between the linear theory and the fully nonlinear model. When the influences of the background surface waves are considered, the corresponding numerical analyses reveal that with the fully nonlinear model the eta(max)(0) near-linearly varies with the wave amplitudes of the surface waves, and the eta(max)(0) has significant dependences on the wave lengths and the wave phases of the surface waves. In addition, the differences between the linear theory and the fully nonlinear model are still obvious, aid these differences are significantly affected by The wave parameters of the background surface waves, such as the wave amplitude, the wave length and the wave phase.

Numerical simulation and validation of ocean waves measured by an Along-Track Interferometric Synthetic Aperture Radar

A new nonlinear integral transform of ocean wave spectra into Along-Track Interferometric Synthetic Aperture Radar (ATI-SAR) image spectra is described. ATI-SAR phase image spectra are calculated for various sea states and radar configurations based on the nonlinear integral transform. The numerical simulations show that the slant range to velocity ratio (R/V), significant wave height to ocean wavelength ratio (H-s/lambda), the baseline (2B) and incident angle (theta) affect ATI-SAR imaging. The ATI-SAR imaging theory is validated by means of Two X-band, HH-polarized ATI-SAR phase images of ocean waves and eight C-band, HH-polarized ATI-SAR phase image spectra of ocean waves. It is shown that ATI-SAR phase image spectra are in agreement with those calculated by forward mapping in situ directional wave spectra collected simultaneously with available ATI-SAR observations. ATI-SAR spectral correlation coefficients between observed and simulated are greater than 0.6 and are not sensitive to the degree of nonlinearity. However, the ATI-SAR phase image spectral turns towards the range direction, even if the real ocean wave direction is 30 degrees. It is also shown that the ATI-SAR imaging mechanism is significantly affected by the degree of velocity bunching nonlinearity, especially for high values of R/V and H-s/lambda.

Factors influencing the climatological mixed layer depth in the South China Sea: numerical simulations

The mixed layer depth (MLD) in the upper ocean is an important physical parameter for describing the upper ocean mixed layer. We analyzed several major factors influencing the climatological mixed layer depth (CMLD), and established a numerical simulation in the South China Sea (SCS) using the Regional Ocean Model System (ROMS) with a high-resolution (1/12A degrees x1/12A degrees) grid nesting method and 50 vertical layers. Several ideal numerical experiments were tested by modifying the existing sea surface boundary conditions. Especially, we analyzed the sensitivity of the results simulated for the CMLD with factors of sea surface wind stress (SSWS), sea surface net heat flux (SSNHF), and the difference between evaporation and precipitation (DEP). The result shows that of the three factors that change the depth of the CMLD, SSWS is in the first place, when ignoring the impact of SSWS, CMLD will change by 26% on average, and its effect is always to deepen the CMLD; the next comes SSNHF (13%) for deepening the CMLD in October to January and shallowing the CMLD in February to September; and the DEP comes in the third (only 2%). Moreover, we analyzed the temporal and spatial characteristics of CMLD and compared the simulation result with the ARGO observational data. The results indicate that ROMS is applicable for studying CMLD in the SCS area.

Numerical investigation of nutrient limitations in the Bohai Sea

River discharges are the important freshwater and nutrient sources for Bohai Sea (BS), and have a profound impact on the local marine environment. In this paper, the annual cycles of nutrient and phytoplankton dynamics in 1980s were reproduced using a coupled biogeochemical-physical model. Based on the validated simulations, the nutrient limitation characters were further investigated by running the model with the riverine nutrient altered, first enriching nitrogen and then phosphorus. It was found that although the riverine N:P ratios in Yellow and Haihe Rivers were much higher than the Redfield number, the nitrogen enrichment was still able to enhance the algae bloom in Laizhou and Bohai Bays. On the other hand, the response of algae growth to phosphorus enrichment was not thus obvious, which suggests that the local phytoplankton dynamics was characterized by the nitrogen limitation. Simulations also show that the nitrogen enrichment is generally accompanied by the phosphorus consumption, so a shift from nitrogen limitation to phosphorus limitation may occur if such a trend continues. (C) 2010 Elsevier Ltd. All rights reserved.

Numerical simulation on the process of saltwater intrusion and its impact on the suspended sediment concentration in the Changjiang (Yangtze) estuary

To study the relationship between sediment transportation and saltwater intrusion in the Changjiang (Yangtze) estuary, a three-dimensional numerical model for temperature, salinity, velocity field, and suspended sediment concentration was established based on the ECOMSED model. Using this model, sediment transportation in the flood season of 2005 was simulated for the Changjiang estuary. A comparison between simulated results and observation data for the tidal level, flow velocity and direction, salinity and suspended sediment concentration indicated that they were consistent in overall. Based on model verification, the simulation of saltwater intrusion and its effect on sediment in the Changjiang estuary was analyzed in detail. The saltwater intrusion in the estuary including the formation, evolution, and disappearance of saltwater wedge and the induced vertical circulation were reproduced, and the crucial impact of the wedge on cohesive and non-cohesive suspended sediment distribution and transportation were successfully simulated. The result shows that near the salinity front, the simulated concentrations of both cohesive and non-cohesive suspended sediment at the surface layer had a strong relationship with the simulated velocity, especially when considering a 1-hour lag. However, in the bottom layer, there was no obvious correlation between them, because the saltwater wedge and its inducing vertical circulation may have resuspended loose sediment on the bed, thus forming a high-concentration area near the bottom even if the velocity near the bottom was very low during the transition phase from flood to ebb.

Numerical study on the velocity structure around tidal fronts in the Yellow Sea

The velocity components across tidal fronts are examined using the Blumberg and Mellor 3-D nonlinear numerical coastal circulation model incorporated with the Mellor and Yamada level 2.5 turbulent closure model based on the reasonable model output of the M-2 tide and density residual currents. In the numerical experiments, upwelling motion appears around all the fronts with different velocity structures, accounting for surface cold water around the fronts. The experiments also suggest that the location and formation of fronts are closely related to topography and tidal mixing, as is the velocity structure around the front.

Descriptive Simulation: Combining Symbolic and Numerical Methods in the Analysis of Chemical Reaction Mechanisms

The Kineticist's Workbench is a computer program currently under development whose purpose is to help chemists understand, analyze, and simplify complex chemical reaction mechanisms. This paper discusses one module of the program that numerically simulates mechanisms and constructs qualitative descriptions of the simulation results. These descriptions are given in terms that are meaningful to the working chemist (e.g., steady states, stable oscillations, and so on); and the descriptions (as well as the data structures used to construct them) are accessible as input to other programs.

Coherent Behavior from Incoherent Knowledge Sources in the Automatic Synthesis of Numerical Computer Programs

A fundamental problem in artificial intelligence is obtaining coherent behavior in rule-based problem solving systems. A good quantitative measure of coherence is time behavior; a system that never, in retrospect, applied a rule needlessly is certainly coherent; a system suffering from combinatorial blowup is certainly behaving incoherently. This report describes a rule-based problem solving system for automatically writing and improving numerical computer programs from specifications. The specifications are in terms of "constraints" among inputs and outputs. The system has solved program synthesis problems involving systems of equations, determining that methods of successive approximation converge, transforming recursion to iteration, and manipulating power series (using differing organizations, control structures, and argument-passing techniques).