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.

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.

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.

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.

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.

Total ionizing dose effect modeling method for CMOS digital-integrated circuit

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.

Modularized and Parametric Modeling Technology for Finite Element Simulations of Underground Engineering under Complicated Geological Conditions

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.

Modeling Analysis of Factors Influencing Wind-Borne Seed Dispersal: A Case Study on Dandelion

A weed is a plant that thrives in areas of human disturbance, such as gardens, fields, pastures, waysides, and waste places where it is not intentionally cultivated. Dispersal affects community dynamics and vegetation response to global change. The process of seed disposal is influenced by wind, which plays a crucial role in determining the distance and probability of seed dispersal. Existing models of seed dispersal consider wind direction but fail to incorporate wind intensity. In this paper, a novel seed disposal model was proposed in this paper, incorporating wind intensity based on relevant references. According to various climatic conditions, including temperate, arid, and tropical regions, three specific regions were selected to establish a wind dispersal model that accurately reflects the density function distribution of dispersal distance. Additionally, dandelions growth is influenced by a multitude of factors, encompassing temperature, humidity, climate, and various environmental variables that necessitate meticulous consideration. Based on Factor Analysis model, which completely considers temperature, precipitation, solar radiation, wind, and land carrying capacity, a conclusion is presented, indicating that the growth of seeds is primarily influenced by plant attributes and climate conditions, with the former exerting a relatively stronger impact. Subsequently, the remaining two plants were chosen based on seed weight, yielding consistent conclusion.

Modeling of multiphase flow in low permeability porous media:Effect of wettability and pore structure properties

Multiphase flow in low permeability porous media is involved in numerous energy and environmental applications.However,a complete description of this process is challenging due to the limited modeling scale and the effects of complex pore structures and wettability.To address this issue,based on the digital rock of low permeability sandstone,a direct numerical simulation is performed considering the interphase drag and boundary slip to clarify the microscopic water-oil displacement process.In addition,a dual-porosity pore network model(PNM)is constructed to obtain the water-oil relative permeability of the sample.The displacement efficiency as a recovery process is assessed under different wetting and pore structure properties.Results show that microscopic displacement mechanisms explain the corresponding macroscopic relative permeability.The injected water breaks through the outlet earlier with a large mass flow,while thick oil films exist in rough hydrophobic surfaces and poorly connected pores.The variation of water-oil relative permeability is significant,and residual oil saturation is high in the oil-wet system.The flooding is extensive,and the residual oil is trapped in complex pore networks for hydrophilic pore surfaces;thus,water relative permeability is lower in the water-wet system.While the displacement efficiency is the worst in mixed-wetting systems for poor water connectivity.Microporosity negatively correlates with invading oil volume fraction due to strong capillary resistance,and a large microporosity corresponds to low residual oil saturation.This work provides insights into the water-oil flow from different modeling perspectives and helps to optimize the development plan for enhanced recovery.

Integrated Hydrological Modeling of the Godavari River Basin in Maharashtra Using the SWAT Model: Streamflow Simulation and Analysis

Hydrological modeling plays a crucial role in efficiently managing water resources and understanding the hydrologic behavior of watersheds. This study aims to simulate daily streamflow in the Godavari River Basin in Maharashtra using the Soil and Water Assessment Tool (SWAT). SWAT is a process-based hydrological model used to predict water balance components, sediment levels, and nutrient contamination. In this research, we used integrated remote sensing and GIS data, including Digital Elevation Models (DEM), land use and land cover (LULC) maps, soil maps, and observed precipitation and temperature data, as input for developing the SWAT model to assess surface runoff in this large river basin. The Godavari River Basin under study was divided into 25 sub-basins, comprising 151 hydrological response units categorized by unique land cover, soil, and slope characteristics using the SWAT model. The model was calibrated and validated against observed runoff data for two time periods: 2003-2006 and 2007-2010 respectively. Model performance was assessed using the Nash-Sutcliffe efficiency (NSE) and the coefficient of determination (R2). The results show the effectiveness of the SWAT2012 model, with R2 value of 0.84 during calibration and 0.86 during validation. NSE values also ranged from 0.84 during calibration to 0.85 during validation. These findings enhance our understanding of surface runoff dynamics in the Godavari River Basin under study and highlight the suit-ability of the SWAT model for this region.

Modeling and Validation of Diamagnetic Rotor Levitated by Permanent Magnetics

As an innovative,low-power consuming,and low-stiffness suspension approach,the diamagnetic levitation technique has attracted considerable interest because of its potential applicability in miniaturized mechanical systems.The foundation of a diamagnetic levitation system is mathematical modeling,which is essential for operating performance optimization and stability prediction.However,few studies on systematic mathematical modeling have been reported.In this study,a systematic mathematical model for a disc-shaped diamagnetically levitated rotor on a permanent magnet array is proposed.Based on the proposed model,the magnetic field distribution characteristics,diamagnetic levitation force characteristics(i.e.,levitation height and stiffness),and optimized theoretical conditions for realizing stable levitation are determined.Experiments are conducted to verify the feasibility of the proposed mathematical model.Theoretical predictions and experimental results indicate that increasing the levitation height enlarges the stable region.Moreover,with a further increase in the rotor radius,the stable regions of the rotor gradually diminish and even vanish.Thus,when the levitation height is fixed,a moderate rotor radius permits stable levitation.This study proposes a mathematical modeling method for a diamagnetic levitation system that has potential applications in miniaturized mechanical systems.

