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.     

An approach to aid INS using time-differenced GPS carrier phase (TDCP) measurements

This paper presents a simple and effective approach that incorporates single-frequency, L1 time-differenced GPS carrier phase (TDCP) measurements without the need of ambiguity resolution techniques and the complexity to accommodate the delayed-state terms. Static trial results are included to illustrate the stochastic characteristics and effectiveness of the TDCP measurements in controlling position error growth. The formulation of the TDCP observation model is also described in a 17-state tightly-coupled GPS/INS iterative, extended Kalman filter (IEKF) approach. Preliminary land vehicle trial results are also presented to illustrate the effectiveness of the TDCP which provides sub-meter positional accuracies when operating for more than 10 min.     

朝鲜半岛古元古代摩天岭群的碎屑锆石U-Pb年龄及其地质意义

克拉通内部活动带变质沉积建造蕴藏了早期超大陆形成中地体聚-散动力学过程的重要信息,碎屑锆石年代学则是揭示这种信息的最佳工具之一。针对华北(中-朝)克拉通最具代表性的古元古代活动带(辽-吉-朝活动带)中巨量沉积建造精细表征方面的空白区域,本文开展了朝鲜半岛检德地区摩天岭群的碎屑锆石年代学研究。SIMS U-Pb定年结果显示,北大川组变火山碎屑岩样品的锆石年龄谱记录了约2.53Ga的主要年龄峰值和约2.19Ga的最年轻峰值,说明北大川组开始沉积的时代应该晚于2.2Ga。南大川组变沉积岩样品碎屑锆石记录的最年轻峰值年龄为约2.02Ga,说明南大川组的沉积时代应该不早于2.02Ga;其锆石年龄谱记录了一个主要年龄峰值区域(2.02~2.18Ga)、一个次要年龄峰值(约2.46Ga)和一些2.78~3.30 Ga的古老年龄,指示其物质源区主要为同区古元古代地质体和邻近地体上的太古宙变质基底。区域对比表明,摩天岭群与中国境内辽吉活动带中的辽河群及集安群/老岭群具有基本一致的碎屑锆石年龄谱,说明它们不仅是同一时代的产物,而且分享类似的物质源区。这些时代一致、物源相近的沉积建造可以为甄别古元古代辽-吉-朝活动带的大地构造属性提供重要参考。与沉积年龄相近的碎屑锆石占据主导的碎屑锆石年龄谱反映了汇聚板块边缘各类盆地的典型特征,而相悖于大部分碎屑锆石年龄远大于地层沉积年龄的裂谷型盆地。     

朝鲜半岛平南盆地中元古代岩浆事件

朝鲜平南盆地翁津地区发育中元古代黄海群和同时期(称之为瓮津期)花岗岩,花岗岩体侵入于黄海群。本文采用锆石原位微区U-Pb定年技术,对黄海群中的酸性火山岩及花岗岩进行了年龄测试。获得的数据表明,黄海群中下部层位及上部层位的酸性火山岩分别在1235±5Ma和1203±7Ma喷发,由此说明黄海群的沉积时代应为中元古代,而不是传统上认为的古元古代;两个翁津期花岗岩体(翁津和黄衣山岩体)的侵位年龄分别为1251±22Ma和1248±13Ma,为中元古代花岗质岩浆活动的产物。上述1251~1203Ma年龄的获得,表明朝鲜半岛发育中元古代岩浆作用,从而明确朝鲜黄海裂谷与华北东缘裂谷在时间上具有同期性,同时也表明中国华北与朝鲜在中元古代具有类似的发展历史。     

Thermal fronts and cross-frontal heat flux in the southern Huanghai Sea and the East China Sea

Synoptic features in/around thermal fronts and cross-frontal heat fluxes in the southern Huanghai./Yellow Sea and East China Sea （HES） were examined using the data collected from four airborne expendable bathythermograph surveys with horizontal approxmately 35 km and vertical 1 m（from the surface to 400 m deep） spacings. Since the fronts are strongly affected by HES current system, the synoptic thermal features in/around them represent the interaction of currents with surrounding water masses. These features can not be obtained from climatological data. The identified thermal features are listed as follows ： （ 1 ） multiple boundaries of cold water, asymmetric thermocline intrusion, locally-split front by homogeneous water of approxmately 18 ℃, and mergence of the front by the Taiwan Warm Current in/around summertime southern Cheju - Changjiang/Yangtze front and Tsushima front; （2） springtime frontal eddy-like feature around Tsushima front; （3） year-round cyclonic meandering and summertime temperature-inversion at the bottom of the surface mixed layer in Cheju - Tsushima front; and （4） multistructure of Kuroshio front. In the Kuroshio front the mean variance of vertical temperature gradient is an order of degree smaller than that in other HES fronts. The southern Cheju- Changjiang front and Cheju -Tsushima front are connected with each other in the summer with comparable cross-frontal temperature gradient. However, cross-frontal heat flux and lateral eddy diffusivity are stronger in the southern Cheju - Changjiang front. The cross-frontal heat exchange is the largest in the mixing zone between the modified Huanghai Sea bottom cold water and the Tsushima Warm Current, which is attributable to enhanced thermocline intrusions.     

