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.

Three-Dimension Mathematical Model of Total Phosphor in the Reservoir and Application

Taking the transport of total phosphor pollutants in the Beijing Miyun reservoir for example,we have obtained three dimensional distributing regularity of total phosphor pollutants by the calculation of the linear interpolation value of each point between horizontal layers. The credibility analysis in allusion to this method was carried out and the programming scheme for realizing this method was set forth.

Model-based Predictive Control for Spatially-distributed Systems Using Dimensional Reduction Models

In this paper, a low-dimensional multiple-input and multiple-output （MIMO） model predictive control （MPC） configuration is presented for partial differential equation （PDE） unknown spatially-distributed systems （SDSs）. First, the dimension reduction with principal component analysis （PCA） is used to transform the high-dimensional spatio-temporal data into a low-dimensional time domain. The MPC strategy is proposed based on the online correction low-dimensional models, where the state of the system at a previous time is used to correct the output of low-dimensional models. Sufficient conditions for closed-loop stability are presented and proven. Simulations demonstrate the accuracy and efficiency of the proposed methodologies.

Simulation of two-dimensional interior ballistics model of solid propellant electrothermal-chem ical launch with discharge rod plasma generator

Instead of the capillary plasma generator(CPG),a discharge rod plasma generator(DRPG)is used in the30 mm electrothermal-chemical(ETC)gun to improve the ignition uniformity of the solid propellant.An axisymmetric two-dimensional interior ballistics model of the solid propellant ETC gun(2D-IB-SPETCG)is presented to describe the process of the ETC launch.Both calculated pressure and projectile muzzle velocity accord well with the experimental results.The feasibility of the 2D-IB-SPETCG model is proved.Depending on the experimental data and initial parameters,detailed distribution of the ballistics parameters can be simulated.With the distribution of pressure and temperature of the gas phase and the propellant,the influence of plasma during the ignition process can be analyzed.Because of the radial flowing plasma,the propellant in the area of the DRPG is ignited within 0.01 ms,while all propellant in the chamber is ignited within 0.09 ms.The radial ignition delay time is much less than the axial delay time.During the ignition process,the radial pressure difference is less than 5 MPa at the place 0.025 m away from the breech.The radial ignition uniformity is proved.The temperature of the gas increases from several thousand K(conventional ignition)to several ten thousand K(plasma ignition).Compare the distribution of the density and temperature of the gas,we know that low density and high temperature gas appears near the exits of the DRPG,while high density and low temperature gas appears at the wall near the breech.The simulation of the 2D-IB-SPETCG model is an effective way to investigate the interior ballistics process of the ETC launch.The 2D-IB-SPETC model can be used for prediction and improvement of experiments.

Seismic data reconstruction based on low dimensional manifold model

Seismic data reconstruction is an essential and yet fundamental step in seismic data processing workflow,which is of profound significance to improve migration imaging quality,multiple suppression effect,and seismic inversion accuracy.Regularization methods play a central role in solving the underdetermined inverse problem of seismic data reconstruction.In this paper,a novel regularization approach is proposed,the low dimensional manifold model(LDMM),for reconstructing the missing seismic data.Our work relies on the fact that seismic patches always occupy a low dimensional manifold.Specifically,we exploit the dimension of the seismic patches manifold as a regularization term in the reconstruction problem,and reconstruct the missing seismic data by enforcing low dimensionality on this manifold.The crucial procedure of the proposed method is to solve the dimension of the patches manifold.Toward this,we adopt an efficient dimensionality calculation method based on low-rank approximation,which provides a reliable safeguard to enforce the constraints in the reconstruction process.Numerical experiments performed on synthetic and field seismic data demonstrate that,compared with the curvelet-based sparsity-promoting L1-norm minimization method and the multichannel singular spectrum analysis method,the proposed method obtains state-of-the-art reconstruction results.

Modeling and validating three dimensional human normal cervix and cervical cancer tissues in vitro

The use of three dimensional in vitro systems in cancer research is a promising path for developing effective anticancer therapies.The aim of this study was to engineer a functional 3-D in vitro model of normal and cancerous cervical tissue.Normal epithelial and immortalized cervical epithelial carcinoma cell lines were used to construct 3-D artificial normal cervical and cervical cancerous tissues.De-epidermised dermis（DED） was used as a scaffold for both models.Morphological analyses were conducted by using hematoxylin and eosin staining and characteristics of the models were studied by analyzing the expression of different structural cytokeratins and differential protein marker MAX dimerisation protein 1（Mad1） using immunohistochemical technique.Haematoxylin and eosin staining results showed that normal cervical tissue had multi epithelial layers while cancerous cervical tissue showed dysplastic changes.Immunohistochemistry staining revealed that for normal cervix model cytokeratin 10 was expressed in the upper stratified layer of the epithelium while cytokeratin 5 was expressed mainly in the middle and basal layer.Cytokeratin 19 was weakly expressed in a few basal cells.Cervical cancer model showed cytokeratin 19 expression in different epithelial layers and weak or no expression for cytokeratin 5 and cytokeratin 10.Madl expression was detected in some suprabasal cells.The 3-D in vitro models showed stratified epithelial layers and expressed the same types and patterns of differentiation marker proteins as seen in corresponding in vivo tissue in either normal cervical or cervical cancerous tissue.These findings imply that they can serve as functional normal and cervical cancer models.

