Atmosphere-Ocean Coupling Schemes in a One-Dimensional Climate Model

In this paper, the coupling schemes of atmosphere-ocean climate models are discussed with one-dimensional advection equations. The convergence and stability for synchronous and asynchronous schemes are demonstrated and compared.Conclusions inferred from the analysis are given below. The synchronous scheme as well as the asynchronous-implicit scheme in this model are stable for arbitrary integrating time intervals. The asynchronous explicit scheme is unstable under certain conditions, which depend upon advection velocities and heat exchange parameters in the atmosphere and oceans. With both synchronous and asynchronous stable schemes the discrete solutions converge to their unique exact ones. Advections in the atmosphere and ocean accelerate the rate of convergence of the asynchronous-implicit scheme. It is suggusted that the asynchronous-implicit coupling scheme is a stable and efficient method for most climatic simulations.

OPERATIONAL FORECAST OF RAINFALL INDUCED BY LANDFALLING TROPICAL CYCLONES ALONG GUANGDONG COAST

Following previous studies of the rainfall forecast in Shenzhen owing to landfalling tropical cyclones(TCs),a nonparametric statistical scheme based on the classification of the landfalling TCs is applied to analyze and forecast the rainfall induced by landfalling TCs in the coastal area of Guangdong province,China.All the TCs landfalling with the distance less than 700 kilometers to the 8 coastal stations in Guangdong province during 1950—2013 are categorized according to their landfalling position and intensity.The daily rainfall records of all the 8 meteorological stations are obtained and analyzed.The maximum daily rainfall and the maximum 3 days’accumulated rainfall at the 8 coastal stations induced by each category of TCs during the TC landfall period(a couple of days before and after TC landfalling time)from 1950 to 2013 are computed by the percentile estimation and illustrated by boxplots.These boxplots can be used to estimate the rainfall induced by landfalling TC of the same category in the future.The statistical boxplot scheme is further coupled with the model outputs from the European Centre for Medium-Range Weather Forecasts(ECMWF)to predict the rainfall induced by landfalling TCs along the coastal area.The TCs landfalling in south China from 2014 to 2017 and the corresponding rainfall at the 8 stations area are used to evaluate the performance of these boxplots and coupled boxplots schemes.Results show that the statistical boxplots scheme and coupled boxplots scheme can perform better than ECMWF model in the operational rainfall forecast along the coastal area in south China.

Dynamic Damage Detection and Assessment on the Bogie Frame of High-Speed Train Based on the Coupled Elastic and Multibody Dynamics Scheme

The dynamic cumulative damage of rigid-flexible coupling model of high-speed train with flexible bogie frame is performed by using the coupled scheme of elastic and multibody dynamics theories.The motion equations of the present problem are firstly established by integrating the finite element method and floating frame of reference approach based on the virtual power principle and D'Alembert principle.The process of condensing the elastic DOFs of the obtained finite element model involving the incorporation of the substructure technique and sparse approximate inverse method is tentatively carried out.Then,the motion equations are further solved by virtue of the generalized α method and the Jacobian-free Newton-Krylov technologies.And the superiority of coupled scheme is proven by comparing with the traditional approach.Finally,besides the dynamic behaviors of the considered vehicle model,the time-variations of stresses on the elastic bogie frame's dangerous nodes and the distributions of stresses of bogie frame at some specified moments are synchronously calculated and analyzed.More importantly,the real-time and time-varying cumulative damages of some typical nodes on bogie frame are investigated.

ANALYSIS ON A NUMERICAL SCHEME WITH SECOND-ORDER TIME ACCURACY FOR NONLINEAR DIFFUSION EQUATIONS

A nonlinear fully implicit finite difference scheme with second-order time evolution for nonlinear diffusion problem is studied.The scheme is constructed with two-layer coupled discretization(TLCD)at each time step.It does not stir numerical oscillation,while permits large time step length,and produces more accurate numerical solutions than the other two well-known second-order time evolution nonlinear schemes,the Crank-Nicolson(CN)scheme and the backward difference formula second-order(BDF2)scheme.By developing a new reasoning technique,we overcome the difficulties caused by the coupled nonlinear discrete diffusion operators at different time layers,and prove rigorously the TLCD scheme is uniquely solvable,unconditionally stable,and has second-order convergence in both s-pace and time.Numerical tests verify the theoretical results,and illustrate its superiority over the CN and BDF2 schemes.

A Coupled Discrete Unified Gas-Kinetic Scheme for Convection Heat Transfer in Porous Media

In this paper,the discrete unified gas-kinetic scheme(DUGKS)is extended to the convection heat transfer in porous media at representative elementary volume(REV)scale,where the changes of velocity and temperature fields are described by two kinetic equations.The effects from the porous medium are incorporated into the method by including the porosity into the equilibrium distribution function,and adding a resistance force in the kinetic equation for the velocity field.The proposed method is systematically validated by several canonical cases,including the mixed convection in porous channel,the natural convection in porous cavity,and the natural convection in a cavity partially filled with porous media.The numerical results are in good agreement with the benchmark solutions and the available experimental data.It is also shown that the coupled DUGKS yields a second-order accuracy in both temporal and spatial spaces.

