Improvement of the Semi-Lagrangian Advection Scheme in the GRAPES Model:Theoretical Analysis and Idealized Tests

ABSTRACT The Global/Regional Assimilation and PrEdiction System （GRAPES） is the newgeneration numerical weather predic- tion （NWP） system developed by the China Meteorological Administration. It is a fully compressible non-hydrostatical global/regional unified model that uses a traditional semi-Lagrangian advection scheme with cubic Lagrangian interpola tion （referred to as the SL_CL scheme）. The SL_CL scheme has been used in many operational NWP models, but there are still some deficiencies, such as the damping effects due to the interpolation and the relatively low accuracy. Based on Reich＇s semi-Lagrangian advection scheme （referred to as the R2007 scheme）, the Re_R2007 scheme that uses the low- and high-order B-spline function for interpolation at the departure point, is developed in this paper. One- and two-dimensional idealized tests in the rectangular coordinate system with uniform grid cells were conducted to compare the Re..R2007 scheme and the SL_CL scheme. The numerical results showed that： （1） the damping effects were remarkably reduced with the Re_R2007 scheme; and （2） the normalized errors of the Re_R2007 scheme were about 7.5 and 3 times smaller than those of the SL_CL scheme in one- and two-dimensional tests, respectively, indicating the higher accuracy of the Re..R2007 scheme. Furthermore, two solid-body rotation tests were conducted in the latitude-longitude spherical coordinate system with non uniform grid cells, which also verified the Re_R2007 scheme＇s advantages. Finally, in comparison with other global advection schemes, the Re_R2007 scheme was competitive in terms of accuracy and flow independence. An encouraging possibility for the application of the Re_R2007 scheme to the GRAPES model is provided.

A three-dimensional semi-implicit unstructured grid finite volume ocean model

A new three-dimensional semi-implicit finite-volume ocean model has been developed for simulating the coastal ocean circulation, which is based on the staggered C-unstructured non-orthogonal grid in the hor- izontal direction and z-level grid in the vertical direction. The three-dimensional model is discretized by the semi-implicit finite-volume method, in that the free-surface and the vertical diffusion are semi-implicit, thereby removing stability limitations associated with the surface gravity wave and vertical diffusion terms. The remaining terms in the momentum equations are discretized explicitly by an integral method. The partial cell method is used for resolving topography, which enables the model to better represent irregular topography. The model has been tested against analytical cases for wind and tidal oscillation circulation, and is applied to simulating the tidal flow in the Bohal Sea. The results are in good agreement both with the analytical solutions and measurement results.

GA-Based Model Predictive Control of Semi-Active Landing Gear

Semi-active landing gear can provide good performance of both landing impact and taxi situation, and has the ability for adapting to various ground conditions and operational conditions. A kind of Nonlinear Model Predictive Control algorithm （NMPC） for semi-active landing gears is developed in this paper. The NMPC algorithm uses Genetic Algorithm （GA） as the optimization technique and chooses damping performance of landing gear at touch down to be the optimization object. The valve＇s rate and magnitude limitations are also considered in the controller＇s design. A simulation model is built for the semi-active landing gear＇s damping process at touchdown. Drop tests are carried out on an experimental passive landing gear systerm to validate the parameters of the simulation model. The result of numerical simulation shows that the isolation of impact load at touchdown can be significantly improved compared to other control algorithms. The strongly nonlinear dynamics of semi-active landing gear coupled with control valve＇s rate and magnitude limitations are handled well with the proposed controller.

The Application of Flux-Form Semi-Lagrangian Transport Scheme in a Spectral Atmosphere Model

