非饱和土壤水运动模型参数优选方法研究

针对数值反演过程中参数优选方法适用性不明确的问题,研究最小化目标函数、带有自适应差分演化Metropolis的马尔可夫链蒙特卡罗方法(MCMC-DREAM)对数值反演效果的影响,为探索高效的参数优选方法提供参考。数值算例反演结果表明:最小化目标函数方法计算复杂度较低,但对参数初值敏感,适合于对目标区域土壤有较为深入了解的场合使用;MCMC-DREAM对参数初值不敏感,但计算复杂度较高,适合于先验信息有限的场合使用。两种参数优选方法都存在"异参同效"现象,先验信息与敏感性分析有助于克服该问题,提高数值反演解决实际问题的能力。

MATLAB-Based Simulation of Buoyancy-Driven Underwater Glider Motion

The mass configuration of the buoyancy-driven underwater glider is decomposed and defined. The coupling between the glider body and its internal masses is addressed using the energy law. A glider motion model is established, and the corresponding simulation program is derived using MATLAB. The characteristics of the glider motion are explored using this program. The simula- tion results show that the basic characteristic of a buoyancy-driven underwater glider is the periodic alternation of downward and upward motions. The glider's spiral motion can be applied to missions in restricted regions. The glider's horizontal velocity, gliding depth and its motion radius in spiral motion can be changed to meet different application purposes by using different glider parameter designs. The simulation also shows that the model is appropriate and the program has strong simulation functions.

Dynamics model of underwater robot motion control in 6 degrees of freedom

In order to analyze underwater robot control system dynamics features, a system 6-DOF dynamics model was founded. Underwater robot linear and nonlinear hydrodynamics were analyzed by Taylor series, based on general motion equation. Special control system motion equation was deduced by cluster of inertial items and non-inertial items. For program convenience, motion equation matrix format was presented. Experimental principles of screw propellers, rudders and wings were discussed. Experimental data least-square curve fitting, interpolation and their corresponding traditional equation helped us to obtain the whole system dynamic response procedure. A series of simulation experiments show that the dynamics model is correct and reliable. The model can provide theory proof for analyzing underwater robot motion control system physics characters and provide a mathematic model for traditional control method.

Ground movement mechanism in tectonic stress metal mines with steep structure planes

When mining metal mines with steep structure planes by the caving method,there is a mechanical model in which the horizontal stress on the rock mass is simplified as a column before surface subsidence.The model is used to deduce critical support load and limiting column length for a given horizontal stress and support pressure.Considering the impact of the column effect,a method is proposed to determine the movement of the ground and caving area in a mine.After surface subsidence,the horizontal stress on a surrounding rock mass can be simplified to a cantilever beam mechanical model.Expressions for its bending fracture length are deduced,and a method is given to determine its stability.On this basis,an explanation for the large ground movement and subsidence scope was given.A case study shows that the damage effect of column and cantilever beam is significant for ground movement in metal-ore mine,and an appropriate correction value should be applied when designing for its angle of ground movements.

CFD-Based Numerical Simulation of Water Film Flash Evaporation with a New Flash Evaporation Model

In this study, a new mass model involving superheat, initial temperature, liquid height, evaporator diameter, and flashing time is established to describe the flash evaporation process of water film. Of 469 sets of flash experimental data from three previous researches, 305 sets were applied to optimize parameters, and the other 164 sets were used to verify the practicability of the model. The results showed that the mean relative error between the literature data and the model values was less than 16.3%, and the model statistics proved that the model was well-posed. Then, the kinetic model was obtained using the time derivative of the new mass model. Computational fluid dynamics simulation of water film flash evaporation was studied based on a user-defined function program of the new evaporation kinetic model. The new kinetic model shows more consistency with the experimental phenomena in terms of evaporated mass and temperature compared with the evaporation–condensation model in Fluent software and Gopalakrishna's model. This new kinetic model can be extended to describe the flash process of water solution under other conditions.

Probability modeling for robustness of multivariate LQG designing based on ship lateral motion

