雾化进样石墨炉原子吸收光谱法的研究 Ⅰ.直接雾化进样装置及其性能

介绍一种自行设计和加工的直接雾化石墨炉进样装置的结构、工作原理和操作方法。该装置由微机控制与石墨炉加热程序同步工作,可以实现自动控制和自动测量。对雾化进样石墨炉分析的仪器参数和实验条件进行了研究和选择。结果表明,该装置与石墨炉结合具有自动化程度高、重现性好、样品利用率高和分析速度快等特点。     

Ice‐cover variability on shallow lakes at high latitudes: model simulations and observations

A one‐dimensional thermodynamic model for simulating lake‐ice phenology is presented and evaluated. The model can be driven with observed daily or hourly atmospheric forcing of air temperature, relative humidity, wind speed, cloud amount and snowfall. In addition to computing the energy balance components, key model output includes the temperature profile at an arbitrary number of levels within the ice/snow (or the water temperature if there is no ice) and ice thickness (clear ice and snow‐ice) on a daily basis, as well as freeze‐up and break‐up dates. The lake‐ice model is used to simulate ice‐growth processes on shallow lakes in arctic, sub‐arctic, and high‐boreal forest environments. Model output is compared with field and remote sensing observations gathered over several ice seasons. Simulated ice thickness, including snow‐ice formation, compares favourably with field measurements. Ice‐on and ice‐off dates are also well simulated when compared with field and satellite observations, with a mean absolute difference of 2 days. Model simulations and observations illustrate the key role that snow cover plays on the seasonal evolution of ice thickness and the timing of spring break‐up. It is also shown that lake morphometry, depth in particular, is a determinant of ice‐off dates for shallow lakes at high latitudes. Copyright © 2003 John Wiley & Sons, Ltd.     

Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

Magnitude-frequency relationship of coastal sand delivery by a southern California stream

Coastal sand delivery by a stream in southern California is estimated based on a numerical model which stimulates unsteady flow, sediment transport, and the associated channel adjustments for a stream-delta system. An average annual sediment yield of 51,400 m³/yr is estimated for the San Dieguito River, which drains a semiarid watershed controlled by dams. Of the total sand delivery by the stream, 20.5 percent is contributed from floods greater than the 100-year flood; 17.6 percent from those between the 50- and 100-year events; 28.4 percent from those between the 25- and 50 year floods; and 33.5 percent from those smaller than the 25-year flood.     

Freely floating-body simulation by a 2D fully nonlinear numerical wave tank

Nonlinear interactions between large waves and freely floating bodies are investigated by a 2D fully nonlinear numerical wave tank (NWT). The fully nonlinear 2D NWT is developed based on the potential theory, MEL/material-node time-marching approach, and boundary element method (BEM). A robust and stable 4th-order Runge–Kutta fully updated time-integration scheme is used with regriding (every time step) and smoothing (every five steps). A special φ_n-η type numerical beach on the free surface is developed to minimize wave reflection from end-wall and wave maker. The acceleration-potential formulation and direct mode-decomposition method are used for calculating the time derivative of velocity potential. The indirect mode-decomposition method is also independently developed for cross-checking. The present fully nonlinear simulations for a 2D freely floating barge are compared with the corresponding linear results, Nojiri and Murayama’s (Trans. West-Jpn. Soc. Nav. Archit. 51 (1975)) experimental results, and Tanizawa and Minami’s (Abstract for the 6th Symposium on Nonlinear and Free-surface Flow, 1998) fully nonlinear simulation results. It is shown that the fully nonlinear results converge to the corresponding linear results as incident wave heights decrease. A noticeable discrepancy between linear and fully nonlinear simulations is observed near the resonance area, where the second and third harmonic sway forces are even bigger than the first harmonic component causing highly nonlinear features in sway time series. The surprisingly large second harmonic heave forces in short waves are also successfully reproduced. The fully updated time-marching scheme is found to be much more robust than the frozen-coefficient method in fully nonlinear simulations with floating bodies. To compare the role of free-surface and body-surface nonlinearities, the body-nonlinear-only case with linearized free-surface condition was separately developed and simulated.     

A comparison of several wave spectra for the random wave diffraction by a semi-infinite breakwater

