Implementation of a particle-in-cell method for the energy solver in 3D spherical geodynamic modeling

The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially in the presence of sharp thermal gradients,such as when modeling subducting slabs and rising plumes.This phenomenon prohibits the correct representation of thermal evolution and may cause incorrect implications of geodynamic processes.After examining several approaches for removing these numerical oscillations,we show that the Lagrangian method provides an ideal way to solve this problem.In this study,we propose a particle-in-cell method as a strategy for improving the solution to the energy equation and demonstrate its effectiveness in both one-dimensional and three-dimensional thermal problems,as well as in a global spherical simulation with data assimilation.We have implemented this method in the open-source finite-element code CitcomS,which features a spherical coordinate system,distributed memory parallel computing,and data assimilation algorithms.

Three-dimensional forward modeling of DC resistivity using the aggregation-based algebraic multigrid method

To speed up three-dimensional （3D） DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method （AGMG）. We first discretize the differential equation of the secondary potential field with mixed boundary conditions by using a seven-point finite-difference method to obtain a large sparse system of linear equations. Then, we introduce the theory behind the pairwise aggregation algorithms for AGMG and use the conjugate-gradient method with the V-cycle AGMG preconditioner （AGMG-CG） to solve the linear equations. We use typical geoelectrical models to test the proposed AGMG-CG method and compare the results with analytical solutions and the 3DDCXH algorithm for 3D DC modeling （3DDCXH）. In addition, we apply the AGMG-CG method to different grid sizes and geoelectrical models and compare it to different iterative methods, such as ILU-BICGSTAB, ILU-GCR, and SSOR-CG. The AGMG-CG method yields nearly linearly decreasing errors, whereas the number of iterations increases slowly with increasing grid size. The AGMG-CG method is precise and converges fast, and thus can improve the computational efficiency in forward modeling of three-dimensional DC resistivity.

Modeling constitutive relationship of 6013 aluminum alloy during hot plane strain compression based on Kriging method

Hot plane strain compression tests of 6013 aluminum alloy were conducted within the temperature range of 613?773 K and the strain rate range of 0.001?10 s?1. Based on the corrected experimental data with temperature compensation, Kriging method is selected to model the constitutive relationship among flow stress, temperature, strain rate and strain. The predictability and reliability of the constructed Kriging model are evaluated by statistical measures, comparative analysis and leave-one-out cross-validation (LOO-CV). The accuracy of Kriging model is validated by the R-value of 0.999 and the AARE of 0.478%. Meanwhile, its superiority has been demonstrated while comparing with the improved Arrhenius-type model. Furthermore, the generalization capability of Kriging model is identified by LOO-CV with 25 times of testing. It is indicated that Kriging method is competent to develop accurate model for describing the hot deformation behavior and predicting the flow stress even beyond the experimental conditions in hot compression tests.

Phase-field modeling of dendritic growth under forced flow based on adaptive finite element method

A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.

Forward modeling of marine DC resistivity method for a layered anisotropic earth

Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell＇s equations for an anisotropic model, we introduce scalar potentials based on the divergence-free characteristic of the electric and magnetic （EM） fields. We then continue the EM fields down into the deep earth and upward into the seawater and couple them at the ocean bottom to the transmitting source. By studying both the DC apparent resistivity curves and their polar plots, we can resolve the anisotropy of the ocean bottom. Forward modeling of a high-resistivity thin layer in an anisotropic half-space demonstrates that the marine DC resistivity method in shallow water is very sensitive to the resistive reservoir but is not influenced by airwaves. As such, it is very suitable for oil and gas exploration in shallowwater areas but, to date, most modeling algorithms for studying marine DC resistivity are based on isotropic models. In this paper, we investigate one-dimensional anisotropic forward modeling for marine DC resistivity method, prove the algorithm to have high accuracy, and thus provide a theoretical basis for 2D and 3D forward modeling.

Nonlinear Auto-Companding Method for Behavioral Modeling of Switched-Current Circuits

