Modern medicine is reliant on various medical imaging technologies for non-invasively observing patients’anatomy.However,the interpretation of medical images can be highly subjective and dependent on the expertise of...Modern medicine is reliant on various medical imaging technologies for non-invasively observing patients’anatomy.However,the interpretation of medical images can be highly subjective and dependent on the expertise of clinicians.Moreover,some potentially useful quantitative information in medical images,especially that which is not visible to the naked eye,is often ignored during clinical practice.In contrast,radiomics performs high-throughput feature extraction from medical images,which enables quantitative analysis of medical images and prediction of various clinical endpoints.Studies have reported that radiomics exhibits promising performance in diagnosis and predicting treatment responses and prognosis,demonstrating its potential to be a non-invasive auxiliary tool for personalized medicine.However,radiomics remains in a developmental phase as numerous technical challenges have yet to be solved,especially in feature engineering and statistical modeling.In this review,we introduce the current utility of radiomics by summarizing research on its application in the diagnosis,prognosis,and prediction of treatment responses in patients with cancer.We focus on machine learning approaches,for feature extraction and selection during feature engineering and for imbalanced datasets and multi-modality fusion during statistical modeling.Furthermore,we introduce the stability,reproducibility,and interpretability of features,and the generalizability and interpretability of models.Finally,we offer possible solutions to current challenges in radiomics research.展开更多
We have proposed a methodology to assess the robustness of underground tunnels against potential failure.This involves developing vulnerability functions for various qualities of rock mass and static loading intensiti...We have proposed a methodology to assess the robustness of underground tunnels against potential failure.This involves developing vulnerability functions for various qualities of rock mass and static loading intensities.To account for these variations,we utilized a Monte Carlo Simulation(MCS)technique coupled with the finite difference code FLAC^(3D),to conduct two thousand seven hundred numerical simulations of a horseshoe tunnel located within a rock mass with different geological strength index system(GSIs)and subjected to different states of static loading.To quantify the severity of damage within the rock mass,we selected one stress-based(brittle shear ratio(BSR))and one strain-based failure criterion(plastic damage index(PDI)).Based on these criteria,we then developed fragility curves.Additionally,we used mathematical approximation techniques to produce vulnerability functions that relate the probabilities of various damage states to loading intensities for different quality classes of blocky rock mass.The results indicated that the fragility curves we obtained could accurately depict the evolution of the inner and outer shell damage around the tunnel.Therefore,we have provided engineers with a tool that can predict levels of damages associated with different failure mechanisms based on variations in rock mass quality and in situ stress state.Our method is a numerically developed,multi-variate approach that can aid engineers in making informed decisions about the robustness of underground tunnels.展开更多
Machine learning(ML)provides a new surrogate method for investigating groundwater flow dynamics in unsaturated soils.Traditional pure data-driven methods(e.g.deep neural network,DNN)can provide rapid predictions,but t...Machine learning(ML)provides a new surrogate method for investigating groundwater flow dynamics in unsaturated soils.Traditional pure data-driven methods(e.g.deep neural network,DNN)can provide rapid predictions,but they do require sufficient on-site data for accurate training,and lack interpretability to the physical processes within the data.In this paper,we provide a physics and equalityconstrained artificial neural network(PECANN),to derive unsaturated infiltration solutions with a small amount of initial and boundary data.PECANN takes the physics-informed neural network(PINN)as a foundation,encodes the unsaturated infiltration physical laws(i.e.Richards equation,RE)into the loss function,and uses the augmented Lagrangian method to constrain the learning process of the solutions of RE by adding stronger penalty for the initial and boundary conditions.Four unsaturated infiltration cases are designed to test the training performance of PECANN,i.e.one-dimensional(1D)steady-state unsaturated infiltration,1D transient-state infiltration,two-dimensional(2D)transient-state infiltration,and 1D coupled unsaturated infiltration and deformation.The predicted results of PECANN are compared with the finite difference solutions or analytical solutions.The results indicate that PECANN can accurately capture the variations of pressure head during the unsaturated infiltration,and present higher precision and robustness than DNN and PINN.It is also revealed that PECANN can achieve the same accuracy as the finite difference method with fewer initial and boundary training data.Additionally,we investigate the effect of the hyperparameters of PECANN on solving RE problem.PECANN provides an effective tool for simulating unsaturated infiltration.展开更多
Natural slopes usually display complicated exposed rock surfaces that are characterized by complex and substantial terrain undulation and ubiquitous undesirable phenomena such as vegetation cover and rockfalls.This st...Natural slopes usually display complicated exposed rock surfaces that are characterized by complex and substantial terrain undulation and ubiquitous undesirable phenomena such as vegetation cover and rockfalls.This study presents a systematic outcrop research of fracture pattern variations in a complicated rock slope,and the qualitative and quantitative study of the complex phenomena impact on threedimensional(3D)discrete fracture network(DFN)modeling.As the studies of the outcrop fracture pattern have been so far focused on local variations,thus,we put forward a statistical analysis of global variations.The entire outcrop is partitioned into several subzones,and the subzone-scale variability of fracture geometric properties is analyzed(including the orientation,the density,and the trace length).The results reveal significant variations in fracture characteristics(such as the concentrative degree,the average orientation,the density,and the trace length)among different subzones.Moreover,the density of fracture sets,which is approximately parallel to the slope surface,exhibits a notably higher value compared to other fracture sets across all subzones.To improve the accuracy of the DFN modeling,the effects of three common phenomena resulting from vegetation and rockfalls are qualitatively analyzed and the corresponding quantitative data processing solutions are proposed.Subsequently,the 3D fracture geometric parameters are determined for different areas of the high-steep rock slope in terms of the subzone dimensions.The results show significant variations in the same set of 3D fracture parameters across different regions with density differing by up to tenfold and mean trace length exhibiting differences of 3e4 times.The study results present precise geological structural information,improve modeling accuracy,and provide practical solutions for addressing complex outcrop issues.展开更多
Global images of auroras obtained by cameras on spacecraft are a key tool for studying the near-Earth environment.However,the cameras are sensitive not only to auroral emissions produced by precipitating particles,but...Global images of auroras obtained by cameras on spacecraft are a key tool for studying the near-Earth environment.However,the cameras are sensitive not only to auroral emissions produced by precipitating particles,but also to dayglow emissions produced by photoelectrons induced by sunlight.Nightglow emissions and scattered sunlight can contribute to the background signal.To fully utilize such images in space science,background contamination must be removed to isolate the auroral signal.Here we outline a data-driven approach to modeling the background intensity in multiple images by formulating linear inverse problems based on B-splines and spherical harmonics.The approach is robust,flexible,and iteratively deselects outliers,such as auroral emissions.The final model is smooth across the terminator and accounts for slow temporal variations and large-scale asymmetries in the dayglow.We demonstrate the model by using the three far ultraviolet cameras on the Imager for Magnetopause-to-Aurora Global Exploration(IMAGE)mission.The method can be applied to historical missions and is relevant for upcoming missions,such as the Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)mission.展开更多
Rock fragmentation plays a critical role in rock avalanches,yet conventional approaches such as classical granular flow models or the bonded particle model have limitations in accurately characterizing the progressive...Rock fragmentation plays a critical role in rock avalanches,yet conventional approaches such as classical granular flow models or the bonded particle model have limitations in accurately characterizing the progressive disintegration and kinematics of multi-deformable rock blocks during rockslides.The present study proposes a discrete-continuous numerical model,based on a cohesive zone model,to explicitly incorporate the progressive fragmentation and intricate interparticle interactions inherent in rockslides.Breakable rock granular assemblies are released along an inclined plane and flow onto a horizontal plane.The numerical scenarios are established to incorporate variations in slope angle,initial height,friction coefficient,and particle number.The evolutions of fragmentation,kinematic,runout and depositional characteristics are quantitatively analyzed and compared with experimental and field data.A positive linear relationship between the equivalent friction coefficient and the apparent friction coefficient is identified.In general,the granular mass predominantly exhibits characteristics of a dense granular flow,with the Savage number exhibiting a decreasing trend as the volume of mass increases.The process of particle breakage gradually occurs in a bottom-up manner,leading to a significant increase in the angular velocities of the rock blocks with increasing depth.The simulation results reproduce the field observations of inverse grading and source stratigraphy preservation in the deposit.We propose a disintegration index that incorporates factors such as drop height,rock mass volume,and rock strength.Our findings demonstrate a consistent linear relationship between this index and the fragmentation degree in all tested scenarios.展开更多
Ethylene glycol(EG)plays a pivotal role as a primary raw material in the polyester industry,and the syngas-to-EG route has become a significant technical route in production.The carbon monoxide(CO)gas-phase catalytic ...Ethylene glycol(EG)plays a pivotal role as a primary raw material in the polyester industry,and the syngas-to-EG route has become a significant technical route in production.The carbon monoxide(CO)gas-phase catalytic coupling to synthesize dimethyl oxalate(DMO)is a crucial process in the syngas-to-EG route,whereby the composition of the reactor outlet exerts influence on the ultimate quality of the EG product and the energy consumption during the subsequent separation process.However,measuring product quality in real time or establishing accurate dynamic mechanism models is challenging.To effectively model the DMO synthesis process,this study proposes a hybrid modeling strategy that integrates process mechanisms and data-driven approaches.The CO gas-phase catalytic coupling mechanism model is developed based on intrinsic kinetics and material balance,while a long short-term memory(LSTM)neural network is employed to predict the macroscopic reaction rate by leveraging temporal relationships derived from archived measurements.The proposed model is trained semi-supervised to accommodate limited-label data scenarios,leveraging historical data.By integrating these predictions with the mechanism model,the hybrid modeling approach provides reliable and interpretable forecasts of mass fractions.Empirical investigations unequivocally validate the superiority of the proposed hybrid modeling approach over conventional data-driven models(DDMs)and other hybrid modeling techniques.展开更多
We used field experimental data to evaluate the ability of the agricultural production system model (APSIM) to simulate soil P availability,maize biomass and grain yield in response to P fertilizer applications on a f...We used field experimental data to evaluate the ability of the agricultural production system model (APSIM) to simulate soil P availability,maize biomass and grain yield in response to P fertilizer applications on a fluvo-aquic soil in the North China Plain.Crop and soil data from a 2-year experiment with three P fertilizer application rates(0,75 and 300 kg P_(2)O_(5) ha^(–1)) were used to calibrate the model.Sensitivity analysis was carried out to investigate the influence of APSIM SoilP parameters on the simulated P availability in soil and maize growth.Crop and soil P parameters were then derived by matching or relating the simulation results to observed crop biomass,yield,P uptake and Olsen-P in soil.The re-parameterized model was further validated against 2 years of independent data at the same sites.The re-parameterized model enabled good simulation of the maize leaf area index (LAI),biomass,grain yield,P uptake,and grain P content in response to different levels of P additions against both the calibration and validation datasets.Our results showed that APSIM needs to be re-parameterized for simulation of maize LAI dynamics through modification of leaf size curve and a reduction in the rate of leaf senescence for modern staygreen maize cultivars in China.The P concentration limits (maximum and minimum P concentrations in organs)at different stages also need to be adjusted.Our results further showed a curvilinear relationship between the measured Olsen-P concentration and simulated labile P content,which could facilitate the initialization of APSIM P pools in the NCP with Olsen-P measurements in future studies.It remains difficult to parameterize the APSIM SoilP module due to the conceptual nature of the pools and simplified conceptualization of key P transformation processes.A fundamental understanding still needs to be developed for modelling and predicting the fate of applied P fertilizers in soils with contrasting physical and chemical characteristics.展开更多
The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling ...The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling efficiency of underground engineering,a modularized and parametric modeling cloud server is developed by using Python codes.The basic framework of the cloud server is as follows:input the modeling parameters into the web platform,implement Rhino software and FLAC3D software to model and run simulations in the cloud server,and return the simulation results to the web platform.The modeling program can automatically generate instructions that can run the modeling process in Rhino based on the input modeling parameters.The main modules of the modeling program include modeling the 3D geological structures,the underground engineering structures,and the supporting structures as well as meshing the geometric models.In particular,various cross-sections of underground caverns are crafted as parametricmodules in themodeling program.Themodularized and parametric modeling program is used for a finite element simulation of the underground powerhouse of the Shuangjiangkou Hydropower Station.This complicatedmodel is rapidly generated for the simulation,and the simulation results are reasonable.Thus,this modularized and parametric modeling program is applicable for three-dimensional finite element simulations and analyses.展开更多
We apply stochastic seismic inversion and Bayesian facies classification for porosity modeling and igneous rock identification in the presalt interval of the Santos Basin. This integration of seismic and well-derived ...We apply stochastic seismic inversion and Bayesian facies classification for porosity modeling and igneous rock identification in the presalt interval of the Santos Basin. This integration of seismic and well-derived information enhances reservoir characterization. Stochastic inversion and Bayesian classification are powerful tools because they permit addressing the uncertainties in the model. We used the ES-MDA algorithm to achieve the realizations equivalent to the percentiles P10, P50, and P90 of acoustic impedance, a novel method for acoustic inversion in presalt. The facies were divided into five: reservoir 1,reservoir 2, tight carbonates, clayey rocks, and igneous rocks. To deal with the overlaps in acoustic impedance values of facies, we included geological information using a priori probability, indicating that structural highs are reservoir-dominated. To illustrate our approach, we conducted porosity modeling using facies-related rock-physics models for rock-physics inversion in an area with a well drilled in a coquina bank and evaluated the thickness and extension of an igneous intrusion near the carbonate-salt interface. The modeled porosity and the classified seismic facies are in good agreement with the ones observed in the wells. Notably, the coquinas bank presents an improvement in the porosity towards the top. The a priori probability model was crucial for limiting the clayey rocks to the structural lows. In Well B, the hit rate of the igneous rock in the three scenarios is higher than 60%, showing an excellent thickness-prediction capability.展开更多
The influences of particle size,shape,and catalyst distribution on the reactivity and hydrocarbon product selectivity of a cobalt-based catalyst for Fischer-Tropsch synthesis were investigated in the present work.A se...The influences of particle size,shape,and catalyst distribution on the reactivity and hydrocarbon product selectivity of a cobalt-based catalyst for Fischer-Tropsch synthesis were investigated in the present work.A self-consistent kinetic model for Fischer-Tropsch reaction proposed here was found to correlate experimental data well and hence was used to describe the consumption rates of reactants and formation rates of hydrocarbon products.The perturbed-chain statistical associating fluid theory equation of state was used to describe vapor-liquid equilibrium behavior associated with Fischer-Tropsch reaction.Local interaction between intraparticle diffusion and Fischer-Tropsch reaction was investigated in detail.Results showed that in order to avoid the adverse influence of intraparticle diffusional limitations on catalyst reactivity and product selectivity,the use of small particles is necessary.Large eggshell spherical particles are shown to keep the original catalyst reactivity and enhance the selectivity of heavy hydrocarbon products.The suitable layer thickness for a spherical particle with a diameter of 2 mm is nearly 0.15 mm.With the same outer diameter of 2 mm,the catalyst reactivity and heavy product selectivity of hollow cylindrical particles with a layer thickness of 0.25 mm are found to be larger than eggshell spherical particles.From the viewpoint of catalytic performance,hollow cylindrical particles are a better choice for industrial applications.展开更多
The global populationhas beenandwill continue to be severely impacted by theCOVID-19 epidemic.The primary objective of this research is to demonstrate the future impact of COVID-19 on those who suffer from other fatal...The global populationhas beenandwill continue to be severely impacted by theCOVID-19 epidemic.The primary objective of this research is to demonstrate the future impact of COVID-19 on those who suffer from other fatal conditions such as cancer,heart disease,and diabetes.Here,using ordinary differential equations(ODEs),two mathematical models are developed to explain the association between COVID-19 and cancer and between COVID-19 and diabetes and heart disease.After that,we highlight the stability assessments that can be applied to these models.Sensitivity analysis is used to examine how changes in certain factors impact different aspects of disease.The sensitivity analysis showed that many people are still nervous about seeing a doctor due to COVID-19,which could result in a dramatic increase in the diagnosis of various ailments in the years to come.The correlation between diabetes and cardiovascular illness is also illustrated graphically.The effects of smoking and obesity are also found to be significant in disease compartments.Model fitting is also provided for interpreting the relationship between real data and the results of thiswork.Diabetic people,in particular,need tomonitor their health conditions closely and practice heart health maintenance.People with heart diseases should undergo regular checks so that they can protect themselves from diabetes and take some precautions including suitable diets.The main purpose of this study is to emphasize the importance of regular checks,to warn people about the effects of COVID-19(including avoiding healthcare centers and doctors because of the spread of infectious diseases)and to indicate the importance of family history of cancer,heart diseases and diabetes.The provision of the recommendations requires an increase in public consciousness.展开更多
Visible light communication(VLC)has attracted much attention in the research of sixthgeneration(6G)systems.Furthermore,channel modeling is the foundation for designing efficient and robust VLC systems.In this paper,we...Visible light communication(VLC)has attracted much attention in the research of sixthgeneration(6G)systems.Furthermore,channel modeling is the foundation for designing efficient and robust VLC systems.In this paper,we present extensive VLC channel measurement campaigns in indoor environments,i.e.,an office and a corridor.Based on the measured data,the large-scale fading characteristics and multipath-related characteristics,including omnidirectional optical path loss(OPL),K-factor,power angular spectrum(PAS),angle spread(AS),and clustering characteristics,are analyzed and modeled through a statistical method.Based on the extracted statistics of the above-mentioned channel characteristics,we propose a statistical spatial channel model(SSCM)capable of modeling multipath in the spatial domain.Furthermore,the simulated statistics of the proposed model are compared with the measured statistics.For instance,in the office,the simulated path loss exponent(PLE)and the measured PLE are 1.96and 1.97,respectively.And,the simulated medians of AS and measured medians of AS are 25.94°and 24.84°,respectively.Generally,the fact that the simulated results fit well with measured results has demonstrated the accuracy of our SSCM.展开更多
Simulating the total ionizing dose(TID)of an electrical system using transistor-level models can be difficult and expensive,particularly for digital-integrated circuits(ICs).In this study,a method for modeling TID eff...Simulating the total ionizing dose(TID)of an electrical system using transistor-level models can be difficult and expensive,particularly for digital-integrated circuits(ICs).In this study,a method for modeling TID effects in complementary metaloxide semiconductor(CMOS)digital ICs based on the input/output buffer information specification(IBIS)was proposed.The digital IC was first divided into three parts based on its internal structure:the input buffer,output buffer,and functional area.Each of these three parts was separately modeled.Using the IBIS model,the transistor V-I characteristic curves of the buffers were processed,and the physical parameters were extracted and modeled using VHDL-AMS.In the functional area,logic functions were modeled in VHDL according to the data sheet.A golden digital IC model was developed by combining the input buffer,output buffer,and functional area models.Furthermore,the golden ratio was reconstructed based on TID experimental data,enabling the assessment of TID effects on the threshold voltage,carrier mobility,and time series of the digital IC.TID experiments were conducted using a CMOS non-inverting multiplexer,NC7SZ157,and the results were compared with the simulation results,which showed that the relative errors were less than 2%at each dose point.This confirms the practicality and accuracy of the proposed modeling method.The TID effect model for digital ICs developed using this modeling technique includes both the logical function of the IC and changes in electrical properties and functional degradation impacted by TID,which has potential applications in the design of radiation-hardening tolerance in digital ICs.展开更多
This paper describes numerical simulation of hydraulic fracturing using fracture-based continuum modeling(FBCM)of coupled geomechanical-hydrological processes to evaluate a technique for high-density fracturing and fr...This paper describes numerical simulation of hydraulic fracturing using fracture-based continuum modeling(FBCM)of coupled geomechanical-hydrological processes to evaluate a technique for high-density fracturing and fracture caging.The simulations are innovative because of modeling discrete fractures explicitly in continuum analysis.A key advantage of FBCM is that fracture initiation and propagation are modeled explicitly without changing the domain grid(i.e.no re-meshing).Further,multiple realizations of a preexisting fracture distribution can be analyzed using the same domain grid.The simulated hydraulic fracturing technique consists of pressurizing multiple wells simultaneously:initially without permeating fluids into the rock,to seed fractures uniformly and at high density in the wall rock of the wells;followed by fluid injection to propagate the seeded fracture density hydraulically.FBCM combines the ease of continuum modeling with the potential accuracy of modeling discrete fractures and fracturing explicitly.Fractures are modeled as piecewise planar based on intersections with domain elements;fracture geometry stored as continuum properties is used to calculate parameters needed to model individual fractures;and rock behavior is modeled through tensorial aggregation of the behavior of discrete fractures and unfractured rock.Simulations are presented for previously unfractured rock and for rock with preexisting fractures of horizontal,shallow-dipping,steeply dipping,or vertical orientation.Simulations of a single-well model are used to determine the pattern and spacing for a multiple-well design.The results illustrate high-density fracturing and fracture caging through simultaneous fluid injection in multiple wells:for previously unfractured rock or rock with preexisting shallow-dipping or horizontal fractures,and in situ vertical compressive stress greater than horizontal.If preexisting fractures are steeply dipping or vertical,and considering the same in situ stress condition,well pressurization without fluid permeation appears to be the only practical way to induce new fractures and contain fracturing within the target domain.展开更多
Staggered-grid finite-difference(SGFD)schemes have been widely used in acoustic wave modeling for geophysical problems.Many improved methods are proposed to enhance the accuracy of numerical modeling.However,these met...Staggered-grid finite-difference(SGFD)schemes have been widely used in acoustic wave modeling for geophysical problems.Many improved methods are proposed to enhance the accuracy of numerical modeling.However,these methods are inevitably limited by the maximum Courant-Friedrichs-Lewy(CFL)numbers,making them unstable when modeling with large time sampling intervals or small grid spacings.To solve this problem,we extend a stable SGFD scheme by controlling SGFD dispersion relations and maximizing the maximum CFL numbers.First,to improve modeling stability,we minimize the error between the FD dispersion relation and the exact relation in the given wave-number region,and make the FD dispersion approach a given function outside the given wave-number area,thus breaking the conventional limits of the maximum CFL number.Second,to obtain high modeling accuracy,we use the SGFD scheme based on the Remez algorithm to compute the FD coefficients.In addition,the hybrid absorbing boundary condition is adopted to suppress boundary reflections and we find a suitable weighting coefficient for the proposed scheme.Theoretical derivation and numerical modeling demonstrate that the proposed scheme can maintain high accuracy in the modeling process and the value of the maximum CFL number of the proposed scheme can exceed that of the conventional SGFD scheme when adopting a small maximum effective wavenumber,indicating that the proposed scheme improves stability during the modeling.展开更多
Deterministic compartment models(CMs)and stochastic models,including stochastic CMs and agent-based models,are widely utilized in epidemic modeling.However,the relationship between CMs and their corresponding stochast...Deterministic compartment models(CMs)and stochastic models,including stochastic CMs and agent-based models,are widely utilized in epidemic modeling.However,the relationship between CMs and their corresponding stochastic models is not well understood.The present study aimed to address this gap by conducting a comparative study using the susceptible,exposed,infectious,and recovered(SEIR)model and its extended CMs from the coronavirus disease 2019 modeling literature.We demonstrated the equivalence of the numerical solution of CMs using the Euler scheme and their stochastic counterparts through theoretical analysis and simulations.Based on this equivalence,we proposed an efficient model calibration method that could replicate the exact solution of CMs in the corresponding stochastic models through parameter adjustment.The advancement in calibration techniques enhanced the accuracy of stochastic modeling in capturing the dynamics of epidemics.However,it should be noted that discrete-time stochastic models cannot perfectly reproduce the exact solution of continuous-time CMs.Additionally,we proposed a new stochastic compartment and agent mixed model as an alternative to agent-based models for large-scale population simulations with a limited number of agents.This model offered a balance between computational efficiency and accuracy.The results of this research contributed to the comparison and unification of deterministic CMs and stochastic models in epidemic modeling.Furthermore,the results had implications for the development of hybrid models that integrated the strengths of both frameworks.Overall,the present study has provided valuable epidemic modeling techniques and their practical applications for understanding and controlling the spread of infectious diseases.展开更多
This research aims to optimize the utilization of long-term sea level data from the TOPEX/Poseidon,Jason1,Jason2,and Jason3 altimetry missions for tidal modeling.We generate a time series of along-track observations a...This research aims to optimize the utilization of long-term sea level data from the TOPEX/Poseidon,Jason1,Jason2,and Jason3 altimetry missions for tidal modeling.We generate a time series of along-track observations and apply a developed method to produce tidal models with specific tidal constituents for each location.Our tidal modeling methodology follows an iterative process:partitioning sea surface height(SSH)observations into analysis/training and prediction/validation parts and ultimately identi-fying the set of tidal constituents that provide the best predictions at each time series location.The study focuses on developing 1256 time series along the altimetry tracks over the Baltic Sea,each with its own set of tidal constituents.Verification of the developed tidal models against the sSH observations within the prediction/validation part reveals mean absolute error(MAE)values ranging from 0.0334 m to 0.1349 m,with an average MAE of 0.089 m.The same validation process is conducted on the FES2014 and EOT20 global tidal models,demonstrating that our tidal model,referred to as BT23(short for Baltic Tide 2023),outperforms both models with an average MAE improvement of 0.0417 m and 0.0346 m,respectively.In addition to providing details on the development of the time series and the tidal modeling procedure,we offer the 1256 along-track time series and their associated tidal models as supplementary materials.We encourage the satellite altimetry community to utilize these resources for further research and applications.展开更多
Implanted neural probes can detect weak discharges of neurons in the brain by piercing soft brain tissue,thus as important tools for brain science research,as well as diagnosis and treatment of brain diseases.However,...Implanted neural probes can detect weak discharges of neurons in the brain by piercing soft brain tissue,thus as important tools for brain science research,as well as diagnosis and treatment of brain diseases.However,the rigid neural probes,such as Utah arrays,Michigan probes,and metal microfilament electrodes,are mechanically unmatched with brain tissue and are prone to rejection and glial scarring after implantation,which leads to a significant degradation in the signal quality with the implantation time.In recent years,flexible neural electrodes are rapidly developed with less damage to biological tissues,excellent biocompatibility,and mechanical compliance to alleviate scarring.Among them,the mechanical modeling is important for the optimization of the structure and the implantation process.In this review,the theoretical calculation of the flexible neural probes is firstly summarized with the processes of buckling,insertion,and relative interaction with soft brain tissue for flexible probes from outside to inside.Then,the corresponding mechanical simulation methods are organized considering multiple impact factors to realize minimally invasive implantation.Finally,the technical difficulties and future trends of mechanical modeling are discussed for the next-generation flexible neural probes,which is critical to realize low-invasiveness and long-term coexistence in vivo.展开更多
The human digital twin(HDT)emerges as a promising human-centric technology in Industry 5.0,but challenges remain in human modeling and simulation.Digital human modeling(DHM)provides solutions for modeling and simulati...The human digital twin(HDT)emerges as a promising human-centric technology in Industry 5.0,but challenges remain in human modeling and simulation.Digital human modeling(DHM)provides solutions for modeling and simulating human physical and cognitive aspects to support ergonomic analysis.However,it has limitations in real-time data usage,personalized services,and timely interaction.The emerging HDT concept offers new possibilities by integrating multi-source data and artificial intelligence for continuous monitoring and assessment.Hence,this paper reviews the evolution from DHM to HDT and proposes a unified HDT framework from a human factors perspective.The framework comprises the physical twin,the virtual twin,and the linkage between these two.The virtual twin integrates human modeling and AI engines to enable model-data-hybrid-enabled simulation.HDT can potentially upgrade traditional ergonomic methods to intelligent services through real-time analysis,timely feedback,and bidirectional interactions.Finally,the future perspectives of HDT for industrial applications as well as technical and social challenges are discussed.In general,this study outlines a human factors perspective on HDT for the first time,which is useful for cross-disciplinary research and human factors innovation to enhance the development of HDT in industry.展开更多
基金supported in part by the National Natural Science Foundation of China(82072019)the Shenzhen Basic Research Program(JCYJ20210324130209023)+5 种基金the Shenzhen-Hong Kong-Macao S&T Program(Category C)(SGDX20201103095002019)the Mainland-Hong Kong Joint Funding Scheme(MHKJFS)(MHP/005/20),the Project of Strategic Importance Fund(P0035421)the Projects of RISA(P0043001)from the Hong Kong Polytechnic University,the Natural Science Foundation of Jiangsu Province(BK20201441)the Provincial and Ministry Co-constructed Project of Henan Province Medical Science and Technology Research(SBGJ202103038,SBGJ202102056)the Henan Province Key R&D and Promotion Project(Science and Technology Research)(222102310015)the Natural Science Foundation of Henan Province(222300420575),and the Henan Province Science and Technology Research(222102310322).
文摘Modern medicine is reliant on various medical imaging technologies for non-invasively observing patients’anatomy.However,the interpretation of medical images can be highly subjective and dependent on the expertise of clinicians.Moreover,some potentially useful quantitative information in medical images,especially that which is not visible to the naked eye,is often ignored during clinical practice.In contrast,radiomics performs high-throughput feature extraction from medical images,which enables quantitative analysis of medical images and prediction of various clinical endpoints.Studies have reported that radiomics exhibits promising performance in diagnosis and predicting treatment responses and prognosis,demonstrating its potential to be a non-invasive auxiliary tool for personalized medicine.However,radiomics remains in a developmental phase as numerous technical challenges have yet to be solved,especially in feature engineering and statistical modeling.In this review,we introduce the current utility of radiomics by summarizing research on its application in the diagnosis,prognosis,and prediction of treatment responses in patients with cancer.We focus on machine learning approaches,for feature extraction and selection during feature engineering and for imbalanced datasets and multi-modality fusion during statistical modeling.Furthermore,we introduce the stability,reproducibility,and interpretability of features,and the generalizability and interpretability of models.Finally,we offer possible solutions to current challenges in radiomics research.
基金funding received by a grant from the Natural Sciences and Engineering Research Council of Canada(NSERC)(Grant No.CRDPJ 469057e14).
文摘We have proposed a methodology to assess the robustness of underground tunnels against potential failure.This involves developing vulnerability functions for various qualities of rock mass and static loading intensities.To account for these variations,we utilized a Monte Carlo Simulation(MCS)technique coupled with the finite difference code FLAC^(3D),to conduct two thousand seven hundred numerical simulations of a horseshoe tunnel located within a rock mass with different geological strength index system(GSIs)and subjected to different states of static loading.To quantify the severity of damage within the rock mass,we selected one stress-based(brittle shear ratio(BSR))and one strain-based failure criterion(plastic damage index(PDI)).Based on these criteria,we then developed fragility curves.Additionally,we used mathematical approximation techniques to produce vulnerability functions that relate the probabilities of various damage states to loading intensities for different quality classes of blocky rock mass.The results indicated that the fragility curves we obtained could accurately depict the evolution of the inner and outer shell damage around the tunnel.Therefore,we have provided engineers with a tool that can predict levels of damages associated with different failure mechanisms based on variations in rock mass quality and in situ stress state.Our method is a numerically developed,multi-variate approach that can aid engineers in making informed decisions about the robustness of underground tunnels.
基金funding support from the science and technology innovation Program of Hunan Province(Grant No.2023RC1017)Hunan Provincial Postgraduate Research and Innovation Project(Grant No.CX20220109)National Natural Science Foundation of China Youth Fund(Grant No.52208378).
文摘Machine learning(ML)provides a new surrogate method for investigating groundwater flow dynamics in unsaturated soils.Traditional pure data-driven methods(e.g.deep neural network,DNN)can provide rapid predictions,but they do require sufficient on-site data for accurate training,and lack interpretability to the physical processes within the data.In this paper,we provide a physics and equalityconstrained artificial neural network(PECANN),to derive unsaturated infiltration solutions with a small amount of initial and boundary data.PECANN takes the physics-informed neural network(PINN)as a foundation,encodes the unsaturated infiltration physical laws(i.e.Richards equation,RE)into the loss function,and uses the augmented Lagrangian method to constrain the learning process of the solutions of RE by adding stronger penalty for the initial and boundary conditions.Four unsaturated infiltration cases are designed to test the training performance of PECANN,i.e.one-dimensional(1D)steady-state unsaturated infiltration,1D transient-state infiltration,two-dimensional(2D)transient-state infiltration,and 1D coupled unsaturated infiltration and deformation.The predicted results of PECANN are compared with the finite difference solutions or analytical solutions.The results indicate that PECANN can accurately capture the variations of pressure head during the unsaturated infiltration,and present higher precision and robustness than DNN and PINN.It is also revealed that PECANN can achieve the same accuracy as the finite difference method with fewer initial and boundary training data.Additionally,we investigate the effect of the hyperparameters of PECANN on solving RE problem.PECANN provides an effective tool for simulating unsaturated infiltration.
基金supported by the National Key Research and Development Program of China(Grant No.2022YFC3080200)the National Natural Science Foundation of China(Grant No.42022053)the China Postdoctoral Science Foundation(Grant No.2023M731264).
文摘Natural slopes usually display complicated exposed rock surfaces that are characterized by complex and substantial terrain undulation and ubiquitous undesirable phenomena such as vegetation cover and rockfalls.This study presents a systematic outcrop research of fracture pattern variations in a complicated rock slope,and the qualitative and quantitative study of the complex phenomena impact on threedimensional(3D)discrete fracture network(DFN)modeling.As the studies of the outcrop fracture pattern have been so far focused on local variations,thus,we put forward a statistical analysis of global variations.The entire outcrop is partitioned into several subzones,and the subzone-scale variability of fracture geometric properties is analyzed(including the orientation,the density,and the trace length).The results reveal significant variations in fracture characteristics(such as the concentrative degree,the average orientation,the density,and the trace length)among different subzones.Moreover,the density of fracture sets,which is approximately parallel to the slope surface,exhibits a notably higher value compared to other fracture sets across all subzones.To improve the accuracy of the DFN modeling,the effects of three common phenomena resulting from vegetation and rockfalls are qualitatively analyzed and the corresponding quantitative data processing solutions are proposed.Subsequently,the 3D fracture geometric parameters are determined for different areas of the high-steep rock slope in terms of the subzone dimensions.The results show significant variations in the same set of 3D fracture parameters across different regions with density differing by up to tenfold and mean trace length exhibiting differences of 3e4 times.The study results present precise geological structural information,improve modeling accuracy,and provide practical solutions for addressing complex outcrop issues.
基金supported by the Research Council of Norway under contracts 223252/F50 and 300844/F50the Trond Mohn Foundation。
文摘Global images of auroras obtained by cameras on spacecraft are a key tool for studying the near-Earth environment.However,the cameras are sensitive not only to auroral emissions produced by precipitating particles,but also to dayglow emissions produced by photoelectrons induced by sunlight.Nightglow emissions and scattered sunlight can contribute to the background signal.To fully utilize such images in space science,background contamination must be removed to isolate the auroral signal.Here we outline a data-driven approach to modeling the background intensity in multiple images by formulating linear inverse problems based on B-splines and spherical harmonics.The approach is robust,flexible,and iteratively deselects outliers,such as auroral emissions.The final model is smooth across the terminator and accounts for slow temporal variations and large-scale asymmetries in the dayglow.We demonstrate the model by using the three far ultraviolet cameras on the Imager for Magnetopause-to-Aurora Global Exploration(IMAGE)mission.The method can be applied to historical missions and is relevant for upcoming missions,such as the Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)mission.
基金support from the National Key R&D plan(Grant No.2022YFC3004303)the National Natural Science Foundation of China(Grant No.42107161)+3 种基金the State Key Laboratory of Hydroscience and Hydraulic Engineering(Grant No.2021-KY-04)the Open Research Fund Program of State Key Laboratory of Hydroscience and Engineering(sklhse-2023-C-01)the Open Research Fund Program of Key Laboratory of the Hydrosphere of the Ministry of Water Resources(mklhs-2023-04)the China Three Gorges Corporation(XLD/2117).
文摘Rock fragmentation plays a critical role in rock avalanches,yet conventional approaches such as classical granular flow models or the bonded particle model have limitations in accurately characterizing the progressive disintegration and kinematics of multi-deformable rock blocks during rockslides.The present study proposes a discrete-continuous numerical model,based on a cohesive zone model,to explicitly incorporate the progressive fragmentation and intricate interparticle interactions inherent in rockslides.Breakable rock granular assemblies are released along an inclined plane and flow onto a horizontal plane.The numerical scenarios are established to incorporate variations in slope angle,initial height,friction coefficient,and particle number.The evolutions of fragmentation,kinematic,runout and depositional characteristics are quantitatively analyzed and compared with experimental and field data.A positive linear relationship between the equivalent friction coefficient and the apparent friction coefficient is identified.In general,the granular mass predominantly exhibits characteristics of a dense granular flow,with the Savage number exhibiting a decreasing trend as the volume of mass increases.The process of particle breakage gradually occurs in a bottom-up manner,leading to a significant increase in the angular velocities of the rock blocks with increasing depth.The simulation results reproduce the field observations of inverse grading and source stratigraphy preservation in the deposit.We propose a disintegration index that incorporates factors such as drop height,rock mass volume,and rock strength.Our findings demonstrate a consistent linear relationship between this index and the fragmentation degree in all tested scenarios.
基金supported in part by the National Key Research and Development Program of China(2022YFB3305300)the National Natural Science Foundation of China(62173178).
文摘Ethylene glycol(EG)plays a pivotal role as a primary raw material in the polyester industry,and the syngas-to-EG route has become a significant technical route in production.The carbon monoxide(CO)gas-phase catalytic coupling to synthesize dimethyl oxalate(DMO)is a crucial process in the syngas-to-EG route,whereby the composition of the reactor outlet exerts influence on the ultimate quality of the EG product and the energy consumption during the subsequent separation process.However,measuring product quality in real time or establishing accurate dynamic mechanism models is challenging.To effectively model the DMO synthesis process,this study proposes a hybrid modeling strategy that integrates process mechanisms and data-driven approaches.The CO gas-phase catalytic coupling mechanism model is developed based on intrinsic kinetics and material balance,while a long short-term memory(LSTM)neural network is employed to predict the macroscopic reaction rate by leveraging temporal relationships derived from archived measurements.The proposed model is trained semi-supervised to accommodate limited-label data scenarios,leveraging historical data.By integrating these predictions with the mechanism model,the hybrid modeling approach provides reliable and interpretable forecasts of mass fractions.Empirical investigations unequivocally validate the superiority of the proposed hybrid modeling approach over conventional data-driven models(DDMs)and other hybrid modeling techniques.
基金funded by the National Natural Science Program of China(2022YFD1900300)the China Scholarship Council(CSC)through the CSC-CSIRO(Commonwealth Scientific and Industrial Research Organisation)Joint Ph D Program,the Zhumadian Major Scientific and Technological Innovation Project,China(170109564016)the Huanghuai University Scientific Research Foundation,China(502310020017)。
文摘We used field experimental data to evaluate the ability of the agricultural production system model (APSIM) to simulate soil P availability,maize biomass and grain yield in response to P fertilizer applications on a fluvo-aquic soil in the North China Plain.Crop and soil data from a 2-year experiment with three P fertilizer application rates(0,75 and 300 kg P_(2)O_(5) ha^(–1)) were used to calibrate the model.Sensitivity analysis was carried out to investigate the influence of APSIM SoilP parameters on the simulated P availability in soil and maize growth.Crop and soil P parameters were then derived by matching or relating the simulation results to observed crop biomass,yield,P uptake and Olsen-P in soil.The re-parameterized model was further validated against 2 years of independent data at the same sites.The re-parameterized model enabled good simulation of the maize leaf area index (LAI),biomass,grain yield,P uptake,and grain P content in response to different levels of P additions against both the calibration and validation datasets.Our results showed that APSIM needs to be re-parameterized for simulation of maize LAI dynamics through modification of leaf size curve and a reduction in the rate of leaf senescence for modern staygreen maize cultivars in China.The P concentration limits (maximum and minimum P concentrations in organs)at different stages also need to be adjusted.Our results further showed a curvilinear relationship between the measured Olsen-P concentration and simulated labile P content,which could facilitate the initialization of APSIM P pools in the NCP with Olsen-P measurements in future studies.It remains difficult to parameterize the APSIM SoilP module due to the conceptual nature of the pools and simplified conceptualization of key P transformation processes.A fundamental understanding still needs to be developed for modelling and predicting the fate of applied P fertilizers in soils with contrasting physical and chemical characteristics.
基金The Construction S&T Project of the Department of Transportation of Sichuan Province(Grant No.2023A02)the National Natural Science Foundation of China(No.52109135).
文摘The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling efficiency of underground engineering,a modularized and parametric modeling cloud server is developed by using Python codes.The basic framework of the cloud server is as follows:input the modeling parameters into the web platform,implement Rhino software and FLAC3D software to model and run simulations in the cloud server,and return the simulation results to the web platform.The modeling program can automatically generate instructions that can run the modeling process in Rhino based on the input modeling parameters.The main modules of the modeling program include modeling the 3D geological structures,the underground engineering structures,and the supporting structures as well as meshing the geometric models.In particular,various cross-sections of underground caverns are crafted as parametricmodules in themodeling program.Themodularized and parametric modeling program is used for a finite element simulation of the underground powerhouse of the Shuangjiangkou Hydropower Station.This complicatedmodel is rapidly generated for the simulation,and the simulation results are reasonable.Thus,this modularized and parametric modeling program is applicable for three-dimensional finite element simulations and analyses.
基金Equinor for financing the R&D projectthe Institute of Science and Technology of Petroleum Geophysics of Brazil for supporting this research。
文摘We apply stochastic seismic inversion and Bayesian facies classification for porosity modeling and igneous rock identification in the presalt interval of the Santos Basin. This integration of seismic and well-derived information enhances reservoir characterization. Stochastic inversion and Bayesian classification are powerful tools because they permit addressing the uncertainties in the model. We used the ES-MDA algorithm to achieve the realizations equivalent to the percentiles P10, P50, and P90 of acoustic impedance, a novel method for acoustic inversion in presalt. The facies were divided into five: reservoir 1,reservoir 2, tight carbonates, clayey rocks, and igneous rocks. To deal with the overlaps in acoustic impedance values of facies, we included geological information using a priori probability, indicating that structural highs are reservoir-dominated. To illustrate our approach, we conducted porosity modeling using facies-related rock-physics models for rock-physics inversion in an area with a well drilled in a coquina bank and evaluated the thickness and extension of an igneous intrusion near the carbonate-salt interface. The modeled porosity and the classified seismic facies are in good agreement with the ones observed in the wells. Notably, the coquinas bank presents an improvement in the porosity towards the top. The a priori probability model was crucial for limiting the clayey rocks to the structural lows. In Well B, the hit rate of the igneous rock in the three scenarios is higher than 60%, showing an excellent thickness-prediction capability.
基金supported by the National Natural Science Foundation of China(21908234)the National Key Research&Development Program of China(2020YFB0606404)+1 种基金the Inner Mongolia Science and Technology Agency Program(2019CG058)Shanxi Province Natural Science Foundation(202103021223063).
文摘The influences of particle size,shape,and catalyst distribution on the reactivity and hydrocarbon product selectivity of a cobalt-based catalyst for Fischer-Tropsch synthesis were investigated in the present work.A self-consistent kinetic model for Fischer-Tropsch reaction proposed here was found to correlate experimental data well and hence was used to describe the consumption rates of reactants and formation rates of hydrocarbon products.The perturbed-chain statistical associating fluid theory equation of state was used to describe vapor-liquid equilibrium behavior associated with Fischer-Tropsch reaction.Local interaction between intraparticle diffusion and Fischer-Tropsch reaction was investigated in detail.Results showed that in order to avoid the adverse influence of intraparticle diffusional limitations on catalyst reactivity and product selectivity,the use of small particles is necessary.Large eggshell spherical particles are shown to keep the original catalyst reactivity and enhance the selectivity of heavy hydrocarbon products.The suitable layer thickness for a spherical particle with a diameter of 2 mm is nearly 0.15 mm.With the same outer diameter of 2 mm,the catalyst reactivity and heavy product selectivity of hollow cylindrical particles with a layer thickness of 0.25 mm are found to be larger than eggshell spherical particles.From the viewpoint of catalytic performance,hollow cylindrical particles are a better choice for industrial applications.
文摘The global populationhas beenandwill continue to be severely impacted by theCOVID-19 epidemic.The primary objective of this research is to demonstrate the future impact of COVID-19 on those who suffer from other fatal conditions such as cancer,heart disease,and diabetes.Here,using ordinary differential equations(ODEs),two mathematical models are developed to explain the association between COVID-19 and cancer and between COVID-19 and diabetes and heart disease.After that,we highlight the stability assessments that can be applied to these models.Sensitivity analysis is used to examine how changes in certain factors impact different aspects of disease.The sensitivity analysis showed that many people are still nervous about seeing a doctor due to COVID-19,which could result in a dramatic increase in the diagnosis of various ailments in the years to come.The correlation between diabetes and cardiovascular illness is also illustrated graphically.The effects of smoking and obesity are also found to be significant in disease compartments.Model fitting is also provided for interpreting the relationship between real data and the results of thiswork.Diabetic people,in particular,need tomonitor their health conditions closely and practice heart health maintenance.People with heart diseases should undergo regular checks so that they can protect themselves from diabetes and take some precautions including suitable diets.The main purpose of this study is to emphasize the importance of regular checks,to warn people about the effects of COVID-19(including avoiding healthcare centers and doctors because of the spread of infectious diseases)and to indicate the importance of family history of cancer,heart diseases and diabetes.The provision of the recommendations requires an increase in public consciousness.
基金supported by the National Science Fund for Distinguished Young Scholars(No.61925102)the National Natural Science Foundation of China(No.62201086,92167202,62201087,62101069)BUPT-CMCC Joint Innovation Center,and State Key Laboratory of IPOC(BUPT)(No.IPOC2023ZT02),China。
文摘Visible light communication(VLC)has attracted much attention in the research of sixthgeneration(6G)systems.Furthermore,channel modeling is the foundation for designing efficient and robust VLC systems.In this paper,we present extensive VLC channel measurement campaigns in indoor environments,i.e.,an office and a corridor.Based on the measured data,the large-scale fading characteristics and multipath-related characteristics,including omnidirectional optical path loss(OPL),K-factor,power angular spectrum(PAS),angle spread(AS),and clustering characteristics,are analyzed and modeled through a statistical method.Based on the extracted statistics of the above-mentioned channel characteristics,we propose a statistical spatial channel model(SSCM)capable of modeling multipath in the spatial domain.Furthermore,the simulated statistics of the proposed model are compared with the measured statistics.For instance,in the office,the simulated path loss exponent(PLE)and the measured PLE are 1.96and 1.97,respectively.And,the simulated medians of AS and measured medians of AS are 25.94°and 24.84°,respectively.Generally,the fact that the simulated results fit well with measured results has demonstrated the accuracy of our SSCM.
基金This work was supported by the special fund of the State Key Laboratory of Intense Pulsed Radiation Simulation and Effect(No.SKLIPR2011).
文摘Simulating the total ionizing dose(TID)of an electrical system using transistor-level models can be difficult and expensive,particularly for digital-integrated circuits(ICs).In this study,a method for modeling TID effects in complementary metaloxide semiconductor(CMOS)digital ICs based on the input/output buffer information specification(IBIS)was proposed.The digital IC was first divided into three parts based on its internal structure:the input buffer,output buffer,and functional area.Each of these three parts was separately modeled.Using the IBIS model,the transistor V-I characteristic curves of the buffers were processed,and the physical parameters were extracted and modeled using VHDL-AMS.In the functional area,logic functions were modeled in VHDL according to the data sheet.A golden digital IC model was developed by combining the input buffer,output buffer,and functional area models.Furthermore,the golden ratio was reconstructed based on TID experimental data,enabling the assessment of TID effects on the threshold voltage,carrier mobility,and time series of the digital IC.TID experiments were conducted using a CMOS non-inverting multiplexer,NC7SZ157,and the results were compared with the simulation results,which showed that the relative errors were less than 2%at each dose point.This confirms the practicality and accuracy of the proposed modeling method.The TID effect model for digital ICs developed using this modeling technique includes both the logical function of the IC and changes in electrical properties and functional degradation impacted by TID,which has potential applications in the design of radiation-hardening tolerance in digital ICs.
文摘This paper describes numerical simulation of hydraulic fracturing using fracture-based continuum modeling(FBCM)of coupled geomechanical-hydrological processes to evaluate a technique for high-density fracturing and fracture caging.The simulations are innovative because of modeling discrete fractures explicitly in continuum analysis.A key advantage of FBCM is that fracture initiation and propagation are modeled explicitly without changing the domain grid(i.e.no re-meshing).Further,multiple realizations of a preexisting fracture distribution can be analyzed using the same domain grid.The simulated hydraulic fracturing technique consists of pressurizing multiple wells simultaneously:initially without permeating fluids into the rock,to seed fractures uniformly and at high density in the wall rock of the wells;followed by fluid injection to propagate the seeded fracture density hydraulically.FBCM combines the ease of continuum modeling with the potential accuracy of modeling discrete fractures and fracturing explicitly.Fractures are modeled as piecewise planar based on intersections with domain elements;fracture geometry stored as continuum properties is used to calculate parameters needed to model individual fractures;and rock behavior is modeled through tensorial aggregation of the behavior of discrete fractures and unfractured rock.Simulations are presented for previously unfractured rock and for rock with preexisting fractures of horizontal,shallow-dipping,steeply dipping,or vertical orientation.Simulations of a single-well model are used to determine the pattern and spacing for a multiple-well design.The results illustrate high-density fracturing and fracture caging through simultaneous fluid injection in multiple wells:for previously unfractured rock or rock with preexisting shallow-dipping or horizontal fractures,and in situ vertical compressive stress greater than horizontal.If preexisting fractures are steeply dipping or vertical,and considering the same in situ stress condition,well pressurization without fluid permeation appears to be the only practical way to induce new fractures and contain fracturing within the target domain.
基金This research is supported by the National Natural Science Foundation of China(NSFC)under contract no.42274147.
文摘Staggered-grid finite-difference(SGFD)schemes have been widely used in acoustic wave modeling for geophysical problems.Many improved methods are proposed to enhance the accuracy of numerical modeling.However,these methods are inevitably limited by the maximum Courant-Friedrichs-Lewy(CFL)numbers,making them unstable when modeling with large time sampling intervals or small grid spacings.To solve this problem,we extend a stable SGFD scheme by controlling SGFD dispersion relations and maximizing the maximum CFL numbers.First,to improve modeling stability,we minimize the error between the FD dispersion relation and the exact relation in the given wave-number region,and make the FD dispersion approach a given function outside the given wave-number area,thus breaking the conventional limits of the maximum CFL number.Second,to obtain high modeling accuracy,we use the SGFD scheme based on the Remez algorithm to compute the FD coefficients.In addition,the hybrid absorbing boundary condition is adopted to suppress boundary reflections and we find a suitable weighting coefficient for the proposed scheme.Theoretical derivation and numerical modeling demonstrate that the proposed scheme can maintain high accuracy in the modeling process and the value of the maximum CFL number of the proposed scheme can exceed that of the conventional SGFD scheme when adopting a small maximum effective wavenumber,indicating that the proposed scheme improves stability during the modeling.
基金supported by the National Natural Science Foundation of China(Grant Nos.82173620 to Yang Zhao and 82041024 to Feng Chen)partially supported by the Bill&Melinda Gates Foundation(Grant No.INV-006371 to Feng Chen)Priority Academic Program Development of Jiangsu Higher Education Institutions.
文摘Deterministic compartment models(CMs)and stochastic models,including stochastic CMs and agent-based models,are widely utilized in epidemic modeling.However,the relationship between CMs and their corresponding stochastic models is not well understood.The present study aimed to address this gap by conducting a comparative study using the susceptible,exposed,infectious,and recovered(SEIR)model and its extended CMs from the coronavirus disease 2019 modeling literature.We demonstrated the equivalence of the numerical solution of CMs using the Euler scheme and their stochastic counterparts through theoretical analysis and simulations.Based on this equivalence,we proposed an efficient model calibration method that could replicate the exact solution of CMs in the corresponding stochastic models through parameter adjustment.The advancement in calibration techniques enhanced the accuracy of stochastic modeling in capturing the dynamics of epidemics.However,it should be noted that discrete-time stochastic models cannot perfectly reproduce the exact solution of continuous-time CMs.Additionally,we proposed a new stochastic compartment and agent mixed model as an alternative to agent-based models for large-scale population simulations with a limited number of agents.This model offered a balance between computational efficiency and accuracy.The results of this research contributed to the comparison and unification of deterministic CMs and stochastic models in epidemic modeling.Furthermore,the results had implications for the development of hybrid models that integrated the strengths of both frameworks.Overall,the present study has provided valuable epidemic modeling techniques and their practical applications for understanding and controlling the spread of infectious diseases.
文摘This research aims to optimize the utilization of long-term sea level data from the TOPEX/Poseidon,Jason1,Jason2,and Jason3 altimetry missions for tidal modeling.We generate a time series of along-track observations and apply a developed method to produce tidal models with specific tidal constituents for each location.Our tidal modeling methodology follows an iterative process:partitioning sea surface height(SSH)observations into analysis/training and prediction/validation parts and ultimately identi-fying the set of tidal constituents that provide the best predictions at each time series location.The study focuses on developing 1256 time series along the altimetry tracks over the Baltic Sea,each with its own set of tidal constituents.Verification of the developed tidal models against the sSH observations within the prediction/validation part reveals mean absolute error(MAE)values ranging from 0.0334 m to 0.1349 m,with an average MAE of 0.089 m.The same validation process is conducted on the FES2014 and EOT20 global tidal models,demonstrating that our tidal model,referred to as BT23(short for Baltic Tide 2023),outperforms both models with an average MAE improvement of 0.0417 m and 0.0346 m,respectively.In addition to providing details on the development of the time series and the tidal modeling procedure,we offer the 1256 along-track time series and their associated tidal models as supplementary materials.We encourage the satellite altimetry community to utilize these resources for further research and applications.
基金support received from the National Natural Science Foundation of China(GrantNos.62204204 and 52175148)Science and Technology Innovation 2030-Major Project(Grant No.2022ZD0208601)+1 种基金Shanghai Sailing Program(Grant No.21YF1451000)Presidential Foundation of CAEP(Grant No.YZJJZQ2022001).
文摘Implanted neural probes can detect weak discharges of neurons in the brain by piercing soft brain tissue,thus as important tools for brain science research,as well as diagnosis and treatment of brain diseases.However,the rigid neural probes,such as Utah arrays,Michigan probes,and metal microfilament electrodes,are mechanically unmatched with brain tissue and are prone to rejection and glial scarring after implantation,which leads to a significant degradation in the signal quality with the implantation time.In recent years,flexible neural electrodes are rapidly developed with less damage to biological tissues,excellent biocompatibility,and mechanical compliance to alleviate scarring.Among them,the mechanical modeling is important for the optimization of the structure and the implantation process.In this review,the theoretical calculation of the flexible neural probes is firstly summarized with the processes of buckling,insertion,and relative interaction with soft brain tissue for flexible probes from outside to inside.Then,the corresponding mechanical simulation methods are organized considering multiple impact factors to realize minimally invasive implantation.Finally,the technical difficulties and future trends of mechanical modeling are discussed for the next-generation flexible neural probes,which is critical to realize low-invasiveness and long-term coexistence in vivo.
基金Supported by National Natural Science Foundation of China(Grant No.72071179)ZJU-Sunon Joint Research Center of Smart Furniture,Zhejiang University,China.
文摘The human digital twin(HDT)emerges as a promising human-centric technology in Industry 5.0,but challenges remain in human modeling and simulation.Digital human modeling(DHM)provides solutions for modeling and simulating human physical and cognitive aspects to support ergonomic analysis.However,it has limitations in real-time data usage,personalized services,and timely interaction.The emerging HDT concept offers new possibilities by integrating multi-source data and artificial intelligence for continuous monitoring and assessment.Hence,this paper reviews the evolution from DHM to HDT and proposes a unified HDT framework from a human factors perspective.The framework comprises the physical twin,the virtual twin,and the linkage between these two.The virtual twin integrates human modeling and AI engines to enable model-data-hybrid-enabled simulation.HDT can potentially upgrade traditional ergonomic methods to intelligent services through real-time analysis,timely feedback,and bidirectional interactions.Finally,the future perspectives of HDT for industrial applications as well as technical and social challenges are discussed.In general,this study outlines a human factors perspective on HDT for the first time,which is useful for cross-disciplinary research and human factors innovation to enhance the development of HDT in industry.