A comparison of the diffraction of multidirectional random waves using several selected wave spectrum models is presented in this paper. Six wave spectrum models, Bretschneider, Pierson–Moskowitz, ISSC, ITTC, Mitsuyasu, and JONSWAP spectrum, are considered. A discrete form for each of the given spectrum models is used to specify the incident wave conditions. Analytical solutions based on both the Fresnel integrals and polynomial approximations of the Fresnel integrals and numerical solutions of a boundary integral approach have been used to obtain the two-dimensional wave diffraction by a semi-infinite breakwater at uniform water depth. The diffraction of random waves is based on the cumulative superposition of linear diffraction solution. The results of predicted random wave diffraction for each of the given spectrum models are compared with those of the published physical model presented by Briggs et al. [1995. Wave diffraction around breakwater. Journal of Waterway, Port, Coastal and Ocean Engineering—ASCE 121(1), 23–35]. Reasonable agreement is obtained in all cases. The effect of the directional spreading function is also examined from the results of the random wave diffraction. Based on these comparisons, the present model for the analysis of various wave spectra is found to be an accurate and efficient tool for predicting the random wave field around a semi-infinite breakwater or inside a harbor of arbitrary geometry in practical applications.     

Sedimentology of a subtidal,tide-dominated sand body in the Yellow Sea,southwest Korea

Odanam Satoe, a subtidal, tide-dominated sand body in the Yellow Sea, Korea, is linear in plan and asymmetrical in cross-section. It consists of fine- to medium-grained, well-sorted subangular sand. Bedforms consist of high-amplitude (1–2 m) sandwaves on the lower flanks of the gentler-sloping bar surface, and medium-amplitude (0.5-1 m) sandwaves on the sand body trough adjoining the steeper face, the bar crest and shallower parts of the gently sloping bar surface. Bedforms are absent on the relatively steeper bar surface, which is characterized by 2° slopes. Bedform orientation on the gentler slope is oblique by 30° to the bar crest, parallel to the sand-body crest on the crest itself, and opposite to the steeper sand-body face in the trough below the steeper slope of the bar.Bottom current velocity data show that tidal currents are semi-rotary with a flood time—velocity asymmetry over the gentler slope, and ebb time—velocity asymmetry over the steeper slope during most of the tidal cycle. Tidal-current flow parallels bar elongation over the steeper slope, whereas over the gentler slope, tidal-current flow is directed at 30° to the bar crest and changes to normal to the crest one hour prior to low tide. Bedform orientation mapped with side-scan sonar shows agreement with these flow directions.Sand dispersal around the sand body is controlled by time—velocity asymmetry and partial rotary flow directions of tidal currents. This circulation causes not only a trapezoidal mode of grain dispersal, but also westerly migration of the sand body documented from comparative bathymetric surveys in 1964 and 1980.     

A general method for calculating three-dimensional laminar and turbulent boundary layers on ship hulls

A general method for representing the flow properties in the three-dimensional boundary layers around ship hulls of arbitrary shape is described. It makes use of an efficient two-point finite-diffirence schem to solve the boundary-layer equations and includes an algebraic eddy-viscosity representaion of the Reynolds-stress ternsor. The numericzal method contains novel and desirable features and allows the calculation of flows in which the circumferential velocity component contains regions of flow reversal across the boundary layer. The inviscid pressure distribution is determined with the Douglas-Neumann method which, if necessary, can conveniently allow for the boundary-layer displacement surface. To allow its application to ships, and particularly to those with double-elliptic and flat-bottomed hulls, a non-orthogonal coordinate system has been developed and is shown to be economical, precise and comparatively easy to use. Present calculations relate to zero Froude number but they can be extended to include the effects of a water wave and local regions of flow separation which may stem from bulbous-bow geometries.     

Fully nonlinear numerical wave tank (NWT) simulations and wave run-up prediction around 3-D structures

A finite-difference scheme and a modified marker-and-cell (MAC) algorithm have been developed to investigate the interactions of fully nonlinear waves with two- or three-dimensional structures of arbitrary shape. The Navier–Stokes (NS) and continuity equations are solved in the computational domain and the boundary values are updated at each time step by the finite-difference time-marching scheme in the framework of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique developed for two fluid layers.To demonstrate the capability and accuracy of the present method, the numerical simulation of backstep flows with free-surface, and the numerical tests of the MDF technique with limit functions are conducted. The 3D program was then applied to nonlinear wave interactions with conical gravity platforms of circular and octagonal cross-sections. The numerical prediction of maximum wave run-up on arctic structures is compared with the prediction of the Shore Protection Manual (SPM) method and those of linear and second-order diffraction analyses based on potential theory and boundary element method (BEM). Through this comparison, the effects of non-linearity and viscosity on wave loading and run-up are discussed.