A wavelet collocation method with nonlinear auto companding is proposed for behavioral modeling of switched current circuits.The companding function is automatically constructed according to the initial error distribution obtained through approximating the input output function of the SI circuit by conventional wavelet collocation method.In practical applications,the proposed method is a general purpose approach,by which both the small signal effect and the large signal effect are modeled in a unified formulation to ease the process of modeling and simulation.Compared with the published modeling approaches,the proposed nonlinear auto companding method works more efficiently not only in controlling the error distribution but also in reducing the modeling errors.To demonstrate the promising features of the proposed method,several SI circuits are employed as examples to be modeled and simulated.

Modeling Rice rgMT as a Plant Metallothionein-Like Protein by the Distance Geometry and Homology Methods

Rice metallothionein-like protein (rgMT) shows characteristics of a three-section pattern composed of two highly conserved cysteine rich (CR) domains in the terminals and a spacer without cysteine (cys) residues in the center of the molecule. In this paper, the two CR domains and the spacer region were modeled by the distance geometry and homology methods separately. For the CR domains, twenty random models were generated for each cys combination based on the constraint conditions of CXC (C represents cys, X represents any amino acid other than cys), and CXXC motifs and a metal-sulfur chelating cluster. Four models for the N-terminal and two for C-terminal CR domain containing metal chelating structures formed by different combinations of cys were selected from 900 possible conformations. The GOR method was used to predict the secondary structure of the spacer region and its model was built by the homology method. After three parts of the protein were modeled, they were connected to form a three-dimensional structure model of rgMT. The whole conformation showed that rgMT could form two independent metal-sulfur chelating structures connected by a spacer peptide, without a structural or energy barrier for them to form two independent metal-chelating clusters just as mammalian metallothionein (MT) proteins. As all plant metallothionein-like (MT-L) proteins have the same primary structural characteristic, two CR domains connected by a spacer region, and many have the same cys arrangement pattern as rgMT, the three-dimensional structure model of rgMT will provide an important reference for the structural study of other plant MT-L proteins.

Precision of meshfree methods and application to forward modeling of two-dimensional electromagnetic sources

Meshfree method offers high accuracy and computational capability and constructs the shape function without relying on predefined elements. We comparatively analyze the global weak form meshfree methods, such as element-free Galerkin method （EFGM）, the point interpolation method （PIM）, and the radial point interpolation method （RPIM）. Taking two dimensional Poisson equation as an example, we discuss the support-domain dimensionless size, the field nodes, and background element settings with respect to their effect on calculation accuracy of the meshfree method. RPIM and EFGM are applied to controlled- source two-dimensional electromagnetic modeling with fixed shape parameters. The accuracy of boundary conditions imposed directly and by a penalty function are discussed in the case of forward modeling of two-dimensional magnetotellurics in a homogeneous medium model. The coupling algorithm of EFG-PIM and EFG-RPIM are generated by integrating the PIM or RPIM and EFGM. The results of the numerical modeling suggest the following. First, the proposed meshfree method and corresponding coupled methods are well-suited for electromagnetic numerical modeling. The accuracy of the algorithm is the highest when the support-domain dimensionless size is 1.0 and the distribution of field nodes is consistent with the nodes of background elements. Second, the accuracy of PIM and RPIM are lower than that of EFGM for the Poisson equation but higher than EFGM for the homogeneous medium MT response. Third, RPIM overcomes the matrix inversion problem of PIM and has a wider selection of support-domain dimensionless sizes as compared to RPIM.

3-D Geological Modeling-Concept,Methods and Key Techniques

3-D geological modeling plays an increasingly important role in Petroleum Geology, Mining Geology and Engineering Geology. The complexity of geological conditions requires different modeling methods in different situations. This paper summarizes the general concept of geological modeling; compares the characteristics of borehole-based modeling, cross-section based modeling and multi- source interactive modeling; analyses key techniques in 3-D geological modeling; and highlights the main difficulties and directions of future studies.

Wall-modeling for large-eddy simulation of flows around an axisymmetric body using the diffuse-interface immersed boundary method

A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equation is integrated along the wall-normal direction to link the tangential component of the effective body force for the IB method to the wall shear stress predicted by the wall model;(ii) a set of Lagrangian points near the wall are introduced to compute the normal component of the effective body force for the IB method by reconstructing the normal component of the velocity. This novel method will be a classical direct-forcing IB method if the grid is fine enough to resolve the flow near the wall. The method is used to simulate the flows around the DARPA SUBOFF model. The results obtained are well comparable to the measured experimental data and wall-resolved LES results.