The robustness of LQG designing for latitudinal movement of ship is mainly discussed, when its hydrodynamic parameters fluctuate around criterion value at random on the proportional distributing. When a given ship state at the speed of 18 kn and the course of 45° under Rank 5 state of sea, and the hydrodynamic parameters of the ship fluctuate at random on the proportional distributing with a range of ± 10% , ± 20% , ± 30% , the robustness of muhivariate LQG designing for ship is analyzed with applying the probability modeling of relative controlling effect, The result of simulating shows that when the hydrodynamic parameters of ship fluctuates the relative controlling effect of the LQG designing submit to normal distribution and the mean value of relative controlling effect has no remarkable changes comparing to that without perturbation of hydrodynamic parameter,

A Modified Logvinovich Model For Hydrodynamic Loads on an Asymmetric Wedge Entering Water with a Roll Motion

The water entry problem of an asymmetric wedge with roll motion was analyzed by the method of a modified Logvinovich model (MLM). The MLM is a kind of analytical model based on the Wagner method, which linearizes the free surface condition and body boundary condition. The difference is that the MLM applies a nonlinear Bernoulli equation to obtain pressure distribution, which has been proven to be helpful to enhance the accuracy of hydrodynamic loads. The Wagner condition in this paper was generalized to solve the problem of the water entry of a wedge body with rotational velocity. The comparison of wet width between the MLM and a fully nonlinear numerical approach was given, and they agree well with each other. The effect of angular velocity on the hydrodynamic loads of a wedge body was investigated.

Fractal Measures of Drainage Network to Investigate Surface Deformation from Remote Sensing Data: A Paradigm from Hindukush (NE-Afghanistan)

This approach represents the relative susceptibility of the topography of the earth to active deformation by means of geometrical distinctiveness of the river networks. This investigation employs the fractal analysis of drainage system extracted from ASTER Global Digital Elevation Model (GDEM-30m resolution). The objective is to mark active structures and to pinpoint the areas robustly influenced by neotectonics. This approach was examined in the Hindukush, NE-Afghanistan. This region is frequently affected by deadly earthquakes and the modern fault activities and deformation are driven by the collision between the northward-moving Indian subcontinent and Eurasia. This attempt is based on the fact that drainage system is strained to linearize due to neotectonic deformation. Hence, the low fractal dimensions of the Kabul, Panjsher, Laghman, Andarab, Alingar and Kocha Rivers are credited to active tectonics. A comprehensive textural examination is conducted to probe the linearization, heterogeneity and connectivity of the drainage patterns. The aspects for these natural textures are computed by using the fractal dimension (FD), lacunarity (LA) and succolarity (SA) approach. All these methods are naturally interrelated, i.e. objects with similar FD can be further differentiated with LA and/or SA analysis. The maps of FD, LA and SA values are generated by using a sliding window of 50 arc seconds by 50 arc seconds (50" × 50"). Afterwards, the maps are interpreted in terms of regional susceptibility to neotectonics. This method is useful to pinpoint numerous zones where the drainage system is highly controlled by Hindukush active structures. In the North-Northeast of the Kabul block, we recognized active tectonic blocks. The region comprising, Kabul, Panjsher, Andrab, Alingar and Badakhshan is more susceptible to damaging events. This investigation concludes that the fractal analysis of the river networks is a bonus tool to localize areas vulnerable to deadly incidents influencing the Earth’s topography and consequently intimidate human lives.

Effect of approach run velocity on the optimal performance of the triple jump

Purpose： The purpose of this study was to determine the effect of horizontal and vertical velocities at the landing of the last step of approach run on the performance and optimal phase ratio of the triple jump. Methods： Three-dimensional kinematic data of 13 elite male triple jumpers were obtained during a competition. Computer simulations were performed using a biomechanical model of the triple jump to determine the longest actual distance using the optimal phase ratio with altered horizontal and vertical velocities at the landing of the last step of approach run. Results： The actual distance obtained using the optimal phase ratio significantly increased as the horizontal velocity at the landing of the last step of approach run increased （p = 0.001） and the corresponding downward vertical velocity decreased （p = 0.001）. Increasing horizontal velocity at the landing of the last step of approach run decreased optimal hop percentage and increased optimal jump percentage （p = 0.001）, while decreasing corresponding downward vertical velocity increased optimal hop percentage and decreased optimal jump percentage （p = 0.001）. Conclusion： The effects of the velocities at the landing of the last step of approach run on the optimal phase ratio were generally small and did not qualitatively alter optimal techniques.

