Three-dimensional magnetotellurics modeling using edgebased finite-element unstructured meshes

Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficiently simulating 3D MT responses. The electric field curl-curl equation in the frequency domain was used to deduce the H （curl） variation weak form of the MT forward problem, the Galerkin rule was used to derive a linear finite-element equation on the linear-edge tetrahedroid space, and, finally, a BI-CGSTAB solver was used to estimate the unknown electric fields. A local mesh refinement technique in the neighbor of the measuring MT stations was used to greatly improve the accuracies of the numerical solutions. Four synthetic models validated the powerful performance of our algorithms. We believe that our method will effectively contribute to processing more complex MT studies.

Three-dimensional land FD-CSEM forward modeling using edge finite-element method

A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.

Indirect tension test of epoxy asphalt mixtureusing microstructural finite-element model

A finite-element model of the thermosetting epoxy asphalt mixture（EAM） microstructure is developed to simulate the indirect tension test（IDT）.Image techniques are used to capture the EAM microstructure which is divided into two phases：aggregates and mastic.A viscoelastic constitutive relationship,which is obtained from the results of a creep test,is used to represent the mastic phase at intermittent temperatures.Model simulation results of the stiffness modulus in IDT compare favorably with experimental data.Different loading directions and velocities are employed in order to account for their influence on the modulus and the localized stress of the microstructure model.It is pointed out that the modulus is not consistent when the loading direction changes since the heterogeneous distribution of the mixture internal structure,and the loading velocity affects the localized stress as a result of the viscoelasticity of the mastic.The study results can provide a theoretical basis for the finite-element method,which can be extended to the numerical simulations of asphalt mixture micromechanical behavior.

Perfectly matched layer-absorbing boundary condition for finite-element time-domain modeling of elastic wave equations

The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.

Locating Error Considering Dimensional Errors Modeling for Multistation Manufacturing System

Multistation machining process is widely applied in contemporary manufacturing environment. Modeling of variation propagation in multistation machining process is one of the most important research scenarios. Due to the existence of multiple variation streams, it is challenging to model and analyze variation propagation in a multi-station system. Current approaches to error modeling for multistation machining process are not explicit enough for error control and ensuring final product quality. In this paper, a mathematic model to depict the part dimensional variation of the complex multistation manufacturing process is formulated. A linear state space dimensional error propagation equation is established through kinematics analysis of the influence of locating parameter variations and locating datum variations on dimensional errors, so the dimensional error accumulation and transformation within the multistation process are quantitatively described. A systematic procedure to build the model is presented, which enhances the way to determine the variation sources in complex machining systems. A simple two-dimensional example is used to illustrate the proposed procedures. Finally, an industrial case of multistation machining part in a manufacturing shop is given to testify the validation and practicability of the method. The proposed analytical model is essential to quality control and improvement for multistation systems in machining quality forecasting and design optimization.

Three-Dimensional TIN Algorithm for Digital Terrain Modeling

The problem of taking an unorganized point cloud in 3D space and fitting a polyhedral surface to those points is both important and difficult. Aiming at increasing applications of full three dimensional digital terrain surface modeling, a new algorithm for the automatic generation of three dimensional triangulated irregular network from a point cloud is pro- posed. Based on the local topological consistency test, a combined algorithm of constrained 3D Delaunay triangulation and region-growing is extended to ensure topologically correct reconstruction. This paper also introduced an efficient neighbor- ing triangle location method by making full use of the surface normal information. Experimental results prove that this algo- rithm can efficiently obtain the most reasonable reconstructed mesh surface with arbitrary topology, wherein the automati- cally reconstructed surface has only small topological difference from the true surface. This algorithm has potential applica- tions to virtual environments, computer vision, and so on.

Finite-element 3D modeling of stress patterns around a dipping fault

Stresses in a block around a dipping fracture simulating a damage zone of a fault are reconstructed by finite-element modeling. A fracture corresponding to a fault of different lengths, with its plane dipping at different angles, is assumed to follow a lithological interface and to experience either compression or shear. The stress associated with the destruction shows an asymmetrical pattern with different distances from the highest stress sites to the fault plane in the hanging and foot walls. As the dip angle decreases,the high-stress zone becomes wider in the hanging wall but its width changes negligibly in the foot wall.The length of the simulated fault and the deformation type affect only the magnitude of maximum stress,which remains asymmetrical relative to the fault plane. The Lh/Lfratio, where Lhand Lfare the widths of high-stress zones in the hanging and foot walls of the fault, respectively, is inversely proportional to the fault plane dip. The arithmetic mean of this ratio over different fault lengths in fractures subject to compression changes from 0.29 at a dip of 80°to 1.67 at 30°. In the case of shift displacement, ratios are increasing to 1.2 and 2.94, respectively.Usually they consider vertical fault planes and symmetry in a damage zone of faults. Following that assumption may cause errors in reconstructions of stress and fault patterns in areas of complex structural setting. According geological data, we know the structures are different and asymmetric in hanging and foot walls of fault. Thus, it is important to quantify zones of that asymmetry. The modeling results have to be taken into account in studies of natural faults, especially for practical applications in seismic risk mapping, engineering geology, hydrogeology, and tectonics.

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.

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.

Testing Rating Scale Unidimensionality Using the Principal Component Analysis (PCA)/<i>t</i>-Test Protocol with the Rasch Model: The Primacy of Theory over Statistics