Research on dynamic updating of three dimensional geological modeling based on the OO-Solid model

The dynamic updating of the model included: the change of space border,addi- tion and reduction of spatial component (disappearing,dividing and merging),the change of the topological relationship and synchronous dynamic updating of database.Firstly, arming at the deficiency of OO-Solid model in the aspect of dynamic updating,modeling primitives of OO-Solid model were modified.And then the algorithms of dynamic updating of 3D geological model with the node data,line data or surface data change were dis- cussed.The core algorithms was done by establishing space index,following the way of facing the object from bottom to top,namely the dynamic updating from the node to arc, and then to polygon,then to the face of the component and finally to the geological object. The research has important theoretical and practical values in the field of three dimen- sional geological modeling and is significant in the field of mineral resources.

A 3-DIMENSIONAL DATA MODEL FOR VISUALIZING CLOVERLEAF JUNCTION IN A CITY MODEL

The research work has been seldom done about cloverleaf junction expression in a 3-dimensional city model (3DCM). The main reason is that the cloverleaf junction is often in a complex and enormous construction. Its main body is bestraddle in air,and has aerial intersections between its parts. This complex feature made cloverleaf junction quite different from buildings and terrain, therefore, it is difficult to express this kind of spatial objects in the same way as for buildings and terrain. In this paper,authors analyze spatial characteristics of cloverleaf junction, propose an all-constraint points TIN algorithm to partition cloverleaf junction road surface, and develop a method to visualize cloverleaf junction road surface using TIN. In order to manage cloverleaf junction data efficiently, the authors also analyzed the mechanism of 3DCM data management, extended BLOB type in relational database, and combined R-tree index to manage 3D spatial data. Based on this extension, an appropriate data

Three-dimensional numerical modeling of gravity anomalies based on Poisson equation in spacewavenumber mixed domain

In gravity-anomaly-based prospecting, the computational and memory requirements for practical numerical modeling are potentially enormous. Achieving an efficient and precise inversion for gravity anomaly imaging over large-scale and complex terrain requires additional methods. To this end, we have proposed a new topography-capable By performing a two-dimensional Fourier transform in the horizontal directions, threedimensional partial differential equations in the spatial domain were transformed into a group of independent, one-dimensional differential equations engaged with different wave numbers. These independent differential equations are highly parallel across different wave numbers. differential equations with different wave numbers, and the efficiency of solving fixedbandwidth linear equations was further improved by a chasing method. In a synthetic test, a prism model was used to verify the accuracy and reliability of the proposed algorithm by comparing the numerical solution with the analytical solution. We studied the computational precision and efficiency with and without topography using different Fourier transform methods. The results showed that the Guass-FFT method has higher numerical precision, while the standard FFT method is superior, in terms of computation time, for inversion and quantitative interpretation under complicated terrain.

Risk Assessment of Coal Mine Water Hazard Based on Three-Dimensional Model of Geology and Underground Space

Mine safety have top-five disasters,which including the water,gas,fire,dust and geological dynamic disaster.The coal mine water disaster is one of the important factors which restricted the development of China’s coal production.It is showed by statistics that 60%of mine accidents are affected by groundwater,which not only result in the production losses,casualties and a variety of

A NEW UNSTEADY THREE DIMENSIONAL MODEL FOR MACROMOLECULAR TRANSPORT AND WATER FILTRATION ACROSS THE ARTERIAL WALL

A new unsteady three-dimensional convective-diffusive mathematical model for the transportation of macromolecules and water across the arterial wall was proposed . After the formation of leaky junctions due to the mitosis of endothelial cell of the arterial wall, the macromolecular transport happens surrounding the leaky cells. The arterial wall was divided into four layers: the endothelial layer, the subendothelial intima, the internal elastic lamina and the media for the convenience of research. The time-dependent concentration growth,the effect of the shape of endothelial cell and the effect of physiological parameters were analyzed. The analytical solution of velocity field and pressure field of water flow across the arterial wall were obtained; and concentration distribution of three macromolecules ; LDL,HRP and Albumin, were calculated with numerical simulation method. The new theory predicts, the maximum and distribution areas of time dependent concentration with round shape endothelial cell are both larger than that with ellipse-shape endothelial cell. The model also predicts the concentration growth is much alike that of a two-dimensional model and it shows that the concentration reaches its peak at the leaky junction where atherosclerotic formation frequently occurs and falls down rapidly in a limited area beginning from its earlier time growth to the state when macromolecular transfer approaches steadily. These predictions of the new model are in agreement with the experimental observation for the growth and concentration distribution of LDL and Albumin.

