The equivalent source(ES)method in the spherical coordinate system has been widely applied to processing,reduction,field modeling,and geophysical and geological interpretation of satellite magnetic anomaly data.Howeve...The equivalent source(ES)method in the spherical coordinate system has been widely applied to processing,reduction,field modeling,and geophysical and geological interpretation of satellite magnetic anomaly data.However,the inversion for the ES model suffers from nonuniqueness and instability,which remain unresolved.To mitigate these issues,we introduce both the minimum and flattest models into the model objective function as an alternative regularization approach in the spherical ES method.We first present the methods,then analyze the accuracy of forward calculation and test the proposed ES method in this study by using synthetic data.The experimental results from simulation data indicate that our proposed regularization effectively suppresses the Backus effect and mitigates inversion instability in the low-latitude region.Finally,we apply the proposed method to magnetic anomaly data from China Seismo-Electromagnetic Satellite-1(CSES-1)and Macao Science Satellite-1(MSS-1)magnetic measurements over Africa by constructing an ES model of the large-scale lithospheric magnetic field.Compared with existing global lithospheric magnetic field models,our ES model demonstrates good consistency at high altitudes and predicts more stable fields at low altitudes.Furthermore,we derive the reduction to the pole(RTP)magnetic anomaly fields and the apparent susceptibility contrast distribution based on the ES model.The latter correlates well with the regional tectonic framework in Africa and surroundings.展开更多
This study introduces a comprehensive theoretical framework for accurately calculating the electronic band-structure of strained long-wavelength InAs/GaSb type-Ⅱsuperlattices.Utilizing an eight-band k·p Hamilto⁃...This study introduces a comprehensive theoretical framework for accurately calculating the electronic band-structure of strained long-wavelength InAs/GaSb type-Ⅱsuperlattices.Utilizing an eight-band k·p Hamilto⁃nian in conjunction with a scattering matrix method,the model effectively incorporates quantum confinement,strain effects,and interface states.This robust and numerically stable approach achieves exceptional agreement with experimental data,offering a reliable tool for analyzing and engineering the band structure of complex multi⁃layer systems.展开更多
Software systems are vulnerable to security breaches as they expand in complexity and functionality.The confidentiality,integrity,and availability of data are gravely threatened by flaws in a system’s design,implemen...Software systems are vulnerable to security breaches as they expand in complexity and functionality.The confidentiality,integrity,and availability of data are gravely threatened by flaws in a system’s design,implementation,or configuration.To guarantee the durability&robustness of the software,vulnerability identification and fixation have become crucial areas of focus for developers,cybersecurity experts and industries.This paper presents a thorough multi-phase mathematical model for efficient patch management and vulnerability detection.To uniquely model these processes,the model incorporated the notion of the learning phenomenon in describing vulnerability fixation using a logistic learning function.Furthermore,the authors have used numerical methods to approximate the solution of the proposed framework where an analytical solution is difficult to attain.The suggested systematic architecture has been demonstrated through statistical analysis using patch datasets,which offers a solid basis for the research conclusions.According to computational research,learning dynamics improves security response and results in more effective vulnerability management.The suggested model offers a systematic approach to proactive vulnerability mitigation and has important uses in risk assessment,software maintenance,and cybersecurity.This study helps create more robust software systems by increasing patch management effectiveness,which benefits developers,cybersecurity experts,and sectors looking to reduce security threats in a growing digital world.展开更多
We propose a novel workflow for fast forward modeling of well logs in axially symmetric 2D models of the nearwellbore environment.The approach integrates the finite element method with deep residual neural networks to...We propose a novel workflow for fast forward modeling of well logs in axially symmetric 2D models of the nearwellbore environment.The approach integrates the finite element method with deep residual neural networks to achieve exceptional computational efficiency and accuracy.The workflow is demonstrated through the modeling of wireline electromagnetic propagation resistivity logs,where the measured responses exhibit a highly nonlinear relationship with formation properties.The motivation for this research is the need for advanced modeling al-gorithms that are fast enough for use in modern quantitative interpretation tools,where thousands of simulations may be required in iterative inversion processes.The proposed algorithm achieves a remarkable enhancement in performance,being up to 3000 times faster than the finite element method alone when utilizing a GPU.While still ensuring high accuracy,this makes it well-suited for practical applications when reliable payzone assessment is needed in complex environmental scenarios.Furthermore,the algorithm’s efficiency positions it as a promising tool for stochastic Bayesian inversion,facilitating reliable uncertainty quantification in subsurface property estimation.展开更多
In this study,a powerful thermo-hydro-mechanical(THM)coupling solution scheme for saturated poroelastic media involving brittle fracturing is developed.Under the local thermal non-equilibrium(LTNE)assumption,this sche...In this study,a powerful thermo-hydro-mechanical(THM)coupling solution scheme for saturated poroelastic media involving brittle fracturing is developed.Under the local thermal non-equilibrium(LTNE)assumption,this scheme seamlessly combines the material point method(MPM)for accurately tracking solid-phase deformation and heat transport,and the Eulerian finite element method(FEM)for effectively capturing fluid flow and heat advection-diffusion behavior.The proposed approach circumvents the substantial challenges posed by large nonlinear equation systems with the monolithic solution scheme.The staggered solution process strategically separates each physical field through explicit or implicit integration.The characteristic-based method is used to stabilize advection-dominated heat flows for efficient numerical implementation.Furthermore,a fractional step approach is employed to decompose fluid velocity and pressure,thereby suppressing pore pressure oscillation on the linear background grid.The fracturing initiation and propagation are simulated by a rate-dependent phase field model.Through a series of quasi-static and transient simulations,the exceptional performance and promising potential of the proposed model in addressing THM fracturing problems in poro-elastic media is demonstrated.展开更多
In this paper,the multi-agent model about shop logistics is set up. This model has 8 agents: raw materials stock agent,process agent,testing agent,transition agent,production information agent,scheduling agent,process...In this paper,the multi-agent model about shop logistics is set up. This model has 8 agents: raw materials stock agent,process agent,testing agent,transition agent,production information agent,scheduling agent,process agent and stock agent. The scheduling agent has three subagents: manager agent (MA),resource agent (RA) and part agent (PA). MA,PA and RA are communicating equally that guarantees agility of the whole MAS system. The part tasks pass between MA,RA and PA as an integer,which can guarantee the consistency of the data. We use a detailed example about shop logistics scheduling in a semiconductor company to explain the principle. In this example,we use two scheduling strategies: FCFS and SPT. The result data indicates that the average flow time and lingering ratio are changed using different strategy. It is proves that the multi-agent scheduling is useful.展开更多
An intelligent manufacturing system is modeled currently from the viewpoint of manufacturing applications,and the network platform’s influence to manufacturing applications is not considered adequately. However any b...An intelligent manufacturing system is modeled currently from the viewpoint of manufacturing applications,and the network platform’s influence to manufacturing applications is not considered adequately. However any bottleneck in service oriented architecture (SOA) for the manufacturing network can affect the agility of the IT environment. In this paper,to achieve a trade-off between manufacturing resources and network resources,the manufacturing network is modeled with multi-agent,in which two kinds of basic elements,the manufacturing application unit and the network carrier of manufacturing information,are presented. And their main characters are described by colored petri net. The manufacturing application model drives the network platform that inversely provides this application model technology supports. The proposed multi-agent system is demonstrated through an example integration scenario involving production plan,resources management and execution subsystems. And the result suggests that analyzing and designing the system architecture of networked manufacturing should give due attention to the operation system as well as manufacturing applications.展开更多
Dynamic architecture of multi-agent systems （MAS） is important for critical systems. As the existing formal specifications of MAS cannot describe its dynamic architecture, a formal approach using n-calculus is prese...Dynamic architecture of multi-agent systems （MAS） is important for critical systems. As the existing formal specifications of MAS cannot describe its dynamic architecture, a formal approach using n-calculus is presented, which is suited for the describing and analyzing of concurrent MAS with dynamic topology, n-calculus describes the belief-desireintention （BDI） model that represents agent＇s mental states and provides many useful facilities to analyze MAS model such as deadlock, behavior equivalence, and model checking. To illustrate the favorable representation capability of n-calculus, an example of dynamic multi-agent systems in e-commerce is provided. Finally, by using an existing n-calculus supporting tool, MAS model and some key behaviors properties are analyzed and verified.展开更多
In a multi-agent system, each agent must adapt itself to the environment and coordinate with other agents dynamically. TO predict or cooperate with the behavior of oiller agents. An agent should dynamically establish ...In a multi-agent system, each agent must adapt itself to the environment and coordinate with other agents dynamically. TO predict or cooperate with the behavior of oiller agents. An agent should dynamically establish and evolve the cooperative behavior model of itself. In this paper, we represent the behavior model of an agent as a f-mite state machine and propose a new method of dynamically evolving the behavior model of an agent by evolutionary programming.展开更多
With the aid of multi-agent based modeling approach to complex systems, the hierarchy simulation models of carrier-based aircraft catapult launch are developed. Ocean, carrier, aircraft, and atmosphere are treated as ...With the aid of multi-agent based modeling approach to complex systems, the hierarchy simulation models of carrier-based aircraft catapult launch are developed. Ocean, carrier, aircraft, and atmosphere are treated as aggregation agents, the detailed components like catapult, landing gears, and disturbances are considered as meta-agents, which belong to their aggregation agent. Thus, the model with two layers is formed i.e. the aggregation agent layer and the meta-agent layer. The information communication among all agents is described. The meta-agents within one aggregation agent communicate with each other directly by information sharing, but the meta-agents, which belong to different aggregation agents exchange their information through the aggregation layer first, and then perceive it from the sharing environment, that is the aggregation agent. Thus, not only the hierarchy model is built, but also the environment perceived by each agent is specified. Meanwhile, the problem of balancing the independency of agent and the resource consumption brought by real-time communication within multi-agent system （MAS） is resolved. Each agent involved in carrier-based aircraft catapult launch is depicted, with considering the interaction within disturbed atmospheric environment and multiple motion bodies including carrier, aircraft, and landing gears. The models of reactive agents among them are derived based on tensors, and the perceived messages and inner frameworks of each agent are characterized. Finally, some results of a simulation instance are given. The simulation and modeling of dynamic system based on multi-agent system is of benefit to express physical concepts and logical hierarchy clearly and precisely. The system model can easily draw in kinds of other agents to achieve a precise simulation of more complex system. This modeling technique makes the complex integral dynamic equations of multibodies decompose into parallel operations of single agent, and it is convenient to expand, maintain, and reuse the program codes.展开更多
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 secondar...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.展开更多
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 c...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.展开更多
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 gr...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.展开更多
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 mo...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.展开更多
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 distri...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.展开更多
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 th...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.展开更多
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 eleme...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 plays an increasingly important role in Petroleum Geology, Mining Geology and Engineering Geology. The complexity of geological conditions requires different modeling methods in different situa...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.展开更多
Surrogate-Based Optimization(SBO) is becoming increasingly popular since it can remarkably reduce the computational cost for design optimizations based on high-fidelity and expensive numerical analyses. However, for c...Surrogate-Based Optimization(SBO) is becoming increasingly popular since it can remarkably reduce the computational cost for design optimizations based on high-fidelity and expensive numerical analyses. However, for complicated optimization problems with a large design space, many design variables, and strong nonlinearity, SBO converges slowly and shows imperfection in local exploitation. This paper proposes a trust region method within the framework of an SBO process based on the Kriging model. In each refinement cycle, new samples are selected by a certain design of experiment method within a variable design space, which is sequentially updated by the trust region method. A multi-dimensional trust-region radius is proposed to improve the adaptability of the developed methodology. Further, the scale factor and the limit factor of the trust region are studied to evaluate their effects on the optimization process. Thereafter, different SBO methods using error-based exploration, prediction-based exploitation, refinement based on the expected improvement function, a hybrid refinement strategy, and the developed trust-regionbased refinement are utilized in four analytical tests. Further, the developed optimization methodology is employed in the drag minimization of an RAE2822 airfoil. Results indicate that it has better robustness and local exploitation capability in comparison with those of other SBO展开更多
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 equ...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.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.42250103 and 42174090)the Opening Fund of Key Laboratory of Geological Survey and Evaluation of Ministry of Education(Grant No.GLAB2023ZR02)the MOST Special Fund from the State Key Laboratory of Geological Processes and Mineral Resources(Grant No.MSFGPMR2022-4).
文摘The equivalent source(ES)method in the spherical coordinate system has been widely applied to processing,reduction,field modeling,and geophysical and geological interpretation of satellite magnetic anomaly data.However,the inversion for the ES model suffers from nonuniqueness and instability,which remain unresolved.To mitigate these issues,we introduce both the minimum and flattest models into the model objective function as an alternative regularization approach in the spherical ES method.We first present the methods,then analyze the accuracy of forward calculation and test the proposed ES method in this study by using synthetic data.The experimental results from simulation data indicate that our proposed regularization effectively suppresses the Backus effect and mitigates inversion instability in the low-latitude region.Finally,we apply the proposed method to magnetic anomaly data from China Seismo-Electromagnetic Satellite-1(CSES-1)and Macao Science Satellite-1(MSS-1)magnetic measurements over Africa by constructing an ES model of the large-scale lithospheric magnetic field.Compared with existing global lithospheric magnetic field models,our ES model demonstrates good consistency at high altitudes and predicts more stable fields at low altitudes.Furthermore,we derive the reduction to the pole(RTP)magnetic anomaly fields and the apparent susceptibility contrast distribution based on the ES model.The latter correlates well with the regional tectonic framework in Africa and surroundings.
文摘This study introduces a comprehensive theoretical framework for accurately calculating the electronic band-structure of strained long-wavelength InAs/GaSb type-Ⅱsuperlattices.Utilizing an eight-band k·p Hamilto⁃nian in conjunction with a scattering matrix method,the model effectively incorporates quantum confinement,strain effects,and interface states.This robust and numerically stable approach achieves exceptional agreement with experimental data,offering a reliable tool for analyzing and engineering the band structure of complex multi⁃layer systems.
基金supported by grants received by the first author and third author from the Institute of Eminence,Delhi University,Delhi,India,as part of the Faculty Research Program via Ref.No./IoE/2024-25/12/FRP.
文摘Software systems are vulnerable to security breaches as they expand in complexity and functionality.The confidentiality,integrity,and availability of data are gravely threatened by flaws in a system’s design,implementation,or configuration.To guarantee the durability&robustness of the software,vulnerability identification and fixation have become crucial areas of focus for developers,cybersecurity experts and industries.This paper presents a thorough multi-phase mathematical model for efficient patch management and vulnerability detection.To uniquely model these processes,the model incorporated the notion of the learning phenomenon in describing vulnerability fixation using a logistic learning function.Furthermore,the authors have used numerical methods to approximate the solution of the proposed framework where an analytical solution is difficult to attain.The suggested systematic architecture has been demonstrated through statistical analysis using patch datasets,which offers a solid basis for the research conclusions.According to computational research,learning dynamics improves security response and results in more effective vulnerability management.The suggested model offers a systematic approach to proactive vulnerability mitigation and has important uses in risk assessment,software maintenance,and cybersecurity.This study helps create more robust software systems by increasing patch management effectiveness,which benefits developers,cybersecurity experts,and sectors looking to reduce security threats in a growing digital world.
基金financially supported by the Russian federal research project No.FWZZ-2022-0026“Innovative aspects of electro-dynamics in problems of exploration and oilfield geophysics”.
文摘We propose a novel workflow for fast forward modeling of well logs in axially symmetric 2D models of the nearwellbore environment.The approach integrates the finite element method with deep residual neural networks to achieve exceptional computational efficiency and accuracy.The workflow is demonstrated through the modeling of wireline electromagnetic propagation resistivity logs,where the measured responses exhibit a highly nonlinear relationship with formation properties.The motivation for this research is the need for advanced modeling al-gorithms that are fast enough for use in modern quantitative interpretation tools,where thousands of simulations may be required in iterative inversion processes.The proposed algorithm achieves a remarkable enhancement in performance,being up to 3000 times faster than the finite element method alone when utilizing a GPU.While still ensuring high accuracy,this makes it well-suited for practical applications when reliable payzone assessment is needed in complex environmental scenarios.Furthermore,the algorithm’s efficiency positions it as a promising tool for stochastic Bayesian inversion,facilitating reliable uncertainty quantification in subsurface property estimation.
基金supported by National Natural Science Foundation of China(Grant No.42377149)the Research Grants Council of Hong Kong(General Research Fund Project No.17202423).
文摘In this study,a powerful thermo-hydro-mechanical(THM)coupling solution scheme for saturated poroelastic media involving brittle fracturing is developed.Under the local thermal non-equilibrium(LTNE)assumption,this scheme seamlessly combines the material point method(MPM)for accurately tracking solid-phase deformation and heat transport,and the Eulerian finite element method(FEM)for effectively capturing fluid flow and heat advection-diffusion behavior.The proposed approach circumvents the substantial challenges posed by large nonlinear equation systems with the monolithic solution scheme.The staggered solution process strategically separates each physical field through explicit or implicit integration.The characteristic-based method is used to stabilize advection-dominated heat flows for efficient numerical implementation.Furthermore,a fractional step approach is employed to decompose fluid velocity and pressure,thereby suppressing pore pressure oscillation on the linear background grid.The fracturing initiation and propagation are simulated by a rate-dependent phase field model.Through a series of quasi-static and transient simulations,the exceptional performance and promising potential of the proposed model in addressing THM fracturing problems in poro-elastic media is demonstrated.
基金Supported by the Zhejiang Province Science Foundation of China( M703022)
文摘In this paper,the multi-agent model about shop logistics is set up. This model has 8 agents: raw materials stock agent,process agent,testing agent,transition agent,production information agent,scheduling agent,process agent and stock agent. The scheduling agent has three subagents: manager agent (MA),resource agent (RA) and part agent (PA). MA,PA and RA are communicating equally that guarantees agility of the whole MAS system. The part tasks pass between MA,RA and PA as an integer,which can guarantee the consistency of the data. We use a detailed example about shop logistics scheduling in a semiconductor company to explain the principle. In this example,we use two scheduling strategies: FCFS and SPT. The result data indicates that the average flow time and lingering ratio are changed using different strategy. It is proves that the multi-agent scheduling is useful.
文摘An intelligent manufacturing system is modeled currently from the viewpoint of manufacturing applications,and the network platform’s influence to manufacturing applications is not considered adequately. However any bottleneck in service oriented architecture (SOA) for the manufacturing network can affect the agility of the IT environment. In this paper,to achieve a trade-off between manufacturing resources and network resources,the manufacturing network is modeled with multi-agent,in which two kinds of basic elements,the manufacturing application unit and the network carrier of manufacturing information,are presented. And their main characters are described by colored petri net. The manufacturing application model drives the network platform that inversely provides this application model technology supports. The proposed multi-agent system is demonstrated through an example integration scenario involving production plan,resources management and execution subsystems. And the result suggests that analyzing and designing the system architecture of networked manufacturing should give due attention to the operation system as well as manufacturing applications.
基金Project supported by the National High-Technology Research and Development Program of China(Grant No.8632003AA721070)
文摘Dynamic architecture of multi-agent systems （MAS） is important for critical systems. As the existing formal specifications of MAS cannot describe its dynamic architecture, a formal approach using n-calculus is presented, which is suited for the describing and analyzing of concurrent MAS with dynamic topology, n-calculus describes the belief-desireintention （BDI） model that represents agent＇s mental states and provides many useful facilities to analyze MAS model such as deadlock, behavior equivalence, and model checking. To illustrate the favorable representation capability of n-calculus, an example of dynamic multi-agent systems in e-commerce is provided. Finally, by using an existing n-calculus supporting tool, MAS model and some key behaviors properties are analyzed and verified.
文摘In a multi-agent system, each agent must adapt itself to the environment and coordinate with other agents dynamically. TO predict or cooperate with the behavior of oiller agents. An agent should dynamically establish and evolve the cooperative behavior model of itself. In this paper, we represent the behavior model of an agent as a f-mite state machine and propose a new method of dynamically evolving the behavior model of an agent by evolutionary programming.
基金Aeronautical Science Foundation of China (2006ZA51004)
文摘With the aid of multi-agent based modeling approach to complex systems, the hierarchy simulation models of carrier-based aircraft catapult launch are developed. Ocean, carrier, aircraft, and atmosphere are treated as aggregation agents, the detailed components like catapult, landing gears, and disturbances are considered as meta-agents, which belong to their aggregation agent. Thus, the model with two layers is formed i.e. the aggregation agent layer and the meta-agent layer. The information communication among all agents is described. The meta-agents within one aggregation agent communicate with each other directly by information sharing, but the meta-agents, which belong to different aggregation agents exchange their information through the aggregation layer first, and then perceive it from the sharing environment, that is the aggregation agent. Thus, not only the hierarchy model is built, but also the environment perceived by each agent is specified. Meanwhile, the problem of balancing the independency of agent and the resource consumption brought by real-time communication within multi-agent system （MAS） is resolved. Each agent involved in carrier-based aircraft catapult launch is depicted, with considering the interaction within disturbed atmospheric environment and multiple motion bodies including carrier, aircraft, and landing gears. The models of reactive agents among them are derived based on tensors, and the perceived messages and inner frameworks of each agent are characterized. Finally, some results of a simulation instance are given. The simulation and modeling of dynamic system based on multi-agent system is of benefit to express physical concepts and logical hierarchy clearly and precisely. The system model can easily draw in kinds of other agents to achieve a precise simulation of more complex system. This modeling technique makes the complex integral dynamic equations of multibodies decompose into parallel operations of single agent, and it is convenient to expand, maintain, and reuse the program codes.
基金supported by the Natural Science Foundation of China(Nos.41404057,41674077 and 411640034)the Nuclear Energy Development Project of China,and the‘555’Project of Gan Po Excellent People
文摘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.
基金Project(51475156)supported by the National Natural Science Foundation of ChinaProject(2014ZX04002071)supported by the National Key Project of Science and Technology of ChinaProject(GXKFJ14-08)supported by the Opening Foundation of Key Laboratory for Non-Ferrous Metal and Featured Material Processing,Guangxi Zhuang Autonomous Region,China
文摘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.
基金Projects(51161011,11364024)supported by the National Natural Science Foundation of ChinaProject(1204GKCA065)supported by the Key Technology R&D Program of Gansu Province,China+1 种基金Project(201210)supported by the Fundamental Research Funds for the Universities of Gansu Province,ChinaProject(J201304)supported by the Funds for Distinguished Young Scientists of Lanzhou University of Technology,China
文摘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.
基金financially supported by the National Hi-tech Research and Development Program of China(863 Program)(No.2012AA09A20103)
文摘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.
文摘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.
文摘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.
基金supported by the National Nature Science Foundation of China(Grant No.40874055)the Natural Science Foundation of the Hunan Province,China(Grant No.14JJ2012)
文摘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 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.
基金co-supported by the National Natural Science Foundation of China (No. 11502209)the Free Research Projects of the Central University Funding of China (No. 3102015ZY007)
文摘Surrogate-Based Optimization(SBO) is becoming increasingly popular since it can remarkably reduce the computational cost for design optimizations based on high-fidelity and expensive numerical analyses. However, for complicated optimization problems with a large design space, many design variables, and strong nonlinearity, SBO converges slowly and shows imperfection in local exploitation. This paper proposes a trust region method within the framework of an SBO process based on the Kriging model. In each refinement cycle, new samples are selected by a certain design of experiment method within a variable design space, which is sequentially updated by the trust region method. A multi-dimensional trust-region radius is proposed to improve the adaptability of the developed methodology. Further, the scale factor and the limit factor of the trust region are studied to evaluate their effects on the optimization process. Thereafter, different SBO methods using error-based exploration, prediction-based exploitation, refinement based on the expected improvement function, a hybrid refinement strategy, and the developed trust-regionbased refinement are utilized in four analytical tests. Further, the developed optimization methodology is employed in the drag minimization of an RAE2822 airfoil. Results indicate that it has better robustness and local exploitation capability in comparison with those of other SBO
基金Project supported by the National Natural Science Foundation of China(Nos.91752118,11672305,11232011,and 11572331)the Strategic Priority Research Program(No.XDB22040104)the Key Research Program of Frontier Sciences of the Chinese Academy of Sciences(No.QYZDJ-SSWSYS002)
文摘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.