Psychometric theory requires unidimensionality (i.e., scale items should represent a common latent variable). One advocated approach to test unidimensionality within the Rasch model is to identify two item sets from a Principal Component Analysis (PCA) of residuals, estimate separate person measures based on the two item sets, compare the two estimates on a person-by-person basis using t-tests and determine the number of cases that differ significantly at the 0.05-level;if ≤5% of tests are significant, or the lower bound of a binomial 95% confidence interval (CI) of the observed proportion overlaps 5%, then it is suggested that strict unidimensionality can be inferred;otherwise the scale is multidimensional. Given its proposed significance and potential implications, this procedure needs detailed scrutiny. This paper explores the impact of sample size and method of estimating the 95% binomial CI upon conclusions according to recommended conventions. Normal approximation, “exact”, Wilson, Agresti-Coull, and Jeffreys binomial CIs were calculated for observed proportions of 0.06, 0.08 and 0.10 and sample sizes from n= 100 to n= 2500. Lower 95%CI boundaries were inspected regarding coverage of the 5% threshold. Results showed that all binomial 95% CIs included as well as excluded 5% as an effect of sample size for all three investigated proportions, except for the Wilson, Agresti-Coull, and JeffreysCIs, which did not include 5% for any sample size with a 10% observed proportion. The normal approximation CI was most sensitive to sample size. These data illustrate that the PCA/t-test protocol should be used and interpreted as any hypothesis testing procedure and is dependent on sample size as well as binomial CI estimation procedure. The PCA/t-test protocol should not be viewed as a “definite” test of unidimensionality and does not replace an integrated quantitative/qualitative interpretation based on an explicit variable definition in view of the perspective, context and purpose of measurement.

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.

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.

Numerical Investigation of Thermal Behavior of CNC Machine Tool and Its Effects on Dimensional Accuracy of Machined Parts

The dimensional accuracy of machined parts is strongly influenced by the thermal behavior of machine tools (MT). Minimizing this influence represents a key objective for any modern manufacturing industry. Thermally induced positioning error compensation remains the most effective and practical method in this context. However, the efficiency of the compensation process depends on the quality of the model used to predict the thermal errors. The model should consistently reflect the relationships between temperature distribution in the MT structure and thermally induced positioning errors. A judicious choice of the number and location of temperature sensitive points to represent heat distribution is a key factor for robust thermal error modeling. Therefore, in this paper, the temperature sensitive points are selected following a structured thermomechanical analysis carried out to evaluate the effects of various temperature gradients on MT structure deformation intensity. The MT thermal behavior is first modeled using finite element method and validated by various experimentally measured temperature fields using temperature sensors and thermal imaging. MT Thermal behavior validation shows a maximum error of less than 10% when comparing the numerical estimations with the experimental results even under changing operation conditions. The numerical model is used through several series of simulations carried out using varied working condition to explore possible relationships between temperature distribution and thermal deformation characteristics to select the most appropriate temperature sensitive points that will be considered for building an empirical prediction model for thermal errors as function of MT thermal state. Validation tests achieved using an artificial neural network based simplified model confirmed the efficiency of the proposed temperature sensitive points allowing the prediction of the thermally induced errors with an accuracy greater than 90%.

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.

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.

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.

Testing a Four-Dimensional Variational Data Assimilation Method Using an Improved Intermediate Coupled Model for ENSO Analysis and Prediction

A four-dimensional variational （4D-Var） data assimilation method is implemented in an improved intermediate coupled model （ICM） of the tropical Pacific. A twin experiment is designed to evaluate the impact of the 4D-Var data assimilation algorithm on ENSO analysis and prediction based on the ICM. The model error is assumed to arise only from the parameter uncertainty. The ＂observation＂ of the SST anomaly, which is sampled from a ＂truth＂ model simulation that takes default parameter values and has Gaussian noise added, is directly assimilated into the assimilation model with its parameters set erroneously. Results show that 4D-Var effectively reduces the error of ENSO analysis and therefore improves the prediction skill of ENSO events compared with the non-assimilation case. These results provide a promising way for the ICM to achieve better real-time ENSO prediction.

A Novel Three-Dimensional Mouse Embryonic Implantation Model In Vitro

To regenerate three-dimensional endometrium in vitro as a novel model for studying the mechanism of implantation of embryos, the luminal epithelial cells and stromal cells of the rabbit uterus were separated and cultured in vitro. The type Ⅰ mouse tail collagen was used as scaffolding material. The stromal cells were inoculated in the type I mouse tail collagen, and the luminal epithelial cells were inoculated on the type i mouse tall collagen to regenerate the endometrium in vitro. The regenerated endometrium was cultured in DMEM-F/12 media containing 100 nmol L^-1 progesterone, 10 nM β-estradiol, and 10% fetal bovine serum （FBS） for 3 d. The media were then replaced with CZB containing 100 nM progesterone, 10 nmol L-1 β-estradiol, and 10% FBS, and the mouse blastulas were co-cultured with it. The results of scanning electronic micrography showed that the epithelial cells on the surface of the reconstructed endometrium were covered with numerous slender microvilli and some epithelial cells protruded pinopodes. After culturing for 12 h with the mouse blastula, the shedding, attachment, and implantation of the blastula were observed. The blastula can escape from zona pellucida and attach to the three-dimensional endometrium and is then implanted into it. This study showed that the reconstructed three-dimensional endometrium can serve as a robust embryo implantation model in vitro.

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.