Simulation of the extreme waves generated by typhoon Bolaven (1215) in the East China Sea and Yellow Sea

Record-breaking high waves occurred during the passage of the typhoon Bolaven (1215) (TYB) in the East China Sea (ECS) and Yellow Sea (YS) although its intensity did not reach the level of a super typhoon. Winds and directional wave measurements were made using a range of in-situ instruments mounted on an ocean tower and buoys. In order to understand how such high waves with long duration occurred, analyses have been made through measurement and numerical simulations. TYB winds were generated using the TC96 typhoon wind model with the best track data calibrated with the measurements. And then the wind fields were blended with the reanalyzed synoptic-scale wind fields for a wave model. Wave fields were simulated using WAM4.5 with adjustment of C_d for gust of winds and bottom friction for the study area. Thus the accuracy of simulations is considerably enhanced, and the computed results are also in better agreement with measured data than before. It is found that the extremely high waves evolved as a result of the superposition of distant large swells and high wind seas generated by strong winds from the front/right quadrant of the typhoon track. As the typhoon moved at a speed a little slower than the dominant wave group velocity in a consistent direction for two days, the wave growth was significantly enhanced by strong wind input in an extended fetch and non-linear interaction.     

Differential rotation of some HK-Project stars and the butterfly diagrams

We analyze the long-term variability of the chromospheric radiation of 20 stars monitored in the course of the HK-Project at the Mount Wilson Observatory. We apply the modified wavelet algorithm for this set of gapped time series. Besides the mean rotational periods for all these stars, we find reliable changes of the rotational periods from year to year for a few stars. Epochs of slower rotation occur when the activity level of the star is high, and the relationship repeats again during the next maximum of an activity cycle. Such an effect is traced in two stars with activity cycles that are not perfectly regular (but labeled “Good” under the classification in [Baliunas, S.L., Donahue, R.A., Soon, W.H., Horne, J.H., Frazer, J., Woodard-Eklund, L., Bradford, M., Rao, L.M., Wilson, O.C., Zhang, Q. et al., 1995. ApJ 438, 269.]) but the two stars have mean activity levels exceed that of the Sun. The averaged rotational period of HD 115404 is 18.5 days but sometimes the period increases up to 21.5 days. The sign of the differential rotation is the same as the Sun’s, and the value ΔΩ/Ω=-0.14. For the star HD 149661, this ratio is −0.074. Characteristic changes of rotational periods occur over around three years when the amplitude of the rotational modulation is large. These changes can be transformed into latitude-time butterfly diagrams with minimal a priori assumptions. We compare these results with those for the Sun as a star and conclude that epochs when surface inhomogeneities rotate slower are synchronous with the reversal of the global magnetic dipole.     

Down-scaled regional ocean modeling system (ROMS) for high-resolution coastal hydrodynamics in Korea

A down-scaled operational oceanographic system is developed for the coastal waters of Korea using a regional ocean modeling system(ROMS).The operational oceanographic modeling system consists of atmospheric and hydrodynamic models.The hydrodynamic model,ROMS,is coupled with wave,sediment transport,and water quality modules.The system forecasts the predicted results twice a day on a 72 h basis,including sea surface elevation,currents,temperature,salinity,storm surge height,and wave information for the coastal waters of Korea.The predicted results are exported to the web-GIS-based coastal information system for real-time dissemination to the public and validation with real-time monitoring data using visualization technologies.The ROMS is two-way coupled with a simulating waves nearshore model,SWAN,for the hydrodynamics and waves,nested with the meteorological model,WRF,for the atmospheric surface forcing,and externally nested with the eutrophication model,CE-QUAL-ICM,for the water quality.The operational model,ROMS,was calibrated with the tidal surface observed with a tide-gage and verified with current data observed by bottom-mounted ADCP or AWAC near the coastal waters of Korea.To validate the predicted results,we used real-time monitoring data derived from remote buoy system,HF-radar,and geostationary ocean color imager(GOCI).This down-scaled operational coastal forecasting system will be used as a part of the Korea operational oceanographic system(KOOS) with other operational oceanographic systems.