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.展开更多
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.展开更多
Fluid-structure interaction (FSI) problems caused by fluid impact loads are com- monly existent in naval architectures and ocean engineering fields. For instance, the impact loads due to non-linear fluid motion in a l...Fluid-structure interaction (FSI) problems caused by fluid impact loads are com- monly existent in naval architectures and ocean engineering fields. For instance, the impact loads due to non-linear fluid motion in a liquid sloshing tank potentially affect the structural safety of cargo tanks or vessels. The challenges of numerical study on FSI problems involve not only multidisciplinary features, but also accurate description of non-linear free surface. A fully Lagrangian particle-based method , the moving particle semi-implicit and finite element coupled method ( MPS-FEM), is developed to numerically study the FSI problems. Taking into account the advantage of the Lagrangian method for large deformations of both fluid and solid bounda- ties, the MPS method is used to simulate the fluid field while the finite element method(FEM) to calculate the structure field. Besides, the partitioning strategy is employed to couple the MPS and FEM modules. To validate accuracy of the proposed algorithm, a benchmark case is numer- ically investigated. Both the patterns of free surface and the deflections of the elastic structures are in good agreement with the experimental data. Then, the present FSI solver is applied to the comparative study of the mitigating effects of rigid baffles and elastic baffles on the sloshing motions and impact loads.展开更多
In order to establish an adaptive turbo-shaft engine model with high accuracy, a new modeling method based on parameter selection （PS） algorithm and multi-input multi-output recursive reduced least square support ve...In order to establish an adaptive turbo-shaft engine model with high accuracy, a new modeling method based on parameter selection （PS） algorithm and multi-input multi-output recursive reduced least square support vector regression （MRR-LSSVR） machine is proposed. Firstly, the PS algorithm is designed to choose the most reasonable inputs of the adaptive module. During this process, a wrapper criterion based on least square support vector regression （LSSVR） machine is adopted, which can not only reduce computational complexity but also enhance generalization performance. Secondly, with the input variables determined by the PS algorithm, a mapping model of engine parameter estimation is trained off-line using MRR-LSSVR, which has a satisfying accuracy within 5＆. Finally, based on a numerical simulation platform of an integrated helicopter/ turbo-shaft engine system, an adaptive turbo-shaft engine model is developed and tested in a certain flight envelope. Under the condition of single or multiple engine components being degraded, many simulation experiments are carried out, and the simulation results show the effectiveness and validity of the proposed adaptive modeling method.展开更多
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 energ...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.展开更多
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 c...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.展开更多
In this paper, several widely applied fracture criteria were first numerically examined and the crack-tip-region Jntegral criterion was confirmed to be more applicable to predict fracture angle in an elastic-plastic m...In this paper, several widely applied fracture criteria were first numerically examined and the crack-tip-region Jntegral criterion was confirmed to be more applicable to predict fracture angle in an elastic-plastic multiphase material. Then, the crack propagation in an idealized an elastic-plastic finite element method. The variation dendritic two-phase AI-7%Si alloy was modeled using of crack growth driving force with crack extension was also demonstrated. It is found that the crack path is significantly influenced by the presence of α-phase near the crack tip, and the crack growth driving force varies drastically from place to place. Lastly, the simulated fracture path in the two-phase model alloy was compared with the experimentally observed fracture path.展开更多
基金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 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.
文摘Fluid-structure interaction (FSI) problems caused by fluid impact loads are com- monly existent in naval architectures and ocean engineering fields. For instance, the impact loads due to non-linear fluid motion in a liquid sloshing tank potentially affect the structural safety of cargo tanks or vessels. The challenges of numerical study on FSI problems involve not only multidisciplinary features, but also accurate description of non-linear free surface. A fully Lagrangian particle-based method , the moving particle semi-implicit and finite element coupled method ( MPS-FEM), is developed to numerically study the FSI problems. Taking into account the advantage of the Lagrangian method for large deformations of both fluid and solid bounda- ties, the MPS method is used to simulate the fluid field while the finite element method(FEM) to calculate the structure field. Besides, the partitioning strategy is employed to couple the MPS and FEM modules. To validate accuracy of the proposed algorithm, a benchmark case is numer- ically investigated. Both the patterns of free surface and the deflections of the elastic structures are in good agreement with the experimental data. Then, the present FSI solver is applied to the comparative study of the mitigating effects of rigid baffles and elastic baffles on the sloshing motions and impact loads.
基金co-supported by Aeronautical Science Foundation of China (No. 2010ZB52011)Funding of Jiangsu Innovation Program for Graduate Education (No.CXLX11_0213)
文摘In order to establish an adaptive turbo-shaft engine model with high accuracy, a new modeling method based on parameter selection （PS） algorithm and multi-input multi-output recursive reduced least square support vector regression （MRR-LSSVR） machine is proposed. Firstly, the PS algorithm is designed to choose the most reasonable inputs of the adaptive module. During this process, a wrapper criterion based on least square support vector regression （LSSVR） machine is adopted, which can not only reduce computational complexity but also enhance generalization performance. Secondly, with the input variables determined by the PS algorithm, a mapping model of engine parameter estimation is trained off-line using MRR-LSSVR, which has a satisfying accuracy within 5＆. Finally, based on a numerical simulation platform of an integrated helicopter/ turbo-shaft engine system, an adaptive turbo-shaft engine model is developed and tested in a certain flight envelope. Under the condition of single or multiple engine components being degraded, many simulation experiments are carried out, and the simulation results show the effectiveness and validity of the proposed adaptive modeling method.
基金Supported by National Natural Science Foundation of China(Grant No.51605058)Chongqing Research Program of Basic Research and Frontier Technology of China(Grant No.cstc2015jcyjBX0088)+2 种基金Fundamental Research Funds for the Central Universities of China(Grant No.106112016CDJCR021226)Six Talent Peaks Project in Jiangsu Province of China(Grant No.2014-ZBZZ-006)"Excellence Plans-Zijin Star" Foundation of Nanjing University of Science and Technology,China(Grant No.2015-zijin-07)
文摘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.
文摘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.
文摘In this paper, several widely applied fracture criteria were first numerically examined and the crack-tip-region Jntegral criterion was confirmed to be more applicable to predict fracture angle in an elastic-plastic multiphase material. Then, the crack propagation in an idealized an elastic-plastic finite element method. The variation dendritic two-phase AI-7%Si alloy was modeled using of crack growth driving force with crack extension was also demonstrated. It is found that the crack path is significantly influenced by the presence of α-phase near the crack tip, and the crack growth driving force varies drastically from place to place. Lastly, the simulated fracture path in the two-phase model alloy was compared with the experimentally observed fracture path.