Effect of Solving Models on Euler Vectors

In the process of solving Euler vectors based on GNSS horizontal movement field,the number of estimated parameters can affect Euler vector results. This issue is analyzed through theoretical deduction and practical example in this paper. Firstly,the difference between the results of Euler vectors in different solving models is deduced. Meanwhile, based on GNSS horizontal movement field in the Chinese mainland from 2004 to 2007,two common models( RRM and REHSM) are used to discuss the impact of solving models on Euler vectors and the follow-up study. The result shows that the maximum value of the difference in a block's entire rotation can reach 2. 6mm /a,and should not be ignored. Therefore,the results of horizontal movement are different using different kinematic block models,and this should be paid more attention in the analysis of crustal horizontal movement.

Application of Mathematical Model to Assess Sediment Movement Pattern along Nearshore Coastal Water

Movement of sediment load and its pattern of transportation along nearshore coastal water is a very important phenomenon to be assessed for different sector of coastal Engineering. To develop and understand the physical processes responsible for shaping the ongoing evolution of the coast and to develop the management strategies to deal the impact of human activities on the coastal zone and as well as for adapting to the hazards associated with the people living on the coast, knowledge of the mechanism, processes and the pattern of sediment movement in the nearshore zone is of utmost importance. Nearshore zone is a very active area, where a series of dynamic processes occur in response to changing wave climates and sediment budgets. Nowadays mathematical modeling is an attractive alternative and becoming a very viable approach to study the sediment movement pattern with the advanced computational facilities and improved understanding on wave mechanics and sediment transport processes. It is very effective, reliable and also comfortable to study the pattern of sediment transportation including yield, distribution and management of sediment with the help of mathematical model. Validity of forecast in sediment transport depends on both mathematical modeling technique and boundary conditions.

The Development of a Geographic Information System （GIS） to Document Research in an Everglades Physical Model

The Loxahatchee Impoundment Landscape Assessment （LILA） facility is a unique physical model of the Everglades ecosystem. LILA has a closed-loop water delivery system and consists of four 0.08 square kilometer （-8 ha） macrocosms, created to be replicates of one another and of the Everglades landscape. Built in 2003, LILA＇s purpose is to provide scientists with an opportunity to design and implement research concerning Everglades restoration techniques in an accessible, controlled and replicated Everglades environment. Key Everglades habitats were sculpted within LILA： tree islands, ridges, sloughs and alligator holes. Water levels and flows in each macrocosm are controlled independently, so that researchers can study the effects of hydrology on Everglades landscape and ecology. Studies have focused upon measuring survival and growth of native trees planted on the tree islands; measuring surface water and ground water movement and chemistry; studying wading bird feeding and the movement of prey species （crayfish）; and measuring erosion and accretion on tree islands and ridges. We developed a Geographic Information System （GIS） data set to identify, characterize, and spatially reference the features of LILA and document research activities. This development included mapping the boundaries of the landscape features, creating a theoretical Digital Elevation Model （DEM） and describing the research projects being carried out. The creation of this GIS data set enhances the ability to schedule and coordinate research, assist scientists in the visualization and spatial representation of their research, and provide a resource for the storage, analysis and synthesis of valuable scientific information.

Discrete Boltzmann Method with Maxwell-Type Boundary Condition for Slip Flow

The rarefied effect of gas flow in microchannel is significant and cannot be well described by traditional hydrodynamic models. It has been known that discrete Boltzmann model(DBM) has the potential to investigate flows in a relatively wider range of Knudsen number because of its intrinsic kinetic nature inherited from Boltzmann equation.It is crucial to have a proper kinetic boundary condition for DBM to capture the velocity slip and the flow characteristics in the Knudsen layer. In this paper, we present a DBM combined with Maxwell-type boundary condition model for slip flow. The tangential momentum accommodation coefficient is introduced to implement a gas-surface interaction model.Both the velocity slip and the Knudsen layer under various Knudsen numbers and accommodation coefficients can be well described. Two kinds of slip flows, including Couette flow and Poiseuille flow, are simulated to verify the model.To dynamically compare results from different models, the relation between the definition of Knudsen number in hard sphere model and that in BGK model is clarified.