A flux-form semi-Lagrangian transport scheme （FFSL） was implemented in a spectral atmospheric GCM developed and used at IAP/LASG. Idealized numerical experiments show that the scheme is good at shape preserving with less dissipation and dispersion, in comparison with other conventional schemes, hnportantly, FFSL can automatically maintain the positive definition of the transported tracers, which was an underlying problem in the previous spectral composite method （SCM）. To comprehensively investigate the impact of FFSL on GCM results, we conducted sensitive experiments. Three main improvements resulted： first, rainfall simulation in both distribution and intensity was notably improved, which led to an improvement in precipitation frequency. Second, the dry bias in the lower troposphere was significantly reduced compared with SCM simulations. Third, according to the Taylor diagram, the FFSL scheme yields simulations that are superior to those using the SCM： a higher correlation between model output and observation data was achieved with the FFSL scheme, especially for humidity in lower troposphere. However, the moist bias in the middle and upper troposphere was more pronounced with the FFSL scheme. This bias led to an over-simulation of precipitable water in comparison with reanalysis data. Possible explanations, as well as solutions, are discussed herein.

RUIN PROBABILITY IN A SEMI-MARKOV RISK MODEL WITH CONSTANT INTEREST FORCE AND HEAVY-TAILED CLAIMS

In the present paper, we consider a kind of semi-Markov risk model （SMRM） with constant interest force and heavy-tailed claims~ in which the claim rates and sizes are conditionally independent, both fluctuating according to the state of the risk business. First, we derive a matrix integro-differential equation satisfied by the survival probabilities. Second, we analyze the asymptotic behaviors of ruin probabilities in a two-state SMRM with special claim amounts. It is shown that the asymptotic behaviors of ruin probabilities depend only on the state 2 with heavy-tailed claim amounts, not on the state 1 with exponential claim sizes.

Fuzzy Control of Model Travel Tracking for Vehicle Semi-Active Suspension

The control strategy of the model travel tracking for the vehicle suspension sys tem is presented based on analyzing the responses of the vehicle suspension tra vel. A fuzzy control system of vehicle suspension is designed, in which the sus pension travel output of the adaptive LQG control system is taken as the tracking objective. The simulation results prove that the suspension travel and vertical acceleration can be tracked simultaneously with the simple fuzzy controller, and the tracking effect of fuzzy control is better than that of the PID controller.

A Problem of a Semi-Infinite Medium Subjected to Exponential Heating Using a Dual-Phase-Lag Thermoelastic Model

The problem of a semi-infinite medium subjected to thermal shock on its plane boundary is solved in the context of the dual-phase-lag thermoelastic model. The expressions for temperature, displacement and stress are presented. The governing equations are expressed in Laplace transform domain and solved in that domain. The solution of the problem in the physical domain is obtained by using a numerical method for the inversion of the Laplace transforms based on Fourier series expansions. The numerical estimates of the displacement, temperature, stress and strain are obtained for a hypothetical material. The results obtained are presented graphically to show the effect phase-lag of the heat flux and a phase-lag of temperature gradient on displacement, temperature, stress.

Static analysis of ultra-thin beams based on a semi-continuum model

A linear semi-continuum model with discrete atomic layers in the thickness direction was developed to investigate the bending behaviors of ultra-thin beams with nanoscale thickness.The theoretical results show that the deflection of an ultra-thin beam may be enhanced or reduced due to different relaxation coefficients.If the relaxation coefficient is greater/less than one,the deflection of micro/nano-scale structures is enhanced/reduced in comparison with macro-scale structures.So,two opposite types of size-dependent behaviors are observed and they are mainly caused by the relaxation coefficients.Comparisons with the classical continuum model,exact nonlocal stress model and finite element model （FEM） verify the validity of the present semi-continuum model.In particular,an explanation is proposed in the debate whether the bending stiffness of a micro/nano-scale beam should be greater or weaker as compared with the macro-scale structures.The characteristics of bending stiffness are proved to be associated with the relaxation coefficients.

Towards a semi-empirical model of the sea ice thickness based on hyperspectral remote sensing in the Bohai Sea

