Mathematical Constraints in Multiscale Subgrid-Scale Modeling of Nonlinear Systems

To shed light on the subgrid-seale （SGS） modeling methodology of nonlinear systems such as the Navier-Stokes turbulence, we define the concepts of assumption and restriction in the modeling procedure, which are shown by generalized derivation of three general mathematical constraints for different combinations of restrictions. These constraints are verified numerically in a one-dimensional nonlinear advection equation. This study is expected to inspire future research on the SGS modeling methodology of nonlinear systems.

Applying resolved-scale linearly forced isotropic turbulence in rational subgrid-scale modeling

In previous attempts of rational subgrid-scale (SGS) modeling by employing the Kolmogorov equation of filtered (KEF) quantities, it was necessary to assume that the resolved-scale second-order structure function is stationary. Forced isotropic turbulence is often used as a framework for establishing and validating such SGS models based on stationary restrictions, for it generates statistical stationary samples. However, traditional forcing method at low wavenumbers cannot provide an analytic form of forcing term for a complete KEF in physical space, which has been illustrated to be essential in the modeling of such SGS models. Thus, an alternative forcing method giving an analytic forcing term in physical space is needed for rational SGS modeling. Giving an analytic linear driving term in physical space, linearly forced isotropic turbulence should be considered an ideal theoretical framework for rational SGS modeling. In this paper, we demonstrate the feasibility of establishing a rational SGS model with stationary restriction based on linearly forced isotropic turbulence. The performance of this rational SGS model is validated. We, therefore, propose the use of linearly forced isotropic turbulence as a complement to free-decaying isotropic turbulence and low-wavenumber forced isotropic turbulence for SGS model validations.

Recent understanding on the subgrid-scale modeling of large-eddy simulation in physical space

In the last 50 years,the methodology of large-eddy simulation(LES)has been greatly developed,while lots of different subgridscale(SGS)models have appeared.However,the understanding of the procedure of SGS modeling is still not clear.The present contribution aims at reviewing the recent SGS models and,more importantly,expressing our recent understanding on the SGS modeling of LES in physical space.Taking the Kolmogorov equation for filtered quantities(KEF)as an example,it is argued that the KEF alone is not enough to be a closure method.Three physical laws are then introduced to complete this closure procedure and are expected to inspire the future researches of SGS modeling.

A promising approach for quantifying focal stroke modeling and assessing stroke progression:optical resolution photoacoustic microscopy photothrombosis

To investigate the mechanisms underlying the onset and progression of ischemic stroke,some methods have been proposed that can simultaneously monitor and create embolisms in the animal cerebral cortex.However,these methods often require complex systems and the effect of age on cerebral embolism has not been adequately studied,although ischemic stroke is strongly age-related.In this study,we propose an optical-resolution photoacoustic microscopy-based visualized photothrombosis methodology to create and monitor ischemic stroke in mice simultaneously using a 532 nm pulsed laser.We observed the molding process in mice of different ages and presented age-dependent vascular embolism differentiation.Moreover,we integrated optical coherence tomography angiography to investigate age-associated trends in cerebrovascular variability following a stroke.Our imaging data and quantitative analyses underscore the differential cerebrovascular responses to stroke in mice of different ages,thereby highlighting the technique's potential for evaluating cerebrovascular health and unraveling age-related mechanisms involved in ischemic strokes.

Stochastic modeling of subgrid-scale effects on particle motion in forced isotropic turbulence

The subgrid-scale effects on particle motion were investigated in forced isotropic turbulence by DNS and priorLES methods.In the DNS field,the importance of Kolmogorov scaling to preferential accumulation was validated by comparing the radial distribution functions under various particle Stokes numbers.The prior-LES fields were generated by filtering the DNS data.The subgrid-scale Stokes number(StSGS)is a useful tool for determining the effects of subgrid-scale eddies on particle motion.The subgrid-scale eddies tend to accumulate particles with StSGSb 1 and disperse particles with 1 b StSGSb 10.For particles with StSGS?1,the effects of subgrid-scale eddies on particle motion can be neglected.In order to restore the subgrid-scale effects,the Langevin-type stochastic model with optimized parameters was adopted in this study.This model is effective for the particles with StSGS N 1 while has an adverse impact on the particles with StSGSb 1.The results show that the Langevin-type stochastic model tends to smooth the particle distribution in the isotropic turbulence.

Evaluation of Subgrid-scale Models in Large-eddy Simulations of Turbulent Flow in a Centrifugal Pump Impeller