Impact of Lower Boundary Condition of Richards＇ Equation on Water, Energy, and Soil Carbon Based on Coupling Land Surface and Biogeochemical Models

Soil moisture has a significant influence on water, energy, and carbon biogeochemical cycles. A numerical method for solving Richards＇ equation is usually used for simulating soil moisture. Selection of a lower boundary condition for Richards＇ equation will further affect the simulation results for soil moisture, water cycle, energy balance, and carbon biogeochemical processes. In this study, the soil water movement dynamic sub-model of a hydrologically based land surface model, the variable infiltration capacity （VIC） model, was modified using the finite difference method （FDM） to solve a mixed form of Richards＇ equation. In addition, the VIC model was coupled with a terrestrial biogeochemical model, the Carnegie Ames Stanford Approach model of carbon, nitrogen, and phosphorus （CASACNP model）. The no-flux boundary （NB） and free-drainage boundary （FB） were selected to investigate their impacts on simulations of the water, energy, and soil carbon cycles based on the coupling model. The NB and FB had different influences on the water, energy, and soil carbon simulations. The water and energy simulations were more sensitive, while the soil carbon simulation was less sensitive to FB than to NB. Free-drainage boundary could result in lower soil moisture, evaporation, runoff, and heterotrophic respiration and higher surface soil temperature, sensible heat flux, and soil carbon content. The impact of the lower boundary condition on simulation would be greater with an increase in soil permeability. In the silt loam soil case, evaporation, runoff, and soil respiration of FB were nearly 169, 13%, and 1% smaller, respectively, compared to those of NB.

Initial results of China's GNSS-R airborne campaign: soil moisture retrievals

The global navigation satellite system reflectometry （GNSS-R） technique has been proven to be a powerful tool for retrieving geophysical parameters of ocean and land/hydrology processes. The ultimate goal for such GNSS-R applications is to achieve large-scale, all- weather, and full-time mapping using spaceborne platforms. In order to ensure both GNSS-R receiver and algorithm meet the requirements of spaceborne observations, airborne experimental campaigns need to be first carried out for early testing and validation purposes. This paper presents a first comprehensive overview of China＇s airborne GNSS-R campaign conducted on May 30, 2014. There were two objectives for this campaign： （l） to examine the capability of the GNSS-R receiver developed by the National Space Science Center, Chinese Academy of Sciences, for airborne observations and （2） to study algo-rithms for soil moisture and altimetry retrievals. In this paper, initial results of soil moisture retrievals are pre- sented. The left-hand circularly polarized-predominant satellite information was successfully used to retrieval soil moisture over the cropland. The right-hand circularly polarized components of the reflected signals were also received and examined. The GPS-derived soil moisture results, on the one hand, correctly represented the spatial variations of the soil moisture along the tracking of the flight; on the other hand, the results underestimated the ground-truth. Errors from the retrieval model and from the positioning and effects from the vegetation layer and from the atmospheric water vapor were the primary causes of the uncertainties in soil moisture retrievals using the airborne GNSS-R data. This airborne experimental campaign firstly investigate that China has the capability to perform airborne GNSS-R observation using the self-developed receiver, although the receiver developed by the NSSC needs to be further examined for its capability for spaceborne observation. The early findings of this study will provide illustrations for planned future airborne campaigns.

A two-dimensional Richards equation solver based on CVODE for variably saturated soil water movement

A numerical model of two-dimensional soil water movement under surface drip irrigation condition was developed. The physical process of soil water movement is described by 2D Richards equation,and the upper boundary condition is depicted by the improved moving ponded area boundary. The partial differential equation(PDE) is transformed into ordinary differential equations(ODEs) through spatial semi-discretization and numerically solved by an ordinary differential equation solver(CVODE) . The numerical and field experiments indicate the good performance of the model in terms of accuracy and efficiency. The model provides a useful tool for long-term simulation of soil water movement under surface drip irrigation. Also,the model can serve as a general 2D Richards equation solver for variably saturated soil water movement,which is named as TIVS model(Tsinghua Integrated Variably Saturated soil water movement model).