Sea ice thickness is one of the most important input parameters for the prevention and mitigation of sea ice disasters and the prediction of local sea environments and climates. Estimating the sea ice thickness is currently the most important issue in the study of sea ice remote sensing. With the Bohai Sea as the study area, a semiempirical model of the sea ice thickness（SEMSIT） that can be used to estimate the thickness of first-year ice based on existing water depth estimation models and hyperspectral remote sensing data according to an optical radiative transfer process in sea ice is proposed. In the model, the absorption and scattering properties of sea ice in different bands（spectral dimension information） are utilized. An integrated attenuation coefficient at the pixel level is estimated using the height of the reflectance peak at 1 088 nm. In addition, the surface reflectance of sea ice at the pixel level is estimated using the 1 550–1 750 nm band reflectance. The model is used to estimate the sea ice thickness with Hyperion images. The first validation results suggest that the proposed model and parameterization scheme can effectively reduce the estimation error associated with the sea ice thickness that is caused by temporal and spatial heterogeneities in the integrated attenuation coefficient and sea ice surface. A practical semi-empirical model and parameterization scheme that may be feasible for the sea ice thickness estimation using hyperspectral remote sensing data are potentially provided.

Analytical models for the penetration of semi-infinite targets by rigid,deformable and erosive long rods

A theoretical study is presented herein on the pen- etration of a semi-infinite target by a spherical-headed long rod for Yp 〉 S, where Yp is the penetrator strength and S is the static target resistance. For Yp 〉 S, depending upon initial impact velocity, there exist three types of penetration, namely, penetration by a rigid long rod, penetration by a deforming non-erosive long rod and penetration by an erosive long rod. If the impact velocity of the penetrator is higher than the hydrodynamic velocity （VH）, it will penetrate the target in an erosive mode; if the impact velocity lies between the hydrodynamic velocity （VH） and the rigid body velocity （VR）, it will penetrate the target in a deformable mode; if the impact velocity is less than the rigid body velocity （VR）, it will penetrate the target in a rigid mode. The critical conditions for the transition among these three penetration modes are proposed. It is demonstrated that the present model predictions correlate well with the experimental observations in terms of depth of penetration （DOP） and the critical transition conditions.

Dependence Model Selection for Semi-Competing Risks Data

We consider the model selection problem of the dependency between the?terminal event and the non-terminal event under semi-competing risks data. When the relationship between the two events is unspecified, the inference on the non-terminal event is not identifiable. We cannot make inference on the non-terminal event without extra assumptions. Thus, an association model for?semi-competing risks data is necessary, and it is important to select an appropriate dependence model for a data set. We construct the likelihood function for semi-competing risks data to select an appropriate dependence model. From?simulation studies, it shows the performance of the proposed approach is well. Finally, we apply our method to a bone marrow transplant data set.

Weather Yield Model for the Semi Tropical Region (Pakistan)

Weather models are essential tools for checking of the effect of the weather elements in terms of their effect on the production of the crop. This research is an attempt to see the effect of only two variables i.e., temperature and rainfall for the division Faisalabad (semitropical region of Pakistan).The model fitted is of the linear form:the values of a,b, c have been found. The expected yield has been calculated by using the aridity indices (X1 and X2 ) and the result in the form of coefficient of determination R2 has been found equal to 0.166. The significance of the regression coefficient has been tested, which shows that the contribution to the yield from aridity index at germination and that at ripening is significant.The wheat yields are the results of a wide variety of variables, most of which show varying degree of relationship with one another, some positive and some negative in terms of output. These variables may be technology, fertilizers, pesticides, epidemics, kinds of seeds used, market price of crop and the area under cultivation etc, which can be the source of variation in the wheat yield. Since rainfall during germination and temperature at the ripening periods are the necessary factors for the yield of wheat, for this purpose these parameters have been studied in order to their contribution.

Comparative Study among Different Semi-Empirical Models for Soil Salinity Prediction in an Arid Environment Using OLI Landsat-8 Data