The study on three-dimensional numerical model and fronts of the Jiulong Estuary and the Xiamen Bay

Applying the methods of on-site observation and dynamic model, the research on the fronts at the Jiulong Estuary has been carried out, during which spatial and temporal distribution, dynamic characteristics and formation mechanism of salinity fronts are analyzed and discussed. The research shows that the estuarine fronts mainly lie in the area from the Jiyu Islet to the Haimen Island, outside of Yuweizai to Hulishan cross-section, the near coast of Yuweizai and the south of the Songyu-Gulangyu Channel. The fronts in the former two regions are formed directly by plume, while the one near the coast of Yuweizai is a tidal intrusion front caused by flood current and the one at the south of the Songyu-Gulangyu Channel is the result of current shear transformation. Under normal circumstances, fresh water of the Jiulong River mainly influences the inside of the Xiamen Bay, and when it is in typhoon seasons, plume front can affect the Taiwan Strait and has an effect on the biogeochemical Drocesses in the strait.

Three-Dimensional Model Test for Port Engineering with Multi-Directional Waves

-The necessity of using irregular waves, especially multi- directional waves to conduct three-dimensional model tests for port engineering and the test method are described in this paper through an example of model test for a port. The test results show that a deep navigation channel has a large effect on the waves in front of the breakwater near the port entrance and on the wave condition in the port.

Numerical Simulation for Two-Phase Water Hammer Flows in Pipe by Quasi-Two-Dimensional Model

The features of a quasi-two-dimensional( quasi-2D) model for simulating two-phase water hammer flows with vaporous cavity in a pipe are investigated. The quasi-2D model with discrete vaporous cavity in the pipe is proposed in this paper. This model uses the quasi-2D model for pure liquid zone and one-dimensional( 1D) discrete vapor cavity model for vaporous cavity zone. The quasi-2D model solves two-dimensional equations for both axial and radial velocities and 1D equations for both pressure head and discharge by the method of characteristics. The 1D discrete vapor cavity model is used to simulate the vaporous cavity occurred when the pressure in the local pipe is lower than the vapor pressure of the liquid. The proposed model is used to simulate two-phase water flows caused by the rapid downstream valve closure in a reservoir-pipe-valve system.The results obtained by the proposed model are compared with those by the corresponding 1D model and the experimental ones provided by the literature,respectively. The comparison shows that the maximum pressure heads simulated by the proposed model are more accurate than those by the corresponding 1D model.

Performance assessment of two-dimensional hydraulic models for generation of flood inundation maps in mountain river basins

Hydraulic models for the generation of flood inundation maps are not commonly applied in mountain river basins because of the difficulty in modeling the hydraulic behavior and the complex topography. This paper presents a comparative analysis of the performance of four twodimensional hydraulic models (HEC-RAS 2D, Iber 2D, Flood Modeller 2D, and PCSWMM 2D) with respect to the generation of flood inundation maps. The study area covers a 5-km reach of the Santa B-arbara River located in the Ecuadorian Andes, at 2330 masl, in Gualaceo. The model's performance was evaluated based on the water surface elevation and flood extent, in terms of the mean absolute difference and measure of fit. The analysis revealed that, for a given case, Iber 2D has the best performance in simulating the water level and inundation for flood events with 20- and 50-year return periods, respectively, followed by Flood Modeller 2D, HEC-RAS 2D, and PCSWMM 2D in terms of their performance. Grid resolution, the way in which hydraulic structures are mimicked, the model code, and the default value of the parameters are considered the main sources of prediction uncertainty.

Development and application of a two-dimensional water quality model for the Daqinghe River Mouth of the Dianchi Lake

