Oil shale is characterized by a dense structure,low proportion of pores and fissures,and low permeability.Pore-fracture systems serve as crucial channels for shale oil migration,directly influencing the production eff...Oil shale is characterized by a dense structure,low proportion of pores and fissures,and low permeability.Pore-fracture systems serve as crucial channels for shale oil migration,directly influencing the production efficiency of shale oil resources.Effectively stimulating oil shale reservoirs remains a challenging and active research topic.This investigation employed shale specimens obtained from the Longmaxi Formation.Scanning electron microscopy,fluid injection experiments,and fluid-structure interaction simulations were used to comprehensively analyze structural changes and fluid flow behavior under high temperatures from microscopic to macroscopic scales.Experimental results indicate that the temperature has little effect on the structure and permeability of shale before 300℃.However,there are two threshold temperatures within the range of 300 to 600℃that have significant effects on the structure and permeability of oil shale.The first threshold temperature is between 300 and 400℃,which causes the oil shale porosity,pore-fracture ratio,and permeability begin to increase.This is manifested by the decrease in micropores and mesopores,the increase in macropores,and the formation of a large number of isolated pores and fissures within the shale.The permeability increases but not significantly.The second threshold temperature is between 500 and 600℃,which increases the permeability of oil shale significantly.During this stage,micropores and mesopores are further reduced,and macropores are significantly enlarged.A large number of connected and penetrated pores and fissures are formed.More numerous and thicker streamlines appear inside the oil shale.The experimental results demonstrate that high temperatures significantly alter the microstructure and permeability of oil shale.At the same time,the experimental results can provide a reference for the research of in-situ heating techniques in oil shale reservoir transformation.展开更多
The oceanic mixed layer in the Southern Ocean is characterized by numerous fronts due to the stirring of freshwater influxes arising from ice melting.The interaction of these fronts with winds modulates the evolution ...The oceanic mixed layer in the Southern Ocean is characterized by numerous fronts due to the stirring of freshwater influxes arising from ice melting.The interaction of these fronts with winds modulates the evolution of the mixed layer and affects atmosphere−ocean energy exchanges.However,the underlying mechanism behind the wind-front interaction remains obscure due to a lack of three-dimensional observations of the ocean,particularly in terms of velocities.To address this issue,this study investigates the dynamics of fronts within the mixed layer during a storm by employing a subset of the global submesoscale-permitting simulation,Northeast Weddell Sea Pre-SWOT Level-4 Hourly MITgcm LLC4320 Native Grid 2km Oceanographic Dataset(ROAM_MIZ).We first compare the ROAM_MIZ data to glider data to assess the performance of the model simulation and find that the ROAM_MIZ can,to a large degree,capture sub-mesoscale features within a mixed layer.Subsequent analyses based on a subset of ROAM_MIZ show that lateral density gradients within the mixed layer rapidly decrease during high winds associated with the storm.Down-front winds accelerate this process as the Ekman buoyancy transport responsible for enhancing the instability of the fronts is primarily dominated by horizontal baroclinic components.After the storm,the fronts strengthen again in the presence of weaker winds due to the frontogenesis by the larger-scale strain.Moreover,the non-geostrophic turbulence induces a modification of the relative vorticity,affecting the instability within the mixed layer.These findings offer valuable guidance for the deployment of observational instruments and subsequent analysis,as well as deepen the understanding of air−sea interactions in the Southern Ocean.展开更多
A two-dimensional Reynolds-averaged Navier-Stokes solver is applied to analyze the aerodynamic behavior of the Shock/Boundary-Layer interaction of rocket with a boosted The K-ε turbulence model and a finite volume m...A two-dimensional Reynolds-averaged Navier-Stokes solver is applied to analyze the aerodynamic behavior of the Shock/Boundary-Layer interaction of rocket with a boosted The K-ε turbulence model and a finite volume method in a unstructured body-fitted curvilinear coordinates have been used. The results indicate that the separation and the reattachment occur in the Boundary-Layer of the main rocket because of the shock interaction. The shape of the booster nose effects the flow field obviously. In the case of the hemisphere booster nose the pressure has complicate distributions and the separation is very clear. The distance between the booster and main rocket has the evident effect on the flow field. If the distance is smaller the pressure coefficient is bigger the separation zone even the separation bubble occurs.展开更多
The efficiency and mechanism of an active control device ＂＇Spark Jet＂ and its application in shock-induced separation control are studied using large-eddy simulation in this paper. The base flow is the interaction ...The efficiency and mechanism of an active control device ＂＇Spark Jet＂ and its application in shock-induced separation control are studied using large-eddy simulation in this paper. The base flow is the interaction of an oblique shock-wave generated by 8° wedge and a spatially-developing Ma = 2.3 turbulent boundary layer. The Reynolds number based on the incoming flow property and the boundary layer displacement thickness at the impinging point without shock-wave is 20000. The detailed numerical approaches were presented. The inflow turbulence was generated using the digital filter method to avoid artificial temporal or streamwise periodicity. The , merical results including velocity profile, Reynolds stress profile, skin friction, and wall pressure were sys- tematically validated against the available wind tunnel particle image velocimetry （PIV） measure- ments of the same flow condition. Further study on the control of flow separation due to the strong shock-viscous interaction using an active control actuator ＂＇Spark Jet＇＂ was conducted. The single-pulsed characteristic of the device was obtained and compared with the experiment. Both instantaneous and time-averaged flow fields have shown that the jet flow issuing from the actuator cavity enhances the flow mixing inside the boundary layer, making the boundary layer more resis- tant to flow separation. Skin friction coefficient distribution shows that the separation bubble length is reduced by about 35% with control exerted.展开更多
The interactions between konjac glucomannan（KGM） and soy protein isolate （SPI） were studied with the method of molecular dynamics simulation. Part representative structures segments of KGM and SPI were used as mod...The interactions between konjac glucomannan（KGM） and soy protein isolate （SPI） were studied with the method of molecular dynamics simulation. Part representative structures segments of KGM and SPI were used as mode, and the force-field was FF03. The stability and sites of KGM/SPI interactions in water were researched at 363 K with the following results： the potential energy （EPOT） of the mixed gel dropped, while that of single KGM gel increased. The surface area （SA） of KGM in the mixed system was decreased to 401.41 from 1 267.54 Az, and that of SPI to 484.94 from 1 943.28 A2. The sum potential energy of KGM and soy protein in the mixed system was decreased to -13 402.41 from -5 768.56 kcal mol^-1. The variations of two parameters showed that the stability of compound gel KGM/SPI was improved, which was consistent with the previous studies. The sites of interactions in the mixed gel were the -OH groups on C（2） in KGM mannose and glucose, and the amide linkage group on Histidine, Asparagine and Leucine in SPI. The hydrogen bond was formed directly or indirectly by the bridge of waters.展开更多
In this paper, the authors present airflow field characteristics of human upper airway and soft palate movement attitude during breathing. On the basis of the data taken from the spiral computerized tomography images ...In this paper, the authors present airflow field characteristics of human upper airway and soft palate movement attitude during breathing. On the basis of the data taken from the spiral computerized tomography images of a healthy person and a patient with Obstructive Sleep Apnea-Hypopnea Syndrome （OSAHS）, three-dimensional models of upper airway cavity and soft palate are reconstructed by the method of surface rendering. Numerical simulation is performed for airflow in the upper airway and displacement of soft palate by fluid-structure interaction analysis. The reconstructed threedimensional models precisely preserve the original configuration of upper airways and soft palate. The results of the pressure and velocity distributions in the airflow field are quantitatively determined, and the displacement of soft palate is presented. Pressure gradients of airway are lower for the healthy person and the airflow distribution is quite uniform in the case of free breathing. However, the OSAHS patient remarkably escalates both the pressure and velocity in the upper airway, and causes higher displacement of the soft palate. The present study is useful in revealing pathogenesis and quantitative mutual relationship between configuration and function of the upper airway as well as in diagnosingdiseases related to anatomical structure and function of the upper airway.展开更多
The interactions between konjac glucomannan and carrageenan were studied with the method of molecular dynamics simulation. Part representative structure segments of KGM and two unit structures of κ-carrageenan （Fig...The interactions between konjac glucomannan and carrageenan were studied with the method of molecular dynamics simulation. Part representative structure segments of KGM and two unit structures of κ-carrageenan （Fig. 2） were used as mode, and the force-field was AMBER2. The stability and sites of konjac glucomannan/carrageenan interactions in water were researched at 373 K with the following results： the potential energy （EPOT） of the mixed gel was dropped, while those of single-konjac glucomannan gel and single carrageenan were increased. The surface area （SA） of KGM in the mixed system was decreased to 1002.2A^°^2, and that of carrageenan to 800.9 A^°^2. The variations of two parameters showed that the stability of compound gel konjac glucomannan/carrageenan was improved, which is consistent with the previous studies. The sites of interactions in the mixed gel were the -OH groups on C（2）, C（4） and C（6）, the acetyl group in KGM mannose, and the -OH group on C（6） in carrageenan. The hydrogen bond was formed directly or indirectly by the bridge of waters.展开更多
This article presents a computer simulation of stress distribution around tunnels and interaction between tunnels using an elasto-plastic model. A finite element method using ANSYS software has been used for the analy...This article presents a computer simulation of stress distribution around tunnels and interaction between tunnels using an elasto-plastic model. A finite element method using ANSYS software has been used for the analyses of one and two tunnels at different overburden depths with different separating distances between the tunnels. The results of numerical analyses indicate that stress distribution and stress concentration around the tunnels vary with the overburden depths. It is found that the coefficients of stress concentration for elasto-plastic medium are smaller than those for elastic one by 1.9%. Furthermore, the interaction between the two tunnels rapidly decreases with the increase of separation distance between them. In addition, for quantitatively describing the interaction between the two tunnels, a critical separation distance is introduced. The critical separation distances between the two tunnels at different overburden depths are 8 m, 12 m, and 14 m respectively. This fact is very important and essential for the design of mining tunnels and to ensure safety in tunnel engineering.展开更多
Three-dimensional unsteady Navier-Stokes equations are numerically solved to simulate the aerodynamic interaction of rotor, canard and horizontal tail in hover based on moving chimera grid. The variations of unsteady ...Three-dimensional unsteady Navier-Stokes equations are numerically solved to simulate the aerodynamic interaction of rotor, canard and horizontal tail in hover based on moving chimera grid. The variations of unsteady aerodynamic forces and moments of the canard and horizontal tail with respect to the rotor azimuth are analyzed with the deflection angle set at 0° and 50°, respectively. The pressure map of aerodynamic surfaces and velocity vector distribution of flow field are investigated to get better understanding of the unsteady aerodynamic interaction. The result shows that the canard and horizontal tail present different characteristics under the downwash of the rotor. The canard produces much vertical force loss with low amplitude fluctuation. Contrarily, the horizontal tail, which is within the flow field induced by the down wash of the rotor, produces only less vertical force loss, but the amplitudes of the lift and pitching moment are larger, implying that a potential deflection angle scheme in hover is 50° for the canard and 0° for the horizontal tail.展开更多
Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics h...Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range. By decomposing the subgrid energy transfer and nonlinear interaction into ‘forward’ and ‘backward’ groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.展开更多
An effective computational method is developed for dynamic analysis offluid-structure interaction problems involving large-amplitude sloshing of the fluid andlarge-displacement motion of the structure. The structure i...An effective computational method is developed for dynamic analysis offluid-structure interaction problems involving large-amplitude sloshing of the fluid andlarge-displacement motion of the structure. The structure is modeled as a rigid container supportedby a system consisting of springs and dashpots. The motion of the fluid is decomposed into twoparts: the large-displacement motion with the container and the large-amplitude sloshing relative tothe container. The former is conveniently dealt with by defining a container-fixed noninertiallocal frame, while the latter is easily handled by adopting an ALE kinematical description. Thisleads to an easy and accurate treatment of both the fluid-structure interface and the fluid freesurface without producing excessive distortion of the computational mesh. The coupling between thefluid and the structure is accomplished through the coupling matrices that can be easilyestablished. Two numerical examples, including a TLD-structure system and a simplified liquid-loadedvehicle system, are presented to demonstrate the effectiveness and reliability of the proposedmethod. The present work can also be applied to simulate fluid-structure problems incorporatingmultibody systems and several fluid domains.展开更多
A rate dependent crystal plasticity constitutive model considering self and latent hardening in finite element analysis was developed to simulate rolling textures of pure aluminum. By changing the assignment of orient...A rate dependent crystal plasticity constitutive model considering self and latent hardening in finite element analysis was developed to simulate rolling textures of pure aluminum. By changing the assignment of orientations to finite elements, i.e. assigning the same set of orientations to all elements or different orientations to different elements, the influences of grain interaction on the formation of rolling textures were numerically simulated with this kind of crystal plasticity finite element model. The simulation results reveal that the grains without considering grain interaction rotate faster than those considering grain interaction, and the rotation of grain boundary is slowed down due to the grain interaction. For a good simulation more elements should be assigned to one grain, in which the effects of both the boundary and interior parts of grain contribute to the formation of rolling textures.展开更多
In this work grand canonical Monte Carlo simulations were performed to study gas separation in three pairs of isoreticular metal-organic frameworks (IRMOFs) with and without catenation at room temperature.Mixture comp...In this work grand canonical Monte Carlo simulations were performed to study gas separation in three pairs of isoreticular metal-organic frameworks (IRMOFs) with and without catenation at room temperature.Mixture composed of CO2 and H2 was selected as the model system to separate.The results show that CO2 selectivity in catenated MOFs with multi-porous frameworks is much higher than their non-catenated counterparts.The simulations also show that the electrostatic interactions are very important for the selectivity,and the contributions of different electrostatic interactions are different,depending on pore size,pressure and mixture composition.In fact,changing the electrostatic interactions can even qualitatively change the adsorption behavior.A general conclusion is that the electrostatic interactions between adsorbate molecules and the framework atoms play a dominant role at low pressures,and these interactions in catenated MOFs have much more pronounced effects than those in their non-catenated counterparts,while the electrostatic interactions between adsorbate molecules become evident with increasing pressure,and eventually dominant.展开更多
Micromagnetic simulations have been performed to obtain the dynamic susceptibility spectra of 4×4 cobalt nanowire arrays with different spatial configurations and geometries. The susceptibility spectra of isolate...Micromagnetic simulations have been performed to obtain the dynamic susceptibility spectra of 4×4 cobalt nanowire arrays with different spatial configurations and geometries. The susceptibility spectra of isolated wires have also been simulated for comparison purposes. It is found that the susceptibility spectrum of nanowire array bears a lot of similarities to that of an isolated wire, such as the occurrences of the edge mode and the bulk resonance mode. The simulation results also reveal that the susceptibility spectrum of nanowire array behaves like that of single isolated wire as the interwire distance grows to an extent, which is believed due to the decrease of magnetostatic interaction among nanowires, and can be further confirmed by the static magnetic hysteresis simulations. In comparison with single nanowire, magnetostatic interaction may increase or decrease the resonance frequencies of nanowire arrays assuming a certain interwire distance when the length of array increases. Our simulation results are also analysed by employing the Kittel equation and recent theoretical studies.展开更多
A numerical simulation of the interaction between laminar flow with low Reynolds number and a highly flexible elastic sheet is presented. The mathematical model for the simulation includes a three-dimensional finitevo...A numerical simulation of the interaction between laminar flow with low Reynolds number and a highly flexible elastic sheet is presented. The mathematical model for the simulation includes a three-dimensional finitevolume based fluid solver for incompressible viscous flow and a combined finite-discrete element method for the three-dimensional deformation of solid. An immersed boundary method is used to couple the simulation of fluid and solid. It is implemented through a set of immersed boundary points scattered on the solid surface. These points provide a deformable solid wall boundary for the fluid by adding body force to Navier-Stokes equations. The force from the fluid is also obtained for each point and then applied on the boundary nodes of the solid. The vortex-induced vibration of the highly flexible elastic sheet is simulated with the established mathematical model. The simulated results for both swing pattern and oscillation frequency of the elastic sheet in low Reynolds number flow agree well with experimental data.展开更多
To improve the accuracy and interactivity of soft tissue delormatlon simulation, a new plate spring model based on physics is proposed. The model is parameterized and thus can be adapted to simulate different organs. ...To improve the accuracy and interactivity of soft tissue delormatlon simulation, a new plate spring model based on physics is proposed. The model is parameterized and thus can be adapted to simulate different organs. Different soft tissues are modeled by changing the width, number of pieces, thickness, and length of a single plate spring. In this paper, the structural design, calcula- tion of soft tissue deformation and real-time feedback operations of our system are also introduced. To evaluate the feasibility of the system and validate the model, an experimental system of haptic in- teraction, in which users can use virtual hands to pull virtual brain tissues, is built using PHANTOM OMNI devices. Experimental results show that the proposed system is stable, accurate and promising for modeling instantaneous soft tissue deformation.展开更多
Modulated high frequency （HF） heating of the ionosphere provides a feasible means of artificially generating ex- tremely low frequency （ELF）/very low frequency （VLF） whistler waves, which can leak into the inner...Modulated high frequency （HF） heating of the ionosphere provides a feasible means of artificially generating ex- tremely low frequency （ELF）/very low frequency （VLF） whistler waves, which can leak into the inner magnetosphere and contribute to resonant interactions with high energy electrons. Combining the ray tracing method and test particle simulations, we evaluate the effects of energetic electron resonant scattering driven by the discrete, multi-frequency arti- ficially generated ELF/VLF waves. The simulation results indicate a stochastic behavior of electrons and a linear profile of pitch angle and kinetic energy variations averaged over all test electrons. These features are similar to those associated with single-frequency waves. The computed local diffusion coefficients show that, although the momentum diffusion of relativistic electrons due to artificial ELF/VLF whistlers with a nominal amplitude of ~ 1 pT is minor, the pitch angle scattering can be notably efficient at low pitch angles near the loss cone, which supports the feasibility of artificial triggering of multi-frequency ELF/VLF whistler waves for the removal of high energy electrons from the magnetosphere. We also investigate the dependences of diffusion coefficients on the frequency interval （△f） of the discrete, multi-frequency waves. We find that there is a threshold value of Af for which the net diffusion coefficient of multi-frequency whistlers is inversely proportional to △f （proportional to the frequency components Nw） when △f is below the threshold value but it remains unchanged with increasing Af when △f is larger than the threshold value. This is explained as being due to the fact that the resonant scattering effect of broadband waves is the sum of the effects of each frequency in the ＇effective frequency band＇. Our results suggest that the modulation frequency of HF heating of the ionosphere can be appropriately selected with reasonable frequency intervals so that better performance of controlled precipitation of high energy electrons in the plasmasphere by artificial ELF/VLF whistler waves can be achieved.展开更多
High-order accurate schemes are employed to numerically simulate the interaction of a supersonic jet and a co-directional supersonic inflow. A double backward-facing step model is proposed to investigate the interacti...High-order accurate schemes are employed to numerically simulate the interaction of a supersonic jet and a co-directional supersonic inflow. A double backward-facing step model is proposed to investigate the interaction between the jet shear layer and the supersonic inflow shear layer. It is found that due to the interaction of the shear layer, a secondary jet is injected into the recirculation zone at the intersection of the two shear layers. The secondary jet produced by the interaction of the two shear layers has a periodicity because of shear layers interaction. The distinction in the shape of double backward-facing steps will induce changes in the period of the secondary jet. The analysis and discussion of the periodicity of the secondary jet are mainly focused in this letter.展开更多
Brazil is the world leader in sugarcane production and the largest sugar exporter. Developing new varieties is one of the main factors that contribute to yield increase. In order to select the best genotypes, during t...Brazil is the world leader in sugarcane production and the largest sugar exporter. Developing new varieties is one of the main factors that contribute to yield increase. In order to select the best genotypes, during the final selection stage, varieties are tested in different environments (locations and years), and breeders need to estimate the phenotypic performance for main traits such as tons of cane yield per hectare (TCH) considering the genotype × environment interaction (GEI) effect. Geneticists and biometricians have used different methods and there is no clear consensus of the best method. In this study, we present a comparison of three methods, viz. Eberhart-Russel (ER), additive main effects and multiplicative interaction (AMMI) and mixed model (REML/BLUP), in a simulation study performed in the R computing environment to verify the effectiveness of each method in detecting GEI, and assess the particularities of each method from a statistical standpoint. In total, 63 cases representing different conditions were simulated, generating more than 34 million data points for analysis by each of the three methods. The results show that each method detects GEI differently in a different way, and each has some limitations. All three methods detected GEI effectively, but the mixed model showed higher sensitivity. When applying the GEI analysis, firstly it is important to verify the assumptions inherent in each method and these limitations should be taken into account when choosing the method to be used.展开更多
The finite element method was used to solve fluid dynamic interaction problems between the crust and mantle of the Earth. To consider different mechanical behaviours, the lithosphere consisting of the crust and upper ...The finite element method was used to solve fluid dynamic interaction problems between the crust and mantle of the Earth. To consider different mechanical behaviours, the lithosphere consisting of the crust and upper mantle was simulated as fluid-saturated porous rocks, while the upper aesthenospheric part of the mantle was simulated as viscous fluids. Since the whole lithosphere was computationally simulated, the dynamic interaction between the crust and the upper mantle was appropriately considered. In particular, the mixing of mantle fluids and crustal fluids was simulated in the corresponding computational model. The related computational simulation results from an example problem demonstrate that the mantle fluids can flow into the crust and mix with the crustal fluids due to the resulting convective flows in the crust-mantle system. Likewise, the crustal fluids can also flow into the upper mantle and mix with the mantle fluids. This kind of fluids mixing and exchange is very important to the better understanding of the governing processes that control the ore body formation and mineralization in the upper crust of the Earth.展开更多
基金supported by the Chongqing Natural Science Foundation of Chongqing,China(No.CSTB2022NSCQ-MSX0333)the Science and Technology Research Program of Chongqing Municipal Education Commission(Grant No.KJZD-K202401205)+1 种基金Chongqing Three Gorges University Graduate Research and Innovation Project Funding(No.YJSKY24045)Chongqing Engineering Research Center of Disaster Prevention&Control for Banks and Structures in Three Gorges Reservoir Area(No.SXAPGC24YB14,No.SXAPGC24YB03,No.SXAPGC24YB12)。
文摘Oil shale is characterized by a dense structure,low proportion of pores and fissures,and low permeability.Pore-fracture systems serve as crucial channels for shale oil migration,directly influencing the production efficiency of shale oil resources.Effectively stimulating oil shale reservoirs remains a challenging and active research topic.This investigation employed shale specimens obtained from the Longmaxi Formation.Scanning electron microscopy,fluid injection experiments,and fluid-structure interaction simulations were used to comprehensively analyze structural changes and fluid flow behavior under high temperatures from microscopic to macroscopic scales.Experimental results indicate that the temperature has little effect on the structure and permeability of shale before 300℃.However,there are two threshold temperatures within the range of 300 to 600℃that have significant effects on the structure and permeability of oil shale.The first threshold temperature is between 300 and 400℃,which causes the oil shale porosity,pore-fracture ratio,and permeability begin to increase.This is manifested by the decrease in micropores and mesopores,the increase in macropores,and the formation of a large number of isolated pores and fissures within the shale.The permeability increases but not significantly.The second threshold temperature is between 500 and 600℃,which increases the permeability of oil shale significantly.During this stage,micropores and mesopores are further reduced,and macropores are significantly enlarged.A large number of connected and penetrated pores and fissures are formed.More numerous and thicker streamlines appear inside the oil shale.The experimental results demonstrate that high temperatures significantly alter the microstructure and permeability of oil shale.At the same time,the experimental results can provide a reference for the research of in-situ heating techniques in oil shale reservoir transformation.
基金supported by the National Natural Science Foundation of China (Grant Nos. 42406241,42325604,42227901)the Ministry of Science and Technology of China (No. 2021YFC2803304)+2 种基金the Program of Shanghai Academic/Technology Research Leader (22XD1403600)supported by the Swedish Research Council (Nos. 2020–03190 and 2024-04209)the Swedish Research Council for the Environment, Agricultural Sciences and Spatial Planning (No. 202400375)
文摘The oceanic mixed layer in the Southern Ocean is characterized by numerous fronts due to the stirring of freshwater influxes arising from ice melting.The interaction of these fronts with winds modulates the evolution of the mixed layer and affects atmosphere−ocean energy exchanges.However,the underlying mechanism behind the wind-front interaction remains obscure due to a lack of three-dimensional observations of the ocean,particularly in terms of velocities.To address this issue,this study investigates the dynamics of fronts within the mixed layer during a storm by employing a subset of the global submesoscale-permitting simulation,Northeast Weddell Sea Pre-SWOT Level-4 Hourly MITgcm LLC4320 Native Grid 2km Oceanographic Dataset(ROAM_MIZ).We first compare the ROAM_MIZ data to glider data to assess the performance of the model simulation and find that the ROAM_MIZ can,to a large degree,capture sub-mesoscale features within a mixed layer.Subsequent analyses based on a subset of ROAM_MIZ show that lateral density gradients within the mixed layer rapidly decrease during high winds associated with the storm.Down-front winds accelerate this process as the Ekman buoyancy transport responsible for enhancing the instability of the fronts is primarily dominated by horizontal baroclinic components.After the storm,the fronts strengthen again in the presence of weaker winds due to the frontogenesis by the larger-scale strain.Moreover,the non-geostrophic turbulence induces a modification of the relative vorticity,affecting the instability within the mixed layer.These findings offer valuable guidance for the deployment of observational instruments and subsequent analysis,as well as deepen the understanding of air−sea interactions in the Southern Ocean.
文摘A two-dimensional Reynolds-averaged Navier-Stokes solver is applied to analyze the aerodynamic behavior of the Shock/Boundary-Layer interaction of rocket with a boosted The K-ε turbulence model and a finite volume method in a unstructured body-fitted curvilinear coordinates have been used. The results indicate that the separation and the reattachment occur in the Boundary-Layer of the main rocket because of the shock interaction. The shape of the booster nose effects the flow field obviously. In the case of the hemisphere booster nose the pressure has complicate distributions and the separation is very clear. The distance between the booster and main rocket has the evident effect on the flow field. If the distance is smaller the pressure coefficient is bigger the separation zone even the separation bubble occurs.
基金supported by the National Natural Science Foundation of China(Nos.11302012,51420105008,51476004,11572025 and 51136003)the National Basic Research Program of China(No.2012CB720205)The computational time for the present study was provided by the UK Turbulence Consortium(EPSRC grant EP/L000261/1)
文摘The efficiency and mechanism of an active control device ＂＇Spark Jet＂ and its application in shock-induced separation control are studied using large-eddy simulation in this paper. The base flow is the interaction of an oblique shock-wave generated by 8° wedge and a spatially-developing Ma = 2.3 turbulent boundary layer. The Reynolds number based on the incoming flow property and the boundary layer displacement thickness at the impinging point without shock-wave is 20000. The detailed numerical approaches were presented. The inflow turbulence was generated using the digital filter method to avoid artificial temporal or streamwise periodicity. The , merical results including velocity profile, Reynolds stress profile, skin friction, and wall pressure were sys- tematically validated against the available wind tunnel particle image velocimetry （PIV） measure- ments of the same flow condition. Further study on the control of flow separation due to the strong shock-viscous interaction using an active control actuator ＂＇Spark Jet＇＂ was conducted. The single-pulsed characteristic of the device was obtained and compared with the experiment. Both instantaneous and time-averaged flow fields have shown that the jet flow issuing from the actuator cavity enhances the flow mixing inside the boundary layer, making the boundary layer more resis- tant to flow separation. Skin friction coefficient distribution shows that the separation bubble length is reduced by about 35% with control exerted.
基金supported by the National Natural Science Foundation of China(30871749,30901004)
文摘The interactions between konjac glucomannan（KGM） and soy protein isolate （SPI） were studied with the method of molecular dynamics simulation. Part representative structures segments of KGM and SPI were used as mode, and the force-field was FF03. The stability and sites of KGM/SPI interactions in water were researched at 363 K with the following results： the potential energy （EPOT） of the mixed gel dropped, while that of single KGM gel increased. The surface area （SA） of KGM in the mixed system was decreased to 401.41 from 1 267.54 Az, and that of SPI to 484.94 from 1 943.28 A2. The sum potential energy of KGM and soy protein in the mixed system was decreased to -13 402.41 from -5 768.56 kcal mol^-1. The variations of two parameters showed that the stability of compound gel KGM/SPI was improved, which was consistent with the previous studies. The sites of interactions in the mixed gel were the -OH groups on C（2） in KGM mannose and glucose, and the amide linkage group on Histidine, Asparagine and Leucine in SPI. The hydrogen bond was formed directly or indirectly by the bridge of waters.
基金The project supported by the National Natural Science Foundation of China(10672036,10472025 and 10421002)the Natural Science Foundation of Liaoning Province(20032109)
文摘In this paper, the authors present airflow field characteristics of human upper airway and soft palate movement attitude during breathing. On the basis of the data taken from the spiral computerized tomography images of a healthy person and a patient with Obstructive Sleep Apnea-Hypopnea Syndrome （OSAHS）, three-dimensional models of upper airway cavity and soft palate are reconstructed by the method of surface rendering. Numerical simulation is performed for airflow in the upper airway and displacement of soft palate by fluid-structure interaction analysis. The reconstructed threedimensional models precisely preserve the original configuration of upper airways and soft palate. The results of the pressure and velocity distributions in the airflow field are quantitatively determined, and the displacement of soft palate is presented. Pressure gradients of airway are lower for the healthy person and the airflow distribution is quite uniform in the case of free breathing. However, the OSAHS patient remarkably escalates both the pressure and velocity in the upper airway, and causes higher displacement of the soft palate. The present study is useful in revealing pathogenesis and quantitative mutual relationship between configuration and function of the upper airway as well as in diagnosingdiseases related to anatomical structure and function of the upper airway.
基金supported by the National Natural Science Foundation of China (30371009, 30471218)
文摘The interactions between konjac glucomannan and carrageenan were studied with the method of molecular dynamics simulation. Part representative structure segments of KGM and two unit structures of κ-carrageenan （Fig. 2） were used as mode, and the force-field was AMBER2. The stability and sites of konjac glucomannan/carrageenan interactions in water were researched at 373 K with the following results： the potential energy （EPOT） of the mixed gel was dropped, while those of single-konjac glucomannan gel and single carrageenan were increased. The surface area （SA） of KGM in the mixed system was decreased to 1002.2A^°^2, and that of carrageenan to 800.9 A^°^2. The variations of two parameters showed that the stability of compound gel konjac glucomannan/carrageenan was improved, which is consistent with the previous studies. The sites of interactions in the mixed gel were the -OH groups on C（2）, C（4） and C（6）, the acetyl group in KGM mannose, and the -OH group on C（6） in carrageenan. The hydrogen bond was formed directly or indirectly by the bridge of waters.
文摘This article presents a computer simulation of stress distribution around tunnels and interaction between tunnels using an elasto-plastic model. A finite element method using ANSYS software has been used for the analyses of one and two tunnels at different overburden depths with different separating distances between the tunnels. The results of numerical analyses indicate that stress distribution and stress concentration around the tunnels vary with the overburden depths. It is found that the coefficients of stress concentration for elasto-plastic medium are smaller than those for elastic one by 1.9%. Furthermore, the interaction between the two tunnels rapidly decreases with the increase of separation distance between them. In addition, for quantitatively describing the interaction between the two tunnels, a critical separation distance is introduced. The critical separation distances between the two tunnels at different overburden depths are 8 m, 12 m, and 14 m respectively. This fact is very important and essential for the design of mining tunnels and to ensure safety in tunnel engineering.
文摘Three-dimensional unsteady Navier-Stokes equations are numerically solved to simulate the aerodynamic interaction of rotor, canard and horizontal tail in hover based on moving chimera grid. The variations of unsteady aerodynamic forces and moments of the canard and horizontal tail with respect to the rotor azimuth are analyzed with the deflection angle set at 0° and 50°, respectively. The pressure map of aerodynamic surfaces and velocity vector distribution of flow field are investigated to get better understanding of the unsteady aerodynamic interaction. The result shows that the canard and horizontal tail present different characteristics under the downwash of the rotor. The canard produces much vertical force loss with low amplitude fluctuation. Contrarily, the horizontal tail, which is within the flow field induced by the down wash of the rotor, produces only less vertical force loss, but the amplitudes of the lift and pitching moment are larger, implying that a potential deflection angle scheme in hover is 50° for the canard and 0° for the horizontal tail.
文摘Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range. By decomposing the subgrid energy transfer and nonlinear interaction into ‘forward’ and ‘backward’ groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.
基金This project is supported by National 863 Hi-Tech Project Foundation (No. 2002AA411030).
文摘An effective computational method is developed for dynamic analysis offluid-structure interaction problems involving large-amplitude sloshing of the fluid andlarge-displacement motion of the structure. The structure is modeled as a rigid container supportedby a system consisting of springs and dashpots. The motion of the fluid is decomposed into twoparts: the large-displacement motion with the container and the large-amplitude sloshing relative tothe container. The former is conveniently dealt with by defining a container-fixed noninertiallocal frame, while the latter is easily handled by adopting an ALE kinematical description. Thisleads to an easy and accurate treatment of both the fluid-structure interface and the fluid freesurface without producing excessive distortion of the computational mesh. The coupling between thefluid and the structure is accomplished through the coupling matrices that can be easilyestablished. Two numerical examples, including a TLD-structure system and a simplified liquid-loadedvehicle system, are presented to demonstrate the effectiveness and reliability of the proposedmethod. The present work can also be applied to simulate fluid-structure problems incorporatingmultibody systems and several fluid domains.
基金Projects(50230310 ,50301016) supported by the National Natural Science Foundation of China project(2004053304)supported by the Doctor Program Foundation of the Ministry of Education of China project(2005CB623706) supported by the State KeyFundamental Research and Development Programof China
文摘A rate dependent crystal plasticity constitutive model considering self and latent hardening in finite element analysis was developed to simulate rolling textures of pure aluminum. By changing the assignment of orientations to finite elements, i.e. assigning the same set of orientations to all elements or different orientations to different elements, the influences of grain interaction on the formation of rolling textures were numerically simulated with this kind of crystal plasticity finite element model. The simulation results reveal that the grains without considering grain interaction rotate faster than those considering grain interaction, and the rotation of grain boundary is slowed down due to the grain interaction. For a good simulation more elements should be assigned to one grain, in which the effects of both the boundary and interior parts of grain contribute to the formation of rolling textures.
基金Supported by the National Natural Science Foundation of China （20725622, 20706002, and 20876006）, Beijing Nova Program （2008B15） and the Dutch STW/CW Separation Technology Program （700.56.655-DPC.6243）.
文摘In this work grand canonical Monte Carlo simulations were performed to study gas separation in three pairs of isoreticular metal-organic frameworks (IRMOFs) with and without catenation at room temperature.Mixture composed of CO2 and H2 was selected as the model system to separate.The results show that CO2 selectivity in catenated MOFs with multi-porous frameworks is much higher than their non-catenated counterparts.The simulations also show that the electrostatic interactions are very important for the selectivity,and the contributions of different electrostatic interactions are different,depending on pore size,pressure and mixture composition.In fact,changing the electrostatic interactions can even qualitatively change the adsorption behavior.A general conclusion is that the electrostatic interactions between adsorbate molecules and the framework atoms play a dominant role at low pressures,and these interactions in catenated MOFs have much more pronounced effects than those in their non-catenated counterparts,while the electrostatic interactions between adsorbate molecules become evident with increasing pressure,and eventually dominant.
基金Project supported by the National Natural Science Foundation of China(Grant No.60701016)the National Natural Science Foundation of China-the Royal Society of United Kingdom International Jointed Project(Grant No.60911130130)
文摘Micromagnetic simulations have been performed to obtain the dynamic susceptibility spectra of 4×4 cobalt nanowire arrays with different spatial configurations and geometries. The susceptibility spectra of isolated wires have also been simulated for comparison purposes. It is found that the susceptibility spectrum of nanowire array bears a lot of similarities to that of an isolated wire, such as the occurrences of the edge mode and the bulk resonance mode. The simulation results also reveal that the susceptibility spectrum of nanowire array behaves like that of single isolated wire as the interwire distance grows to an extent, which is believed due to the decrease of magnetostatic interaction among nanowires, and can be further confirmed by the static magnetic hysteresis simulations. In comparison with single nanowire, magnetostatic interaction may increase or decrease the resonance frequencies of nanowire arrays assuming a certain interwire distance when the length of array increases. Our simulation results are also analysed by employing the Kittel equation and recent theoretical studies.
基金Supported by Marie Curie International Incoming Fellowship (No. PIIF-GA-2009-253453)
文摘A numerical simulation of the interaction between laminar flow with low Reynolds number and a highly flexible elastic sheet is presented. The mathematical model for the simulation includes a three-dimensional finitevolume based fluid solver for incompressible viscous flow and a combined finite-discrete element method for the three-dimensional deformation of solid. An immersed boundary method is used to couple the simulation of fluid and solid. It is implemented through a set of immersed boundary points scattered on the solid surface. These points provide a deformable solid wall boundary for the fluid by adding body force to Navier-Stokes equations. The force from the fluid is also obtained for each point and then applied on the boundary nodes of the solid. The vortex-induced vibration of the highly flexible elastic sheet is simulated with the established mathematical model. The simulated results for both swing pattern and oscillation frequency of the elastic sheet in low Reynolds number flow agree well with experimental data.
基金Supported by the National High Technology Research and Development Programme of China(No.2013AA010803,2009AA01Z311,2009AA01Z314)the National Natural Science Foundation of China(No.61304205,61203316,61272379,61103086,41301037)+3 种基金the Natural Science Foundation of Jiangsu Province(BK20141002)the Open Funding Project of State Key Laboratory of Virtual Reality Technology and Systems,Beihang University,Jiangsu Ordinary University Science Research Project(No.13KJB120007)Innovation and Entrepreneurship Training Project of College Students(No.201410300153,201410300165)the Excellent Undergraduate Paper(design)Supporting Project of NUIST
文摘To improve the accuracy and interactivity of soft tissue delormatlon simulation, a new plate spring model based on physics is proposed. The model is parameterized and thus can be adapted to simulate different organs. Different soft tissues are modeled by changing the width, number of pieces, thickness, and length of a single plate spring. In this paper, the structural design, calcula- tion of soft tissue deformation and real-time feedback operations of our system are also introduced. To evaluate the feasibility of the system and validate the model, an experimental system of haptic in- teraction, in which users can use virtual hands to pull virtual brain tissues, is built using PHANTOM OMNI devices. Experimental results show that the proposed system is stable, accurate and promising for modeling instantaneous soft tissue deformation.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.41204120 and 41304130)the Fundamental Research Funds for the Central Universities(Grant No.2042014kf0251)
文摘Modulated high frequency （HF） heating of the ionosphere provides a feasible means of artificially generating ex- tremely low frequency （ELF）/very low frequency （VLF） whistler waves, which can leak into the inner magnetosphere and contribute to resonant interactions with high energy electrons. Combining the ray tracing method and test particle simulations, we evaluate the effects of energetic electron resonant scattering driven by the discrete, multi-frequency arti- ficially generated ELF/VLF waves. The simulation results indicate a stochastic behavior of electrons and a linear profile of pitch angle and kinetic energy variations averaged over all test electrons. These features are similar to those associated with single-frequency waves. The computed local diffusion coefficients show that, although the momentum diffusion of relativistic electrons due to artificial ELF/VLF whistlers with a nominal amplitude of ~ 1 pT is minor, the pitch angle scattering can be notably efficient at low pitch angles near the loss cone, which supports the feasibility of artificial triggering of multi-frequency ELF/VLF whistler waves for the removal of high energy electrons from the magnetosphere. We also investigate the dependences of diffusion coefficients on the frequency interval （△f） of the discrete, multi-frequency waves. We find that there is a threshold value of Af for which the net diffusion coefficient of multi-frequency whistlers is inversely proportional to △f （proportional to the frequency components Nw） when △f is below the threshold value but it remains unchanged with increasing Af when △f is larger than the threshold value. This is explained as being due to the fact that the resonant scattering effect of broadband waves is the sum of the effects of each frequency in the ＇effective frequency band＇. Our results suggest that the modulation frequency of HF heating of the ionosphere can be appropriately selected with reasonable frequency intervals so that better performance of controlled precipitation of high energy electrons in the plasmasphere by artificial ELF/VLF whistler waves can be achieved.
基金supported by the National Key Research and Development Program of China(Grant 2016YFA0401201)the National Natural Science Foundation of China(Grants 11872066,11472281,11727901,and 11532014)。
文摘High-order accurate schemes are employed to numerically simulate the interaction of a supersonic jet and a co-directional supersonic inflow. A double backward-facing step model is proposed to investigate the interaction between the jet shear layer and the supersonic inflow shear layer. It is found that due to the interaction of the shear layer, a secondary jet is injected into the recirculation zone at the intersection of the two shear layers. The secondary jet produced by the interaction of the two shear layers has a periodicity because of shear layers interaction. The distinction in the shape of double backward-facing steps will induce changes in the period of the secondary jet. The analysis and discussion of the periodicity of the secondary jet are mainly focused in this letter.
文摘Brazil is the world leader in sugarcane production and the largest sugar exporter. Developing new varieties is one of the main factors that contribute to yield increase. In order to select the best genotypes, during the final selection stage, varieties are tested in different environments (locations and years), and breeders need to estimate the phenotypic performance for main traits such as tons of cane yield per hectare (TCH) considering the genotype × environment interaction (GEI) effect. Geneticists and biometricians have used different methods and there is no clear consensus of the best method. In this study, we present a comparison of three methods, viz. Eberhart-Russel (ER), additive main effects and multiplicative interaction (AMMI) and mixed model (REML/BLUP), in a simulation study performed in the R computing environment to verify the effectiveness of each method in detecting GEI, and assess the particularities of each method from a statistical standpoint. In total, 63 cases representing different conditions were simulated, generating more than 34 million data points for analysis by each of the three methods. The results show that each method detects GEI differently in a different way, and each has some limitations. All three methods detected GEI effectively, but the mixed model showed higher sensitivity. When applying the GEI analysis, firstly it is important to verify the assumptions inherent in each method and these limitations should be taken into account when choosing the method to be used.
基金Project(10872219) supported by the National Natural Science Foundation of China
文摘The finite element method was used to solve fluid dynamic interaction problems between the crust and mantle of the Earth. To consider different mechanical behaviours, the lithosphere consisting of the crust and upper mantle was simulated as fluid-saturated porous rocks, while the upper aesthenospheric part of the mantle was simulated as viscous fluids. Since the whole lithosphere was computationally simulated, the dynamic interaction between the crust and the upper mantle was appropriately considered. In particular, the mixing of mantle fluids and crustal fluids was simulated in the corresponding computational model. The related computational simulation results from an example problem demonstrate that the mantle fluids can flow into the crust and mix with the crustal fluids due to the resulting convective flows in the crust-mantle system. Likewise, the crustal fluids can also flow into the upper mantle and mix with the mantle fluids. This kind of fluids mixing and exchange is very important to the better understanding of the governing processes that control the ore body formation and mineralization in the upper crust of the Earth.