WT-FCTGN:A wavelet-enhanced fully connected time-gated neural network for complex noisy traffic flow modeling

Accurate forecasting of traffic flow provides a powerful traffic decision-making basis for an intelligent transportation system. However, the traffic data's complexity and fluctuation, as well as the noise produced during collecting information and summarizing original data of traffic flow, cause large errors in the traffic flow forecasting results. This article suggests a solution to the above mentioned issues and proposes a fully connected time-gated neural network based on wavelet reconstruction(WT-FCTGN). To eliminate the potential noise and strengthen the potential traffic trend in the data, we adopt the methods of wavelet reconstruction and periodic data introduction to preprocess the data. The model introduces fully connected time-series blocks to model all the information including time sequence information and fluctuation information in the flow of traffic, and establishes the time gate block to comprehend the periodic characteristics of the flow of traffic and predict its flow. The performance of the WT-FCTGN model is validated on the public Pe MS data set. The experimental results show that the WT-FCTGN model has higher accuracy, and its mean absolute error(MAE), mean absolute percentage error(MAPE) and root mean square error(RMSE) are obviously lower than those of the other algorithms. The robust experimental results prove that the WT-FCTGN model has good anti-noise ability.

Two-temperature modeling of lamellar cathode arc

A three-dimensional, two-temperature(2T) model of a lamellar cathode arc is constructed,drawing upon the conservation equations for mass, momentum, electron energy, and heavy particle energy, in addition to Maxwell's equations. The model aims to elucidate how the physical properties of electrons and heavy particles affect heat transfer and fluid flow in a lamellar cathode arc. This is achieved by solving and comparing the fields of electron temperature,heavy particle temperature, fluid flow, current density, and Lorentz force distribution under varying welding currents. The results show that the guiding effect of the lamellar cathode on current density, the inertial drag effect of moving arc, and the attraction effect of Lorentz force at the lamellar cathode tip primarily govern the distribution of the arc's physical fields. The guiding effect localizes the current density to the front end of the lamellar cathode, particularly where the discharge gap is minimal. Both the inertial drag effect and the attraction effect of Lorentz force direct arc flow toward its periphery. Under the influence of the aforementioned factors, the physical fields of the lamellar cathode arc undergo expansion and shift counter to the arc's direction of motion. A reduction in welding current substantially weakens the guiding effect,causing the arc's physical fields to deviate further in the direction opposite to the arc motion. In comparison with a cylindrical cathode arc, the physical fields of the lamellar cathode arc are markedly expanded, leading to a reduction in current density, electron temperature, heavy particle temperature, cathode jet flow velocity, and Lorentz force.

Combining field data and modeling to better understand maize growth response to phosphorus(P) fertilizer application and soil P dynamics in calcareous soils

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.

A Novel Framework to Construct S-Box Quantum Circuits Using System Modeling: Application to 4-Bit S-Boxes

Quantum computers accelerate many algorithms based on the superposition principle of quantum mechanics.The Grover algorithm provides significant performance to malicious users attacking symmetric key systems.Since the performance of attacks using quantum computers depends on the efficiency of the quantum circuit of the encryption algorithms,research research on the implementation of quantum circuits is essential.This paper presents a new framework to construct quantum circuits of substitution boxes(S-boxes)using system modeling.We model the quantum circuits of S-boxes using two layers:Toffoli and linear layers.We generate vector spaces based on the values of qubits used in the linear layers and apply them to find quantum circuits.The framework finds the circuit bymatching elements of vector spaces generated fromthe input and output of a given S-box,using the forward search or themeet-in-the-middle strategy.We developed a tool to apply this framework to 4-bit S-boxes.While the 4-bit S-box quantum circuit construction tool LIGHTER-R only finds circuits that can be implemented with four qubits,the proposed tool achieves the circuits with five qubits.The proposed tool can find quantum circuits of 4-bit odd permutations based on the controlled NOT,NOT,and Toffoli gates,whereas LIGHTER-R is unable to perform this task in the same environment.We expect this technique to become a critical step toward optimizing S-box quantum circuits.

A function-based behavioral modeling method for air combat simulation

ion by proposing multiple levels of cascaded hierarchi-cal structures from the perspective of function,namely,the func-tional decision tree.This method is developed to represent behavioral modeling of air combat systems,and its metamodel,execution mechanism,and code generation can provide a sound basis for function-based behavioral modeling.As a proof of con-cept,an air combat simulation is developed to validate this method and the results show that the fighter Alpha built using the proposed framework provides better performance than that using default scripts.