Modeling Method for Flexible Energy Behaviors in CNC Machining Systems

CNC machining systems are inevitably confronted with frequent changes in energy behaviors because they are widely used to perform various machining tasks. It is a challenge to understand and analyze the flexible energy behaviors in CNC machining systems. A method to model flexible energy behaviors in CNC machining systems based on hierarchical objected-oriented Petri net(HOONet) is proposed. The structure of the HOONet is constructed of a high-level model and detail models. The former is used to model operational states for CNC machining systems, and the latter is used to analyze the component models for operational states. The machining parameters having great impacts on energy behaviors in CNC machining systems are declared with the data dictionary in HOONet models. A case study based on a CNC lathe is presented to demonstrate the proposed modeling method. The results show that it is effective for modeling flexible energy behaviors and providing a fine-grained description to quantitatively analyze the energy consumption of CNC machining systems.

Numerical Method for Modeling the Constitutive Relationship of Sand under Different Stress Paths

A numerical method was used in order to establish the constitutive relationship of sands under different stress paths, Firstly, based on the numerical method modeling the constitutive law of sands, the elastoplastic constitutive relationship of sand was established for three paths： the constant proportion of principle stress path, the conventional triaxial compression （CTC） path, and the p=constant （TC） path. The yield lines of plastic volumetric strain and plastic generalized shear strain were given. Through visualization, the three dimensional surface of the stress-strain relationship in the whole stress field （p, q） obtained under the three paths was plotted. Also, by comparing the stress-strain surfaces and yield locus of the three stress paths, the differences were found to be obvious, which demonstrates that the influence of the stress paths on constitutive law was not neglected. The numerical modeling method overcame the difficulty of finding an analytical expression for plastic potential. The results simulated the experimental data with an accuracy of 90% on average, so the constitutive model established in this paper provides an effective constitutive equation for this kind of engineering, reflecting the effect of practical stress paths that occur in sands.

A Finite Volume Method with Unstructured Triangular Grids for Numerical Modeling of Tidal Current

The finite volume method （FVM） has many advantages in 2-D shallow water numerical simulation. In this study, the finite volume method is used with unstructured triangular grids to simulate the tidal currents. The Roe scheme is applied in the calculation of the intercell numerical flux, and the MUSCL method is introduced to improve its accuracy. The time integral is a two-step scheme of forecast and revision. For the verification of the present method, the Stoker＇s problem is calculated and the result is compared with the mathematically analytic solutions. The comparison indicates that the method is feasible. A sea area of a port is used as an example to test the method established here. The result shows that the present computational method is satisfactory, and it could be applied to the engineering fields.

Microstructure modeling and virtual test of asphalt mixture based on three-dimensional discrete element method

The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete element method.A randomly generating algorithm was proposed to capture the three-dimensional irregular shape of coarse aggregate.And then,modeling algorithm and method for graded aggregates were built.Based on the combination of modeling of coarse aggregates,asphalt mastic and air voids,three-dimensional virtual sample of asphalt mixture was modeled by using PFC^(3D).Virtual tests for penetration test of aggregate and uniaxial creep test of asphalt mixture were built and conducted by using PFC^(3D).By comparison of the testing results between virtual tests and actual laboratory tests,the validity of the microstructure modeling and virtual test built in this study was verified.Additionally,compared with laboratory test,the virtual test is easier to conduct and has less variability.It is proved that microstructure modeling and virtual test based on three-dimensional discrete element method is a promising way to conduct research of asphalt mixture.

Infinitesimal dividing modeling method for dual suppliers inventory model with random lead times