Sharp interface direct forcing immersed boundary methods： A summary of some algorithms and applications

Body-fitted mesh generation has long been the bottleneck of simulating fluid flows involving complex geometries. Immersed boundary methods are non-boundary-conforming methods that have gained great popularity in the last two decades for their simplicity and flexibility, as well as their non-compromised accuracy. This paper presents a summary of some numerical algori- thms along the line of sharp interface direct forcing approaches and their applications in some practical problems. The algorithms include basic Navier-Stokes solvers, immersed boundary setup procedures, treatments of stationary and moving immersed bounda- ries, and fluid-structure coupling schemes. Applications of these algorithms in particulate flows, flow-induced vibrations, biofluid dynamics, and free-surface hydrodynamics are demonstrated. Some concluding remarks are made, including several future research directions that can further expand the application regime of immersed boundary methods.

Incorporation of TOPMODEL into land surface model SSiB and numerically testing the effects of the corporation at basin scale

In order to examine and analyze the effects of integration of land surface models with TOPMODEL and different approaches for the integration on the model simulation of water and energy balances,the coupled models have been developed,which incorporate TOPMODEL into the Simplified Biosphere Model(SSiB) with different approaches(one divides a basin into a number of zones according to the distribution of topographic index,and the other only divides the basin into saturated and unsaturated zones).The coupled models are able to(but SSiB is not able to) take into account the impacts of topography variation and vertical heterogeneity of soil saturated hydraulic conductivity on horizontal distribution of soil moisture and in turn on water and energy balances within the basin(or a grid cell).By using the coupled models and SSiB model itself,the daily hydrological components such as runoffs are simulated and final results are analyzed carefully.Simulated daily results of hydrological components produced by both SSiB and coupled models show that(i) There is significant difference between results of soil wetness,its vertical distribution and seasonal variation,water and energy balance,and daily runoff in the basin predicted by SSiB and by the coupled models.The land surface model currently used such as SSiB is likely to misrepresent real feature of water and energy balances in the basin.(ii) Compared with the results for the basin predicted by SSiB,the coupled models predict more strong vertical and seasonal changes in soil wetness,higher evaporation and lower runoff,and improve the base flow simulation obviously.(iii) Comparing the results for the basin predicted by two coupled models with different integration approach and SSiB one by one,the results of daily runoffs and soil wetness predicted by the two coupled models are quite similar.It means,for the coupled models,the approach by dividing a region being considered into more subzones may have limited effects on improving runoff simulation results.The scheme only to divide the region into saturated and unsaturated zones may be a convenient and effective scheme.But then,if the results from the two coupled models for the basin are carefully compared,the simulated results by the coupled model with dividing the basin into more subzones will show higher evaporation and surface runoff but lower subsurface flow,lower total runoff,and lower ground water level averaged for five years.

CFD-DEM simulation of the hole cleaning process in a deviated well drilling： The effects of particle shape

We investigate the effect of particle shape on the transportation mechanism in well-drilling using a three-dimensional model that couples computational fluid dynamics （CFD） with the discrete element method （DEM）. This numerical method allows us to incorporate the fluid-particle interactions （drag force, contact force, Saffman lift force, Magnus lift force, buoyancy force） using momentum exchange and the non-Newtonian behavior of the fluid. The interactions of particle-particle, particle-wall, and particle-drill pipe are taken into account with the Hertz-Mindlin model. We compare the transport of spheres with non-spherical particles （non-smooth sphere, disc, and cubic） constructed via the multi- sphere method for a range of fluid inlet velocities and drill pipe inclination angles. The simulations are carried out for laboratory-scale drilling configurations. Our results demonstrate good agreement with published experimental data. We evaluate the fluid-particle flow patterns, the particle velocities, and the particle concentration profiles. The results reveal that particle sphericity plays a major role in the fluid-solid interaction. The traditional assumption of an ideal spherical particle may cause inaccurate results.

NUMERICAL SIMULATION OF TURBULENT SPOTS IN INCLINED OPEN-CHANNEL FLOW

The generation and evolution of turbulent spots in the open-channel flow are simulated numerically by using the Navier-Stokes equations. An effective numerical method with high accuracy and high resolution is developed. The fourth-order time splitting methods with high accuracy is proposed. Three-dimensional coupling difference methods are presented for the spatial discretization of the Poisson equation of pressure and Hemholtz equations of velocity, therefore, the fourth-order three-dimensional coupling central difference schemes are constituted. The fourth-order explicit upwind-biased compact difference schemes are designed to overcome the difficulty for the general higher-order central difference scheme which is inadaptable in the boundary neighborhood. The iterative algorithm and overall time marching is used to enhance efficiency. The method is applied in the numerical simulation of turbulent spots at various complex boundary conditions and flow domains. The generation and the developing process of turbulent spots are given, and the basic characteristics of turbulent spots are shown by simulating the evolution of the wall pulse in inclined open-channel flow.