Arctic sea ice is an important component of the global climate system and has experienced rapid changes during in the past few decades,the prediction of which is a significant application for climate models.In this st...Arctic sea ice is an important component of the global climate system and has experienced rapid changes during in the past few decades,the prediction of which is a significant application for climate models.In this study,a Localized Error Subspace Transform Kalman Filter is employed in a coupled climate system model(the Flexible Global Ocean–Atmosphere–Land System Model,version f3-L(FGOALS-f3-L))to assimilate sea-ice concentration(SIC)and sea-ice thickness(SIT)data for melting-season ice predictions.The scheme is applied through the following steps:(1)initialization for generating initial ensembles;(2)analysis for assimilating observed data;(3)adoption for dividing ice states into five thickness categories;(4)forecast for evolving the model;(5)resampling for updating model uncertainties.Several experiments were conducted to examine its results and impacts.Compared with the control experiment,the continuous assimilation experiments(CTNs)indicate assimilations improve model SICs and SITs persistently and generate realistic initials.Assimilating SIC+SIT data better corrects overestimated model SITs spatially than when only assimilating SIC data.The continuous assimilation restart experiments indicate the initials from the CTNs correct the overestimated marginal SICs and overall SITs remarkably well,as well as the cold biases in the oceanic and atmospheric models.The initials with SIC+SIT assimilated show more reasonable spatial improvements.Nevertheless,the SICs in the central Arctic undergo abnormal summer reductions,which is probably because overestimated SITs are reduced in the initials but the strong seasonal cycle(summer melting)biases are unchanged.Therefore,since systematic biases are complicated in a coupled system,for FGOALS-f3-L to make better ice predictions,oceanic and atmospheric assimilations are expected required.展开更多
This study evaluates the 1995-2020 global ocean-sea ice simulation using the unstructured-mesh model for prediction across scales(MPAS)-ocean/sea ice model within energy exascale earth system model(E3SM)version 2.1(E3...This study evaluates the 1995-2020 global ocean-sea ice simulation using the unstructured-mesh model for prediction across scales(MPAS)-ocean/sea ice model within energy exascale earth system model(E3SM)version 2.1(E3SMv2-MPAS)at 60 km to 10 km resolution.Multi-source observational data are utilized to validate sea surface temperature/salinity,sea ice,three-dimensional thermal-saline structures,mixed layer depth,ocean heat content,and sea surface height.Key results show the following:(1)E3SMv2-MPAS captures seasonal-to-decadal variability in surface fields and sea ice,but shows systematic biases in sea surface temperature of western boundary currents(inadequate eddy parameterization)and Arctic sea surface salinity(misrepresented freshwater fluxes and mixing processes).(2)The model robustly represents three-dimensional climate variability,yet underestimates mixed layer depth in key regions(Antarctic Circumpolar Current and North Atlantic),revealing deficiencies in extreme mixing.(3)Ocean heat content distributions are well-simulated.(4)Sea surface height spatial patterns and interannual variability are accurately reproduced.This work identifies critical refinements for unstructured-mesh models:mesoscale eddy parameterization,polar ocean-sea ice coupling,and multi-scale energy processes,advancing high-resolution climate model development and laying the groundwork for improved ocean forecasting systems.展开更多
A coupled discrete-continuum simulation incorporating a 3D aspect and non-circular particles was performed to analyze soil-pile interactions during pile penetration in sand.A self-developed non-circular particle numer...A coupled discrete-continuum simulation incorporating a 3D aspect and non-circular particles was performed to analyze soil-pile interactions during pile penetration in sand.A self-developed non-circular particle numerical simulation program was used which considered sand near the pile as interacted particles using a discrete element method;the sand away from the pile was simulated as a continuous medium exhibiting linear elastic behaviors.The domain analyzed was divided into two zones.Contact forces at the interface between the two zones were obtained from a discrete zone and applied to the continuum boundaries as nodal forces,while the interface velocities were obtained from the continuum zone and applied to the discrete boundaries.We show that the coupled discrete-continuum simulation can give a microscopic description of the pile penetration process without losing the discrete nature of the zone concerned,and may significantly improve computational efficiency.展开更多
The datasets for the tier-1 Scenario Model Intercomparison Project(ScenarioMIP)experiments from the Chinese Academy of Sciences(CAS)Flexible Global Ocean-Atmosphere-Land System model,finite-volume version 3(CAS FGOALS...The datasets for the tier-1 Scenario Model Intercomparison Project(ScenarioMIP)experiments from the Chinese Academy of Sciences(CAS)Flexible Global Ocean-Atmosphere-Land System model,finite-volume version 3(CAS FGOALS-f3-L)are described in this study.ScenarioMIP is one of the core MIP experiments in phase 6 of the Coupled Model Intercomparison Project(CMIP6).Considering future CO2,CH4,N2O and other gases’concentrations,as well as land use,the design of ScenarioMIP involves eight pathways,including two tiers(tier-1 and tier-2)of priority.Tier-1 includes four combined Shared Socioeconomic Pathways(SSPs)with radiative forcing,i.e.,SSP1-2.6,SSP2-4.5,SSP3-7.0 and SSP5-8.5,in which the globally averaged radiative forcing at the top of the atmosphere around the year 2100 is approximately 2.6,4.5,7.0 and 8.5 W m−2,respectively.This study provides an introduction to the ScenarioMIP datasets of this model,such as their storage location,sizes,variables,etc.Preliminary analysis indicates that surface air temperatures will increase by about 1.89℃,3.07℃,4.06℃ and 5.17℃ by around 2100 under these four scenarios,respectively.Meanwhile,some other key climate variables,such as sea-ice extension,precipitation,heat content,and sea level rise,also show significant long-term trends associated with the radiative forcing increases.These datasets will help us understand how the climate will change under different anthropogenic and radiative forcings.展开更多
This paper points out that a successful numerical simulation is to construct a correct conceptional model which is very dose to the natural condition. A new model, two dimensional coupled three dimensional model (2D -...This paper points out that a successful numerical simulation is to construct a correct conceptional model which is very dose to the natural condition. A new model, two dimensional coupled three dimensional model (2D -3D ) is presented in the Present paper,which is the most suitable one for the dual - structured - aquifer system. An example of Wenyinghu area is shown.By using the 2D-3D model, a satisfied result of the simulated area is achieved.展开更多
The three-member historical simulations by the Chinese Academy of Sciences Flexible Global Ocean–Atmosphere–Land System model,version f3-L(CAS FGOALS-f3-L),which is contributing to phase 6 of the Coupled Model Inter...The three-member historical simulations by the Chinese Academy of Sciences Flexible Global Ocean–Atmosphere–Land System model,version f3-L(CAS FGOALS-f3-L),which is contributing to phase 6 of the Coupled Model Intercomparison Project(CMIP6),are described in this study.The details of the CAS FGOALS-f3-L model,experiment settings and output datasets are briefly introduced.The datasets include monthly and daily outputs from the atmospheric,oceanic,land and sea-ice component models of CAS FGOALS-f3-L,and all these data have been published online in the Earth System Grid Federation(ESGF,https://esgf-node.llnl.gov/projects/cmip6/).The three ensembles are initialized from the 600th,650th and 700th model year of the preindustrial experiment(piControl)and forced by the same historical forcing provided by CMIP6 from 1850 to 2014.The performance of the coupled model is validated in comparison with some recent observed atmospheric and oceanic datasets.It is shown that CAS FGOALS-f3-L is able to reproduce the main features of the modern climate,including the climatology of air surface temperature and precipitation,the long-term changes in global mean surface air temperature,ocean heat content and sea surface steric height,and the horizontal and vertical distribution of temperature in the ocean and atmosphere.Meanwhile,like other state-of-the-art coupled GCMs,there are still some obvious biases in the historical simulations,which are also illustrated.This paper can help users to better understand the advantages and biases of the model and the datasets。展开更多
A combined beam model representing the periodicity of the microstructure and micro deformation of 3D woven composites is developed for predicting mechanical properties. The model considers the effects of off axial ten...A combined beam model representing the periodicity of the microstructure and micro deformation of 3D woven composites is developed for predicting mechanical properties. The model considers the effects of off axial tension/compression and bending/shearing couplings as well as the mutual reactions of fiber yarns. The method determining microstructure by using woven parameters is described for a typical 3D woven composite material. An analytical cell, constructed by a minimum periodic section of yarn and interlayer matrix, is adopted. Micro stresses in the cell under in-plane tensile loading are obtained by using the proposed beam model and macro modulus is then obtained by the averaging method. Material tests and a 2D micro FEM analysis are made to evaluate this model. Analyses reveal that micro stress caused by tensile/bending coupling effect is not negligible in the stress analysis.展开更多
A regional atmosphere-ocean coupled model, RegCM3-POM, was developed by coupling the regional climate model （RegCM3） with the Princeton Ocean Model （POM）. The performance of RegCM3-POM in simulating a persistent s...A regional atmosphere-ocean coupled model, RegCM3-POM, was developed by coupling the regional climate model （RegCM3） with the Princeton Ocean Model （POM）. The performance of RegCM3-POM in simulating a persistent snow storm over southern China and the impact of the Madden Julian oscillation （MJO） on this persistent snow storm were investigated. Compared with the stand-alone RegCM3, the coupled model performed better at reproducing the spatial-temporal evolution and intensity of the precipitation episodes. The power spectral analysis indicated that the coupled model successfully captured the dominant period between 30 and 60 days in the precipitation field, leading to a notable improvement in simulating the magnitude of intraseasonal precipitation variation, and further in enhancing the intensity of the simulated precipitation. These improvements were mainly due to the well-simulated low-frequency oscillation center and its eastward propagation characteristics in each MJO phase by RegCM3-POM, which improved the simulations of MJO-related low-frequency vertical motions, water vapor transport, and the deep inversion layer that can directly influence the precipitation event and that further improved the simulated MJOprecipitation relationship. Analysis of the phase relationship between convection and SST indicated that RegCM3-POM exhibits a near-quadrature relation between the simulated convection and SST anomalies, which was consistent with the observations. However, such a near-quadrature relation was not as significant when the stand-alone RegCM3 was used. This difference indicated that the inherent coupled feedback process between the ocean and atmosphere in RegCM3-POM played an important part in reproducing the features of the MJO that accompanied the snow storm.展开更多
The pneumatic muscle actuator(PMA)has many advantages,such as good flexibility,high power/weight ratio,but its nonlinearity makes it difficult to build a static mathematical model with high precision.A new method is p...The pneumatic muscle actuator(PMA)has many advantages,such as good flexibility,high power/weight ratio,but its nonlinearity makes it difficult to build a static mathematical model with high precision.A new method is proposed to establish the model of PMA.The concept of hybrid elastic modulus which is related to the static characteristic of PMA is put forward,and the energy conservation law is used to achieve the expression of the hybrid elastic modulus,which can be fitted out based on experimental data,and the model of PMA can be derived from this expression.At the same time,a 3-DOF parallel mechanism(a new bionic shoulder joint)driven by five PMAs is designed.This bionic shoulder joint adopts the structure of two antagonistic PMAs actualizing a rotation control and three PMAs controlling another two rotations to get better rotation characteristics.The kinematic and dynamic characteristics of the mechanism are analyzed and a new static model of PMA is used to control it.Experimental results demonstrate the effectiveness of this new static model.展开更多
In most studies of tunnel boring machine(TBM)tunnelling, the groundwater pressure was not considered, or was simplified and exerted on the boundary of lining structure. Meanwhile, the leakage, which mainly occurs in t...In most studies of tunnel boring machine(TBM)tunnelling, the groundwater pressure was not considered, or was simplified and exerted on the boundary of lining structure. Meanwhile, the leakage, which mainly occurs in the segment joints, was often ignored in the relevant studies of TBM tunnelling. Additionally, the geological models in these studies were simplified to different extents, and mostly were simplified as homogenous bodies. Considering the deficiencies above, a 3D refined model of the surrounding rock of a tunnel is firstly established using NURBS-TIN-BRe P hybrid data structure in this paper. Then the seepage field of the surrounding rock considering the leakage in the segment joints is simulated. Finally, the stability of TBM water diversion tunnel is studied coupled with the seepage simulation, to analyze the stress-strain conditions, the axial force and the bending moment of tunnel segment considering the leakage in the segment joints. The results illustrate that the maximum radial displacement, the minimum principal stress, the maximum principal stress and the axial force of segment lining considering the seepage effect are all larger than those disregarding the seepage effect.展开更多
Following the High-Resolution Model Intercomparison Project(HighResMIP)Tier 2 protocol under the Coupled Model Intercomparison Project Phase 6(CMIP6),three numerical experiments are conducted with the Chinese Academy ...Following the High-Resolution Model Intercomparison Project(HighResMIP)Tier 2 protocol under the Coupled Model Intercomparison Project Phase 6(CMIP6),three numerical experiments are conducted with the Chinese Academy of Sciences Flexible Global Ocean-Atmosphere-Land System Model,version f3-H(CAS FGOALS-f3-H),and a 101-year(1950–2050)global high-resolution simulation dataset is presented in this study.The basic configuration of the FGOALSf3-H model and numerical experiments design are briefly described,and then the historical simulation is validated.Forced by observed radiative agents from 1950 to 2014,the coupled model essentially reproduces the observed long-term trends of temperature,precipitation,and sea ice extent,as well as the large-scale pattern of temperature and precipitation.With an approximate 0.25°horizontal resolution in the atmosphere and 0.1°in the ocean,the coupled models also simulate energetic western boundary currents and the Antarctic Circulation Current(ACC),reasonable characteristics of extreme precipitation,and realistic frontal scale air-sea interaction.The dataset and supporting detailed information have been published in the Earth System Grid Federation.展开更多
The shear failure of intact rock under thermo-mechanical(TM)coupling conditions is common,such as in enhanced geothermal mining and deep mine construction.Under the effect of a continuous engineering disturbance,shear...The shear failure of intact rock under thermo-mechanical(TM)coupling conditions is common,such as in enhanced geothermal mining and deep mine construction.Under the effect of a continuous engineering disturbance,shear-formed fractures are prone to secondary instability,posing a severe threat to deep engineering.Although numerous studies regarding three-dimensional(3D)morphologies of fracture surfaces have been conducted,the understanding of shear-formed fractures under TM coupling conditions is limited.In this study,direct shear tests of intact granite under various TM coupling conditions were conducted,followed by 3D laser scanning tests of shear-formed fractures.Test results demonstrated that the peak shear strength of intact granite is positively correlated with the normal stress,whereas it is negatively correlated with the temperature.The internal friction angle and cohesion of intact granite significantly decrease with an increase in the temperature.The anisotropy,roughness value,and height of the asperities on the fracture surfaces are reduced as the normal stress increases,whereas their variation trends are the opposite as the temperature increases.The macroscopic failure mode of intact granite under TM coupling conditions is dominated by mixed tensileeshear and shear failures.As the normal stress increases,intragranular fractures are developed ranging from a local to a global distribution,and the macroscopic failure mode of intact granite changes from mixed tensileeshear to shear failure.Finally,3D morphological characteristics of the asperities on the shear-formed fracture surfaces were analyzed,and a quadrangular pyramid conceptual model representing these asperities was proposed and sufficiently verified.展开更多
As one of the most important ways to reduce the greenhouse gas emission,carbon dioxide（CO2）enhanced gas recovery（CO2-EGR） is attractive since the gas recovery can be enhanced simultaneously with CO2sequestration.B...As one of the most important ways to reduce the greenhouse gas emission,carbon dioxide（CO2）enhanced gas recovery（CO2-EGR） is attractive since the gas recovery can be enhanced simultaneously with CO2sequestration.Based on the existing equation of state（EOS） module of TOUGH2 MP,extEOS7C is developed to calculate the phase partition of H2O-CO2-CH4-NaCl mixtures accurately with consideration of dissolved NaCI and brine properties at high pressure and temperature conditions.Verifications show that it can be applied up to the pressure of 100 MPa and temperature of 150℃.The module was implemented in the linked simulator TOUGH2MP-FLAC3 D for the coupled hydro-mechanical simulations.A simplified three-dimensional（3D）1/4 model（2.2 km×1 km×1 km） which consists of the whole reservoir,caprock and baserock was generated based on the geological conditions of a gas field in the North German Basin.The simulation results show that,under an injection rate of 200,000 t/yr and production rate of 200,000 sm3/d,CO2breakthrough occurred in the case with the initial reservoir pressure of 5 MPa but did not occur in the case of 42 MPa.Under low pressure conditions,the pressure driven horizontal transport is the dominant process;while under high pressure conditions,the density driven vertical flow is dominant.Under the considered conditions,the CO2-EGR caused only small pressure changes.The largest pore pressure increase（2 MPa） and uplift（7 mm） occurred at the caprock bottom induced by only CO2injection.The caprock had still the primary stress state and its integrity was not affected.The formation water salinity and temperature variations of ±20℃ had small influences on the CO2-EGR process.In order to slow down the breakthrough,it is suggested that CO2-EGR should be carried out before the reservoir pressure drops below the critical pressure of CO2.展开更多
Tension Leg Platform(TLP)in deepwater oil and gas field development usually consists of a hull,tendons,and top tension risers(TTRs).To maintain its top tension,each TTR is connected with a tensioner system to the hull...Tension Leg Platform(TLP)in deepwater oil and gas field development usually consists of a hull,tendons,and top tension risers(TTRs).To maintain its top tension,each TTR is connected with a tensioner system to the hull.Owing to the complicated configuration of the tensioners,the hull and TTRs form a strong coupled system.Traditionally,some simplified tensioner models are applied to analyze the TLP structures.There is a large discrepancy between their analysis results and the actual mechanism behaviors of a tensioner.It is very necessary to develop a more detailed tensioner model to consider the coupling effects between TLP and TTRs.In the present study,a fully coupled TLP hull-TTR system for hydrodynamic numerical simulation is established.A specific hydraulic pneumatic tensioner is modeled by considering 4 cylinders.The production TTR model is stacked up by specific riser joints.The simulation is also extended to analyze an array of TTRs.Different regular and irregular waves are considered.The behaviors of different cylinders are presented.The results show that it is important to consider the specific configurations of the tensioner and TTRs,which may lead to obviously different response behaviors,compared with those from a simplified model.展开更多
Acoustic-elastic coupled media is often encountered in most marine explorations, and accurate simulation of acoustic-elastic coupled media is of great significance. At present, the study of acoustic-elastic coupled me...Acoustic-elastic coupled media is often encountered in most marine explorations, and accurate simulation of acoustic-elastic coupled media is of great significance. At present, the study of acoustic-elastic coupled media still assumes that the solid of the acoustic-elastic coupled media is isotropic, but this assumption is not in accordance with the actual situation. In this paper, we derive the solid media of acoustic-elastic coupled media from isotropic media to anisotropic media, and propose an acoustic-elastic coupled medium based ontransverse isotropic media with vertical symmetric axes(VTI) to improve the accuracy of forward modeling. Based on the relationship between the Thomsen parameter and the coefficient matrix of the anisotropic elastic wave equation, we transform the Thomson parameter into a velocity model with anisotropic properties. We use a staggered grid finite difference method to simulate the propagation of a wavefield in a three-dimensional acoustic-elastic coupled media. We obtain the snapshots of the wave field when the solid of the acoustic-elastic coupled media is an isotropic medium and a VTI media. When the solid of the acoustic-elastic coupled media is considered VTI media, we can observe the qP wave and qS wave that cannot be observed in the isotropic medium from the wave field snapshot. We can also find that the seismic records obtained by the method we use are more realistic. The algorithm proposed in this paper is of great significance for high-precision ocean numerical simulation.展开更多
Typically,seat or floor acceleration is used to evaluate the ride comfort of a high-speed train.However,the dynamic performance of the human body significantly differs from that of the floor.Therefore,using the car bo...Typically,seat or floor acceleration is used to evaluate the ride comfort of a high-speed train.However,the dynamic performance of the human body significantly differs from that of the floor.Therefore,using the car body floor and seat accelerations to calculate the ride comfort index of a high-speed train may not reflect the true feelings of passengers.In this study,a 3D human-seat-vehicle-track coupling model was established to investigate the ride comfort of highspeed train passengers.The seated human model,which considers the longitudinal,lateral,vertical,pitching,yawing,and rolling motions,comprises the head,upper torso,lower torso,pelvis,thighs,and shanks.The model parameters were determined using multi-axis excitation measurement data based on a genetic algorithm.Subsequently,the applicability of the small-angle assumption and natural modes of the human model is analyzed.Using the coupling system model,the vibration characteristics of the human-seat interaction surface were analyzed.The ride comfort of the high-speed train and human body dynamic performance were analyzed under normal conditions,track geometric irregularities and train meeting conditions.The results showed that the passenger seats in the front and rear rows adjacent to the window had a higher acceleration value than the others.The human backrest and seat pad connection points have higher vibration amplitudes than the car body floor in the human-sensitive frequency range,indicating that using the acceleration values on the floor may underestimate the discomfort of passengers.The ride comfort of high-speed trains diminishes in the presence of track geometric irregularities and when trains pass each other.When the excitation frequency of track geometry irregularities approached the natural frequency of the human-seat-vehicle system,ride comfort in high-speed trains decreased significantly.Moreover,using seat acceleration to evaluate passenger ride comfort overlooks the vibration characteristics of the human body.The transient aerodynamic force generated when the train meets can cause a larger car body roll and lateral motion at 2 Hz,which,in turn,decreases the passenger ride comfort.This study presents a detailed human-seat-vehicle-track coupling system that can reflect a passenger’s dynamic performance under complex operating conditions.展开更多
This study proposes a three-dimensional(3D)coupled magneto-electro-elastic problem for the static analysis of multilayered plates embedding piezomagnetic and piezoelectric layers by considering both sensor and actuato...This study proposes a three-dimensional(3D)coupled magneto-electro-elastic problem for the static analysis of multilayered plates embedding piezomagnetic and piezoelectric layers by considering both sensor and actuator configurations.The 3D governing equations for the magneto-electro-elastic static behavior of plates are explicitly show that are made by the three 3D equilibrium equations,the 3D divergence equation for magnetic induction,and the 3D divergence equation for the electric displacement.The proposed solution involves the exponential matrix in the thickness direction and primary variables’harmonic forms in the in-plane ones.A closed-form solution is performed considering simply-supported boundary conditions.Interlaminar continuity conditions are imposed for displacements,magnetic potential,electric potential,transverse shear/normal stresses,transverse normal magnetic induction and transverse normal electric displacement.Therefore,a layerwise approach is adopted.The results section is composed of an assessment part,where the present model is compared to past 3D electro-elastic or magneto-elastic formulations and a new benchmark part.Benchmarks consider sensor and actuator plate configurations for the fully coupled magneto-electro-elastic cases for different thickness ratios.Tabular and graphical results are presented for displacements,stresses,magnetic potential,electric potential,transverse normal magnetic induction and transverse normal electric displacement.For each presented benchmark,magneto-electro-elastic coupling and thickness and material layer effects are discussed in depth.展开更多
基金jointly funded by the National Natural Science Foundation of China(NSFC)[grant number 42130608]the China Postdoctoral Science Foundation[grant number 2024M753169]。
文摘Arctic sea ice is an important component of the global climate system and has experienced rapid changes during in the past few decades,the prediction of which is a significant application for climate models.In this study,a Localized Error Subspace Transform Kalman Filter is employed in a coupled climate system model(the Flexible Global Ocean–Atmosphere–Land System Model,version f3-L(FGOALS-f3-L))to assimilate sea-ice concentration(SIC)and sea-ice thickness(SIT)data for melting-season ice predictions.The scheme is applied through the following steps:(1)initialization for generating initial ensembles;(2)analysis for assimilating observed data;(3)adoption for dividing ice states into five thickness categories;(4)forecast for evolving the model;(5)resampling for updating model uncertainties.Several experiments were conducted to examine its results and impacts.Compared with the control experiment,the continuous assimilation experiments(CTNs)indicate assimilations improve model SICs and SITs persistently and generate realistic initials.Assimilating SIC+SIT data better corrects overestimated model SITs spatially than when only assimilating SIC data.The continuous assimilation restart experiments indicate the initials from the CTNs correct the overestimated marginal SICs and overall SITs remarkably well,as well as the cold biases in the oceanic and atmospheric models.The initials with SIC+SIT assimilated show more reasonable spatial improvements.Nevertheless,the SICs in the central Arctic undergo abnormal summer reductions,which is probably because overestimated SITs are reduced in the initials but the strong seasonal cycle(summer melting)biases are unchanged.Therefore,since systematic biases are complicated in a coupled system,for FGOALS-f3-L to make better ice predictions,oceanic and atmospheric assimilations are expected required.
基金The National Key R&D Program of China under contract No.2021YFC3101503the Science and Technology Innovation Program of Hunan Province under contract No.2022RC3070+1 种基金the National Natural Science Foundation of China under contract Nos 42305176 and 42276205the Hunan Provincial Natural Science Foundation of China under contract No.2023JJ10053.
文摘This study evaluates the 1995-2020 global ocean-sea ice simulation using the unstructured-mesh model for prediction across scales(MPAS)-ocean/sea ice model within energy exascale earth system model(E3SM)version 2.1(E3SMv2-MPAS)at 60 km to 10 km resolution.Multi-source observational data are utilized to validate sea surface temperature/salinity,sea ice,three-dimensional thermal-saline structures,mixed layer depth,ocean heat content,and sea surface height.Key results show the following:(1)E3SMv2-MPAS captures seasonal-to-decadal variability in surface fields and sea ice,but shows systematic biases in sea surface temperature of western boundary currents(inadequate eddy parameterization)and Arctic sea surface salinity(misrepresented freshwater fluxes and mixing processes).(2)The model robustly represents three-dimensional climate variability,yet underestimates mixed layer depth in key regions(Antarctic Circumpolar Current and North Atlantic),revealing deficiencies in extreme mixing.(3)Ocean heat content distributions are well-simulated.(4)Sea surface height spatial patterns and interannual variability are accurately reproduced.This work identifies critical refinements for unstructured-mesh models:mesoscale eddy parameterization,polar ocean-sea ice coupling,and multi-scale energy processes,advancing high-resolution climate model development and laying the groundwork for improved ocean forecasting systems.
基金Project (No.90815008) supported by the National Natural Science Foundation of China
文摘A coupled discrete-continuum simulation incorporating a 3D aspect and non-circular particles was performed to analyze soil-pile interactions during pile penetration in sand.A self-developed non-circular particle numerical simulation program was used which considered sand near the pile as interacted particles using a discrete element method;the sand away from the pile was simulated as a continuous medium exhibiting linear elastic behaviors.The domain analyzed was divided into two zones.Contact forces at the interface between the two zones were obtained from a discrete zone and applied to the continuum boundaries as nodal forces,while the interface velocities were obtained from the continuum zone and applied to the discrete boundaries.We show that the coupled discrete-continuum simulation can give a microscopic description of the pile penetration process without losing the discrete nature of the zone concerned,and may significantly improve computational efficiency.
基金supported by the Strategic Priority Research Program of the Chinese Academy of Sciences(Grant Nos.XDA19060102 and XDB42000000)the National Natural Science Foundation of China(Grants Nos.41530426 and 91958201)。
文摘The datasets for the tier-1 Scenario Model Intercomparison Project(ScenarioMIP)experiments from the Chinese Academy of Sciences(CAS)Flexible Global Ocean-Atmosphere-Land System model,finite-volume version 3(CAS FGOALS-f3-L)are described in this study.ScenarioMIP is one of the core MIP experiments in phase 6 of the Coupled Model Intercomparison Project(CMIP6).Considering future CO2,CH4,N2O and other gases’concentrations,as well as land use,the design of ScenarioMIP involves eight pathways,including two tiers(tier-1 and tier-2)of priority.Tier-1 includes four combined Shared Socioeconomic Pathways(SSPs)with radiative forcing,i.e.,SSP1-2.6,SSP2-4.5,SSP3-7.0 and SSP5-8.5,in which the globally averaged radiative forcing at the top of the atmosphere around the year 2100 is approximately 2.6,4.5,7.0 and 8.5 W m−2,respectively.This study provides an introduction to the ScenarioMIP datasets of this model,such as their storage location,sizes,variables,etc.Preliminary analysis indicates that surface air temperatures will increase by about 1.89℃,3.07℃,4.06℃ and 5.17℃ by around 2100 under these four scenarios,respectively.Meanwhile,some other key climate variables,such as sea-ice extension,precipitation,heat content,and sea level rise,also show significant long-term trends associated with the radiative forcing increases.These datasets will help us understand how the climate will change under different anthropogenic and radiative forcings.
文摘This paper points out that a successful numerical simulation is to construct a correct conceptional model which is very dose to the natural condition. A new model, two dimensional coupled three dimensional model (2D -3D ) is presented in the Present paper,which is the most suitable one for the dual - structured - aquifer system. An example of Wenyinghu area is shown.By using the 2D-3D model, a satisfied result of the simulated area is achieved.
基金This study is jointly supported by the Strategic Priority Research Program of Chinese Academy of Sciences(Grant Nos.XDA19060102 and XDB42010400)the Natural Science Foundation of China(Grant Nos.41530426,91958201 and 41931183).
文摘The three-member historical simulations by the Chinese Academy of Sciences Flexible Global Ocean–Atmosphere–Land System model,version f3-L(CAS FGOALS-f3-L),which is contributing to phase 6 of the Coupled Model Intercomparison Project(CMIP6),are described in this study.The details of the CAS FGOALS-f3-L model,experiment settings and output datasets are briefly introduced.The datasets include monthly and daily outputs from the atmospheric,oceanic,land and sea-ice component models of CAS FGOALS-f3-L,and all these data have been published online in the Earth System Grid Federation(ESGF,https://esgf-node.llnl.gov/projects/cmip6/).The three ensembles are initialized from the 600th,650th and 700th model year of the preindustrial experiment(piControl)and forced by the same historical forcing provided by CMIP6 from 1850 to 2014.The performance of the coupled model is validated in comparison with some recent observed atmospheric and oceanic datasets.It is shown that CAS FGOALS-f3-L is able to reproduce the main features of the modern climate,including the climatology of air surface temperature and precipitation,the long-term changes in global mean surface air temperature,ocean heat content and sea surface steric height,and the horizontal and vertical distribution of temperature in the ocean and atmosphere.Meanwhile,like other state-of-the-art coupled GCMs,there are still some obvious biases in the historical simulations,which are also illustrated.This paper can help users to better understand the advantages and biases of the model and the datasets。
文摘A combined beam model representing the periodicity of the microstructure and micro deformation of 3D woven composites is developed for predicting mechanical properties. The model considers the effects of off axial tension/compression and bending/shearing couplings as well as the mutual reactions of fiber yarns. The method determining microstructure by using woven parameters is described for a typical 3D woven composite material. An analytical cell, constructed by a minimum periodic section of yarn and interlayer matrix, is adopted. Micro stresses in the cell under in-plane tensile loading are obtained by using the proposed beam model and macro modulus is then obtained by the averaging method. Material tests and a 2D micro FEM analysis are made to evaluate this model. Analyses reveal that micro stress caused by tensile/bending coupling effect is not negligible in the stress analysis.
基金supported by grants from the National Basic Research Program of China (973 Program,Grant Nos.2012CB955901 and 2011CB952002)the National Science and Technology Support Program of China (Grant No.2009BAC51B03)the National Natural Science Foundation of China (Grant Nos. 41105044 and 41105045)
文摘A regional atmosphere-ocean coupled model, RegCM3-POM, was developed by coupling the regional climate model （RegCM3） with the Princeton Ocean Model （POM）. The performance of RegCM3-POM in simulating a persistent snow storm over southern China and the impact of the Madden Julian oscillation （MJO） on this persistent snow storm were investigated. Compared with the stand-alone RegCM3, the coupled model performed better at reproducing the spatial-temporal evolution and intensity of the precipitation episodes. The power spectral analysis indicated that the coupled model successfully captured the dominant period between 30 and 60 days in the precipitation field, leading to a notable improvement in simulating the magnitude of intraseasonal precipitation variation, and further in enhancing the intensity of the simulated precipitation. These improvements were mainly due to the well-simulated low-frequency oscillation center and its eastward propagation characteristics in each MJO phase by RegCM3-POM, which improved the simulations of MJO-related low-frequency vertical motions, water vapor transport, and the deep inversion layer that can directly influence the precipitation event and that further improved the simulated MJOprecipitation relationship. Analysis of the phase relationship between convection and SST indicated that RegCM3-POM exhibits a near-quadrature relation between the simulated convection and SST anomalies, which was consistent with the observations. However, such a near-quadrature relation was not as significant when the stand-alone RegCM3 was used. This difference indicated that the inherent coupled feedback process between the ocean and atmosphere in RegCM3-POM played an important part in reproducing the features of the MJO that accompanied the snow storm.
基金supported by the National Natural Science Foundation of China(No. 51405229)the Natural Science Foundation of Jiangsu Province of China(Nos. BK20151470,BK20171416)
文摘The pneumatic muscle actuator(PMA)has many advantages,such as good flexibility,high power/weight ratio,but its nonlinearity makes it difficult to build a static mathematical model with high precision.A new method is proposed to establish the model of PMA.The concept of hybrid elastic modulus which is related to the static characteristic of PMA is put forward,and the energy conservation law is used to achieve the expression of the hybrid elastic modulus,which can be fitted out based on experimental data,and the model of PMA can be derived from this expression.At the same time,a 3-DOF parallel mechanism(a new bionic shoulder joint)driven by five PMAs is designed.This bionic shoulder joint adopts the structure of two antagonistic PMAs actualizing a rotation control and three PMAs controlling another two rotations to get better rotation characteristics.The kinematic and dynamic characteristics of the mechanism are analyzed and a new static model of PMA is used to control it.Experimental results demonstrate the effectiveness of this new static model.
基金Supported by the Foundation for Innovation Research Groups of the National Natural Science Foundation of China(No.51321065)Tianjin Research Program of Application Foundation and Advanced Technology(No.12JCZDJC29200)Tianjin Natural Science Foundation(No.13JCYBJC19500)
文摘In most studies of tunnel boring machine(TBM)tunnelling, the groundwater pressure was not considered, or was simplified and exerted on the boundary of lining structure. Meanwhile, the leakage, which mainly occurs in the segment joints, was often ignored in the relevant studies of TBM tunnelling. Additionally, the geological models in these studies were simplified to different extents, and mostly were simplified as homogenous bodies. Considering the deficiencies above, a 3D refined model of the surrounding rock of a tunnel is firstly established using NURBS-TIN-BRe P hybrid data structure in this paper. Then the seepage field of the surrounding rock considering the leakage in the segment joints is simulated. Finally, the stability of TBM water diversion tunnel is studied coupled with the seepage simulation, to analyze the stress-strain conditions, the axial force and the bending moment of tunnel segment considering the leakage in the segment joints. The results illustrate that the maximum radial displacement, the minimum principal stress, the maximum principal stress and the axial force of segment lining considering the seepage effect are all larger than those disregarding the seepage effect.
基金jointly supported by the Strategic Priority Research Program of the Chinese Academy of Sciences(Grant Nos.XDA19060102 and XDB42000000)National Natural Science Foundation of China(Grant Nos.91958201 and 42130608)+1 种基金the National Key Research and Development Program of China(Grant No.2020YFA0608800)supported by the National Key Scientific and Technological Infrastructure project“Earth System Numerical Simulation Facility”(EarthLab)。
文摘Following the High-Resolution Model Intercomparison Project(HighResMIP)Tier 2 protocol under the Coupled Model Intercomparison Project Phase 6(CMIP6),three numerical experiments are conducted with the Chinese Academy of Sciences Flexible Global Ocean-Atmosphere-Land System Model,version f3-H(CAS FGOALS-f3-H),and a 101-year(1950–2050)global high-resolution simulation dataset is presented in this study.The basic configuration of the FGOALSf3-H model and numerical experiments design are briefly described,and then the historical simulation is validated.Forced by observed radiative agents from 1950 to 2014,the coupled model essentially reproduces the observed long-term trends of temperature,precipitation,and sea ice extent,as well as the large-scale pattern of temperature and precipitation.With an approximate 0.25°horizontal resolution in the atmosphere and 0.1°in the ocean,the coupled models also simulate energetic western boundary currents and the Antarctic Circulation Current(ACC),reasonable characteristics of extreme precipitation,and realistic frontal scale air-sea interaction.The dataset and supporting detailed information have been published in the Earth System Grid Federation.
基金supported by the National Natural Science Foundation of China(Grant No.51974173)the Natural Science Foundation of Shandong Province,China(Grant No.ZR2020QD122).
文摘The shear failure of intact rock under thermo-mechanical(TM)coupling conditions is common,such as in enhanced geothermal mining and deep mine construction.Under the effect of a continuous engineering disturbance,shear-formed fractures are prone to secondary instability,posing a severe threat to deep engineering.Although numerous studies regarding three-dimensional(3D)morphologies of fracture surfaces have been conducted,the understanding of shear-formed fractures under TM coupling conditions is limited.In this study,direct shear tests of intact granite under various TM coupling conditions were conducted,followed by 3D laser scanning tests of shear-formed fractures.Test results demonstrated that the peak shear strength of intact granite is positively correlated with the normal stress,whereas it is negatively correlated with the temperature.The internal friction angle and cohesion of intact granite significantly decrease with an increase in the temperature.The anisotropy,roughness value,and height of the asperities on the fracture surfaces are reduced as the normal stress increases,whereas their variation trends are the opposite as the temperature increases.The macroscopic failure mode of intact granite under TM coupling conditions is dominated by mixed tensileeshear and shear failures.As the normal stress increases,intragranular fractures are developed ranging from a local to a global distribution,and the macroscopic failure mode of intact granite changes from mixed tensileeshear to shear failure.Finally,3D morphological characteristics of the asperities on the shear-formed fracture surfaces were analyzed,and a quadrangular pyramid conceptual model representing these asperities was proposed and sufficiently verified.
基金funded by the National Natural Science Foundation of China(Grant No.NSFC51374147)the German Society for Petroleum and Coal Science and Technology(Grant No.DGMK680-4)
文摘As one of the most important ways to reduce the greenhouse gas emission,carbon dioxide（CO2）enhanced gas recovery（CO2-EGR） is attractive since the gas recovery can be enhanced simultaneously with CO2sequestration.Based on the existing equation of state（EOS） module of TOUGH2 MP,extEOS7C is developed to calculate the phase partition of H2O-CO2-CH4-NaCl mixtures accurately with consideration of dissolved NaCI and brine properties at high pressure and temperature conditions.Verifications show that it can be applied up to the pressure of 100 MPa and temperature of 150℃.The module was implemented in the linked simulator TOUGH2MP-FLAC3 D for the coupled hydro-mechanical simulations.A simplified three-dimensional（3D）1/4 model（2.2 km×1 km×1 km） which consists of the whole reservoir,caprock and baserock was generated based on the geological conditions of a gas field in the North German Basin.The simulation results show that,under an injection rate of 200,000 t/yr and production rate of 200,000 sm3/d,CO2breakthrough occurred in the case with the initial reservoir pressure of 5 MPa but did not occur in the case of 42 MPa.Under low pressure conditions,the pressure driven horizontal transport is the dominant process;while under high pressure conditions,the density driven vertical flow is dominant.Under the considered conditions,the CO2-EGR caused only small pressure changes.The largest pore pressure increase（2 MPa） and uplift（7 mm） occurred at the caprock bottom induced by only CO2injection.The caprock had still the primary stress state and its integrity was not affected.The formation water salinity and temperature variations of ±20℃ had small influences on the CO2-EGR process.In order to slow down the breakthrough,it is suggested that CO2-EGR should be carried out before the reservoir pressure drops below the critical pressure of CO2.
基金The research was financially supported by the National Natural Science Foundation of China for Youth(Grant No.51609169)Guangxi Science and Technology Major Project(Grant No.Guike AA17292007)+2 种基金the National Key R&D Program of China(Grant No.2018YFC0310502)National Natural Science Foundation of China(Grant No.51779173)China Scholarship Council(CSC).
文摘Tension Leg Platform(TLP)in deepwater oil and gas field development usually consists of a hull,tendons,and top tension risers(TTRs).To maintain its top tension,each TTR is connected with a tensioner system to the hull.Owing to the complicated configuration of the tensioners,the hull and TTRs form a strong coupled system.Traditionally,some simplified tensioner models are applied to analyze the TLP structures.There is a large discrepancy between their analysis results and the actual mechanism behaviors of a tensioner.It is very necessary to develop a more detailed tensioner model to consider the coupling effects between TLP and TTRs.In the present study,a fully coupled TLP hull-TTR system for hydrodynamic numerical simulation is established.A specific hydraulic pneumatic tensioner is modeled by considering 4 cylinders.The production TTR model is stacked up by specific riser joints.The simulation is also extended to analyze an array of TTRs.Different regular and irregular waves are considered.The behaviors of different cylinders are presented.The results show that it is important to consider the specific configurations of the tensioner and TTRs,which may lead to obviously different response behaviors,compared with those from a simplified model.
基金Supported by Major Project of National Science and Technology of China(No.2016ZX05026-002-003)National Natural Science Foundation of China(No.41374108)
文摘Acoustic-elastic coupled media is often encountered in most marine explorations, and accurate simulation of acoustic-elastic coupled media is of great significance. At present, the study of acoustic-elastic coupled media still assumes that the solid of the acoustic-elastic coupled media is isotropic, but this assumption is not in accordance with the actual situation. In this paper, we derive the solid media of acoustic-elastic coupled media from isotropic media to anisotropic media, and propose an acoustic-elastic coupled medium based ontransverse isotropic media with vertical symmetric axes(VTI) to improve the accuracy of forward modeling. Based on the relationship between the Thomsen parameter and the coefficient matrix of the anisotropic elastic wave equation, we transform the Thomson parameter into a velocity model with anisotropic properties. We use a staggered grid finite difference method to simulate the propagation of a wavefield in a three-dimensional acoustic-elastic coupled media. We obtain the snapshots of the wave field when the solid of the acoustic-elastic coupled media is an isotropic medium and a VTI media. When the solid of the acoustic-elastic coupled media is considered VTI media, we can observe the qP wave and qS wave that cannot be observed in the isotropic medium from the wave field snapshot. We can also find that the seismic records obtained by the method we use are more realistic. The algorithm proposed in this paper is of great significance for high-precision ocean numerical simulation.
基金Supported by National Natural Science Foundation of China(Grant No.U1934203)Research and Development Project of Science and Technology of China Railway Corporation(Grant No.P2023T002)。
文摘Typically,seat or floor acceleration is used to evaluate the ride comfort of a high-speed train.However,the dynamic performance of the human body significantly differs from that of the floor.Therefore,using the car body floor and seat accelerations to calculate the ride comfort index of a high-speed train may not reflect the true feelings of passengers.In this study,a 3D human-seat-vehicle-track coupling model was established to investigate the ride comfort of highspeed train passengers.The seated human model,which considers the longitudinal,lateral,vertical,pitching,yawing,and rolling motions,comprises the head,upper torso,lower torso,pelvis,thighs,and shanks.The model parameters were determined using multi-axis excitation measurement data based on a genetic algorithm.Subsequently,the applicability of the small-angle assumption and natural modes of the human model is analyzed.Using the coupling system model,the vibration characteristics of the human-seat interaction surface were analyzed.The ride comfort of the high-speed train and human body dynamic performance were analyzed under normal conditions,track geometric irregularities and train meeting conditions.The results showed that the passenger seats in the front and rear rows adjacent to the window had a higher acceleration value than the others.The human backrest and seat pad connection points have higher vibration amplitudes than the car body floor in the human-sensitive frequency range,indicating that using the acceleration values on the floor may underestimate the discomfort of passengers.The ride comfort of high-speed trains diminishes in the presence of track geometric irregularities and when trains pass each other.When the excitation frequency of track geometry irregularities approached the natural frequency of the human-seat-vehicle system,ride comfort in high-speed trains decreased significantly.Moreover,using seat acceleration to evaluate passenger ride comfort overlooks the vibration characteristics of the human body.The transient aerodynamic force generated when the train meets can cause a larger car body roll and lateral motion at 2 Hz,which,in turn,decreases the passenger ride comfort.This study presents a detailed human-seat-vehicle-track coupling system that can reflect a passenger’s dynamic performance under complex operating conditions.
文摘This study proposes a three-dimensional(3D)coupled magneto-electro-elastic problem for the static analysis of multilayered plates embedding piezomagnetic and piezoelectric layers by considering both sensor and actuator configurations.The 3D governing equations for the magneto-electro-elastic static behavior of plates are explicitly show that are made by the three 3D equilibrium equations,the 3D divergence equation for magnetic induction,and the 3D divergence equation for the electric displacement.The proposed solution involves the exponential matrix in the thickness direction and primary variables’harmonic forms in the in-plane ones.A closed-form solution is performed considering simply-supported boundary conditions.Interlaminar continuity conditions are imposed for displacements,magnetic potential,electric potential,transverse shear/normal stresses,transverse normal magnetic induction and transverse normal electric displacement.Therefore,a layerwise approach is adopted.The results section is composed of an assessment part,where the present model is compared to past 3D electro-elastic or magneto-elastic formulations and a new benchmark part.Benchmarks consider sensor and actuator plate configurations for the fully coupled magneto-electro-elastic cases for different thickness ratios.Tabular and graphical results are presented for displacements,stresses,magnetic potential,electric potential,transverse normal magnetic induction and transverse normal electric displacement.For each presented benchmark,magneto-electro-elastic coupling and thickness and material layer effects are discussed in depth.