As one of the basic inventory cost models, the （Q, τ）inventory cost model of dual suppliers with random procurement lead time is mostly formulated by using the concepts of ＂effective lead time＂ and ＂lead time demand＂, which may lead to an imprecise inventory cost. Through the real-time statistic of the inventory quantities, this paper considers the precise （Q, τ） inventory cost model of dual supplier procurement by using an infinitesimal dividing method. The traditional modeling method of the inventory cost for dual supplier procurement includes complex procedures. To reduce the complexity effectively, the presented method investigates the statistics properties in real-time of the inventory quantities with the application of the infinitesimal dividing method. It is proved that the optimal holding and shortage costs of dual supplier procurement are less than those of single supplier procurement respectively. With the assumption that both suppliers have the same distribution of lead times, the convexity of the cost function per unit time is proved. So the optimal solution can be easily obtained by applying the classical convex optimization methods. The numerical examples are given to verify the main conclusions.

Hybrid LEAP modeling method for long-term energy demand forecasting of regions with limited statistical data

An accurate long-term energy demand forecasting is essential for energy planning and policy making. However, due to the immature energy data collecting and statistical methods, the available data are usually limited in many regions. In this paper, on the basis of comprehensive literature review, we proposed a hybrid model based on the long-range alternative energy planning (LEAP) model to improve the accuracy of energy demand forecasting in these regions. By taking Hunan province, China as a typical case, the proposed hybrid model was applied to estimating the possible future energy demand and energy-saving potentials in different sectors. The structure of LEAP model was estimated by Sankey energy flow, and Leslie matrix and autoregressive integrated moving average (ARIMA) models were used to predict the population, industrial structure and transportation turnover, respectively. Monte-Carlo method was employed to evaluate the uncertainty of forecasted results. The results showed that the hybrid model combined with scenario analysis provided a relatively accurate forecast for the long-term energy demand in regions with limited statistical data, and the average standard error of probabilistic distribution in 2030 energy demand was as low as 0.15. The prediction results could provide supportive references to identify energy-saving potentials and energy development pathways.

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.

Optimization of a precise integration method for seismic modeling based on graphic processing unit

General purpose graphic processing unit （GPU） calculation technology is gradually widely used in various fields. Its mode of single instruction, multiple threads is capable of seismic numerical simulation which has a huge quantity of data and calculation steps. In this study, we introduce a GPU-based parallel calculation method of a precise integration method （PIM） for seismic forward modeling. Compared with CPU single-core calculation, GPU parallel calculating perfectly keeps the features of PIM, which has small bandwidth, high accuracy and capability of modeling complex substructures, and GPU calculation brings high computational efficiency, which means that high-performing GPU parallel calculation can make seismic forward modeling closer to real seismic records.

A novel method of rapidly modeling optical properties of actual photonic crystal fibres

The flexible structure of photonic crystal fibre not only offers novel optical properties but also brings some difficulties in keeping the fibre structure in the fabrication process which inevitably cause the optical properties of the resulting fibre to deviate from the designed properties. Therefore, a method of evaluating the optical properties of the actual fibre is necessary for the purpose of application. Up to now, the methods employed to measure the properties of the actual photonic crystal fibre often require long fibre samples or complex expensive equipments. To our knowledge, there are few studies of modeling an actual photonic crystal fibre and evaluating its properties rapidly. In this paper, a novel method, based on the combination model of digital image processing and the finite element method, is proposed to rapidly model the optical properties of the actual photonic crystal fibre. Two kinds of photonic crystal fibres made by Crystal Fiber A/S are modeled. It is confirmed from numerical results that the proposed method is simple, rapid and accurate for evaluating the optical properties of the actual photonic crystal fibre without requiring complex equipment.

Noise-induced phase transition in the Vicsek model through eigen microstate methodology

This paper presents a comprehensive framework for analyzing phase transitions in collective models such as theVicsek model under various noise types. The Vicsek model, focusing on understanding the collective behaviors of socialanimals, is known due to its discontinuous phase transitions under vector noise. However, its behavior under scalar noiseremains less conclusive. Renowned for its efficacy in the analysis of complex systems under both equilibrium and nonequilibriumstates, the eigen microstate method is employed here for a quantitative examination of the phase transitions inthe Vicsek model under both vector and scalar noises. The study finds that the Vicsek model exhibits discontinuous phasetransitions regardless of noise type. Furthermore, the dichotomy method is utilized to identify the critical points for thesephase transitions. A significant finding is the observed increase in the critical point for discontinuous phase transitions withescalation of population density.