Salt-affected soils, caused by natural or human activities, are a common environmental hazard in semi-arid and arid landscapes. Excess salts in soils affect plant growth and production, soil and water quality and, therefore, increase soil erosion and land degradation. This research investigates the performance of five different semi-empirical predictive models for soil salinity spatial distribution mapping in arid environment using OLI sensor image data. This is the first attempt to test remote sensing based semi-empirical salinity predictive models in this area: the Kingdom of Bahrain. To achieve our objectives, OLI data were standardized from the atmosphere interferences, the sensor radiometric drift, and the topographic and geometric distortions. Then, the five semi-empirical predictive models based on the Normalized Difference Salinity Index (NDSI), the Salinity Index-ASTER (SI-ASTER), the Salinity Index-1 (SI-1), the Soil Salinity and Sodicity Index-1 and Index-2 (SSSI-1 and SSSI-2), developed for slight and moderate salinity in agricultural land, were implemented and applied to OLI image data. For validation purposes, a fieldwork was organized and different important spots-locations representing different salinity levels were visited, photographed, and localized using an accurate GPS (σ ≤ ±30 cm). Based on this a priori knowledge of the soil salinity, six validation sites were selected to reflect non-saline, low, moderate, high and extreme salinity classes, descriptive statistics extracted from polygons and/or transects over these sites were used. The obtained results showed that the models based on NDSI, SI-1 and SI-ASTER all failed to detect salinity bounds for both extreme salinity (Sabkhah) and non-saline conditions. In Fact, NDSI and SI-ASTER gave respectively only 35% dS/m and 25% dS/m in extreme salinity validation site, while SI-1 and SI-ASTER indicated 38% dS/m and 39% dS/m in non-saline validation site. Therefore, these three models were deemed inadequate for the study site. However, both SSSI-1 and SSSI-2 allowed a detection of the previous salinity bounds and furthermore described similarly and correctly the urban-vegetation areas and the open-land areas. Their predicted EC is around 10% dS/m for non-saline urban soil, about 25% dS/m for low salinity urban-vegetation soil, approximately 30% to 75% dS/m, respectively, for moderate to high salinity soils. SSSI-2 based semi-empirical salinity models was able to differentiate the high salinity versus extreme salinity in areas where both exist and was very accurate to highlight the pure salt where SSSI-1 has reach saturation for both salinity classes. In conclusion, reliable salinity map was produced using the model based on SSSI-2 and OLI sensor data that allows a better characterization of the soil salinity problem in an Arid Environment.

Semi-Automatic Modeling Technique of Torque Converter Flow Passage

The modeling technique of hydrodynamic torque converter flow passage was investigated. The semi-automatic modeling technique of torque converter flow passage was proposed. The flow passage model of each converter wheel is considered as a revolution entity sliced by two curved surfaces. In order to generate the revolution entity, a new approximation method, condition optimum arc approximation, was proposed. The method was used to approximate the meridional streamlines of the inner and outer wall. As a result, the three-dimensional revolution entity can be conveniently generated. In order to create slice surfaces, the central stream surface of flow passage was approximated with a quadric surface. The normal vector of the quadric surface and the thickness/thickness-function of bade were used to calculate the discrete point coordinates of blade surfaces. Via the rotation transformation to the coordinates, the discrete point coordinates of slice surfaces were obtained. A parameterized program code used for the hydrodynamic torque converter design and semi-automatic modeling was developed. Modeling errors were calculated and analyzed. The flow passage model was generated in several minutes with the help of the program code, Auto CAD and Solidworks software. Finally, the model was inputted into Gambit, and the pre-processing task used for the numerical simulation of torque converter flow field was successfully completed. The investigation results show that the semi-automatic modeling not only can ensure the accuracy of modeling, but also librates the research and design workers of torque converter from the time-consuming modeling work, which paves the way for the numerical simulation of the complex flow field of the hydrodynamic torque converter.

Mechanistic Drifting Forecast Model for A Small Semi-Submersible Drifter Under Tide–Wind–Wave Conditions