The current research of large eddy simulation （LES） of turbulent flow in pumps mainly concentrates in applying conventional subgrid-scale （SGS） model to simulate turbulent flow, which aims at obtaining the flow field in pump. The selection of SGS model is usually not considered seriously, so the accuracy and efficiency of the simulation cannot be ensured. Three SGS models including Smagorinsky-Lilly model, dynamic Smagorinsky model and dynamic mixed model are comparably studied by using the commercial CFD code Fluent combined with its user define function. The simulations are performed for the turbulent flow in a centrifugal pump impeller. The simulation results indicate that the mean flows predicted by the three SGS models agree well with the experimental data obtained from the test that detailed measurements of the flow inside the rotating passages of a six-bladed shrouded centrifugal pump impeller performed using particle image velocimetry （PIV） and laser Doppler velocimetry （LDV）. The comparable results show that dynamic mixed model gives the most accurate results for mean flow in the centrifugal pump impeller. The SGS stress of dynamic mixed model is decompose into the scale similar part and the eddy viscous part. The scale similar part of SGS stress plays a significant role in high curvature regions, such as the leading edge and training edge of pump blade. It is also found that the dynamic mixed model is more adaptive to compute turbulence in the pump impeller. The research results presented is useful to improve the computational accuracy and efficiency of LES for centrifugal pumps, and provide important reference for carrying out simulation in similar fluid machineries.

Predictive modeling for postoperative delirium in elderly patients with abdominal malignancies using synthetic minority oversampling technique

BACKGROUND Postoperative delirium,particularly prevalent in elderly patients after abdominal cancer surgery,presents significant challenges in clinical management.AIM To develop a synthetic minority oversampling technique(SMOTE)-based model for predicting postoperative delirium in elderly abdominal cancer patients.METHODS In this retrospective cohort study,we analyzed data from 611 elderly patients who underwent abdominal malignant tumor surgery at our hospital between September 2020 and October 2022.The incidence of postoperative delirium was recorded for 7 d post-surgery.Patients were divided into delirium and non-delirium groups based on the occurrence of postoperative delirium or not.A multivariate logistic regression model was used to identify risk factors and develop a predictive model for postoperative delirium.The SMOTE technique was applied to enhance the model by oversampling the delirium cases.The model’s predictive accuracy was then validated.RESULTS In our study involving 611 elderly patients with abdominal malignant tumors,multivariate logistic regression analysis identified significant risk factors for postoperative delirium.These included the Charlson comorbidity index,American Society of Anesthesiologists classification,history of cerebrovascular disease,surgical duration,perioperative blood transfusion,and postoperative pain score.The incidence rate of postoperative delirium in our study was 22.91%.The original predictive model(P1)exhibited an area under the receiver operating characteristic curve of 0.862.In comparison,the SMOTE-based logistic early warning model(P2),which utilized the SMOTE oversampling algorithm,showed a slightly lower but comparable area under the curve of 0.856,suggesting no significant difference in performance between the two predictive approaches.CONCLUSION This study confirms that the SMOTE-enhanced predictive model for postoperative delirium in elderly abdominal tumor patients shows performance equivalent to that of traditional methods,effectively addressing data imbalance.

Volume-averaged modeling of multiphase solidification with equiaxed crystal sedimentation in a steel ingot

Macrosegregation is a critical factor that limits the mechanical properties of materials.The impact of equiaxed crystal sedimentation on macrosegregation has been extensively studied,as it plays a significant role in determining the distribution of alloying elements and impurities within a material.To improve macrosegregation in steel connecting shafts,a multiphase solidification model that couples melt flow,heat transfer,microstructure evolution,and solute transport was established based on the volume-averaged Eulerian-Eulerian approach.In this model,the effects of liquid phase,equiaxed crystals,columnar dendrites,and columnar-to-equiaxed transition(CET)during solidification and evolution of microstructure can be considered simultaneously.The sedimentation of equiaxed crystals contributes to negative macrosegregation,where regions between columnar dendrites and equiaxed crystals undergo significant A-type positive macrosegregation due to the CET.Additionally,noticeable positive macrosegregation occurs in the area of final solidification in the ingot.The improvement in macrosegregation is beneficial for enhancing the mechanical properties of connecting shafts.To mitigate the thermal convection of molten steel resulting from excessive superheating,reducing the superheating during casting without employing external fields or altering the design of the ingot mold is indeed an effective approach to control macrosegregation.

Predicting the probability distribution of Martian rocks mechanical property based on microscale rock mechanical experiments and accurate grain-based modeling

The exploration of Mars would heavily rely on Martian rocks mechanics and engineering technology.As the mechanical property of Martian rocks is uncertain,it is of utmost importance to predict the probability distribution of Martian rocks mechanical property for the success of Mars exploration.In this paper,a fast and accurate probability distribution method for predicting the macroscale elastic modulus of Martian rocks was proposed by integrating the microscale rock mechanical experiments(micro-RME),accurate grain-based modeling(AGBM)and upscaling methods based on reliability principles.Firstly,the microstructure of NWA12564 Martian sample and elastic modulus of each mineral were obtained by micro-RME with TESCAN integrated mineral analyzer(TIMA)and nanoindentation.The best probability distribution function of the minerals was determined by Kolmogorov-Smirnov(K-S)test.Secondly,based on best distribution function of each mineral,the Monte Carlo simulations(MCS)and upscaling methods were implemented to obtain the probability distribution of upscaled elastic modulus.Thirdly,the correlation between the upscaled elastic modulus and macroscale elastic modulus obtained by AGBM was established.The accurate probability distribution of the macroscale elastic modulus was obtained by this correlation relationship.The proposed method can predict the probability distribution of Martian rocks mechanical property with any size and shape samples.

