Flow boiling in open microchannels offers highly efficient heat transfer performance and has attracted increasing attention in the fields of heat transfer and thermalmanagement of electronic devices in recent years.Ho...Flow boiling in open microchannels offers highly efficient heat transfer performance and has attracted increasing attention in the fields of heat transfer and thermalmanagement of electronic devices in recent years.However,the continuous rise in power density of electronic components imposesmore stringent requirements on the heat transfer capability of microchannel flow boiling.HFE-7100,a dielectric coolant with favorable thermophysical properties,has become a focal point of research for enhancing flow boiling performance in open microchannels.The flow boiling heat transfer performance ofHFE-7100 was investigated in this study by fabricating micro-nano composite structures on the bottom surface of open microchannels using laser ablation technology.Based on visualization results,a comparative analysis was conducted on the bubble dynamics and flow pattern characteristics of HFE-7100 flow boiling in micronano structured open microchannels(MNSOMC)and smooth-surface open microchannels(SSOMC),to elucidate the enhancement mechanism of micro-nano structures on flow boiling heat transfer in open microchannels.The results indicate that the surface structures and strong wettability of MNSOMC accelerated bubble nucleation and departure.Moreover,bubbles in the channel tended to coalesce along the flow direction,forming elongated slug bubbles with high aspect ratios,which enabled efficient thin film evaporation in conjunction with intense nucleate boiling,thereby significantly enhancing flow boiling heat transfer.Under the experimental conditions of this study,the maximum enhancements in the heat transfer coefficient(HTC)and critical heat flux(CHF)of HFE-7100 inMNSOMC were 33.4%and 133.1%,respectively,with the CHF reaching up to 1542.3 kW⋅m^(−2).Furthermore,due to the superior wettability and capillary wicking capability of the micro-nano composite structures,the significant enhancement in flow boiling heat transfer was achieved without incurring a noticeable pressure drop penalty.展开更多
Poly(_(L)-lactide)(PLLA),a leading biodegradable polyester,has demonstrated potential as a sustainable alternative,owing to its excellent biodegradability and rigidity.However,their slow crystallization kinetics and p...Poly(_(L)-lactide)(PLLA),a leading biodegradable polyester,has demonstrated potential as a sustainable alternative,owing to its excellent biodegradability and rigidity.However,their slow crystallization kinetics and poor heat resistance limit their application scope.Recent advances have highlighted that the combination of extensional flow and thermal fields can achieve toughness–stiffness balance,high transparency,and good heat resistance.However,the effect of extensional flow on the post-non-isothermal crystallization of PLLA during heating and the resulting crystalline texture remains unclear.In this study,PLLA with a heterogeneous amorphous structure and oriented polymorph was prepared by extensional flow.The effect of heterogeneous amorphous structures on non-isothermal crystallization kinetics during the heating process was studied by thermal analysis,polarized optical microscopy,infrared spectroscopy,and ex situ/in situ X-ray characterization.These results clearly illustrate that extensional flow enhances the formation of oriented crystalline structures,accelerates non-isothermal crystallization,and modulates the polymorphic composition of PLLA.Moreover,an unexpected dual cold-crystallization behavior is identified in ordered PLLA samples upon extensional flow,which is from the extensional flow-induced heterogeneous amorphous phase into α' phase(low-temperature peak)and the pristine amorphous phase intoαphase(high-temperature peak).The extensional flow primarily promotes the formation of the more perfectαandα'phases,but has a negative effect on the final content ofαphase formed after cold crystallization andα'-to-αphase transformation.The findings of this work advance the understanding of PLLA non-isothermal crystallization after extensional flow and offer valuable guidance for high-performance PLLA upon heat treatment in practical processing.展开更多
The sluggish electrochemical catalytic activity of the graphite felt electrodes for anode reaction is still a barrier for achieving high-performance vanadium redox flow battery(VRFB).It is significant to leverage the ...The sluggish electrochemical catalytic activity of the graphite felt electrodes for anode reaction is still a barrier for achieving high-performance vanadium redox flow battery(VRFB).It is significant to leverage the exceptional conductivity,excellent electrocatalytic activity,and structural tunability of MXene to address this issue.Herein,this work introduces nitrogen atoms to modulate the carbon layer structure of Ti_(3)C_(2)T_(x)MXene,inducing a reconfiguration of the local electronic structure,which enhances the anode interface activity and thereby improves the performance of VRFB.Ti_(3)C_(2)T_(x)exhibits high conductivity,excellent hydrophilicity,and a large specific surface area,providing excellent interface characteristics for V^(3+)/V^(2+)redox reaction.Moreover,interlayer treatment to modulate the mesoporous structure of MXene further increases the reactive surface area.Importantly,doping nitrogen atoms at carbon layer induces lattice distortions in Ti_(3)C_(2)T_(x),which enhances the charge transfer processes of the V^(3+)/V^(2+)redox reaction.The catalysis mechanism is also validated through density functional theory.Furthermore,the modified graphite felt electrode,as the anode of VRFB,relieves a higher energy efficiency of 68%at 250 mA cm^(-2),while the pristine electrode cannot operate at this current density.In addition,at 150 mA cm^(-2),the modified battery maintains energy efficiency at 75%without degradation after 500 cycles.This study utilizes rational atomic-level engineering for effective structural modulation to significantly enhance the catalytic activity of electrode reaction,offering a unique perspective for developing high-performance MXene electrocatalysts of VRFB.展开更多
An experimental system for forming a rotating paraboloid shaped shallow water with a free surface was conducted to study coherent vortex structures in a differentially rotating quasi two dimensional zonal flow.Flow...An experimental system for forming a rotating paraboloid shaped shallow water with a free surface was conducted to study coherent vortex structures in a differentially rotating quasi two dimensional zonal flow.Flow visualization and laser light scattering techniques were used to obtain the information of spatial flow patterns.Experimental results show that the coexistence of Coriolis effect and strong shear in latitudinal zones may lead to formation of coherent vortices.Power spectra analysis and photographs which were taken in a reference frame rotating with the observed vortices also justified the emergence,drift and evolution of persistent vortices on the large scale.Locked vortex state manifests the cyclone and anticyclone asymmetry.展开更多
Inlet recirculation is proved as an effective way for centrifugal compressor surge margin extension,and is successively used in some engineering applications.Unfortunately its working mechanism is still not being well...Inlet recirculation is proved as an effective way for centrifugal compressor surge margin extension,and is successively used in some engineering applications.Unfortunately its working mechanism is still not being well understood,which leads to redesigning of inlet recirculation mostly by experience.Also,most study about inlet recirculation is steady to date.It is necessary to study surge margin extension mechanism about inlet recirculation.To expose the mechanism in detail,steady and unsteady numerical simulations were performed on a centrifugal compressor with and without inlet recirculation.The results showed that,with inlet recirculation,the inlet axial velocity is augmented,relative Mach number around blade tip leading edge area is significantly reduced and so is the flow angle.As the flow angle decreased,the incidence angle reduced which greatly improves the flow field inside the impeller.Moreover,inlet recirculation changes the blade loading around blade tip and restrains the flow separation on the blade suction side at the leading edge area.The unsteady results of static pressure around blade surface,entropy at inlet crossflow section and vorticity distributions at near tip span surface indicated that,at near stall condition,strong fluctuation exists in the vicinity of tip area due to the interaction between tip leakage flow and core flow.By inlet recirculation these strong flow fluctuations are eliminated so the flow stability is greatly enhanced.All these improvements mentioned above are the reason for inlet recirculation delays compressor stall.This research reveals the surge margin extension reason of inlet recirculation from an unsteady flow viewpoint and provides important reference for inlet recirculation structure design.展开更多
The theoretical research on the propulsive principle of aquatic animal becomes more important and attracted more researchers to make efforts on it. In the present study, a computational fluid dynamic (CFD) simulatio...The theoretical research on the propulsive principle of aquatic animal becomes more important and attracted more researchers to make efforts on it. In the present study, a computational fluid dynamic (CFD) simulation of a three-dimensional traveling-wave undulations body of tuna has been developed to investigate the fluid flow features and vorticity structures around this body when moving in a straight line. The undulation only takes place in the posterior half of the fish, and the tuna-tail is considered as a lunate fin oscillating with the mode combined swaying with yawing. A Reynolds-averaged Navier-Stokes (RANS) equation is developed, employing a control-volume method and a k-omega SST turbulent model; meanwhile an unstructured tetrahedral grid, which is generated for the three-dimensional geometry, is used based on the deformation of the hind parts of the body and corresponding movement of the tail. We calculated the hydrodynamic performance of tuna-like body when a tuna swims in a uniform velocity, and compared the input power coefficient, output power coefficient and propulsive efficiency of the oscillating tuna-tail with or without body vortex shedding. Additionally, the load distribution on the body, flow features and vorticity structures around the body were demonstrated. The effect of interaction between the body-generated vortices and the tail-generated vorticity on the hydrodynamic performance can be obtained.展开更多
Implementation of an opposing jet in design of a hypersonic blunt body significantly modifies the external flowfield and yields a considerable reduction in the aerodynamic drag.This study aims to investigate the effec...Implementation of an opposing jet in design of a hypersonic blunt body significantly modifies the external flowfield and yields a considerable reduction in the aerodynamic drag.This study aims to investigate the effects of flowfield modeling parameters of injection and freestream on the flow structure and aerodynamics of a blunt body with an opposing jet in hypersonic flow.Reynolds-Averaged Navier-Stokes(RANS)equations with a Shear Stress Transport(SST)turbulence model are employed to simulate the intricate jet flow interaction.Through utilizing a Non-Intrusive Polynomial Chaos(NIPC)method to construct surrogates,a functional relation is established between input modeling parameters and output flowfield and aerodynamic quantities in concern.Sobol indices in sensitivity analysis are introduced to represent the relative contribution of each parameter.It is found that variations in modeling parameters produce large variations in the flow structure and aerodynamics.The jet-to-freestream total-pressure ratio,jet Mach number,and freestream Mach number are the major contributors to variation in surface pressure,demonstrating an evident location-dependent behavior.The penetration length of injection,reattachment angle of the shear layer,and aerodynamic drag are also most sensitive to the three crucial parameters above.In comparison,the contributions of freestream temperature,freestream density,and jet total temperature are nearly negligible.展开更多
According to the analysis of the turbulent intensity level around the high-speed train, the maximum turbulent intensity ranges from 0.2 to 0.5 which belongs to high turbulent flow. The flow field distribution law was ...According to the analysis of the turbulent intensity level around the high-speed train, the maximum turbulent intensity ranges from 0.2 to 0.5 which belongs to high turbulent flow. The flow field distribution law was studied and eight types of flow regions were proposed. They are high pressure with air stagnant region, pressure decreasing with air accelerating region, low pressure with high air flow velocity region I, turbulent region, steady flow region, low pressure with high air flow velocity region II,pressure increasing with air decelerating region and wake region. The analysis of the vortex structure around the train shows that the vortex is mainly induced by structures with complex mutation and large curvature change. The head and rear of train, the underbody structure, the carriage connection section and the wake region are the main vortex generating sources while the train body with even cross-section has rare vortexes. The wake structure development law studied lays foundation for the train drag reduction.展开更多
A series of numerical analyses have been performed to investigate the flow structures in a narrow confined channel with 12 staggered circular impingement holes and one bigger exit hole. The flow enters the channel thr...A series of numerical analyses have been performed to investigate the flow structures in a narrow confined channel with 12 staggered circular impingement holes and one bigger exit hole. The flow enters the channel through the impingement holes and exits through the far end outlet. The flow fields corresponding to two jet Reynolds numbers (25000 and 65000) and three channel con- figurations with different ratios of the channel height to the impingement hole diameter (Zr 1, 3, 5) are analyzed by solving the Reynolds averaged Navier-Stokes equations with the realizable k-e turbulence model. Detailed flow field information including the secondary flow, the interaction between the jets and the cross flow, and flow distribution along the channel has been obtained. Comparisons between the numerical and experimental results of the flow fields at the four planes along the channel are performed to validate the numerical method. The calculated impingement pattern, high velocity flow distribution, low velocity separation region and vortices are in good agreement with the experimental data, implying the validity and effectiveness of the employed numerical approach for analyzing relevant flow field.展开更多
Work on dynamic topology optimization of engineering structures for vibration suppression has mainly addressed the maximization of eigenfrequencies and gaps between consecutive eigenfrequencies of free vibration, mini...Work on dynamic topology optimization of engineering structures for vibration suppression has mainly addressed the maximization of eigenfrequencies and gaps between consecutive eigenfrequencies of free vibration, minimization of the dynamic compliance subject to forced vibration, and minimization of the structural frequency response. A dynamic topology optimization method of bi-material plate structures is presented based on power flow analysis. Topology optimization problems formulated directly with the design objective of minimizing the power flow response are dealt with. In comparison to the displacement or velocity response, the power flow response takes not only the amplitude of force and velocity into account, but also the phase relationship of the two vector quantities. The complex expression of power flow response is derived based on time-harmonic external mechanical loading and Rayleigh damping. The mathematical formulation of topology optimization is established based on power flow response and bi-material solid isotropic material with penalization(SIMP) model. Computational optimization procedure is developed by using adjoint design sensitivity analysis and the method of moving asymptotes(MMA). Several numerical examples are presented for bi-material plate structures with different loading frequencies, which verify the feasibility and effectiveness of this method. Additionally, optimum results between topological design of minimum power flow response and minimum dynamic compliance are compared, showing that the present method has strong adaptability for structural dynamic topology optimization problems. The proposed research provides a more accurate and effective approach for dynamic topology optimization of vibrating structures.展开更多
In a Mach 3.8 wind tunnel, both instantaneous and time-averaged flow structures of different scales around a blunt double-cone with or without supersonic film cooling were visualized via nano-tracer planar laser scatt...In a Mach 3.8 wind tunnel, both instantaneous and time-averaged flow structures of different scales around a blunt double-cone with or without supersonic film cooling were visualized via nano-tracer planar laser scattering (NPLS), which has a high spatiotemporal resolution. Three experimental cases with different injection mass flux rates were carried out. Many typical flow structures were clearly shown, such as shock waves, expansion fans, shear layers, mixing layers, and turbulent boundary layers. The analysis of two NPLS images with an interval of 5 us revealed the temporal evolution characteristics of flow structures. With matched pressures, the laminar length of the mixing layer was longer than that in the case with a larger mass flux rate, but the full covered region was shorter. Structures like K-H (Kelvin-Helmholtz) vortices were clearly seen in both flows. Without injection, the flow was similar to the supersonic flow over a backward- facing step, and the structures were relatively simpler, and there was a longer laminar region. Large scale structures such as hairpin vortices were visualized. In addition, the results were compared in part with the schlieren images captured by others under similar conditions.展开更多
This study investigates the heterogeneous structure and its influence on drag coefficient for concurrent-up gas-solid flow. The energy-minimization multi-scale (EMMS) model is modified to simulate the variation of str...This study investigates the heterogeneous structure and its influence on drag coefficient for concurrent-up gas-solid flow. The energy-minimization multi-scale (EMMS) model is modified to simulate the variation of structure parameters with solids concentration, showing the tendency for particles to aggregate to form clusters and for fluid to pass around clusters. The global drag coefficient is resolved into that for the dense phase, for the dilute phase and for the so-called inter-phase, all of which can be obtained from their respective phase-specific structure parameters. The computational results show that the drag coefficients of the different phases are quite different, and the global drag coefficient calculated from the EMMS approach is much lower than that from the correlation of Wen and Yu. The simulation results demonstrate that the EMMS approach can well describe the heterogeneous flow structure, and is very promising for incorporation into the two-fluid model or the discrete particle model as the closure law for drag coefficient.展开更多
Pipe-in-pipe(PIP)structures are widely used in offshore oil and gas pipelines to settle thermal insulation issues.A PIP structure system usually consists of two concentric pipes and one softer layer for thermal insula...Pipe-in-pipe(PIP)structures are widely used in offshore oil and gas pipelines to settle thermal insulation issues.A PIP structure system usually consists of two concentric pipes and one softer layer for thermal insulation consideration.The total response of the system is related to the dynamics of both pipes and the interactions between these two concentric pipes.In the current work,a theoretical model for flow-induced vibrations of a PIP structure system is proposed and analyzed in the presence of an internal axial flow and an external cross flow.The interactions between the two pipes are modeled by a linear distributed damper,a linear distributed spring and a nonlinear distributed spring along the pipe length.The unsteady hydrodynamic forces due to cross flow are modeled by two distributed van der Pol wake oscillators.The nonlinear partial differential equations for the two pipes and the wake are further discretized by the aid of Galerkin’s technique,resulting in a set of ordinary differential equations.These ordinary differential equations are further numeri cally solved by using a fourth-order Runge-Kutta integration algorithm.Phase portraits,bifurcation diagrams,an Argand diagram and oscillation shape diagrams are plotted,showing the existence of a lock-in phenomenon and figure-of-eight trajectory.The PIP system subjected to cross flow displays some interesting dynamical behaviors different from that of a single-pipe structure.展开更多
The structure of wind-sand flow under different total sand transport rates was measured with field vertical anemometer and sand trap on the crest of typical coastal transverse ridge in Changli Gold Coast of Hebei Prov...The structure of wind-sand flow under different total sand transport rates was measured with field vertical anemometer and sand trap on the crest of typical coastal transverse ridge in Changli Gold Coast of Hebei Province, which is one of the most typical coastal aeolian distribution regions in China and famous for the tall and typical coastal transverse ridges. The measurement results show that, on the conditions of approximate wind velocities and same surface materials and environments, some changes happen to the structure of wind-sand flow with the increase of total sand transport rate on the crest of coastal transverse ridge. First, the sand transport rates of layers at different heights in the wind-sand flow increase, with the maximum increase at the height layer of 4-8cm. Second, the ratios of sand trans-port rates of layers at different heights to total sand transport rate decrease at the low height layer (0-4cm), but increase at the high height layer (4-60cm). Third, the distribution of the sand transport rate in the wind-sand flow can be expressed by an exponential function at the height layer of 0-40cm, but it changes from power function model to ex-ponential function model in the whole height layer (0-60cm) and changes into polynomial function model at the height layer of 40-60cm with the increase of total sand transport rate. Those changes have a close relationship with the limit of sand grain size of wind flow transporting and composition of sand grain size in the wind-sand flow.展开更多
The three-dimensional internal flow field of centrifugal pump is complex and variable with design parameters and operation conditions. The post-processing technique named differential amplification method was proposed...The three-dimensional internal flow field of centrifugal pump is complex and variable with design parameters and operation conditions. The post-processing technique named differential amplification method was proposed for the comparison study of different flow structures. The full steady flow fields of an industrial centrifugal pump working on-design and off-design points were numerically investigated by solving Reynolds average Navier-Stokes equations together with a shear-stress transport(SST) k-? turbulence model. And the numerically predicted performance curves of the studied pump agree well with test measurement results. Compared with the flow flied on design point under the help of differential amplification method, the disturbance caused by interaction between blade and volute tongue is very obvious and it extends to the diffuser pipe on the working point with 0.8 times rated flux. While on the point with 1.2 times rated flux, the flow distribution in impeller region is roughly even and it flows more to the bottom section of the diffuser pipe. The above method was proved to be good at displaying the subtle secondary flow structure changes with a higher resolution effect relative to single isolated case observation, which helps the optimization decision-making from multiple design cases.展开更多
The objective of this paper is to address the transient flow structures around a pitching hydrofoil by com- bining physical and numerical studies. In order to predict the dynamic behavior of the flow structure effecti...The objective of this paper is to address the transient flow structures around a pitching hydrofoil by com- bining physical and numerical studies. In order to predict the dynamic behavior of the flow structure effectively, the Lagrangian coherent structures (LCS) defined by the ridges of the finite-time Lyapunov exponent (FTLE) are utilized under the framework of Navier-Stokes flow computations. In the numerical simulations, the k-w shear stress trans- port (SST) turbulence model, coupled with a two-equation F-Reo transition model, is used for the turbulence closure. Results are presented for a NACA66 hydrofoil undergoing slowly and rapidly pitching motions from 0° to 15° then back to 0° at a moderate Reynolds number Re = 7.5 × 105. The results reveal that the transient flow structures can be observed by the LCS method. For the slowly pitching case, it consists of five stages: quasi-steady and laminar, transition from laminar to turbulent, vortex development, large-scale vortex shedding, and reverting to laminar. The observation of LCS and Lagrangian particle tracers elucidates that the trailing edge vortex is nearly attached and stable during the vortex development stage and the interaction between the leading and trailing edge vortex caused by the adverse pres- sure gradient forces the vortexes to shed downstream during the large-scale vortex shedding stage, which corresponds to obvious fluctuations of the hydrodynamic response. For the rapidly pitching case, the inflection is hardly to be observed and the stall is delayed. The vortex formation, interaction, and shedding occurred once instead of being repeated three times, which is responsible for just one fluctuation in the hydrody- namic characteristics. The numerical results also show that the FTLE field has the potential to identify the transient flows, and the LCS can represent the divergence extent of infinite neighboring particles and capture the interface of the vortex region.展开更多
This research investigates the aerodynamic performance and flow characteristics of a delta wing with 65° sweep angle and with coarse axial riblets,and then compares with that of a smooth-surface delta wing.Partic...This research investigates the aerodynamic performance and flow characteristics of a delta wing with 65° sweep angle and with coarse axial riblets,and then compares with that of a smooth-surface delta wing.Particle Image Velocimetry(PIV)were utilized to visualize the flow over the wing at 6 cross-sections upright to the wing surface and parallel to the wing span,as well as 3 longitudinal sections on the leading edge,symmetry plane,and a plane between them at Angles of Attack(AOA)=20°and 30°and Re=1.2×10~5,2.4×10~5,and 3.6×10~5.The effects of the riblets were studied on the vortices diameter,vortex breakdown location,vortices distance from the wing surface,flow lines pattern nearby the wing,circulation distribution,and separation.The results show that the textured model has a positive effect on some of the parameters related to drag reduction and lift increase.The riblets increase the flow momentum near the wing’s upper surface except near the apex.They also increase the flow momentum behind the wing.展开更多
The marine ecosystem of the Jiaozhou Bay has degraded significantly in fisheries productivity and its ecological roles as spawning and nursery ground for many species of commercial importance has been declining in rec...The marine ecosystem of the Jiaozhou Bay has degraded significantly in fisheries productivity and its ecological roles as spawning and nursery ground for many species of commercial importance has been declining in recent years. A mass-balanced trophic model was developed using Ecopath with Ecosim to evaluate the trophic structure of the Jiaozhou Bay for improving ecosystem management. The model were parameterized based on the fisheries survey data in the Jiaozhou Bay in 2011, including 23 species groups and one detritus group according to their ecological roles. The trophic levels of these ecological groups ranged from 1(primary producers and detritus) to4.3(large demersal fishes). The estimated total system throughput was 12 917.10 t/(km^2·a), with 74.59% and25.41% contribution of the total energy flows from phytoplankton and detritus, respectively. Network analyses showed that the overall transfer efficiency of the ecosystem was 14.4%, and the mean transfer efficiency was 14.5%for grazing food chain and 13.9% for detritus food chain. The system omnivory index(SOI), Finn's cycled index(FCI) and connectance index(CI) were relatively low in this area while the total primary production/total respiration(TPP/TR) was high, indicating an immature and unstable status of the Jiaozhou Bay ecosystem. Mixed trophic impact analysis revealed that the cultured shellfish had substantial negative impacts on most functional groups. This study contributed to ecosystem-level evaluation and management planning of the Jiaozhou Bay ecosystem.展开更多
The use of multi-perspective and multi-scalar city networks has gradually developed into a range of critical approaches to understand spatial interactions and linkages. In particular, road linkages represent key chara...The use of multi-perspective and multi-scalar city networks has gradually developed into a range of critical approaches to understand spatial interactions and linkages. In particular, road linkages represent key characteristics of spatial dependence and distance decay, and are of great significance in depicting spatial relationships at the regional scale. Therefore, based on highway passenger flow data between prefecture-level administrative units, this paper attempted to identify the functional structures and regional impacts of city networks in China, and to further explore the spatial organization patterns of the existing functional regions, aiming to deepen our understanding of city network structures and to provide new cognitive perspectives for ongoing research. The research results lead to four key conclusions. First, city networks that are based on highway flows exhibit strong spatial dependence and hierarchical characteristics, to a large extent spatially coupled with the distributions of major megaregions in China. These phenomena are a reflection of spatial relationships at regional scales as well as core-periphery structure. Second, 19 communities that belong to an important type of spatial configuration are identified through community detection algorithm, and we suggest they are correspondingly urban economic regions within urban China. Their spatial metaphors include the administrative region economy, spatial spillover effects of megaregions, and core-periphery structure. Third, each community possesses a specific city network system and exhibits strong spatial dependence and various spatial organization patterns. Regional patterns have emerged as the result of multi-level, dynamic, and networked characteristics. Fourth, adopting a morphology-based perspective, the regional city network systems can be basically divided into monocentric, dual-nuclei, polycentric, and low-level equilibration spatial structures, while most are developing monocentrically.展开更多
基金funded by the National Natural Science Foundation of China(Grant No.52276047)the Open Fund of NationalKey Laboratory of SpacecraftThermal Control(Grant No.NKLST-JJ-202401011)the Beijing Municipal Science&Technology Commission(Grant No.Z231100006123010).
文摘Flow boiling in open microchannels offers highly efficient heat transfer performance and has attracted increasing attention in the fields of heat transfer and thermalmanagement of electronic devices in recent years.However,the continuous rise in power density of electronic components imposesmore stringent requirements on the heat transfer capability of microchannel flow boiling.HFE-7100,a dielectric coolant with favorable thermophysical properties,has become a focal point of research for enhancing flow boiling performance in open microchannels.The flow boiling heat transfer performance ofHFE-7100 was investigated in this study by fabricating micro-nano composite structures on the bottom surface of open microchannels using laser ablation technology.Based on visualization results,a comparative analysis was conducted on the bubble dynamics and flow pattern characteristics of HFE-7100 flow boiling in micronano structured open microchannels(MNSOMC)and smooth-surface open microchannels(SSOMC),to elucidate the enhancement mechanism of micro-nano structures on flow boiling heat transfer in open microchannels.The results indicate that the surface structures and strong wettability of MNSOMC accelerated bubble nucleation and departure.Moreover,bubbles in the channel tended to coalesce along the flow direction,forming elongated slug bubbles with high aspect ratios,which enabled efficient thin film evaporation in conjunction with intense nucleate boiling,thereby significantly enhancing flow boiling heat transfer.Under the experimental conditions of this study,the maximum enhancements in the heat transfer coefficient(HTC)and critical heat flux(CHF)of HFE-7100 inMNSOMC were 33.4%and 133.1%,respectively,with the CHF reaching up to 1542.3 kW⋅m^(−2).Furthermore,due to the superior wettability and capillary wicking capability of the micro-nano composite structures,the significant enhancement in flow boiling heat transfer was achieved without incurring a noticeable pressure drop penalty.
基金supported by the National Natural Science Foundation of China(Nos.U23A20583,52033005,U21A2090,and 52173040)Department of Science and Technology of Sichuan Province(No.2024NSFTD0003)。
文摘Poly(_(L)-lactide)(PLLA),a leading biodegradable polyester,has demonstrated potential as a sustainable alternative,owing to its excellent biodegradability and rigidity.However,their slow crystallization kinetics and poor heat resistance limit their application scope.Recent advances have highlighted that the combination of extensional flow and thermal fields can achieve toughness–stiffness balance,high transparency,and good heat resistance.However,the effect of extensional flow on the post-non-isothermal crystallization of PLLA during heating and the resulting crystalline texture remains unclear.In this study,PLLA with a heterogeneous amorphous structure and oriented polymorph was prepared by extensional flow.The effect of heterogeneous amorphous structures on non-isothermal crystallization kinetics during the heating process was studied by thermal analysis,polarized optical microscopy,infrared spectroscopy,and ex situ/in situ X-ray characterization.These results clearly illustrate that extensional flow enhances the formation of oriented crystalline structures,accelerates non-isothermal crystallization,and modulates the polymorphic composition of PLLA.Moreover,an unexpected dual cold-crystallization behavior is identified in ordered PLLA samples upon extensional flow,which is from the extensional flow-induced heterogeneous amorphous phase into α' phase(low-temperature peak)and the pristine amorphous phase intoαphase(high-temperature peak).The extensional flow primarily promotes the formation of the more perfectαandα'phases,but has a negative effect on the final content ofαphase formed after cold crystallization andα'-to-αphase transformation.The findings of this work advance the understanding of PLLA non-isothermal crystallization after extensional flow and offer valuable guidance for high-performance PLLA upon heat treatment in practical processing.
基金financially supported by the National Natural Science Foundation of China(51872090,51772097)Hebei Natural Science Fund for Distinguished Young Scholar(E2019209433)+3 种基金Youth Talent Program of Hebei Provincial Education Department(BJ2018020)Natural Science Foundation of Hebei Province(E2020209151,E2024209029)National Key R&D Plan Project(2022YFB4200305)Research Projects of China National Petroleum Corporation(2024ZG50,2023DQ03-04)。
文摘The sluggish electrochemical catalytic activity of the graphite felt electrodes for anode reaction is still a barrier for achieving high-performance vanadium redox flow battery(VRFB).It is significant to leverage the exceptional conductivity,excellent electrocatalytic activity,and structural tunability of MXene to address this issue.Herein,this work introduces nitrogen atoms to modulate the carbon layer structure of Ti_(3)C_(2)T_(x)MXene,inducing a reconfiguration of the local electronic structure,which enhances the anode interface activity and thereby improves the performance of VRFB.Ti_(3)C_(2)T_(x)exhibits high conductivity,excellent hydrophilicity,and a large specific surface area,providing excellent interface characteristics for V^(3+)/V^(2+)redox reaction.Moreover,interlayer treatment to modulate the mesoporous structure of MXene further increases the reactive surface area.Importantly,doping nitrogen atoms at carbon layer induces lattice distortions in Ti_(3)C_(2)T_(x),which enhances the charge transfer processes of the V^(3+)/V^(2+)redox reaction.The catalysis mechanism is also validated through density functional theory.Furthermore,the modified graphite felt electrode,as the anode of VRFB,relieves a higher energy efficiency of 68%at 250 mA cm^(-2),while the pristine electrode cannot operate at this current density.In addition,at 150 mA cm^(-2),the modified battery maintains energy efficiency at 75%without degradation after 500 cycles.This study utilizes rational atomic-level engineering for effective structural modulation to significantly enhance the catalytic activity of electrode reaction,offering a unique perspective for developing high-performance MXene electrocatalysts of VRFB.
文摘An experimental system for forming a rotating paraboloid shaped shallow water with a free surface was conducted to study coherent vortex structures in a differentially rotating quasi two dimensional zonal flow.Flow visualization and laser light scattering techniques were used to obtain the information of spatial flow patterns.Experimental results show that the coexistence of Coriolis effect and strong shear in latitudinal zones may lead to formation of coherent vortices.Power spectra analysis and photographs which were taken in a reference frame rotating with the observed vortices also justified the emergence,drift and evolution of persistent vortices on the large scale.Locked vortex state manifests the cyclone and anticyclone asymmetry.
文摘Inlet recirculation is proved as an effective way for centrifugal compressor surge margin extension,and is successively used in some engineering applications.Unfortunately its working mechanism is still not being well understood,which leads to redesigning of inlet recirculation mostly by experience.Also,most study about inlet recirculation is steady to date.It is necessary to study surge margin extension mechanism about inlet recirculation.To expose the mechanism in detail,steady and unsteady numerical simulations were performed on a centrifugal compressor with and without inlet recirculation.The results showed that,with inlet recirculation,the inlet axial velocity is augmented,relative Mach number around blade tip leading edge area is significantly reduced and so is the flow angle.As the flow angle decreased,the incidence angle reduced which greatly improves the flow field inside the impeller.Moreover,inlet recirculation changes the blade loading around blade tip and restrains the flow separation on the blade suction side at the leading edge area.The unsteady results of static pressure around blade surface,entropy at inlet crossflow section and vorticity distributions at near tip span surface indicated that,at near stall condition,strong fluctuation exists in the vicinity of tip area due to the interaction between tip leakage flow and core flow.By inlet recirculation these strong flow fluctuations are eliminated so the flow stability is greatly enhanced.All these improvements mentioned above are the reason for inlet recirculation delays compressor stall.This research reveals the surge margin extension reason of inlet recirculation from an unsteady flow viewpoint and provides important reference for inlet recirculation structure design.
基金supported by the National Natural Science Foundation of China(Grant No. 50579007)the Doctoral Program of Higher Education of China(Grant No. 200802170010)
文摘The theoretical research on the propulsive principle of aquatic animal becomes more important and attracted more researchers to make efforts on it. In the present study, a computational fluid dynamic (CFD) simulation of a three-dimensional traveling-wave undulations body of tuna has been developed to investigate the fluid flow features and vorticity structures around this body when moving in a straight line. The undulation only takes place in the posterior half of the fish, and the tuna-tail is considered as a lunate fin oscillating with the mode combined swaying with yawing. A Reynolds-averaged Navier-Stokes (RANS) equation is developed, employing a control-volume method and a k-omega SST turbulent model; meanwhile an unstructured tetrahedral grid, which is generated for the three-dimensional geometry, is used based on the deformation of the hind parts of the body and corresponding movement of the tail. We calculated the hydrodynamic performance of tuna-like body when a tuna swims in a uniform velocity, and compared the input power coefficient, output power coefficient and propulsive efficiency of the oscillating tuna-tail with or without body vortex shedding. Additionally, the load distribution on the body, flow features and vorticity structures around the body were demonstrated. The effect of interaction between the body-generated vortices and the tail-generated vorticity on the hydrodynamic performance can be obtained.
文摘Implementation of an opposing jet in design of a hypersonic blunt body significantly modifies the external flowfield and yields a considerable reduction in the aerodynamic drag.This study aims to investigate the effects of flowfield modeling parameters of injection and freestream on the flow structure and aerodynamics of a blunt body with an opposing jet in hypersonic flow.Reynolds-Averaged Navier-Stokes(RANS)equations with a Shear Stress Transport(SST)turbulence model are employed to simulate the intricate jet flow interaction.Through utilizing a Non-Intrusive Polynomial Chaos(NIPC)method to construct surrogates,a functional relation is established between input modeling parameters and output flowfield and aerodynamic quantities in concern.Sobol indices in sensitivity analysis are introduced to represent the relative contribution of each parameter.It is found that variations in modeling parameters produce large variations in the flow structure and aerodynamics.The jet-to-freestream total-pressure ratio,jet Mach number,and freestream Mach number are the major contributors to variation in surface pressure,demonstrating an evident location-dependent behavior.The penetration length of injection,reattachment angle of the shear layer,and aerodynamic drag are also most sensitive to the three crucial parameters above.In comparison,the contributions of freestream temperature,freestream density,and jet total temperature are nearly negligible.
基金Project(U1134203)supported by the National Natural Science Foundation of China
文摘According to the analysis of the turbulent intensity level around the high-speed train, the maximum turbulent intensity ranges from 0.2 to 0.5 which belongs to high turbulent flow. The flow field distribution law was studied and eight types of flow regions were proposed. They are high pressure with air stagnant region, pressure decreasing with air accelerating region, low pressure with high air flow velocity region I, turbulent region, steady flow region, low pressure with high air flow velocity region II,pressure increasing with air decelerating region and wake region. The analysis of the vortex structure around the train shows that the vortex is mainly induced by structures with complex mutation and large curvature change. The head and rear of train, the underbody structure, the carriage connection section and the wake region are the main vortex generating sources while the train body with even cross-section has rare vortexes. The wake structure development law studied lays foundation for the train drag reduction.
基金supported by the National Natural Science Foundation of China(No.51206180)the Natural Science Basic Research Plan in Shaanxi Province of China(No.2014JQ7276)
文摘A series of numerical analyses have been performed to investigate the flow structures in a narrow confined channel with 12 staggered circular impingement holes and one bigger exit hole. The flow enters the channel through the impingement holes and exits through the far end outlet. The flow fields corresponding to two jet Reynolds numbers (25000 and 65000) and three channel con- figurations with different ratios of the channel height to the impingement hole diameter (Zr 1, 3, 5) are analyzed by solving the Reynolds averaged Navier-Stokes equations with the realizable k-e turbulence model. Detailed flow field information including the secondary flow, the interaction between the jets and the cross flow, and flow distribution along the channel has been obtained. Comparisons between the numerical and experimental results of the flow fields at the four planes along the channel are performed to validate the numerical method. The calculated impingement pattern, high velocity flow distribution, low velocity separation region and vortices are in good agreement with the experimental data, implying the validity and effectiveness of the employed numerical approach for analyzing relevant flow field.
基金supported by China Armament Pre-research Foundation(Grant No. 51318010402)UK Engineering and Physical Science Research Council (EPSRC), and China Scholarship Council (Grant No.2010611054)
文摘Work on dynamic topology optimization of engineering structures for vibration suppression has mainly addressed the maximization of eigenfrequencies and gaps between consecutive eigenfrequencies of free vibration, minimization of the dynamic compliance subject to forced vibration, and minimization of the structural frequency response. A dynamic topology optimization method of bi-material plate structures is presented based on power flow analysis. Topology optimization problems formulated directly with the design objective of minimizing the power flow response are dealt with. In comparison to the displacement or velocity response, the power flow response takes not only the amplitude of force and velocity into account, but also the phase relationship of the two vector quantities. The complex expression of power flow response is derived based on time-harmonic external mechanical loading and Rayleigh damping. The mathematical formulation of topology optimization is established based on power flow response and bi-material solid isotropic material with penalization(SIMP) model. Computational optimization procedure is developed by using adjoint design sensitivity analysis and the method of moving asymptotes(MMA). Several numerical examples are presented for bi-material plate structures with different loading frequencies, which verify the feasibility and effectiveness of this method. Additionally, optimum results between topological design of minimum power flow response and minimum dynamic compliance are compared, showing that the present method has strong adaptability for structural dynamic topology optimization problems. The proposed research provides a more accurate and effective approach for dynamic topology optimization of vibrating structures.
基金Project supported by the National Basic Research Program of China (Grant No. 2009 CB724100)the National Natural Science Foundation of China (Grant No. 11172326)
文摘In a Mach 3.8 wind tunnel, both instantaneous and time-averaged flow structures of different scales around a blunt double-cone with or without supersonic film cooling were visualized via nano-tracer planar laser scattering (NPLS), which has a high spatiotemporal resolution. Three experimental cases with different injection mass flux rates were carried out. Many typical flow structures were clearly shown, such as shock waves, expansion fans, shear layers, mixing layers, and turbulent boundary layers. The analysis of two NPLS images with an interval of 5 us revealed the temporal evolution characteristics of flow structures. With matched pressures, the laminar length of the mixing layer was longer than that in the case with a larger mass flux rate, but the full covered region was shorter. Structures like K-H (Kelvin-Helmholtz) vortices were clearly seen in both flows. Without injection, the flow was similar to the supersonic flow over a backward- facing step, and the structures were relatively simpler, and there was a longer laminar region. Large scale structures such as hairpin vortices were visualized. In addition, the results were compared in part with the schlieren images captured by others under similar conditions.
基金Supported by the National Key Program for Developing Basic Sciences of China (No. G1999022103) and the National Natural Science Foundation of China (No. 20176059).
文摘This study investigates the heterogeneous structure and its influence on drag coefficient for concurrent-up gas-solid flow. The energy-minimization multi-scale (EMMS) model is modified to simulate the variation of structure parameters with solids concentration, showing the tendency for particles to aggregate to form clusters and for fluid to pass around clusters. The global drag coefficient is resolved into that for the dense phase, for the dilute phase and for the so-called inter-phase, all of which can be obtained from their respective phase-specific structure parameters. The computational results show that the drag coefficients of the different phases are quite different, and the global drag coefficient calculated from the EMMS approach is much lower than that from the correlation of Wen and Yu. The simulation results demonstrate that the EMMS approach can well describe the heterogeneous flow structure, and is very promising for incorporation into the two-fluid model or the discrete particle model as the closure law for drag coefficient.
基金The work was supported by the National Natural Science Foundation of China(Grant 11622216).
文摘Pipe-in-pipe(PIP)structures are widely used in offshore oil and gas pipelines to settle thermal insulation issues.A PIP structure system usually consists of two concentric pipes and one softer layer for thermal insulation consideration.The total response of the system is related to the dynamics of both pipes and the interactions between these two concentric pipes.In the current work,a theoretical model for flow-induced vibrations of a PIP structure system is proposed and analyzed in the presence of an internal axial flow and an external cross flow.The interactions between the two pipes are modeled by a linear distributed damper,a linear distributed spring and a nonlinear distributed spring along the pipe length.The unsteady hydrodynamic forces due to cross flow are modeled by two distributed van der Pol wake oscillators.The nonlinear partial differential equations for the two pipes and the wake are further discretized by the aid of Galerkin’s technique,resulting in a set of ordinary differential equations.These ordinary differential equations are further numeri cally solved by using a fourth-order Runge-Kutta integration algorithm.Phase portraits,bifurcation diagrams,an Argand diagram and oscillation shape diagrams are plotted,showing the existence of a lock-in phenomenon and figure-of-eight trajectory.The PIP system subjected to cross flow displays some interesting dynamical behaviors different from that of a single-pipe structure.
基金Under the auspices of National Natural Science Foundation of China (No 40571019)
文摘The structure of wind-sand flow under different total sand transport rates was measured with field vertical anemometer and sand trap on the crest of typical coastal transverse ridge in Changli Gold Coast of Hebei Province, which is one of the most typical coastal aeolian distribution regions in China and famous for the tall and typical coastal transverse ridges. The measurement results show that, on the conditions of approximate wind velocities and same surface materials and environments, some changes happen to the structure of wind-sand flow with the increase of total sand transport rate on the crest of coastal transverse ridge. First, the sand transport rates of layers at different heights in the wind-sand flow increase, with the maximum increase at the height layer of 4-8cm. Second, the ratios of sand trans-port rates of layers at different heights to total sand transport rate decrease at the low height layer (0-4cm), but increase at the high height layer (4-60cm). Third, the distribution of the sand transport rate in the wind-sand flow can be expressed by an exponential function at the height layer of 0-40cm, but it changes from power function model to ex-ponential function model in the whole height layer (0-60cm) and changes into polynomial function model at the height layer of 40-60cm with the increase of total sand transport rate. Those changes have a close relationship with the limit of sand grain size of wind flow transporting and composition of sand grain size in the wind-sand flow.
基金Project(2014GK3150)supported by Science and Technology Plan of Hunan Province,China
文摘The three-dimensional internal flow field of centrifugal pump is complex and variable with design parameters and operation conditions. The post-processing technique named differential amplification method was proposed for the comparison study of different flow structures. The full steady flow fields of an industrial centrifugal pump working on-design and off-design points were numerically investigated by solving Reynolds average Navier-Stokes equations together with a shear-stress transport(SST) k-? turbulence model. And the numerically predicted performance curves of the studied pump agree well with test measurement results. Compared with the flow flied on design point under the help of differential amplification method, the disturbance caused by interaction between blade and volute tongue is very obvious and it extends to the diffuser pipe on the working point with 0.8 times rated flux. While on the point with 1.2 times rated flux, the flow distribution in impeller region is roughly even and it flows more to the bottom section of the diffuser pipe. The above method was proved to be good at displaying the subtle secondary flow structure changes with a higher resolution effect relative to single isolated case observation, which helps the optimization decision-making from multiple design cases.
基金project was supported by the National Natural Science Foundation of China (Grants 51306020, 11172040)the Natural Science Foundation of Beijing (Grant 3144034)the Excellent Young Scholars Research Fund of Beijing Institute of Technology
文摘The objective of this paper is to address the transient flow structures around a pitching hydrofoil by com- bining physical and numerical studies. In order to predict the dynamic behavior of the flow structure effectively, the Lagrangian coherent structures (LCS) defined by the ridges of the finite-time Lyapunov exponent (FTLE) are utilized under the framework of Navier-Stokes flow computations. In the numerical simulations, the k-w shear stress trans- port (SST) turbulence model, coupled with a two-equation F-Reo transition model, is used for the turbulence closure. Results are presented for a NACA66 hydrofoil undergoing slowly and rapidly pitching motions from 0° to 15° then back to 0° at a moderate Reynolds number Re = 7.5 × 105. The results reveal that the transient flow structures can be observed by the LCS method. For the slowly pitching case, it consists of five stages: quasi-steady and laminar, transition from laminar to turbulent, vortex development, large-scale vortex shedding, and reverting to laminar. The observation of LCS and Lagrangian particle tracers elucidates that the trailing edge vortex is nearly attached and stable during the vortex development stage and the interaction between the leading and trailing edge vortex caused by the adverse pres- sure gradient forces the vortexes to shed downstream during the large-scale vortex shedding stage, which corresponds to obvious fluctuations of the hydrodynamic response. For the rapidly pitching case, the inflection is hardly to be observed and the stall is delayed. The vortex formation, interaction, and shedding occurred once instead of being repeated three times, which is responsible for just one fluctuation in the hydrody- namic characteristics. The numerical results also show that the FTLE field has the potential to identify the transient flows, and the LCS can represent the divergence extent of infinite neighboring particles and capture the interface of the vortex region.
基金supported by the Brain Pool Program through the Korean Federation of Science and Technology Societies (KOFST), which is funded by the Ministry of Science, ICT and Future Planningprovided by the National Research Foundation of Korea (NRF) grant, which is funded by the Korean government (MSIT) (Nos. 2011-0030013, 2018R1A2B2007117 and NRF-2017K1A3A1A30084513)
文摘This research investigates the aerodynamic performance and flow characteristics of a delta wing with 65° sweep angle and with coarse axial riblets,and then compares with that of a smooth-surface delta wing.Particle Image Velocimetry(PIV)were utilized to visualize the flow over the wing at 6 cross-sections upright to the wing surface and parallel to the wing span,as well as 3 longitudinal sections on the leading edge,symmetry plane,and a plane between them at Angles of Attack(AOA)=20°and 30°and Re=1.2×10~5,2.4×10~5,and 3.6×10~5.The effects of the riblets were studied on the vortices diameter,vortex breakdown location,vortices distance from the wing surface,flow lines pattern nearby the wing,circulation distribution,and separation.The results show that the textured model has a positive effect on some of the parameters related to drag reduction and lift increase.The riblets increase the flow momentum near the wing’s upper surface except near the apex.They also increase the flow momentum behind the wing.
基金The Public Science and Technology Research Funds Projects of Ministry of Agriculture under contract No.201303050-02the Scientific and Technological Innovation Project for the Qingdao National Laboratory for Marine Science and Technology under contract No.2015ASKJ02+2 种基金the Fundamental Research Funds for the Central Universities under contract No.201262004the National Natural Science Foundation of China under contract No.41006083the Shandong Provincial Natural Science Foundation,China under contract No.ZR2010DQ026
文摘The marine ecosystem of the Jiaozhou Bay has degraded significantly in fisheries productivity and its ecological roles as spawning and nursery ground for many species of commercial importance has been declining in recent years. A mass-balanced trophic model was developed using Ecopath with Ecosim to evaluate the trophic structure of the Jiaozhou Bay for improving ecosystem management. The model were parameterized based on the fisheries survey data in the Jiaozhou Bay in 2011, including 23 species groups and one detritus group according to their ecological roles. The trophic levels of these ecological groups ranged from 1(primary producers and detritus) to4.3(large demersal fishes). The estimated total system throughput was 12 917.10 t/(km^2·a), with 74.59% and25.41% contribution of the total energy flows from phytoplankton and detritus, respectively. Network analyses showed that the overall transfer efficiency of the ecosystem was 14.4%, and the mean transfer efficiency was 14.5%for grazing food chain and 13.9% for detritus food chain. The system omnivory index(SOI), Finn's cycled index(FCI) and connectance index(CI) were relatively low in this area while the total primary production/total respiration(TPP/TR) was high, indicating an immature and unstable status of the Jiaozhou Bay ecosystem. Mixed trophic impact analysis revealed that the cultured shellfish had substantial negative impacts on most functional groups. This study contributed to ecosystem-level evaluation and management planning of the Jiaozhou Bay ecosystem.
基金National Natural Science Foundation of China,No.41530751,No.41471113,No.41601165
文摘The use of multi-perspective and multi-scalar city networks has gradually developed into a range of critical approaches to understand spatial interactions and linkages. In particular, road linkages represent key characteristics of spatial dependence and distance decay, and are of great significance in depicting spatial relationships at the regional scale. Therefore, based on highway passenger flow data between prefecture-level administrative units, this paper attempted to identify the functional structures and regional impacts of city networks in China, and to further explore the spatial organization patterns of the existing functional regions, aiming to deepen our understanding of city network structures and to provide new cognitive perspectives for ongoing research. The research results lead to four key conclusions. First, city networks that are based on highway flows exhibit strong spatial dependence and hierarchical characteristics, to a large extent spatially coupled with the distributions of major megaregions in China. These phenomena are a reflection of spatial relationships at regional scales as well as core-periphery structure. Second, 19 communities that belong to an important type of spatial configuration are identified through community detection algorithm, and we suggest they are correspondingly urban economic regions within urban China. Their spatial metaphors include the administrative region economy, spatial spillover effects of megaregions, and core-periphery structure. Third, each community possesses a specific city network system and exhibits strong spatial dependence and various spatial organization patterns. Regional patterns have emerged as the result of multi-level, dynamic, and networked characteristics. Fourth, adopting a morphology-based perspective, the regional city network systems can be basically divided into monocentric, dual-nuclei, polycentric, and low-level equilibration spatial structures, while most are developing monocentrically.