Understanding the drifting motion of a small semi-submersible drifter is of vital importance regarding monitoring surface currents and the floating pollutants in coastal regions. This work addresses this issue by establishing a mechanistic drifting forecast model based on kinetic analysis. Taking tide–wind–wave into consideration, the forecast model is validated against in situ drifting experiment in the Radial Sand Ridges. Model results show good performance with respect to the measured drifting features, characterized by migrating back and forth twice a day with daily downwind displacements. Trajectory models are used to evaluate the influence of the individual hydrodynamic forcing. The tidal current is the fundamental dynamic condition in the Radial Sand Ridges and has the greatest impact on the drifting distance. However, it loses its leading position in the field of the daily displacement of the used drifter. The simulations reveal that different hydrodynamic forces dominate the daily displacement of the used drifter at different wind scales. The wave-induced mass transport has the greatest influence on the daily displacement at Beaufort wind scale 5–6; while wind drag contributes mostly at wind scale 2–4.

Construction of Crack Perturbation Model and Forward Semi-analytical Model of Attached Eddy Current Sensor

Fatigue damage monitoring is critical metallic structure health monitoring of aircraft.The sensor should be high sensitive,easy to be integrated into structure and well adaptable for poor working conditions.Therefore,an attached eddy current sensor with flexible plane is put forward and its characteristics are analyzed.By extracting material′s conductivity as the crack features,forward semi-analytical model is established and parameter optimizations are carried out.Crack perturbation model of attached eddy current sensor is constructed,and perturbation voltages of sensing channels under three-dimension structural crack are obtained.To verify the sensor′s performance,monitoring experiment on crack extension is conducted under condition of 3 MHz frequency.The validation experimental results show that perturbation model of 2A12-T4 aluminum alloy agrees well with experiment results,and perturbation model errors of four sensing channels are within 25%.The attached eddy current sensor is capable of testing the crack nondestructively and measuring the crack extension quantitatively with the accuracy of 1mm.

A Study of Aviation Weapon Equipment Maintenance Based on the Semi-Markov Model

Based on the Semi-Markov mathematical description, the multiple states of maintenance processes for aviation weapon equipment are studied. Six kinds of maintenance states are determined and the Semi-Markov model of the maintenance process is given. According to maintenance characteristic, the multiple states maintenance processes are divided into the wait, use and alternate stages. Through using the mathematical model for the different stages, the probability in different states and effective index on different stages are obtained. These results are available to the maintenance practice.

Application of spatially varying storage capacity model for runoff parameterization in semi-arid catchment

This paper introduces the method of designation of water storage capacity for each grid cell within a catchment, which considers topography, vegetation and soil synthetically. For the purpose of hydrological process simulation in semi-arid regions, a spatially varying storage capacity （VSC） model was developed based on the spatial distribution of water storage capacity and the vertical hybrid runoff mechanism. To verify the applicability of the VSC model, both the VSC model and a hybrid runoff model were used to simulate daily runoff processes in the catchment upstream of the Dianzi hydrological station from 1973 to 1979. The results showed that the annual average Nash-Sutcliffe coefficient was 0.80 for the VSC model, and only 0.67 for the hybrid runoff model. The higher annual average Nash-Sutcliffe coefficient of the VSC model means that this hydrological model can better simulate daily runoff processes in semi-arid regions. Furthermore, as a distributed hydrological model, the VSC model can be applied in regional water resource management.

Nondispersive Periodic Solution of a Barotropic Semi-Geostrophic Model

The existence and solution of the non-dispersive periodic solution are achieved concerning nonlinear barotropic Rossby waves of a barotropic semi-geostrophic model, demonstrating the likelihood of the Taylor evolution, together with the related dimensionless a-criterion. Finally, the wave velocity expression is proposed with some diagnostic relations among the wave parameters.

Simulation of soil erosion dynamics in grassland ecosystem of semi-arid China using the CENTURY model

Soil erosion can cause considerable effect on global natural resources and eco-environment. In the paper, the CENTURY model has been used to simulate soil erosion in Xilin Gol Grassland of Inner Mongolia. The results showed before the 1960s, the soil erosion amount was over 2 kg /m2.a in grassland ecosystem in the study area because no trees had been planted. But after the 1960s the mean annual accumulator C lost from soil organic matter due to soil erosion was only 0.3 kg /m2.a in forest ecosystem. So afforestation has exerted notable effect on decreasing soil erosion amount in Xilin Gol Grassland.