Background removal from global auroral images:Data-driven dayglow modeling

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.

Surrogate modeling for unsaturated infiltration via the physics and equality-constrained artificial neural networks

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.

A methodology for damage evaluation of underground tunnels subjected to static loading using numerical modeling

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.

Artificial intelligence-driven radiomics study in cancer:the role of feature engineering and modeling

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.

Modeling time-dependent mechanical behavior of hard rock considering excavation-induced damage and complex 3D stress states

To investigate the long-term stability of deep rocks,a three-dimensional(3D)time-dependent model that accounts for excavation-induced damage and complex stress state is developed.This model comprises three main components:a 3D viscoplastic isotropic constitutive relation that considers excavation damage and complex stress state,a quantitative relationship between critical irreversible deformation and complex stress state,and evolution characteristics of strength parameters.The proposed model is implemented in a self-developed numerical code,i.e.CASRock.The reliability of the model is validated through experiments.It is indicated that the time-dependent fracturing potential index(xTFPI)at a given time during the attenuation creep stage shows a negative correlation with the extent of excavationinduced damage.The time-dependent fracturing process of rock demonstrates a distinct interval effect of the intermediate principal stress,thereby highlighting the 3D stress-dependent characteristic of the model.Finally,the influence of excavation-induced damage and intermediate principal stress on the time-dependent fracturing characteristics of the surrounding rocks around the tunnel is discussed.

Progressive fragmentation of granular assemblies within rockslides: Insights from discrete-continuous numerical modeling

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.

The relationship between compartment models and their stochastic counterparts:A comparative study with examples of the COVID-19 epidemic 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 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.

Discontinuity development patterns and the challenges for 3D discrete fracture network modeling on complicated exposed rock surfaces

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.

Hybrid modeling for carbon monoxide gas-phase catalytic coupling to synthesize dimethyl oxalate process

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.

Exploring Capillary Fringe Flow:Quasilinear Modeling with Kirchhoff Transforms and Gardner Model

Recent studies have underscored the significance of the capillary fringe in hydrological and biochemical processes.Moreover,its role in shallow waters is expected to be considerable.Traditionally,the study of groundwater flow has centered on unsaturated-saturated zones,often overlooking the impact of the capillary fringe.In this study,we introduce a steady-state two-dimensional model that integrates the capillary fringe into a 2-D numerical solution.Our novel approach employs the potential form of the Richards equation,facilitating the determination of boundaries,pressures,and velocities across different ground surface zones.We utilized a two-dimensional Freefem++finite element model to compute the stationary solution.The validation of the model was conducted using experimental data.We employed the OFAT(One_Factor-At-Time)method to identify the most sensitive soil parameters and understand how changes in these parameters may affect the behavior and water dynamics of the capillary fringe.The results emphasize the role of hydraulic conductivity as a key parameter influencing capillary fringe shape and dynamics.Velocity values within the capillary fringe suggest the prevalence of horizontal flow.By variation of the water table level and the incoming flow q0,we have shown the correlation between water table elevation and the upper limit of the capillary fringe.

A Hierarchical Modeling Framework for Electrochemical Behaviors in Lithium-Ion Batteries with Detailed Structures

The accurate representation of lithium plating and aging phenomena has posed a persistent challenge within the battery research community.Empirical evidence underscores the pivotal role of cell structure in influencing aging behaviors and lithium plating within lithium-ion batteries(LIBs).Available lithium-ion plating models often falter in detailed description when integrating the structural intricacies.To address this challenge,this study proposes an innovative hierarchical model that intricately incorporates the layered rolling structure in cells.Notably,our model demonstrates a remarkable capacity to predict the non-uniform distribution of current density and overpotential along the rolling direction of LIBs.Subsequently,we delve into an insightful exploration of the structural factors that influence lithium plating behavior,leveraging the foundation laid by our established model.Furthermore,we easily update the hierarchical model by considering aging factors.This aging model effectively anticipates capacity fatigue and lithium plating tendencies across individual layers of LIBs,all while maintaining computational efficiency.In light of our findings,this model yields novel perspectives on capacity fatigue dynamics and local lithium plating behaviors,offering a substantial advancement compared to existing models.This research paves the way for more efficient and tailored LIB design and operation,with broad implications for energy storage technologies.