Water quality models are important tools to support the optimization of aquatic ecosystem rehabilitation programs and assess their efficiency. Basing on the flow conditions of the Daqinghe River Mouth of the Dianchi Lake, China, a two-dimensional water quality model was developed in the research. The hydrodynamics module was numerically solved by the alternating direction iteration （ADI） method. The parameters of the water quality module were obtained through the in situ experiments and the laboratory analyses that were conducted from 2006 to 2007. The model was calibrated and verified by the observation data in 2007. Among the four modelled key variables, i.e., water level, COD （in CODcr）, NH4＋-N and PO43-P the minimum value of the coefficient of determination （COD） was 0.69, indicating the model performed reasonably well. The developed model was then applied to simulate the water quality changes at a downstream cross-section assuming that the designed restoration programs were implemented. According to the simulated results, the restoration programs could cut down the loads of COD and PO43-P about 15%. Such a load reduction, unfortunately, would have very little effect on the NH4^＋-N removal. Moreover, the water quality at the outlet cross-section would be still in class V （3838-02）, indicating more measures should be taken to further reduce the loads. The study demonstrated the capability of water quality models to support aquatic ecosystem restorations.

Global Two-Dimensional Chemistry Model and Simulation of Atmospheric Chemical Composition

A global two-dimensional zonally averaged chemistry model is developed to study the chemi-cal composition of atmosphere. The region of the model is from 90°S to 90°N and from the ground to the altitude of 20 km with a resolution of 5° x 1 km. The wind field is residual circulation calcu-lated from diabatic rate. 34 species and 104 chemical and photochemical reactions are considered in the model. The sources of CH4, CO and NOx, which are divided into seasonal sources and non-seasonal sources, are parameterized as a function of latitude and time. The chemical composi-tion of atmosphere was simulated with emission level of CH4, CO and NOx in 1990. The results are compared with observations and other model results, showing that the model is successful to simu-late the atmospheric chemical composition and distribution of CH4. Key words Global two-dimensional chemistry model - Atmospheric composition - Emission This work was supported by the State Key Program for basic research “ Climate Dynamics and Cli-mate Prediction Theory” (Pandeng-yu-21).The authors would like to express their thanks to the National Oceanic and Atmospheric Administration (NOAA), Climate Monitoring and Diagnostics Laboratory (CMDL), Carbon Cycle Group for providing the observational data of CO and CH4.

Numerical Simulation for One Dimensional Steady Quasineutral Hybrid Model of Stationary Plasma Thruster

Based on the analysis of the physical mechanism of the Stationary Plasma Thruster （SPT）, an integral equation describing the ion density of the steady SPT and the ion velocity distribution function at an arbitrary axial position of the steady SPT channel are derived. The integral equation is equivalent to the Vlasov equation, but the former is simpler than the latter. A one dimensional steady quasineutral hybrid model is established. In this model, ions are described by the above integral equation, and neutrals and electrons are described by hydrodynamic equations. The transferred equivalency to the differential equation and the integral equation, together with other equations, are solved by an ordinary differential equation （ODE） solver in the Matlab. The numerical simulation results show that under various circumstances, the ion average velocity would be different and needs to be deduced separately.

THREE-DIMENSIONAL MATHEMATICAL MODEL AND ITS APPLICATION IN THE NEIGHBORHOOD OF THE THREE GORGES RESERVOIR DAM IN THE YANGTZE RIVER

The calculation of flow and sediment transport is one of the most important tasks in river engineering. The task is particularly difficult because a number of complex physical phenomena should be accounted for more realistically in a model with a predictive power. Three-dimensional calculations of river flow and suspended sediment transport are performed in this paper with application in the Three Gorges Reservoir in the Yangtze River. A period of 76 years after the dam is built is simulated and the results are compared with laboratory measurements obtained by Tsinghua University whereby the model is verified and calibrated. Generally speaking, the calculated results agree well with the experiments, demonstrating that the present model can be used for flow and sediment transport prediction in rivers.

Numerical Simulation of Saline Intrusion and Purging in Long SeaOutfalls with Three-Dimensional Model

Saline intrusion into marine sewage ouffalls will greatly decrease the efficiency of sewage disposal. In order to investigate the mechanisms of this flow, in this paper, a three-dimensional numerical model based on FVM （Finite Volume Method） is established, The RNG κ-ε model is selected for turbulence modeling. The time-averaged vohtme fraction equations are introduced to simulate the stratification and inteffaeial exchange of sewage and seawater in outfalls. Validity of the established three-dimensional numerical model is evaluated by comparisons of numerical results with experimental data. With this three-dimensional numerical model, the internal flow characteristics in ouffalls for different sewage discharges are simulated. The results indicate that for a low sewage discharge, saline circulates in the outfall due to intrusion and both the inflowing momentum and the inteffaeial turbulent mixing are important mechanisms to extrude the saline. For a high sewage discharge, saline intrusion could be avoided. The inflow momentum is the main mechanism to extrude the saline and the inteffacial turbulent mixing is nut important relatively. Even at a high sewage discharge, the saline wedge would be retained in the main ouffall pipe after the risers are purged. It takes a long time for this saline wedge to be extruded by interracial turbulent mixing.