This paper investigates ruin,capital injection,and dividends for a two-dimensional risk model.The model posits that surplus levels of insurance companies are governed by a perturbed composite Poisson risk model.This m...This paper investigates ruin,capital injection,and dividends for a two-dimensional risk model.The model posits that surplus levels of insurance companies are governed by a perturbed composite Poisson risk model.This model introduces a dependence between the two surplus levels,present in both the associated perturbations and the claims resulting from common shocks.Critical levels of capital injection and dividends are established for each of the two risks.The surplus levels are observed discretely at fixed intervals,guiding decisions on capital injection,dividends,and ruin at these junctures.This study employs a two-dimensional Fourier cosine series expansion method to approximate the finite time expected discounted operating cost until ruin.The ensuing approximation error is also quantified.The validity and accuracy of the method are corroborated through numerical examples.Furthermore,the research delves into the optimal capital allocation problem.展开更多
Rip currents are a significant threat to swimmers worldwide,responsible for numerous drowning incidents each year.In Vietnam,Bai Dai Beach in Cam Ranh Bay,Khanh Hoa Province,has experienced an increase in drowning eve...Rip currents are a significant threat to swimmers worldwide,responsible for numerous drowning incidents each year.In Vietnam,Bai Dai Beach in Cam Ranh Bay,Khanh Hoa Province,has experienced an increase in drowning events due to rip currents in recent years.To address this issue,a comprehensive study was conducted based on developing a depth-averaged 2D hydrodynamic model to simulate rip currents in the Bai Dai-Cam Ranh coast.The HYDIST-2D numerical model was applied to simulate the rip current evolution in space and time for the study area.The results showed that the HYDIST-2D numerical model can accurately predict the location,magnitude,and microstructure of rip currents,including rip current speed,width,and length.The simulation results revealed that the rip current speed is greater during the low tide phase,with an average speed of 0.5 m s^(-1),while during high tide,the rip current speed is lower,around 0.1–0.8 m s^(-1).The width and length of the rip current also vary with the tide phase,with a wider and longer rip current observed during the low tide phase.The results also showed that the rip current speed and microstructure are influenced by the wave features,tide current,and bathymetry of the study area.The present study provides valuable insights into the dynamics of rip currents in the Bai Dai-Cam Ranh coast.The findings can be used to support the management of bathing activities and provide early warnings for potential risks associated with rip currents.展开更多
Systems with quenched disorder possess complex energy landscapes that are challenging to explore under conventional Monte Carlo methods.In this work,we implement an efficient entropy sampling scheme for accurate compu...Systems with quenched disorder possess complex energy landscapes that are challenging to explore under conventional Monte Carlo methods.In this work,we implement an efficient entropy sampling scheme for accurate computation of the entropy function in low-energy regions.The method is applied to the two-dimensional±J random-bond Ising model,where frustration is controlled by the fraction p of ferromagnetic bonds.We investigate the low-temperature paramagnetic–ferromagnetic phase boundary below the multicritical point at T_(N)=0.9530(4),P_(N)=0.89078(8),as well as the zerotemperature ferromagnetic–spin-glass transition.Finite-size scaling analysis reveals that the phase boundary for T<T_(N) exhibits reentrant behavior.By analyzing the evolution of the magnetizationresolved density of states g(E,M)and ground-state spin configurations against increasing frustration,we provide strong evidence that the zero-temperature transition is a mixed-order.Finite-size scaling conducted on the spin-glass side supports the validity of β=0,whereβis the magnetization exponent,with a correlation length exponentν=1.50(8).Our results provide new insights into the nature of the ferromagnetic-to-spin-glass phase transition in an extensively degenerate ground state.展开更多
Vehicle electrification,an important method for reducing carbon emissions from road transport,has been promoted globally.In this study,we analyze how individuals adapt to this transition in transportation and its subs...Vehicle electrification,an important method for reducing carbon emissions from road transport,has been promoted globally.In this study,we analyze how individuals adapt to this transition in transportation and its subsequent impact on urban structure.Considering the varying travel costs associated with electric and fuel vehicles,we analyze the dynamic choices of households concerning house locations and vehicle types in a two-dimensional monocentric city.A spatial equilibrium is developed to model the interactions between urban density,vehicle age and vehicle type.An agent-based microeconomic residential choice model dynamically coupled with a house rent market is developed to analyze household choices of home locations and vehicle energy types,considering vehicle ages and competition for public charging piles.Key findings from our proposed models show that the proportion of electric vehicles(EVs)peaks at over 50%by the end of the first scrappage period,accompanied by more than a 40%increase in commuting distance and time compared to the scenario with only fuel vehicles.Simulation experiments on a theoretical grid indicate that heterogeneity-induced residential segregation can lead to urban sprawl and congestion.Furthermore,households with EVs tend to be located farther from the city center,and an increase in EV ownership contributes to urban expansion.Our study provides insights into how individuals adapt to EV transitions and the resulting impacts on home locations and land use changes.It offers a novel perspective on the dynamic interactions between EV adoption and urban development.展开更多
The measurement of the pairing gap is crucial for investigating the physical properties of superconductors or superfluids.We propose a strategy to measure the pairing gap through the dynamical excitations.With the ran...The measurement of the pairing gap is crucial for investigating the physical properties of superconductors or superfluids.We propose a strategy to measure the pairing gap through the dynamical excitations.With the random phase approximation(RPA),we study the dynamical excitations of a two-dimensional attractive Fermi-Hubbard model by calculating its dynamical structure factor.Two distinct collective modes emerge:a Goldstone phonon mode at transferred momentum q=[0,0]and a roton mode at q=[p,p].The roton mode exhibits a sharp molecular peak in the low-energy regime.Notably,the area under the roton molecular peak scales with the square of the pairing gap,which holds even in three-dimensional and spin-orbit coupled(SOC)optical lattices.This finding suggests an experimental approach to measure the pairing gap in lattice systems by analyzing the dynamical structure factor at q=[p,p].展开更多
The work considers the problem of gas hydrate dissociation in a porous medium using the two-term Forchheimer law,corresponding to high flow rates of reservoir fluids.Such rates can arise during the decomposition of ga...The work considers the problem of gas hydrate dissociation in a porous medium using the two-term Forchheimer law,corresponding to high flow rates of reservoir fluids.Such rates can arise during the decomposition of gas hydrates,since a large amount of gas is released.Intensive emissions of gases from the earth’s interior are observed on the ocean floor.They are also associated with a large number of subvertical geological structures under the ocean floor,coming to the surface in the formof local ring funnels(pockmarks).Many similar objects have also been found on land.Particular interest in this problemis caused by climate threats associated with the release of greenhouse gases.The movement of gas released fromthe hydrate to the breakthrough channel is similar to the gas inflow to the well(without hydrate),which is usually described by a two-term filtration law.In this work,a mathematical model of gas hydrate dissociation with a nonlinear Forchheimer-type law ofmotion is developed.The systemis split in two blocks by physical processes,taking into account the quadratic correction to the velocity in the filtration law.The first block is responsible for the convective transfer of saturation parameters in the model,water,gas and hydrate saturations are taken into account.The second block corresponds to the equation of dissipative piezoconductivity with a different number of thermodynamic degrees of freedom,taking into account heat and mass transfer in a porous medium.The performed splitting allows using explicit-implicit difference schemes when solving problems and avoiding strong refinement of the step in time and space.For numerical modeling,the support operator method is used,which makes it possible to discretize partial differential equations on irregular grids,which allows taking into account the complex geometry and lithology of the reservoir.A difference scheme based on the support operator method is developed,which,due to the mutually consistent approximation of vector analysis operations(divergence and gradient),allows to take into account the various flux laws between adjacent grid cells,including quadratic corrections to the velocity.Based on the developed numerical algorithms and their program implementations,calculations of gas hydrate dissociation are performed both in a reservoir of simple geometric structure and in a heterogeneous reservoir of complex configuration.The results obtained correspond to the physics of the processes under consideration.展开更多
Phase transitions,as one of the most intriguing phenomena in nature,are divided into first-order phase transitions(FOPTs)and continuous ones in current classification.While the latter shows striking phenomena of scali...Phase transitions,as one of the most intriguing phenomena in nature,are divided into first-order phase transitions(FOPTs)and continuous ones in current classification.While the latter shows striking phenomena of scaling and universality,the former has recently also been demonstrated to exhibit scaling and universal behavior within a mesoscopic,coarse-grained Landau-Ginzburg theory.Here we apply this theory to a microscopic model-the paradigmatic Ising model,which undergoes FOPTs between two ordered phases below its critical temperature-and unambiguously demonstrate universal scaling behavior in such FOPTs.These results open the door for extending the theory to other microscopic FOPT systems and experimentally testing them to systematically uncover their scaling and universal behavior.展开更多
Accurate simulation of the horizontal-two-dimension(H2D)focused wave group in deep water requires high accuracy of a numerical model.The two-layer Boussinesq-type model(Liu and Fang,2016;Liu et al.,2018)with the highe...Accurate simulation of the horizontal-two-dimension(H2D)focused wave group in deep water requires high accuracy of a numerical model.The two-layer Boussinesq-type model(Liu and Fang,2016;Liu et al.,2018)with the highest spatial derivative of 2 has high accuracy in both linear and nonlinear properties.Based on the further development of the velocity equations(Liu et al.,2023),the H2D numerical model for water waves is established with the prediction-correction-iteration model in the finite difference method,and a composite fourth-order Adams-Bashforth-Moulton scheme is used for time integration.The wave generation method proposed by Hsiao et al.(2005)is applied and calibrated in this H2D model.The numerical calculations lead to the following three main conclusions:First,compared with the analytical solution of Stokes linear waves,the calculated velocity profiles show higher accuracy by using the improved velocity formulas.Second,the simulations of the focused multidirectional wave group are carried out,and good agreements are found,demonstrating that the present H2D numerical model shows high accuracy in simulating focused multidirectional wave groups,and the effectiveness of the improved velocity formulas is also validated.Furthermore,the velocity profiles throughout the computational domain at the time of maximum wave crest are given.Finally,the FFT method is used to obtain the amplitude with different frequencies for several locations,and the changes of the wavelet energy spectrum at different locations are presented for several cases.展开更多
AIGaN/GaN HEMTs are investigated by numerical simulation from the self-consistent solution of Schr6dinger-Poisson-hydrodynamic （HD） systems. The influences of polarization charge and quantum effects are considered i...AIGaN/GaN HEMTs are investigated by numerical simulation from the self-consistent solution of Schr6dinger-Poisson-hydrodynamic （HD） systems. The influences of polarization charge and quantum effects are considered in this model. Then the two-dimensional conduction band and electron distribution, electron temperature characteristics, Id versus Vd and Id versus Vg, transfer characteristics and transconductance curves are obtained. Corresponding analysis and discussion based on the simulation results are subsequently given.展开更多
Water quality models are important tools to support the optimization of aquatic ecosystem rehabilitation programs and assess their efficiency. Basing on the flow conditions of the Daqinghe River Mouth of the Dianchi L...Water quality models are important tools to support the optimization of aquatic ecosystem rehabilitation programs and assess their efficiency. Basing on the flow conditions of the Daqinghe River Mouth of the Dianchi Lake, China, a two-dimensional water quality model was developed in the research. The hydrodynamics module was numerically solved by the alternating direction iteration （ADI） method. The parameters of the water quality module were obtained through the in situ experiments and the laboratory analyses that were conducted from 2006 to 2007. The model was calibrated and verified by the observation data in 2007. Among the four modelled key variables, i.e., water level, COD （in CODcr）, NH4＋-N and PO43-P the minimum value of the coefficient of determination （COD） was 0.69, indicating the model performed reasonably well. The developed model was then applied to simulate the water quality changes at a downstream cross-section assuming that the designed restoration programs were implemented. According to the simulated results, the restoration programs could cut down the loads of COD and PO43-P about 15%. Such a load reduction, unfortunately, would have very little effect on the NH4^＋-N removal. Moreover, the water quality at the outlet cross-section would be still in class V （3838-02）, indicating more measures should be taken to further reduce the loads. The study demonstrated the capability of water quality models to support aquatic ecosystem restorations.展开更多
Instead of the capillary plasma generator(CPG),a discharge rod plasma generator(DRPG)is used in the30 mm electrothermal-chemical(ETC)gun to improve the ignition uniformity of the solid propellant.An axisymmetric two-d...Instead of the capillary plasma generator(CPG),a discharge rod plasma generator(DRPG)is used in the30 mm electrothermal-chemical(ETC)gun to improve the ignition uniformity of the solid propellant.An axisymmetric two-dimensional interior ballistics model of the solid propellant ETC gun(2D-IB-SPETCG)is presented to describe the process of the ETC launch.Both calculated pressure and projectile muzzle velocity accord well with the experimental results.The feasibility of the 2D-IB-SPETCG model is proved.Depending on the experimental data and initial parameters,detailed distribution of the ballistics parameters can be simulated.With the distribution of pressure and temperature of the gas phase and the propellant,the influence of plasma during the ignition process can be analyzed.Because of the radial flowing plasma,the propellant in the area of the DRPG is ignited within 0.01 ms,while all propellant in the chamber is ignited within 0.09 ms.The radial ignition delay time is much less than the axial delay time.During the ignition process,the radial pressure difference is less than 5 MPa at the place 0.025 m away from the breech.The radial ignition uniformity is proved.The temperature of the gas increases from several thousand K(conventional ignition)to several ten thousand K(plasma ignition).Compare the distribution of the density and temperature of the gas,we know that low density and high temperature gas appears near the exits of the DRPG,while high density and low temperature gas appears at the wall near the breech.The simulation of the 2D-IB-SPETCG model is an effective way to investigate the interior ballistics process of the ETC launch.The 2D-IB-SPETC model can be used for prediction and improvement of experiments.展开更多
Hydraulic models for the generation of flood inundation maps are not commonly applied in mountain river basins because of the difficulty in modeling the hydraulic behavior and the complex topography. This paper presen...Hydraulic models for the generation of flood inundation maps are not commonly applied in mountain river basins because of the difficulty in modeling the hydraulic behavior and the complex topography. This paper presents a comparative analysis of the performance of four twodimensional hydraulic models (HEC-RAS 2D, Iber 2D, Flood Modeller 2D, and PCSWMM 2D) with respect to the generation of flood inundation maps. The study area covers a 5-km reach of the Santa B-arbara River located in the Ecuadorian Andes, at 2330 masl, in Gualaceo. The model's performance was evaluated based on the water surface elevation and flood extent, in terms of the mean absolute difference and measure of fit. The analysis revealed that, for a given case, Iber 2D has the best performance in simulating the water level and inundation for flood events with 20- and 50-year return periods, respectively, followed by Flood Modeller 2D, HEC-RAS 2D, and PCSWMM 2D in terms of their performance. Grid resolution, the way in which hydraulic structures are mimicked, the model code, and the default value of the parameters are considered the main sources of prediction uncertainty.展开更多
The features of a quasi-two-dimensional( quasi-2D) model for simulating two-phase water hammer flows with vaporous cavity in a pipe are investigated. The quasi-2D model with discrete vaporous cavity in the pipe is pro...The features of a quasi-two-dimensional( quasi-2D) model for simulating two-phase water hammer flows with vaporous cavity in a pipe are investigated. The quasi-2D model with discrete vaporous cavity in the pipe is proposed in this paper. This model uses the quasi-2D model for pure liquid zone and one-dimensional( 1D) discrete vapor cavity model for vaporous cavity zone. The quasi-2D model solves two-dimensional equations for both axial and radial velocities and 1D equations for both pressure head and discharge by the method of characteristics. The 1D discrete vapor cavity model is used to simulate the vaporous cavity occurred when the pressure in the local pipe is lower than the vapor pressure of the liquid. The proposed model is used to simulate two-phase water flows caused by the rapid downstream valve closure in a reservoir-pipe-valve system.The results obtained by the proposed model are compared with those by the corresponding 1D model and the experimental ones provided by the literature,respectively. The comparison shows that the maximum pressure heads simulated by the proposed model are more accurate than those by the corresponding 1D model.展开更多
Erosion is an important issue in soil science and is related to many environmental problems,such as soil erosion and sediment transport.Establishing a simulation model suitable for soil erosion prediction is of great ...Erosion is an important issue in soil science and is related to many environmental problems,such as soil erosion and sediment transport.Establishing a simulation model suitable for soil erosion prediction is of great significance not only to accurately predict the process of soil separation by runoff,but also improve the physical model of soil erosion.In this study,we develop a graphic processing unit(GPU)-based numerical model that combines two-dimensional(2D)hydrodynamic and Green-Ampt(G-A)infiltration modelling to simulate soil erosion.A Godunov-type scheme on a uniform and structured square grid is then generated to solve the relevant shallow water equations(SWEs).The highlight of this study is the use of GPU-based acceleration technology to enable numerical models to simulate slope and watershed erosion in an efficient and high-resolution manner.The results show that the hydrodynamic model performs well in simulating soil erosion process.Soil erosion is studied by conducting calculation verification at the slope and basin scales.The first case involves simulating soil erosion process of a slope surface under indoor artificial rainfall conditions from 0 to 1000 s,and there is a good agreement between the simulated values and the measured values for the runoff velocity.The second case is a river basin experiment(Coquet River Basin)that involves watershed erosion.Simulations of the erosion depth change and erosion cumulative amount of the basin during a period of 1-40 h show an elevation difference of erosion at 0.5-3.0 m,especially during the period of 20-30 h.Nine cross sections in the basin are selected for simulation and the results reveal that the depth of erosion change value ranges from-0.86 to-2.79 m and the depth of deposition change value varies from 0.38 to 1.02 m.The findings indicate that the developed GPU-based hydrogeomorphological model can reproduce soil erosion processes.These results are valuable for rainfall runoff and soil erosion predictions on rilled hillslopes and river basins.展开更多
A global two-dimensional zonally averaged chemistry model is developed to study the chemi-cal composition of atmosphere. The region of the model is from 90°S to 90°N and from the ground to the altitude of 20...A global two-dimensional zonally averaged chemistry model is developed to study the chemi-cal composition of atmosphere. The region of the model is from 90°S to 90°N and from the ground to the altitude of 20 km with a resolution of 5° x 1 km. The wind field is residual circulation calcu-lated from diabatic rate. 34 species and 104 chemical and photochemical reactions are considered in the model. The sources of CH4, CO and NOx, which are divided into seasonal sources and non-seasonal sources, are parameterized as a function of latitude and time. The chemical composi-tion of atmosphere was simulated with emission level of CH4, CO and NOx in 1990. The results are compared with observations and other model results, showing that the model is successful to simu-late the atmospheric chemical composition and distribution of CH4. Key words Global two-dimensional chemistry model - Atmospheric composition - Emission This work was supported by the State Key Program for basic research “ Climate Dynamics and Cli-mate Prediction Theory” (Pandeng-yu-21).The authors would like to express their thanks to the National Oceanic and Atmospheric Administration (NOAA), Climate Monitoring and Diagnostics Laboratory (CMDL), Carbon Cycle Group for providing the observational data of CO and CH4.展开更多
River ice is a natural phenomenon in cold regions, influenced by meteorology, geomorphology, and hydraulic conditions. River ice processes involve complex interactions between hydrodynamic, mechanical, and thermal pro...River ice is a natural phenomenon in cold regions, influenced by meteorology, geomorphology, and hydraulic conditions. River ice processes involve complex interactions between hydrodynamic, mechanical, and thermal processes, and they are also influenced by weather and hydrologic conditions. Because natural rivers are serpentine, with bends, narrows, and straight reaches, the commonly-used one-dimensional river ice models and two-dimensional models based on the rectangular Cartesian coordinates are incapable of simulating the physical phenomena accurately. In order to accurately simulate the complicated river geometry and overcome the difficulties of numerical simulation resulting from both complex boundaries and differences between length and width scales, a two-dimensional river ice numerical model based on a boundary-fitted coordinate transformation method was developed. The presented model considers the influence of the frazil ice accumulation under ice cover and the shape of the leading edge of ice cover during the freezing process. The model is capable of determining the velocity field, the distribution of water temperature, the concentration distribution of frazil ice, the transport of floating ice, the progression, stability, and thawing of ice cover, and the transport, accumulation, and erosion of ice under ice cover. A MacCormack scheme was used to solve the equations numerically. The model was validated with field observations from the Hequ Reach of the Yellow River. Comparison of simulation results with field data indicates that the model is capable of simulating the river ice process with high accuracy.展开更多
In this article,a high-order scheme,which is formulated by combining the quadratic finite element method in space with a second-order time discrete scheme,is developed for looking for the numerical solution of a two-d...In this article,a high-order scheme,which is formulated by combining the quadratic finite element method in space with a second-order time discrete scheme,is developed for looking for the numerical solution of a two-dimensional nonlinear time fractional thermal diffusion model.The time Caputo fractional derivative is approximated by using the L2-1formula,the first-order derivative and nonlinear term are discretized by some second-order approximation formulas,and the quadratic finite element is used to approximate the spatial direction.The error accuracy O(h3+t2)is obtained,which is verified by the numerical results.展开更多
On the conditions of low-resolution radar, a parametric model for two-dimensional radar target is described here according to the theory of electromagnetic scattering and the geometrical theory of diffraction. A high ...On the conditions of low-resolution radar, a parametric model for two-dimensional radar target is described here according to the theory of electromagnetic scattering and the geometrical theory of diffraction. A high resolution estimation algorithm to extract the model parameters is also developed by building the relation of the scattering model and Prony model. The analysis of Cramer-Rao bound and simulation show that the method here has better statistical performance. The simulated analysis also indicates that the accurate extraction of the diffraction coefficient of scattering center is restricted by signal to noise ratio, radar center frequency and radar bandwidth.展开更多
Water vapor, cloud, and surface rainfall budgets associated with the landfall of Typhoon Krosa on 6-8 October 2007 are analyzed based on a two-dimensional cloud-resolving model simulation. The model is integrated with...Water vapor, cloud, and surface rainfall budgets associated with the landfall of Typhoon Krosa on 6-8 October 2007 are analyzed based on a two-dimensional cloud-resolving model simulation. The model is integrated with imposed zonally-uniform vertical velocity, zonal wind, horizontal temperature, and vapor advection from NCEP/Global Data Assimilation System （GDAS） data. The simulation data that are validated with observations are examined to study physical causes associated with surface rainfall processes during the landfall. The time- and domain-mean analysis shows that when Krosa approached the eastern coast of China on 6 October, the water vapor convergence over land caused a local atmospheric moistening and a net condensation that further produced surface rainfall and an increase of cloud hydrometeor concentration. Meanwhile, latent heating was balanced by advective cooling and a local atmospheric warming. One day later, the enhancement of net condensation led to an increase of surface rainfall and a local atmospheric drying, while the water vapor convergence weakened as a result of the landfall-induced deprivation of water vapor flux. At the same time, the latent heating is mainly compensated the advective cooling. Further weakening of vapor convergence on 8 October enhanced the local atmospheric drying while the net condensation and associated surface rainfall was maintained. The latent heating is balanced by advective cooling and a local atmospheric cooling.展开更多
A modified two-dimensional Eulerian air quality model was used to simulate both the gaseous and particulate pollutant concentrations during October 21-24, 2004 in the Pearl River Delta （PRD） region, China. The most ...A modified two-dimensional Eulerian air quality model was used to simulate both the gaseous and particulate pollutant concentrations during October 21-24, 2004 in the Pearl River Delta （PRD） region, China. The most significant improvement to the model is the added capability to predict the secondary organic aerosols （SOA） concentrations because of the inclusion of the SOA formation chemistry. The meteorological input data were prepared using the CALMET meteorological model. The concentrations of aerosol-bound species such as NO3^-, NH4^＋, SO4^2-, and SOA were calculated in the fine particle size range （〈2.5 μm）. The results of the two-dimensional model were compared to the measurements at the ground level during the PRD Intensive Monitoring Campaign （IMC）. Overall, there were good agreements between the measured and modeled concentrations of inorganic aerosol components and O3. Both the measured and the modeled results indicated that the maximum hourly O3 concentrations exceeded the China National Air Quality Standard. The predicted 24-h average SOA concentrations were in reasonable agreement with those predicted by the method of minimum OC/EC ratio.展开更多
基金supported by the Shihezi University High-Level Talents Research Startup Project(Project No.RCZK202521)the National Natural Science Foundation of China(Grant Nos.12271066,11871121,12171405)+1 种基金the Chongqing Natural Science Foundation Joint Fund for Innovation and Development Project(Project No.CSTB2024NSCQLZX0085)the Chongqing Normal University Foundation(Grant No.23XLB018).
文摘This paper investigates ruin,capital injection,and dividends for a two-dimensional risk model.The model posits that surplus levels of insurance companies are governed by a perturbed composite Poisson risk model.This model introduces a dependence between the two surplus levels,present in both the associated perturbations and the claims resulting from common shocks.Critical levels of capital injection and dividends are established for each of the two risks.The surplus levels are observed discretely at fixed intervals,guiding decisions on capital injection,dividends,and ruin at these junctures.This study employs a two-dimensional Fourier cosine series expansion method to approximate the finite time expected discounted operating cost until ruin.The ensuing approximation error is also quantified.The validity and accuracy of the method are corroborated through numerical examples.Furthermore,the research delves into the optimal capital allocation problem.
文摘Rip currents are a significant threat to swimmers worldwide,responsible for numerous drowning incidents each year.In Vietnam,Bai Dai Beach in Cam Ranh Bay,Khanh Hoa Province,has experienced an increase in drowning events due to rip currents in recent years.To address this issue,a comprehensive study was conducted based on developing a depth-averaged 2D hydrodynamic model to simulate rip currents in the Bai Dai-Cam Ranh coast.The HYDIST-2D numerical model was applied to simulate the rip current evolution in space and time for the study area.The results showed that the HYDIST-2D numerical model can accurately predict the location,magnitude,and microstructure of rip currents,including rip current speed,width,and length.The simulation results revealed that the rip current speed is greater during the low tide phase,with an average speed of 0.5 m s^(-1),while during high tide,the rip current speed is lower,around 0.1–0.8 m s^(-1).The width and length of the rip current also vary with the tide phase,with a wider and longer rip current observed during the low tide phase.The results also showed that the rip current speed and microstructure are influenced by the wave features,tide current,and bathymetry of the study area.The present study provides valuable insights into the dynamics of rip currents in the Bai Dai-Cam Ranh coast.The findings can be used to support the management of bathing activities and provide early warnings for potential risks associated with rip currents.
基金supported by NKRDPC-2022YFA1402802,NSFC-92165204the Research Grants Council of the HKSAR under Grant Nos.12304020 and 12301723+2 种基金Guangdong Provincial Key Laboratory of Magnetoelectric Physics and Devices under Grant No.2022B1212010008Guangdong Fundamental Research Center for Magnetoelectric Physics under Grant No.2024B0303390001Guangdong Provincial Quantum Science Strategic Initiative under Grant No.GDZX2401010。
文摘Systems with quenched disorder possess complex energy landscapes that are challenging to explore under conventional Monte Carlo methods.In this work,we implement an efficient entropy sampling scheme for accurate computation of the entropy function in low-energy regions.The method is applied to the two-dimensional±J random-bond Ising model,where frustration is controlled by the fraction p of ferromagnetic bonds.We investigate the low-temperature paramagnetic–ferromagnetic phase boundary below the multicritical point at T_(N)=0.9530(4),P_(N)=0.89078(8),as well as the zerotemperature ferromagnetic–spin-glass transition.Finite-size scaling analysis reveals that the phase boundary for T<T_(N) exhibits reentrant behavior.By analyzing the evolution of the magnetizationresolved density of states g(E,M)and ground-state spin configurations against increasing frustration,we provide strong evidence that the zero-temperature transition is a mixed-order.Finite-size scaling conducted on the spin-glass side supports the validity of β=0,whereβis the magnetization exponent,with a correlation length exponentν=1.50(8).Our results provide new insights into the nature of the ferromagnetic-to-spin-glass phase transition in an extensively degenerate ground state.
基金supported by National Natural Science Foundation of China(72288101,72361137002,and 72101018)the Dutch Research Council(NWO Grant 482.22.01).
文摘Vehicle electrification,an important method for reducing carbon emissions from road transport,has been promoted globally.In this study,we analyze how individuals adapt to this transition in transportation and its subsequent impact on urban structure.Considering the varying travel costs associated with electric and fuel vehicles,we analyze the dynamic choices of households concerning house locations and vehicle types in a two-dimensional monocentric city.A spatial equilibrium is developed to model the interactions between urban density,vehicle age and vehicle type.An agent-based microeconomic residential choice model dynamically coupled with a house rent market is developed to analyze household choices of home locations and vehicle energy types,considering vehicle ages and competition for public charging piles.Key findings from our proposed models show that the proportion of electric vehicles(EVs)peaks at over 50%by the end of the first scrappage period,accompanied by more than a 40%increase in commuting distance and time compared to the scenario with only fuel vehicles.Simulation experiments on a theoretical grid indicate that heterogeneity-induced residential segregation can lead to urban sprawl and congestion.Furthermore,households with EVs tend to be located farther from the city center,and an increase in EV ownership contributes to urban expansion.Our study provides insights into how individuals adapt to EV transitions and the resulting impacts on home locations and land use changes.It offers a novel perspective on the dynamic interactions between EV adoption and urban development.
基金supported by the National Natural Science Foundation of China[Grant Nos.U23A2073(P.Z.)and 11547034(H.Z.)].
文摘The measurement of the pairing gap is crucial for investigating the physical properties of superconductors or superfluids.We propose a strategy to measure the pairing gap through the dynamical excitations.With the random phase approximation(RPA),we study the dynamical excitations of a two-dimensional attractive Fermi-Hubbard model by calculating its dynamical structure factor.Two distinct collective modes emerge:a Goldstone phonon mode at transferred momentum q=[0,0]and a roton mode at q=[p,p].The roton mode exhibits a sharp molecular peak in the low-energy regime.Notably,the area under the roton molecular peak scales with the square of the pairing gap,which holds even in three-dimensional and spin-orbit coupled(SOC)optical lattices.This finding suggests an experimental approach to measure the pairing gap in lattice systems by analyzing the dynamical structure factor at q=[p,p].
基金the framework of the state assignment of Keldysh Institute of Applied Mathematics of RAS(Project No.125020701776-0)the Ministry of Education and Science of Russia for IO RAS(Project No.FMWE-2024-0018).
文摘The work considers the problem of gas hydrate dissociation in a porous medium using the two-term Forchheimer law,corresponding to high flow rates of reservoir fluids.Such rates can arise during the decomposition of gas hydrates,since a large amount of gas is released.Intensive emissions of gases from the earth’s interior are observed on the ocean floor.They are also associated with a large number of subvertical geological structures under the ocean floor,coming to the surface in the formof local ring funnels(pockmarks).Many similar objects have also been found on land.Particular interest in this problemis caused by climate threats associated with the release of greenhouse gases.The movement of gas released fromthe hydrate to the breakthrough channel is similar to the gas inflow to the well(without hydrate),which is usually described by a two-term filtration law.In this work,a mathematical model of gas hydrate dissociation with a nonlinear Forchheimer-type law ofmotion is developed.The systemis split in two blocks by physical processes,taking into account the quadratic correction to the velocity in the filtration law.The first block is responsible for the convective transfer of saturation parameters in the model,water,gas and hydrate saturations are taken into account.The second block corresponds to the equation of dissipative piezoconductivity with a different number of thermodynamic degrees of freedom,taking into account heat and mass transfer in a porous medium.The performed splitting allows using explicit-implicit difference schemes when solving problems and avoiding strong refinement of the step in time and space.For numerical modeling,the support operator method is used,which makes it possible to discretize partial differential equations on irregular grids,which allows taking into account the complex geometry and lithology of the reservoir.A difference scheme based on the support operator method is developed,which,due to the mutually consistent approximation of vector analysis operations(divergence and gradient),allows to take into account the various flux laws between adjacent grid cells,including quadratic corrections to the velocity.Based on the developed numerical algorithms and their program implementations,calculations of gas hydrate dissociation are performed both in a reservoir of simple geometric structure and in a heterogeneous reservoir of complex configuration.The results obtained correspond to the physics of the processes under consideration.
基金supported by the National Natural Science Foundation of China(Grant No.12175316).
文摘Phase transitions,as one of the most intriguing phenomena in nature,are divided into first-order phase transitions(FOPTs)and continuous ones in current classification.While the latter shows striking phenomena of scaling and universality,the former has recently also been demonstrated to exhibit scaling and universal behavior within a mesoscopic,coarse-grained Landau-Ginzburg theory.Here we apply this theory to a microscopic model-the paradigmatic Ising model,which undergoes FOPTs between two ordered phases below its critical temperature-and unambiguously demonstrate universal scaling behavior in such FOPTs.These results open the door for extending the theory to other microscopic FOPT systems and experimentally testing them to systematically uncover their scaling and universal behavior.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52171247,51779022,52071057,51709054)the National Key Research and Development Program of China(Grant No.2022YFC3106101)。
文摘Accurate simulation of the horizontal-two-dimension(H2D)focused wave group in deep water requires high accuracy of a numerical model.The two-layer Boussinesq-type model(Liu and Fang,2016;Liu et al.,2018)with the highest spatial derivative of 2 has high accuracy in both linear and nonlinear properties.Based on the further development of the velocity equations(Liu et al.,2023),the H2D numerical model for water waves is established with the prediction-correction-iteration model in the finite difference method,and a composite fourth-order Adams-Bashforth-Moulton scheme is used for time integration.The wave generation method proposed by Hsiao et al.(2005)is applied and calibrated in this H2D model.The numerical calculations lead to the following three main conclusions:First,compared with the analytical solution of Stokes linear waves,the calculated velocity profiles show higher accuracy by using the improved velocity formulas.Second,the simulations of the focused multidirectional wave group are carried out,and good agreements are found,demonstrating that the present H2D numerical model shows high accuracy in simulating focused multidirectional wave groups,and the effectiveness of the improved velocity formulas is also validated.Furthermore,the velocity profiles throughout the computational domain at the time of maximum wave crest are given.Finally,the FFT method is used to obtain the amplitude with different frequencies for several locations,and the changes of the wavelet energy spectrum at different locations are presented for several cases.
文摘AIGaN/GaN HEMTs are investigated by numerical simulation from the self-consistent solution of Schr6dinger-Poisson-hydrodynamic （HD） systems. The influences of polarization charge and quantum effects are considered in this model. Then the two-dimensional conduction band and electron distribution, electron temperature characteristics, Id versus Vd and Id versus Vg, transfer characteristics and transconductance curves are obtained. Corresponding analysis and discussion based on the simulation results are subsequently given.
基金supported by the National Hi-Tech Research and Development Program (863) of China (No.2007AA06A405, 2005AA6010100401)
文摘Water quality models are important tools to support the optimization of aquatic ecosystem rehabilitation programs and assess their efficiency. Basing on the flow conditions of the Daqinghe River Mouth of the Dianchi Lake, China, a two-dimensional water quality model was developed in the research. The hydrodynamics module was numerically solved by the alternating direction iteration （ADI） method. The parameters of the water quality module were obtained through the in situ experiments and the laboratory analyses that were conducted from 2006 to 2007. The model was calibrated and verified by the observation data in 2007. Among the four modelled key variables, i.e., water level, COD （in CODcr）, NH4＋-N and PO43-P the minimum value of the coefficient of determination （COD） was 0.69, indicating the model performed reasonably well. The developed model was then applied to simulate the water quality changes at a downstream cross-section assuming that the designed restoration programs were implemented. According to the simulated results, the restoration programs could cut down the loads of COD and PO43-P about 15%. Such a load reduction, unfortunately, would have very little effect on the NH4^＋-N removal. Moreover, the water quality at the outlet cross-section would be still in class V （3838-02）, indicating more measures should be taken to further reduce the loads. The study demonstrated the capability of water quality models to support aquatic ecosystem restorations.
文摘Instead of the capillary plasma generator(CPG),a discharge rod plasma generator(DRPG)is used in the30 mm electrothermal-chemical(ETC)gun to improve the ignition uniformity of the solid propellant.An axisymmetric two-dimensional interior ballistics model of the solid propellant ETC gun(2D-IB-SPETCG)is presented to describe the process of the ETC launch.Both calculated pressure and projectile muzzle velocity accord well with the experimental results.The feasibility of the 2D-IB-SPETCG model is proved.Depending on the experimental data and initial parameters,detailed distribution of the ballistics parameters can be simulated.With the distribution of pressure and temperature of the gas phase and the propellant,the influence of plasma during the ignition process can be analyzed.Because of the radial flowing plasma,the propellant in the area of the DRPG is ignited within 0.01 ms,while all propellant in the chamber is ignited within 0.09 ms.The radial ignition delay time is much less than the axial delay time.During the ignition process,the radial pressure difference is less than 5 MPa at the place 0.025 m away from the breech.The radial ignition uniformity is proved.The temperature of the gas increases from several thousand K(conventional ignition)to several ten thousand K(plasma ignition).Compare the distribution of the density and temperature of the gas,we know that low density and high temperature gas appears near the exits of the DRPG,while high density and low temperature gas appears at the wall near the breech.The simulation of the 2D-IB-SPETCG model is an effective way to investigate the interior ballistics process of the ETC launch.The 2D-IB-SPETC model can be used for prediction and improvement of experiments.
基金supported by the Research Directorate of the University of Cuenca(DIUC)
文摘Hydraulic models for the generation of flood inundation maps are not commonly applied in mountain river basins because of the difficulty in modeling the hydraulic behavior and the complex topography. This paper presents a comparative analysis of the performance of four twodimensional hydraulic models (HEC-RAS 2D, Iber 2D, Flood Modeller 2D, and PCSWMM 2D) with respect to the generation of flood inundation maps. The study area covers a 5-km reach of the Santa B-arbara River located in the Ecuadorian Andes, at 2330 masl, in Gualaceo. The model's performance was evaluated based on the water surface elevation and flood extent, in terms of the mean absolute difference and measure of fit. The analysis revealed that, for a given case, Iber 2D has the best performance in simulating the water level and inundation for flood events with 20- and 50-year return periods, respectively, followed by Flood Modeller 2D, HEC-RAS 2D, and PCSWMM 2D in terms of their performance. Grid resolution, the way in which hydraulic structures are mimicked, the model code, and the default value of the parameters are considered the main sources of prediction uncertainty.
基金Sponsored by the National Natural Science Foundation of China(Grant No.51208160)the Natural Science Foundation of Heilongjiang Province(Grant No.QC2012C056)
文摘The features of a quasi-two-dimensional( quasi-2D) model for simulating two-phase water hammer flows with vaporous cavity in a pipe are investigated. The quasi-2D model with discrete vaporous cavity in the pipe is proposed in this paper. This model uses the quasi-2D model for pure liquid zone and one-dimensional( 1D) discrete vapor cavity model for vaporous cavity zone. The quasi-2D model solves two-dimensional equations for both axial and radial velocities and 1D equations for both pressure head and discharge by the method of characteristics. The 1D discrete vapor cavity model is used to simulate the vaporous cavity occurred when the pressure in the local pipe is lower than the vapor pressure of the liquid. The proposed model is used to simulate two-phase water flows caused by the rapid downstream valve closure in a reservoir-pipe-valve system.The results obtained by the proposed model are compared with those by the corresponding 1D model and the experimental ones provided by the literature,respectively. The comparison shows that the maximum pressure heads simulated by the proposed model are more accurate than those by the corresponding 1D model.
基金This research was funded by the National Natural Science Foundation of China(52079106,52009104,51609199)the National Key Research and Development Program of China(2016YFC0402704).
文摘Erosion is an important issue in soil science and is related to many environmental problems,such as soil erosion and sediment transport.Establishing a simulation model suitable for soil erosion prediction is of great significance not only to accurately predict the process of soil separation by runoff,but also improve the physical model of soil erosion.In this study,we develop a graphic processing unit(GPU)-based numerical model that combines two-dimensional(2D)hydrodynamic and Green-Ampt(G-A)infiltration modelling to simulate soil erosion.A Godunov-type scheme on a uniform and structured square grid is then generated to solve the relevant shallow water equations(SWEs).The highlight of this study is the use of GPU-based acceleration technology to enable numerical models to simulate slope and watershed erosion in an efficient and high-resolution manner.The results show that the hydrodynamic model performs well in simulating soil erosion process.Soil erosion is studied by conducting calculation verification at the slope and basin scales.The first case involves simulating soil erosion process of a slope surface under indoor artificial rainfall conditions from 0 to 1000 s,and there is a good agreement between the simulated values and the measured values for the runoff velocity.The second case is a river basin experiment(Coquet River Basin)that involves watershed erosion.Simulations of the erosion depth change and erosion cumulative amount of the basin during a period of 1-40 h show an elevation difference of erosion at 0.5-3.0 m,especially during the period of 20-30 h.Nine cross sections in the basin are selected for simulation and the results reveal that the depth of erosion change value ranges from-0.86 to-2.79 m and the depth of deposition change value varies from 0.38 to 1.02 m.The findings indicate that the developed GPU-based hydrogeomorphological model can reproduce soil erosion processes.These results are valuable for rainfall runoff and soil erosion predictions on rilled hillslopes and river basins.
文摘A global two-dimensional zonally averaged chemistry model is developed to study the chemi-cal composition of atmosphere. The region of the model is from 90°S to 90°N and from the ground to the altitude of 20 km with a resolution of 5° x 1 km. The wind field is residual circulation calcu-lated from diabatic rate. 34 species and 104 chemical and photochemical reactions are considered in the model. The sources of CH4, CO and NOx, which are divided into seasonal sources and non-seasonal sources, are parameterized as a function of latitude and time. The chemical composi-tion of atmosphere was simulated with emission level of CH4, CO and NOx in 1990. The results are compared with observations and other model results, showing that the model is successful to simu-late the atmospheric chemical composition and distribution of CH4. Key words Global two-dimensional chemistry model - Atmospheric composition - Emission This work was supported by the State Key Program for basic research “ Climate Dynamics and Cli-mate Prediction Theory” (Pandeng-yu-21).The authors would like to express their thanks to the National Oceanic and Atmospheric Administration (NOAA), Climate Monitoring and Diagnostics Laboratory (CMDL), Carbon Cycle Group for providing the observational data of CO and CH4.
基金supported by the National Natural Science Foundation of China(Grant No.50579030)
文摘River ice is a natural phenomenon in cold regions, influenced by meteorology, geomorphology, and hydraulic conditions. River ice processes involve complex interactions between hydrodynamic, mechanical, and thermal processes, and they are also influenced by weather and hydrologic conditions. Because natural rivers are serpentine, with bends, narrows, and straight reaches, the commonly-used one-dimensional river ice models and two-dimensional models based on the rectangular Cartesian coordinates are incapable of simulating the physical phenomena accurately. In order to accurately simulate the complicated river geometry and overcome the difficulties of numerical simulation resulting from both complex boundaries and differences between length and width scales, a two-dimensional river ice numerical model based on a boundary-fitted coordinate transformation method was developed. The presented model considers the influence of the frazil ice accumulation under ice cover and the shape of the leading edge of ice cover during the freezing process. The model is capable of determining the velocity field, the distribution of water temperature, the concentration distribution of frazil ice, the transport of floating ice, the progression, stability, and thawing of ice cover, and the transport, accumulation, and erosion of ice under ice cover. A MacCormack scheme was used to solve the equations numerically. The model was validated with field observations from the Hequ Reach of the Yellow River. Comparison of simulation results with field data indicates that the model is capable of simulating the river ice process with high accuracy.
基金the National Natural Science Fund(11661058,11761053)Natural Science Fund of Inner Mongolia Autonomous Region(2016MS0102,2017MS0107)+1 种基金Program for Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region(NJYT-17-A07)National Undergraduate Innovative Training Project of Inner Mongolia University(201710126026).
文摘In this article,a high-order scheme,which is formulated by combining the quadratic finite element method in space with a second-order time discrete scheme,is developed for looking for the numerical solution of a two-dimensional nonlinear time fractional thermal diffusion model.The time Caputo fractional derivative is approximated by using the L2-1formula,the first-order derivative and nonlinear term are discretized by some second-order approximation formulas,and the quadratic finite element is used to approximate the spatial direction.The error accuracy O(h3+t2)is obtained,which is verified by the numerical results.
文摘On the conditions of low-resolution radar, a parametric model for two-dimensional radar target is described here according to the theory of electromagnetic scattering and the geometrical theory of diffraction. A high resolution estimation algorithm to extract the model parameters is also developed by building the relation of the scattering model and Prony model. The analysis of Cramer-Rao bound and simulation show that the method here has better statistical performance. The simulated analysis also indicates that the accurate extraction of the diffraction coefficient of scattering center is restricted by signal to noise ratio, radar center frequency and radar bandwidth.
基金supported by the National Natural Science Foundation of China(Grants Nos.40875025,40875030,and 40775033)the Shanghai Natural Science Foundation of China(Grant No.08ZR1422900)
文摘Water vapor, cloud, and surface rainfall budgets associated with the landfall of Typhoon Krosa on 6-8 October 2007 are analyzed based on a two-dimensional cloud-resolving model simulation. The model is integrated with imposed zonally-uniform vertical velocity, zonal wind, horizontal temperature, and vapor advection from NCEP/Global Data Assimilation System （GDAS） data. The simulation data that are validated with observations are examined to study physical causes associated with surface rainfall processes during the landfall. The time- and domain-mean analysis shows that when Krosa approached the eastern coast of China on 6 October, the water vapor convergence over land caused a local atmospheric moistening and a net condensation that further produced surface rainfall and an increase of cloud hydrometeor concentration. Meanwhile, latent heating was balanced by advective cooling and a local atmospheric warming. One day later, the enhancement of net condensation led to an increase of surface rainfall and a local atmospheric drying, while the water vapor convergence weakened as a result of the landfall-induced deprivation of water vapor flux. At the same time, the latent heating is mainly compensated the advective cooling. Further weakening of vapor convergence on 8 October enhanced the local atmospheric drying while the net condensation and associated surface rainfall was maintained. The latent heating is balanced by advective cooling and a local atmospheric cooling.
基金Project supported by the National Natural Science Foundation of China (No. 40375038)the National Basic Research Program of China (No. 2002CB410802, 2002CB410801).
文摘A modified two-dimensional Eulerian air quality model was used to simulate both the gaseous and particulate pollutant concentrations during October 21-24, 2004 in the Pearl River Delta （PRD） region, China. The most significant improvement to the model is the added capability to predict the secondary organic aerosols （SOA） concentrations because of the inclusion of the SOA formation chemistry. The meteorological input data were prepared using the CALMET meteorological model. The concentrations of aerosol-bound species such as NO3^-, NH4^＋, SO4^2-, and SOA were calculated in the fine particle size range （〈2.5 μm）. The results of the two-dimensional model were compared to the measurements at the ground level during the PRD Intensive Monitoring Campaign （IMC）. Overall, there were good agreements between the measured and modeled concentrations of inorganic aerosol components and O3. Both the measured and the modeled results indicated that the maximum hourly O3 concentrations exceeded the China National Air Quality Standard. The predicted 24-h average SOA concentrations were in reasonable agreement with those predicted by the method of minimum OC/EC ratio.
