A one -dimensional time-dependent photochemical model is used to simulate the influence of ion-produced NOx and HOx radicals on the Antarctic ozone depletion in polar night and polar spring at a latitude of 73 degrees...A one -dimensional time-dependent photochemical model is used to simulate the influence of ion-produced NOx and HOx radicals on the Antarctic ozone depletion in polar night and polar spring at a latitude of 73 degrees south.Vertical transport and nitrogen-oxygen (NOx). hydrogen-oxygen (HOx) production by ionic reactions have been introduced into the model.NOx and HOx produced by precipitating ions are transported into the lower stratosphere by vertical motion and have some effects in the development of the Antarctic ozone depletion.From winter through spring the calculated ozone column decreases to 269.4 DU. However, this value is significantly higher than the total ozone observed at several Antarctic ozone stations.展开更多
Achieving the detachment of divertor can help to alleviate excessive heat load and sputtering problems on the target plates,thereby extending the lifetime of divertor components for fusion devices.In order to provide ...Achieving the detachment of divertor can help to alleviate excessive heat load and sputtering problems on the target plates,thereby extending the lifetime of divertor components for fusion devices.In order to provide a fast but relatively reliable prediction of plasma parameters along the flux tube for future device design,a one-dimensional(1D)modeling code for the operating point of impurity seeded detached divertor is developed based on Python language,which is a fluid model based on previous work(Plasma Phys.Control.Fusion 58045013(2016)).The experimental observation of the onset of divertor detachment by neon(Ne)and argon(Ar)seeding in EAST is well reproduced by using the 1D modeling code.The comparison between the 1D modeling and two-dimensional(2D)simulation by the SOLPS-ITER code for CFETR detachment operation with Ne and Ar seeding also shows that they are in good agreement.We also predict the radiative power loss and corresponding impurity concentration requirement for achieving divertor detachment via different impurity seeding under high heating power conditions in EAST and CFETR phase II by using the 1D model.Based on the predictions,the optimized parameter space for divertor detachment operation on EAST and CFETR is also determined.Such a simple but reliable 1D model can provide a reasonable parameter input for a detailed and accurate analysis by 2D or three-dimensional(3D)modeling tools through rapid parameter scanning.展开更多
This study proposes a quasi-one-dimensional model to predict the chemical nonequilibrium flow along the stagnation streamline of hypersonic flow past a blunt body. The model solves reduced equations along the stagnati...This study proposes a quasi-one-dimensional model to predict the chemical nonequilibrium flow along the stagnation streamline of hypersonic flow past a blunt body. The model solves reduced equations along the stagnation streamline and predicts nearly identical results as the numerical solution of the full-field Navier-Stokes equations. The high efficiency of this model makes it useful to investigate the overall quantitative behavior of related physical-chemical phenomena. In this paper two important properties of hypersonic flow, shock stand-off distance and oxygen dissociation, are studied using the quasi-one-dimensional model with the ideal dissociating gas model. It is found that the shock stand-off distance is affected by both chemical and thermal non-equilibrium.The shock stand-off distance will increase when the flow conditions are changed from equilibrium to non-equilibrium, because the average density of the shock-compressed gas will decrease as a result of the increase in translational energy. For oxygen dissociation, the maximum value of its dissociation degree along the stagnation line varies with the flight altitude. It is increased at first and decreased thereafter with the altitude, which is due to the combination effect of the equilibrium shift and chemical non-equilibrium relaxation. The overall variation of the maximum dissociation is then plotted in the speed and altitude coordinates as a reference for engineering application.展开更多
We study the thermodynamic properties of the classical one-dimensional generalized nonlinear Klein-Gordon lattice model(n≥2)by using the cluster variation method with linear response theory.The results of this method...We study the thermodynamic properties of the classical one-dimensional generalized nonlinear Klein-Gordon lattice model(n≥2)by using the cluster variation method with linear response theory.The results of this method are exact in the thermodynamic limit.We present the single-site reduced densityρ^((1))(z),averages such as(z^(2)),<|z^(n)|>,and<(z_(1)-z_(2))^(2)>,the specific heat C_(v),and the static correlation functions.We analyze the scaling behavior of these quantities and obtain the exact scaling powers at the low and high temperatures.Using these results,we gauge the accuracy of the projective truncation approximation for theφ^(4)lattice model.展开更多
We propose an eigen microstate approach(EMA)for analyzing quantum phase transitions in quantum many-body systems,introducing a novel framework that does not require prior knowledge of an order parameter.Using the tran...We propose an eigen microstate approach(EMA)for analyzing quantum phase transitions in quantum many-body systems,introducing a novel framework that does not require prior knowledge of an order parameter.Using the transversefield Ising model(TFIM)as a case study,we demonstrate the effectiveness of EMA by identifying key features of the phase transition through the scaling behavior of eigenvalues and the structure of associated eigen microstates.Our results reveal substantial changes in the ground state of the TFIM as it undergoes a phase transition,as reflected in the behavior of specific componentsξ_(i)^((k))within the eigen microstates.This method is expected to be applicable to other quantum systems where predefining an order parameter is challenging.展开更多
We propose a new section-averaged one-dimensional model for blood flows in deformable arteries.The model is derived from the three-dimensional Navier-Stokes equations,written in cylindrical coordinates,under the“thin...We propose a new section-averaged one-dimensional model for blood flows in deformable arteries.The model is derived from the three-dimensional Navier-Stokes equations,written in cylindrical coordinates,under the“thin-artery”assumption(similar to the“shallow-water”assumption for free surface models).The blood flow/artery interaction is taken into account through suitable boundary conditions.The obtained equations enter the scope of the nonlinear convection-diffusion problems.We show that the resulting model is energetically consistent.The proposed model extends most extant models by adding more scope,depending on an additional viscous term.We compare both models computationally based on an Incomplete Interior Penalty Galerkin(IIPG)method for the parabolic part,and on a Runge Kutta Discontinuous Galerkin(RKDG)method for the hyperbolic part.The time discretization explicit/implicit is based on the well-known Additive Runge-Kutta(ARK)method.Moreover,through a suitable change of variables,by construction,we show that the numerical scheme is well-balanced,i.e.,it preserves exactly still-steady state solutions.To end,we numerically investigate its efficiency through several test cases with a confrontation to an exact solution.展开更多
A Discrete Element Method （DEM） model is developed to study the particle break- age effect on the one-dimensional compression behavior of silica sands. The ＇maximum tensile stress＇ breakage criterion considering m...A Discrete Element Method （DEM） model is developed to study the particle break- age effect on the one-dimensional compression behavior of silica sands. The ＇maximum tensile stress＇ breakage criterion considering multiple contacts is adopted to simulate the crushing of circular particles in the DEM. The model is compared with published experimental results. Com- parison between the compression curves obtained from the numerical and experimental results shows that the proposed method is very effective in studying the compression behavior of silica sands considering particle breakage. The evolution of compression curves at different stress levels is extensively studied using contact force distribution, variation of contact number and particle size distribution curve with loading. It is found that particle breakage has great impact on com- pression behavior of sand, particularly after the yield stress is reached and particle breakage starts. The crushing probability of particles is found to be macroscopically affected by stress level and particle size distribution curve, and microscopically related to the evolutions of contact force and coordination number. Once the soil becomes well-graded and the average coordination number is greater than 4 in two-dimension, the crushing probability of parent particles can reduce by up to 5/6. It is found that the average contact force does not always increase with loading, but increases to a peak value then decreases once the soil becomes more well-graded. It is found through the loading rate sensitivity analysis that the compression behavior of sand samples in the DEM is also affected by the loading rate. Higher yield stresses are obtained at higher loading rates.展开更多
In this paper, on the basis of the heat conduction equation without consideration of the advection and turbulence effects, one-dimensional model for describing surface sea temperature ( T1), bottom sea temperature ( T...In this paper, on the basis of the heat conduction equation without consideration of the advection and turbulence effects, one-dimensional model for describing surface sea temperature ( T1), bottom sea temperature ( Tt ) and the thickness of the upper homogeneous layer ( h ) is developed in terms of the dimensionless temperature θT and depth η and self-simulation function θT - f(η) of vertical temperature profile by means of historical temperature data.The results of trial prediction with our one-dimensional model on T, Th, h , the thickness and gradient of thermocline are satisfactory to some extent.展开更多
Based on a set of fully nonlinear Boussinesq equations up to the order of O（μ^2, ε^3μ^2） （where ε is the ratio of wave amplitude to water depth and ,μ is the ratio of water depth to wave length） a numerical w...Based on a set of fully nonlinear Boussinesq equations up to the order of O（μ^2, ε^3μ^2） （where ε is the ratio of wave amplitude to water depth and ,μ is the ratio of water depth to wave length） a numerical wave model is formulated. The model＇s linear dispersion is acceptably accurate to μ ≌ 1.0, which is confirmed by comparisons between the simulat- ed and measured time series of the regular waves propagating on a submerged bar. The moving shoreline is treated numer- ically by replacing the solid beach with a permeable beach. Run-up of nonbreaking waves is verified against the analytical solution for nonlinear shallow water waves. The inclusion of wave breaking is fulfilled by introducing an eddy term in the momentum equation to serve as the breaking wave force term to dissipate wave energy in the surf zone. The model is applied to cross-shore motions of regular waves including various types of breaking on plane sloping beaches. Comparisons of the model test results comprising spatial distribution of wave height and mean water level with experimental data are presented.展开更多
A one-dimensional discrete Boltzmann model for detonation simulation is presented. Instead of numerical solving Navier-Stokes equations, this model obtains the information of flow field through numerical solving speci...A one-dimensional discrete Boltzmann model for detonation simulation is presented. Instead of numerical solving Navier-Stokes equations, this model obtains the information of flow field through numerical solving specially discretized Boltzmann equation. Several classical benchmarks including Sod shock wave tube, Colella explosion problem,and one-dimensional self-sustainable stable detonation are simulated to validate the new model. Based on the new model,the influence of negative temperature coefficient of reaction rate on detonation is further investigated. It is found that an abnormal detonation with two wave heads periodically appears under negative temperature coefficient condition.The causes of the abnormal detonation are analyzed. One typical cycle of the periodic abnormal detonation and its development process are discussed.展开更多
A new method for driving a One-Dimensional Stratiform Cold （1DSC） cloud model with Weather Research and Fore casting （WRF） model outputs was developed by conducting numerical experiments for a typical large-scale ...A new method for driving a One-Dimensional Stratiform Cold （1DSC） cloud model with Weather Research and Fore casting （WRF） model outputs was developed by conducting numerical experiments for a typical large-scale stratiform rainfall event that took place on 4-5 July 2004 in Changchun, China. Sensitivity test results suggested that, with hydrometeor pro files extracted from the WRF outputs as the initial input, and with continuous updating of soundings and vertical velocities （including downdraft） derived from the WRF model, the new WRF-driven 1DSC modeling system （WRF-1DSC） was able to successfully reproduce both the generation and dissipation processes of the precipitation event. The simulated rainfall intensity showed a time-lag behind that observed, which could have been caused by simulation errors of soundings, vertical velocities and hydrometeor profiles in the WRF output. Taking into consideration the simulated and observed movement path of the precipitation system, a nearby grid point was found to possess more accurate environmental fields in terms of their similarity to those observed in Changchun Station. Using profiles from this nearby grid point, WRF-1DSC was able to repro duce a realistic precipitation pattern. This study demonstrates that 1D cloud-seeding models do indeed have the potential to predict realistic precipitation patterns when properly driven by accurate atmospheric profiles derived from a regional short range forecasting system, This opens a novel and important approach to developing an ensemble-based rain enhancement prediction and operation system under a probabilistic framework concept.展开更多
In sub nanometer carbon nanotubes,water exhibits unique dynamic characteristics,and in the high-frequency region of the infrared spectrum,where the stretching vibrations of the internal oxygen-hydrogen(O-H)bonds are c...In sub nanometer carbon nanotubes,water exhibits unique dynamic characteristics,and in the high-frequency region of the infrared spectrum,where the stretching vibrations of the internal oxygen-hydrogen(O-H)bonds are closely related to the hydrogen bonds(H-bonds)network between water molecules.Therefore,it is crucial to analyze the relationship between these two aspects.In this paper,the infrared spectrum and motion characteristics of the stretching vibrations of the O-H bonds in one-dimensional confined water(1DCW)and bulk water(BW)in(6,6)single-walled carbon nanotubes(SWNT)are studied by molecular dynamics simulations.The results show that the stretching vibrations of the two O-H bonds in 1DCW exhibit different frequencies in the infrared spectrum,while the O-H bonds in BW display two identical main frequency peaks.Further analysis using the spring oscillator model reveals that the difference in the stretching amplitude of the O-H bonds is the main factor causing the change in vibration frequency,where an increase in stretching amplitude leads to a decrease in spring stiffness and,consequently,a lower vibration frequency.A more in-depth study found that the interaction of H-bonds between water molecules is the fundamental cause of the increased stretching amplitude and decreased vibration frequency of the O-H bonds.Finally,by analyzing the motion trajectory of the H atoms,the dynamic differences between 1DCW and BW are clearly revealed.These findings provide a new perspective for understanding the behavior of water molecules at the nanoscale and are of significant importance in advancing the development of infrared spectroscopy detection technology.展开更多
A one-dimensional stationary nonisentropic hydrodynamic model for semiconductor devices with non-constant lattice temperature is studied. This model consists of the equations for the electron density, the electron cur...A one-dimensional stationary nonisentropic hydrodynamic model for semiconductor devices with non-constant lattice temperature is studied. This model consists of the equations for the electron density, the electron current density and electron temperature, coupled with the Poisson equation of the electrostatic potential in a bounded interval supplemented with proper boundary conditions. The existence and uniqueness of a strong subsonic steady-state solution with positive particle density and positive temperature is established. The proof is based on the fixed-point arguments, the Stampacchia truncation methods, and the basic energy estimates.展开更多
Different phenomenological equations based on plasticity, primary creep （as a viscoplastic mechanism）, secondary creep （as another viscoplastic mechanism） and different combinations of these equations are presente...Different phenomenological equations based on plasticity, primary creep （as a viscoplastic mechanism）, secondary creep （as another viscoplastic mechanism） and different combinations of these equations are presented and used to describe the material inelastic deformation in uniaxial test. Agreement of the models with experimental results and with the theoretical concepts and physical realities is the criterion of choosing the most appropriate formulation for uniaxial test. A model is thus proposed in which plastic deformation, primary creep and secondary creep contribute to the inelastic deformation. However, it is believed that the hardening parameter is composed of plastic and primary creep parts. Accordingly, the axial plastic strain in a uniaxial test may no longer be considered as the hardening parameter. Therefore, a proportionality concept is proposed to calculate the plastic contribution of deformation.展开更多
This paper is concerned with a singular limit for the one-dimensional compress- ible radiation hydrodynamics model. The singular limit we consider corresponds to the physical problem of letting the Bouguer number infi...This paper is concerned with a singular limit for the one-dimensional compress- ible radiation hydrodynamics model. The singular limit we consider corresponds to the physical problem of letting the Bouguer number infinite while keeping the Boltzmann number constant. In the case when the corresponding Euler system admits a contact discontinuity wave, Wang and Xie （2011） [12] recently verified this singular limit and proved that the solution of the compressible radiation hydrodynamics model converges to the strong contact 1 discontinuity wave in the L∞-norm away from the discontinuity line at a rate of ε1/4, as the reciprocal of the Bouguer number tends to zero. In this paper, Wang and Xie＇s convergence rate is improved to ε7/8 by introducing a new a priori assumption and some refined energy estimates. Moreover, it is shown that the radiation flux q tends to zero in the L∞-norm away from the discontinuity line, at a convergence rate as the reciprocal of the Bouguer number tends to zero.展开更多
The observation data for 5 d at a station in the South China Sea is presented. After brief anMysis of the wind speed, air temperature from the ship-borne meteorological instruments and temperature and salinity profile...The observation data for 5 d at a station in the South China Sea is presented. After brief anMysis of the wind speed, air temperature from the ship-borne meteorological instruments and temperature and salinity profiles from the CTD （conductivity, temperature, depth recorder） data, the authors find that the CTD casts are too sparse for us to understand the diurnal evolution of the thermal structure in the upper ocean. A one-dimensional （1D） numericM code based on Mellor-Yamada turbulence closure model is used to reconstruct the upper ocean thermal structure, utilizing the atmospheric forcing data from ship-borne weather station. The simulation results show good agreement with the observational data; the significance of breaking waves is also briefly discussed. The evolution of turbulence kinetic energy （TKE） and the contribution from shear production and buoy- ancy production are discussed respectively. Finally, several possible factors which might influence the numerical results are briefly analyzed.展开更多
A modified one-dimensional transient hygrothermal model for multilayer wall was proposed using air humidity ratio and temperature as the driving potentials.The solution for the governing equations was obtained numeric...A modified one-dimensional transient hygrothermal model for multilayer wall was proposed using air humidity ratio and temperature as the driving potentials.The solution for the governing equations was obtained numerically by implementing the finite-difference scheme.To evaluate the accuracy of the model,a test system was built up to measure relative humidity and temperature within a porous wall and compare with the prediction of the model.The prediction results have good agreement with the experimental results.For the interface close to indoor side,the maximum deviation of temperature between calculated and test data is 1.87 K,and the average deviation is 0.95 K;the maximum deviation of relative humidity is 11.4%,and the average deviation is 5.7%.For the interface close to outdoor side,the maximum deviation of temperature between prediction and measurement is 1.78 K,and the average deviation is 1.1 K;the maximum deviation of relative humidity is 9.9%,and the average deviation is 4.2%.展开更多
A stochastic local limited one-dimensional rice-pile model is numerically investigated. The distributions for avalanche sizes have a clear power-law behavior and it displays a simple finite size scaling. We obtain the...A stochastic local limited one-dimensional rice-pile model is numerically investigated. The distributions for avalanche sizes have a clear power-law behavior and it displays a simple finite size scaling. We obtain the avalanche exponents Ts= 1.54±0.10,βs = 2.17±0.10 and TT = 1.80±0.10, βT =1.46 ± 0.10. This self-organized critical model belongs to the same universality class with the Oslo rice-pile model studied by K. Christensen et al. [Phys. Rev. Lett. 77 （1996） 107], a rice-pile model studied by L.A.N. Amaral et al. [Phys. Rev. E 54 （1996） 4512], and a simple deterministic self-organized critical model studied by M.S. Vieira [Phys. Rev. E 61 （2000） 6056].展开更多
This article aims to investigate the transient behavior of a planar direct internal reforming solid oxide fuel cell （DIR-SOFC） comprehensively. A one-dimensional dynamic model of a planar D1R-SOFC is first developed...This article aims to investigate the transient behavior of a planar direct internal reforming solid oxide fuel cell （DIR-SOFC） comprehensively. A one-dimensional dynamic model of a planar D1R-SOFC is first developed based on mass and energy balances, and electrochemical principles. Further, a solution strategy is presented to solve the model, and the International Energy Agency （IEA） benchmark test is used to validate the model. Then, through model-based simulations, the steady-state performance of a co-flow planar DIR-SOFC under specified initial operating conditions and its dynamic response to introduced operating parameter disturbances are studied. The dynamic responses of important SOFC variables, such as cell temperature, current density, and cell voltage are all investigated when the SOFC is subjected to the step-changes in various operating parameters including both the load current and the inlet fuel and air flow rates. The results indicate that the rapid dynamics of the current density and the cell voltage are mainly influenced by the gas composition, particularly the H2 molar fraction in anode gas channels, while their slow dynamics are both dominated by the SOLID （including the PEN and interconnects） temperature. As the load current increases, the SOLID temperature and the maximum SOLID temperature gradient both increase, and thereby, the cell breakdown is apt to occur because of excessive thermal stresses. Changing the inlet fuel flow rate might lead to the change in the anode gas composition and the consequent change in the current density distribution and cell voltage. The inlet air flow rate has a great impact on the cell temperature distribution along the cell, and thus, is a suitable manipulated variable to control the cell temperature.展开更多
In this paper, the coupling schemes of atmosphere-ocean climate models are discussed with one-dimensional advection equations. The convergence and stability for synchronous and asynchronous schemes are demonstrated an...In this paper, the coupling schemes of atmosphere-ocean climate models are discussed with one-dimensional advection equations. The convergence and stability for synchronous and asynchronous schemes are demonstrated and compared.Conclusions inferred from the analysis are given below. The synchronous scheme as well as the asynchronous-implicit scheme in this model are stable for arbitrary integrating time intervals. The asynchronous explicit scheme is unstable under certain conditions, which depend upon advection velocities and heat exchange parameters in the atmosphere and oceans. With both synchronous and asynchronous stable schemes the discrete solutions converge to their unique exact ones. Advections in the atmosphere and ocean accelerate the rate of convergence of the asynchronous-implicit scheme. It is suggusted that the asynchronous-implicit coupling scheme is a stable and efficient method for most climatic simulations.展开更多
文摘A one -dimensional time-dependent photochemical model is used to simulate the influence of ion-produced NOx and HOx radicals on the Antarctic ozone depletion in polar night and polar spring at a latitude of 73 degrees south.Vertical transport and nitrogen-oxygen (NOx). hydrogen-oxygen (HOx) production by ionic reactions have been introduced into the model.NOx and HOx produced by precipitating ions are transported into the lower stratosphere by vertical motion and have some effects in the development of the Antarctic ozone depletion.From winter through spring the calculated ozone column decreases to 269.4 DU. However, this value is significantly higher than the total ozone observed at several Antarctic ozone stations.
基金Project supported by the National Key Research and Development Program of China (Grant No.2022YFE03030001)the National Natural Science Foundation of China (Grant No.12075283)。
文摘Achieving the detachment of divertor can help to alleviate excessive heat load and sputtering problems on the target plates,thereby extending the lifetime of divertor components for fusion devices.In order to provide a fast but relatively reliable prediction of plasma parameters along the flux tube for future device design,a one-dimensional(1D)modeling code for the operating point of impurity seeded detached divertor is developed based on Python language,which is a fluid model based on previous work(Plasma Phys.Control.Fusion 58045013(2016)).The experimental observation of the onset of divertor detachment by neon(Ne)and argon(Ar)seeding in EAST is well reproduced by using the 1D modeling code.The comparison between the 1D modeling and two-dimensional(2D)simulation by the SOLPS-ITER code for CFETR detachment operation with Ne and Ar seeding also shows that they are in good agreement.We also predict the radiative power loss and corresponding impurity concentration requirement for achieving divertor detachment via different impurity seeding under high heating power conditions in EAST and CFETR phase II by using the 1D model.Based on the predictions,the optimized parameter space for divertor detachment operation on EAST and CFETR is also determined.Such a simple but reliable 1D model can provide a reasonable parameter input for a detailed and accurate analysis by 2D or three-dimensional(3D)modeling tools through rapid parameter scanning.
基金supported by the National Natural Science Foundation of China (Nos. 1372325 and 91116013)
文摘This study proposes a quasi-one-dimensional model to predict the chemical nonequilibrium flow along the stagnation streamline of hypersonic flow past a blunt body. The model solves reduced equations along the stagnation streamline and predicts nearly identical results as the numerical solution of the full-field Navier-Stokes equations. The high efficiency of this model makes it useful to investigate the overall quantitative behavior of related physical-chemical phenomena. In this paper two important properties of hypersonic flow, shock stand-off distance and oxygen dissociation, are studied using the quasi-one-dimensional model with the ideal dissociating gas model. It is found that the shock stand-off distance is affected by both chemical and thermal non-equilibrium.The shock stand-off distance will increase when the flow conditions are changed from equilibrium to non-equilibrium, because the average density of the shock-compressed gas will decrease as a result of the increase in translational energy. For oxygen dissociation, the maximum value of its dissociation degree along the stagnation line varies with the flight altitude. It is increased at first and decreased thereafter with the altitude, which is due to the combination effect of the equilibrium shift and chemical non-equilibrium relaxation. The overall variation of the maximum dissociation is then plotted in the speed and altitude coordinates as a reference for engineering application.
基金supported by the National Natural Science Foundation of China(Grant No.11974420).
文摘We study the thermodynamic properties of the classical one-dimensional generalized nonlinear Klein-Gordon lattice model(n≥2)by using the cluster variation method with linear response theory.The results of this method are exact in the thermodynamic limit.We present the single-site reduced densityρ^((1))(z),averages such as(z^(2)),<|z^(n)|>,and<(z_(1)-z_(2))^(2)>,the specific heat C_(v),and the static correlation functions.We analyze the scaling behavior of these quantities and obtain the exact scaling powers at the low and high temperatures.Using these results,we gauge the accuracy of the projective truncation approximation for theφ^(4)lattice model.
基金supported by the National Natural Science Foundation of China(Grant Nos.12475033,12135003,12174194,and 12405032)the National Key Research and Development Program of China(Grant No.2023YFE0109000)+1 种基金supported by the Fundamental Research Funds for the Central Universitiessupport from the China Postdoctoral Science Foundation(Grant No.2023M730299).
文摘We propose an eigen microstate approach(EMA)for analyzing quantum phase transitions in quantum many-body systems,introducing a novel framework that does not require prior knowledge of an order parameter.Using the transversefield Ising model(TFIM)as a case study,we demonstrate the effectiveness of EMA by identifying key features of the phase transition through the scaling behavior of eigenvalues and the structure of associated eigen microstates.Our results reveal substantial changes in the ground state of the TFIM as it undergoes a phase transition,as reflected in the behavior of specific componentsξ_(i)^((k))within the eigen microstates.This method is expected to be applicable to other quantum systems where predefining an order parameter is challenging.
基金supported by the ADEN-MED project(Adaptability to Extreme Events and Natural Risks-Application to the Mediterranean and Djibouti)funded by the Région Sud Provence-Alpes-Côte d’Azur under the AAP MEDCLIMAT“Natural Risks and Food Sovereignty”.
文摘We propose a new section-averaged one-dimensional model for blood flows in deformable arteries.The model is derived from the three-dimensional Navier-Stokes equations,written in cylindrical coordinates,under the“thin-artery”assumption(similar to the“shallow-water”assumption for free surface models).The blood flow/artery interaction is taken into account through suitable boundary conditions.The obtained equations enter the scope of the nonlinear convection-diffusion problems.We show that the resulting model is energetically consistent.The proposed model extends most extant models by adding more scope,depending on an additional viscous term.We compare both models computationally based on an Incomplete Interior Penalty Galerkin(IIPG)method for the parabolic part,and on a Runge Kutta Discontinuous Galerkin(RKDG)method for the hyperbolic part.The time discretization explicit/implicit is based on the well-known Additive Runge-Kutta(ARK)method.Moreover,through a suitable change of variables,by construction,we show that the numerical scheme is well-balanced,i.e.,it preserves exactly still-steady state solutions.To end,we numerically investigate its efficiency through several test cases with a confrontation to an exact solution.
基金Project supported by the National Natural Science Foundation of China(Nos.50909057,51208294 and 41372319)the Innovation Program of Shanghai Municipal Education Commission(No.15ZZ081)the Innovation Project of Shanghai Postgraduate Education(No.20131129)
文摘A Discrete Element Method （DEM） model is developed to study the particle break- age effect on the one-dimensional compression behavior of silica sands. The ＇maximum tensile stress＇ breakage criterion considering multiple contacts is adopted to simulate the crushing of circular particles in the DEM. The model is compared with published experimental results. Com- parison between the compression curves obtained from the numerical and experimental results shows that the proposed method is very effective in studying the compression behavior of silica sands considering particle breakage. The evolution of compression curves at different stress levels is extensively studied using contact force distribution, variation of contact number and particle size distribution curve with loading. It is found that particle breakage has great impact on com- pression behavior of sand, particularly after the yield stress is reached and particle breakage starts. The crushing probability of particles is found to be macroscopically affected by stress level and particle size distribution curve, and microscopically related to the evolutions of contact force and coordination number. Once the soil becomes well-graded and the average coordination number is greater than 4 in two-dimension, the crushing probability of parent particles can reduce by up to 5/6. It is found that the average contact force does not always increase with loading, but increases to a peak value then decreases once the soil becomes more well-graded. It is found through the loading rate sensitivity analysis that the compression behavior of sand samples in the DEM is also affected by the loading rate. Higher yield stresses are obtained at higher loading rates.
文摘In this paper, on the basis of the heat conduction equation without consideration of the advection and turbulence effects, one-dimensional model for describing surface sea temperature ( T1), bottom sea temperature ( Tt ) and the thickness of the upper homogeneous layer ( h ) is developed in terms of the dimensionless temperature θT and depth η and self-simulation function θT - f(η) of vertical temperature profile by means of historical temperature data.The results of trial prediction with our one-dimensional model on T, Th, h , the thickness and gradient of thermocline are satisfactory to some extent.
基金This work was financially supported by the National Natural Science Foundation of China (Grant No.50679010)
文摘Based on a set of fully nonlinear Boussinesq equations up to the order of O（μ^2, ε^3μ^2） （where ε is the ratio of wave amplitude to water depth and ,μ is the ratio of water depth to wave length） a numerical wave model is formulated. The model＇s linear dispersion is acceptably accurate to μ ≌ 1.0, which is confirmed by comparisons between the simulat- ed and measured time series of the regular waves propagating on a submerged bar. The moving shoreline is treated numer- ically by replacing the solid beach with a permeable beach. Run-up of nonbreaking waves is verified against the analytical solution for nonlinear shallow water waves. The inclusion of wave breaking is fulfilled by introducing an eddy term in the momentum equation to serve as the breaking wave force term to dissipate wave energy in the surf zone. The model is applied to cross-shore motions of regular waves including various types of breaking on plane sloping beaches. Comparisons of the model test results comprising spatial distribution of wave height and mean water level with experimental data are presented.
基金Supported by National Natural Science Foundation of China under Grant Nos.11772064,and 11502117CAEP Foundation under Grant No.CX2019033+1 种基金the opening project of State Key Laboratory of Explosion Science and Technology(Beijing Institute of Technology)Science Challenge Project under Grant No.JCKY2016212A501
文摘A one-dimensional discrete Boltzmann model for detonation simulation is presented. Instead of numerical solving Navier-Stokes equations, this model obtains the information of flow field through numerical solving specially discretized Boltzmann equation. Several classical benchmarks including Sod shock wave tube, Colella explosion problem,and one-dimensional self-sustainable stable detonation are simulated to validate the new model. Based on the new model,the influence of negative temperature coefficient of reaction rate on detonation is further investigated. It is found that an abnormal detonation with two wave heads periodically appears under negative temperature coefficient condition.The causes of the abnormal detonation are analyzed. One typical cycle of the periodic abnormal detonation and its development process are discussed.
基金jointly supported by the Main Direction Program of Knowledge Innovation of the Chinese Academy of Sciences(Grant No.KZCX2EW203)the National Key Basic Research Program of China(Grant No.2013CB430105)the National Department of Public Benefit Research Foundation(Grant No.GYHY201006031)
文摘A new method for driving a One-Dimensional Stratiform Cold （1DSC） cloud model with Weather Research and Fore casting （WRF） model outputs was developed by conducting numerical experiments for a typical large-scale stratiform rainfall event that took place on 4-5 July 2004 in Changchun, China. Sensitivity test results suggested that, with hydrometeor pro files extracted from the WRF outputs as the initial input, and with continuous updating of soundings and vertical velocities （including downdraft） derived from the WRF model, the new WRF-driven 1DSC modeling system （WRF-1DSC） was able to successfully reproduce both the generation and dissipation processes of the precipitation event. The simulated rainfall intensity showed a time-lag behind that observed, which could have been caused by simulation errors of soundings, vertical velocities and hydrometeor profiles in the WRF output. Taking into consideration the simulated and observed movement path of the precipitation system, a nearby grid point was found to possess more accurate environmental fields in terms of their similarity to those observed in Changchun Station. Using profiles from this nearby grid point, WRF-1DSC was able to repro duce a realistic precipitation pattern. This study demonstrates that 1D cloud-seeding models do indeed have the potential to predict realistic precipitation patterns when properly driven by accurate atmospheric profiles derived from a regional short range forecasting system, This opens a novel and important approach to developing an ensemble-based rain enhancement prediction and operation system under a probabilistic framework concept.
基金Supported by the Natural Science Foundation of China(51705326,52075339)。
文摘In sub nanometer carbon nanotubes,water exhibits unique dynamic characteristics,and in the high-frequency region of the infrared spectrum,where the stretching vibrations of the internal oxygen-hydrogen(O-H)bonds are closely related to the hydrogen bonds(H-bonds)network between water molecules.Therefore,it is crucial to analyze the relationship between these two aspects.In this paper,the infrared spectrum and motion characteristics of the stretching vibrations of the O-H bonds in one-dimensional confined water(1DCW)and bulk water(BW)in(6,6)single-walled carbon nanotubes(SWNT)are studied by molecular dynamics simulations.The results show that the stretching vibrations of the two O-H bonds in 1DCW exhibit different frequencies in the infrared spectrum,while the O-H bonds in BW display two identical main frequency peaks.Further analysis using the spring oscillator model reveals that the difference in the stretching amplitude of the O-H bonds is the main factor causing the change in vibration frequency,where an increase in stretching amplitude leads to a decrease in spring stiffness and,consequently,a lower vibration frequency.A more in-depth study found that the interaction of H-bonds between water molecules is the fundamental cause of the increased stretching amplitude and decreased vibration frequency of the O-H bonds.Finally,by analyzing the motion trajectory of the H atoms,the dynamic differences between 1DCW and BW are clearly revealed.These findings provide a new perspective for understanding the behavior of water molecules at the nanoscale and are of significant importance in advancing the development of infrared spectroscopy detection technology.
基金the Educational Department of Hubei province(Q200628002)the National Science Foundation of China(10701057)
文摘A one-dimensional stationary nonisentropic hydrodynamic model for semiconductor devices with non-constant lattice temperature is studied. This model consists of the equations for the electron density, the electron current density and electron temperature, coupled with the Poisson equation of the electrostatic potential in a bounded interval supplemented with proper boundary conditions. The existence and uniqueness of a strong subsonic steady-state solution with positive particle density and positive temperature is established. The proof is based on the fixed-point arguments, the Stampacchia truncation methods, and the basic energy estimates.
文摘Different phenomenological equations based on plasticity, primary creep （as a viscoplastic mechanism）, secondary creep （as another viscoplastic mechanism） and different combinations of these equations are presented and used to describe the material inelastic deformation in uniaxial test. Agreement of the models with experimental results and with the theoretical concepts and physical realities is the criterion of choosing the most appropriate formulation for uniaxial test. A model is thus proposed in which plastic deformation, primary creep and secondary creep contribute to the inelastic deformation. However, it is believed that the hardening parameter is composed of plastic and primary creep parts. Accordingly, the axial plastic strain in a uniaxial test may no longer be considered as the hardening parameter. Therefore, a proportionality concept is proposed to calculate the plastic contribution of deformation.
基金supported by the Doctoral Scientific Research Funds of Anhui University(J10113190005)the Tian Yuan Foundation of China(11426031)
文摘This paper is concerned with a singular limit for the one-dimensional compress- ible radiation hydrodynamics model. The singular limit we consider corresponds to the physical problem of letting the Bouguer number infinite while keeping the Boltzmann number constant. In the case when the corresponding Euler system admits a contact discontinuity wave, Wang and Xie （2011） [12] recently verified this singular limit and proved that the solution of the compressible radiation hydrodynamics model converges to the strong contact 1 discontinuity wave in the L∞-norm away from the discontinuity line at a rate of ε1/4, as the reciprocal of the Bouguer number tends to zero. In this paper, Wang and Xie＇s convergence rate is improved to ε7/8 by introducing a new a priori assumption and some refined energy estimates. Moreover, it is shown that the radiation flux q tends to zero in the L∞-norm away from the discontinuity line, at a convergence rate as the reciprocal of the Bouguer number tends to zero.
基金The National Basic Research Program of China under contract Nos 2011CB403501 and 2009CB421201the National Natural Science Foundation of China under contract Nos 41176016 and 41076007
文摘The observation data for 5 d at a station in the South China Sea is presented. After brief anMysis of the wind speed, air temperature from the ship-borne meteorological instruments and temperature and salinity profiles from the CTD （conductivity, temperature, depth recorder） data, the authors find that the CTD casts are too sparse for us to understand the diurnal evolution of the thermal structure in the upper ocean. A one-dimensional （1D） numericM code based on Mellor-Yamada turbulence closure model is used to reconstruct the upper ocean thermal structure, utilizing the atmospheric forcing data from ship-borne weather station. The simulation results show good agreement with the observational data; the significance of breaking waves is also briefly discussed. The evolution of turbulence kinetic energy （TKE） and the contribution from shear production and buoy- ancy production are discussed respectively. Finally, several possible factors which might influence the numerical results are briefly analyzed.
基金Project(51078127) supported by the National Natural Science Foundation of ChinaProject(JJ201109091631) supported by the Foundation for Young Scientists of Jiangxi Education Department, China
文摘A modified one-dimensional transient hygrothermal model for multilayer wall was proposed using air humidity ratio and temperature as the driving potentials.The solution for the governing equations was obtained numerically by implementing the finite-difference scheme.To evaluate the accuracy of the model,a test system was built up to measure relative humidity and temperature within a porous wall and compare with the prediction of the model.The prediction results have good agreement with the experimental results.For the interface close to indoor side,the maximum deviation of temperature between calculated and test data is 1.87 K,and the average deviation is 0.95 K;the maximum deviation of relative humidity is 11.4%,and the average deviation is 5.7%.For the interface close to outdoor side,the maximum deviation of temperature between prediction and measurement is 1.78 K,and the average deviation is 1.1 K;the maximum deviation of relative humidity is 9.9%,and the average deviation is 4.2%.
基金supported by the Science Foundation of Henan University of Science and Technology under Grant Nos.05-032 and 2006QN033
文摘A stochastic local limited one-dimensional rice-pile model is numerically investigated. The distributions for avalanche sizes have a clear power-law behavior and it displays a simple finite size scaling. We obtain the avalanche exponents Ts= 1.54±0.10,βs = 2.17±0.10 and TT = 1.80±0.10, βT =1.46 ± 0.10. This self-organized critical model belongs to the same universality class with the Oslo rice-pile model studied by K. Christensen et al. [Phys. Rev. Lett. 77 （1996） 107], a rice-pile model studied by L.A.N. Amaral et al. [Phys. Rev. E 54 （1996） 4512], and a simple deterministic self-organized critical model studied by M.S. Vieira [Phys. Rev. E 61 （2000） 6056].
基金Supported by the National High Technology Research and Development Program of China (2006AA05Z148)
文摘This article aims to investigate the transient behavior of a planar direct internal reforming solid oxide fuel cell （DIR-SOFC） comprehensively. A one-dimensional dynamic model of a planar D1R-SOFC is first developed based on mass and energy balances, and electrochemical principles. Further, a solution strategy is presented to solve the model, and the International Energy Agency （IEA） benchmark test is used to validate the model. Then, through model-based simulations, the steady-state performance of a co-flow planar DIR-SOFC under specified initial operating conditions and its dynamic response to introduced operating parameter disturbances are studied. The dynamic responses of important SOFC variables, such as cell temperature, current density, and cell voltage are all investigated when the SOFC is subjected to the step-changes in various operating parameters including both the load current and the inlet fuel and air flow rates. The results indicate that the rapid dynamics of the current density and the cell voltage are mainly influenced by the gas composition, particularly the H2 molar fraction in anode gas channels, while their slow dynamics are both dominated by the SOLID （including the PEN and interconnects） temperature. As the load current increases, the SOLID temperature and the maximum SOLID temperature gradient both increase, and thereby, the cell breakdown is apt to occur because of excessive thermal stresses. Changing the inlet fuel flow rate might lead to the change in the anode gas composition and the consequent change in the current density distribution and cell voltage. The inlet air flow rate has a great impact on the cell temperature distribution along the cell, and thus, is a suitable manipulated variable to control the cell temperature.
文摘In this paper, the coupling schemes of atmosphere-ocean climate models are discussed with one-dimensional advection equations. The convergence and stability for synchronous and asynchronous schemes are demonstrated and compared.Conclusions inferred from the analysis are given below. The synchronous scheme as well as the asynchronous-implicit scheme in this model are stable for arbitrary integrating time intervals. The asynchronous explicit scheme is unstable under certain conditions, which depend upon advection velocities and heat exchange parameters in the atmosphere and oceans. With both synchronous and asynchronous stable schemes the discrete solutions converge to their unique exact ones. Advections in the atmosphere and ocean accelerate the rate of convergence of the asynchronous-implicit scheme. It is suggusted that the asynchronous-implicit coupling scheme is a stable and efficient method for most climatic simulations.
