The least-square gridless method was extended to simulate the compressible multi-material flows. The algorithm was accomplished to solve the Arbitrary Lagrange-Euler( ALE) formulation. The local least-square curve fit...The least-square gridless method was extended to simulate the compressible multi-material flows. The algorithm was accomplished to solve the Arbitrary Lagrange-Euler( ALE) formulation. The local least-square curve fits was adopted to approximate the spatial derivatives of a point on the base of the points in its circular support domain,and the basis function was linear. The HLLC( Harten-Lax-van Leer-Contact) scheme was used to calculate the inviscid flux. On the material interfaces,the gridless points were endued with a dual definition corresponding to different materials. The moving velocity of the interface points was updated by solving the Riemann problem. The interface boundary condition was built by using the Ghost Fluid Method( GFM).Computations were performed for several one and two dimensional typical examples. The numerical results show that the interface and the shock wave are well captured,which proves the effectiveness of gridless method in dealing with multi-material flow problems.展开更多
The stress minimization multi-material topology optimization(MMTO)approach has recently attracted significant attention because of its applications in aerospace and mechanical engineering.Nonetheless,the stress minimi...The stress minimization multi-material topology optimization(MMTO)approach has recently attracted significant attention because of its applications in aerospace and mechanical engineering.Nonetheless,the stress minimization MMTO approach may result in stress surpassing the material's tolerance limit,potentially culminating in failure.This research proposes a novel way for imposing stress constraints on each material to regulate their respective stress levels.The fundamental concept is that each material possesses its own interpolation function for the stress model.The maximum von Mises stress for each material can be established with the definition of an upper limit,ensuring that the materials will perform safely and effectively.This aids topological structures in resisting failure and augmenting strength.A multi-physics system including thermoelastic and self-weight loads is concurrently examined alongside stress limitations.The global stress constraint utilizes the p-norm function,and the adjoint method is used to derive sensitivity.This work employs a three-field strategy utilizing density filtering and Heaviside projection functions to mitigate the artificial stress in low density.The technique is assessed through two-dimensional(2D)and three-dimensional(3D)examples,illustrating the influence of stress limits on the compliance minimization under heat and self-weight loads.The optimized results indicate a substantial decrease in the stress levels accompanied by a minor gain in compliance,while maintaining the stress within the specified range for all materials.展开更多
In this paper,a direct arbitrary Lagrangian-Eulerian(ALE)discontinuous Galerkin(DG)scheme is proposed for simulating compressible multi-material flows on the adaptive quadrilateral meshes.Our scheme couples a conserva...In this paper,a direct arbitrary Lagrangian-Eulerian(ALE)discontinuous Galerkin(DG)scheme is proposed for simulating compressible multi-material flows on the adaptive quadrilateral meshes.Our scheme couples a conservative equation related to the volume-fraction model with the Euler equations for describing the dynamics of the fluid mixture.The coupled system is discretized in the reference element and we use a kind of Taylor expansion basis functions to construct the interpolation polynomials of the variables.We show the property that the material derivatives of the basis functions in the DG discretization are equal to zero,with which the scheme is simplified.In addition,the mesh velocity in the ALE framework is obtained by using the adaptive mesh method from[H.Z.Tang and T.Tang,Adaptive mesh methods for one-and two-dimensional hyperbolic conservation laws,SIAMJ.NUMER.ANAL].This adaptivemesh method can automatically concentrate the mesh nodes near the regions with large gradient values and greatly reduces the numerical dissipation near the material interfaces in the simulations.With the help of this adaptivemesh method,the resolution of the solution near the target regions can be greatly improved and the computational efficiency of the simulation is increased.Our scheme can be applied in the simulations for the gas andwatermedia efficiently,and it ismore concise compared to some other methods such as the indirect ALE methods.Several examples including the gas-water flow problemare presented to demonstrate the efficiency of our scheme,and the results show that our scheme can capture the wave structures sharply with high robustness.展开更多
Soft robots, inspired by the flexibility and versatility of biological organisms, have potential in a variety of applications. Recent advancements in magneto-soft robots have demonstrated their abilities to achieve pr...Soft robots, inspired by the flexibility and versatility of biological organisms, have potential in a variety of applications. Recent advancements in magneto-soft robots have demonstrated their abilities to achieve precise remote control through magnetic fields, enabling multi-modal locomotion and complex manipulation tasks. Nonetheless, two main hurdles must be overcome to advance the field: developing a multi-component substrate with embedded magnetic particles to ensure the requisite flexibility and responsiveness, and devising a cost-effective,straightforward method to program three-dimensional distributed magnetic domains without complex processing and expensive machinery. Here, we introduce a cost-effective and simple heat-assisted in-situ integrated molding fabrication method for creating magnetically driven soft robots with three-dimensional programmable magnetic domains. By synthesizing a composite material with neodymium-iron-boron(NdFeB) particles embedded in a polydimethylsiloxane(PDMS) and Ecoflex matrix(PDMS:Ecoflex = 1:2 mass ratio, 50% magnetic particle concentration), we achieved an optimized balance of flexibility, strength, and magnetic responsiveness. The proposed heat-assisted in-situ magnetic domains programming technique,performed at an experimentally optimized temperature of 120℃, resulted in a 2 times magnetization strength(9.5 mT) compared to that at 20℃(4.8 m T), reaching a saturation level comparable to a commercial magnetizer. We demonstrated the versatility of our approach through the fabrication of six kinds of robots, including two kinds of two-dimensional patterned soft robots(2D-PSR), a circular six-pole domain distribution magnetic robot(2D-CSPDMR), a quadrupedal walking magnetic soft robot(QWMSR), an object manipulation robot(OMR), and a hollow thin-walled spherical magneto-soft robot(HTWSMSR). The proposed method provides a practical solution to create highly responsive and adaptable magneto-soft robots.展开更多
Multi-material 3D sand printing has gained significant attention;however,research has mainly focused on materials and mechanisms,with limited exploration of optimizing the sand-laying process through numerical simulat...Multi-material 3D sand printing has gained significant attention;however,research has mainly focused on materials and mechanisms,with limited exploration of optimizing the sand-laying process through numerical simulations.In this study,we investigated the dynamic behavior of sand particles during a vibratory sand-laying process for multi-material additive manufacturing using discrete element simulations.The objective is to enable precise control over the amount and distribution of sand for multi-material printing.In this study,we combined experiments and simulations to calibrate the contact parameters of different sands and establish a relationship between the curing agent content and surface energy of sand particles.A model for the vibratory fall of multimaterial sand was developed to study the motion characteristics of sand particles.This allows for macro-control over the sand spreading flow and high-quality multi-material sand laying.The results show that the flow rate of falling sand increases with decreasing surface energy of the particles,wider spreader openings,and higher vibration frequencies.For silica and chromite sands,when their surface energy ranged from 0.15 to 25 J·m^(2)and0.01-0.03 J·m^(2),respectively,and the sand spreader opening was 6 mm with a vibration frequency of 500 Hz,the sand flow rates of both materials became nearly identical.However,a higher sand paving speed and height increased the scattering of sand particles outside the target area,thereby decreasing the paving quality.The results accomplished in this study enable precise and uniform sand particle deposition and offers guidelines for optimizing sand speed and height,thus expanding the application of multi-material sand 3D printing in complex and high-performance manufacturing.展开更多
The application of multi-material topology optimization affords greater design flexibility compared to traditional single-material methods.However,density-based topology optimization methods encounter three unique cha...The application of multi-material topology optimization affords greater design flexibility compared to traditional single-material methods.However,density-based topology optimization methods encounter three unique challenges when inertial loads become dominant:non-monotonous behavior of the objective function,possible unconstrained characterization of the optimal solution,and parasitic effects.Herein,an improved Guide-Weight approach is introduced,which effectively addresses the structural topology optimization problem when subjected to inertial loads.Smooth and fast convergence of the compliance is achieved by the approach,while also maintaining the effectiveness of the volume constraints.The rational approximation of material properties model and smooth design are utilized to guarantee clear boundaries of the final structure,facilitating its seamless integration into manufacturing processes.The framework provided by the alternating active-phase algorithm is employed to decompose the multi-material topological problem under inertial loading into a set of sub-problems.The optimization of multi-material under inertial loads is accomplished through the effective resolution of these sub-problems using the improved Guide-Weight method.The effectiveness of the proposed approach is demonstrated through numerical examples involving two-phase and multi-phase materials.展开更多
To investigate the influence of Al-Zn-Mg-Cu alloy with as-homogenized and as-rolled initial microstructures on the tensile flow behavior,isothermal tensile tests were conducted on a GLEEBLE-3500 isothermal simulator a...To investigate the influence of Al-Zn-Mg-Cu alloy with as-homogenized and as-rolled initial microstructures on the tensile flow behavior,isothermal tensile tests were conducted on a GLEEBLE-3500 isothermal simulator at temperatures of 380-440℃and strain rates of 0.05-1 s^(−1).The Johnson-Cook model,Hensel-Spittel model,strain-compensated Arrhenius model,and critical fracture strain model were established.Results show that through the evaluation of the models using the correlation coefficient(R)and the average absolute relative error,the strain-compensated Arrhenius model can represent the flow behavior of the alloy more accurately.Shear bands are more pronounced in the as-homogenized specimens,whereas dynamic recrystallization is predominantly observed in as-rolled specimens.Fracture morphology analysis reveals that a mixed fracture mechanism is prevalent in the as-homogenized specimen,whereas a ductile fracture mechanism is predominant in the as-rolled specimen.The processing maps indicate that the unstable region is reduced in the as-rolled specimens compared with that in the as-homogenized specimens.The optimal hot working windows for the as-homogenized and as-rolled specimens are determined as 410-440℃/0.14-1 s^(−1)and 380-400℃/0.05-0.29 s^(−1),respectively.展开更多
Vanadium redox flow batteries(VRFBs)are a means of large-scale energy storage due to their excellent scalability,safety,long cycling life,and decoupled power and energy capacities.However,the slow redox kinetics of va...Vanadium redox flow batteries(VRFBs)are a means of large-scale energy storage due to their excellent scalability,safety,long cycling life,and decoupled power and energy capacities.However,the slow redox kinetics of vanadium species on conventional carbon electrodes remains a major limitation to their performance.We investigated the deposition of carbon black,carbon nanotubes,and electrochemically exfoliated graphene(Exf-Gr)onto thermally-activated carbon paper(ACP)by spray coating to increase the electrode electrocatalytic activity.The modified electrodes were characterized using scanning electron microscopy,X-ray diffraction,Raman spectroscopy,X-ray photoelectron microscopy,and surface area analysis,while their electrochemical properties were evaluated by cyclic voltammetry,electrochemical impedance spectroscopy,and singlecell VRFB testing.Among the modified electrodes,Exf-Gr/ACP had the best performance,achieving a 2.9-fold reduction in charge transfer resistance compared to pristine ACP and delivering 2.5 times the discharge capacity in single-cell tests.This improvement is attributed to Exf-Gr’s high surface area,favorable catalytic activity,and excellent dispersion on the ACP substrate.Surface modification with electrochemically exfoliated graphene is a highly effective strategy for improving the electrode performance in VRFB systems,with significant implications for large-scale energy storage.展开更多
Large-eddy simulation(LES)is conducted to study the statistical properties of mixed-phase turbulence induced by the breaking of bow waves in flow past a partially submerged plate.The simulation is performed using a fi...Large-eddy simulation(LES)is conducted to study the statistical properties of mixed-phase turbulence induced by the breaking of bow waves in flow past a partially submerged plate.The simulation is performed using a finite difference method,with the air-water interface captured by a coupled level-set and volume-of-fluid method.Four cases are conducted to investigate the effects of Froude number on turbulent statistics,including the mean velocity,turbulence kinetic energy,and turbulence mass flux(TMF),which is an additional unclosed term in the Reynolds-averaged momentum equation.The TMF,especially its vertical component,shows a complex behaviour with respect to the Froude number.This property of the TMF imposes high demands on the robustness of the closure model of TMF.The present LES data is further used to examine a closure model of the TMF production term,which shows a high correlation with the data obtained from LES.展开更多
To clarify fluid flow mechanisms and establish effective development conditions in continental shale oil reservoirs,a high-temperature,high-pressure steady-state flow system integrated with nuclear magnetic resonance(...To clarify fluid flow mechanisms and establish effective development conditions in continental shale oil reservoirs,a high-temperature,high-pressure steady-state flow system integrated with nuclear magnetic resonance(NMR)technology has been developed.The apparatus combines sample evacuation,rapid pressurization and saturation,and controlled displacement,enabling systematic investigation of single-phase shale oil flow under representative reservoir conditions.Related experiments allow proper quantification of the activation thresholds and relative contributions of different pore types to flow.A movable fluid index(MFI),defined using dual T_(2) cutoff values,is introduced accordingly and linked to key flow parameters.The results reveal distinct multi-scale characteristics of single-phase shale oil transport,namely micro-scale graded displacement and macro-scale segmented nonlinear behavior.As the injection-production pressure difference increases,flow pathways are activated progressively,beginning with fractures,followed by large and then smaller macropores,leading to a pronounced enhancement in apparent permeability.Although mesopores and micropores contribute little to direct flow,their indirect influence becomes increasingly important,and apparent permeability gradually approaches a stable limit at higher pressure difference.It is also shown that the MFI exhibits a strong negative correlation with the starting pressure gradient and a positive correlation with apparent permeability,providing a rapid and reliable indicator of shale oil flow capacity.Samples containing through-going fractures display consistently higher MFI values and superior flowability compared with those dominated by laminated fractures,highlighting the pivotal role of well-connected fracture networks generated by large-scale hydraulic fracturing in improving shale oil production.展开更多
Investigating the wind-sand flow response regularity in the longitudinal slope sections of desert highways provides a scientific basis for selecting the slope of desert roads.This study uses the Tengger Desert section...Investigating the wind-sand flow response regularity in the longitudinal slope sections of desert highways provides a scientific basis for selecting the slope of desert roads.This study uses the Tengger Desert section of the Wuhai-Maqin Expressway as a case study,employing CFD numerical simulation methods to calculate and analyze the wind-sand flow field distribution characteristics in different longitudinal slope sections.The results show that:(1)Along with the direction of the incoming flow,the windward and leeward slope toes of the embankment are low-wind-speed zones,with the wind speed at the leeward slope toe being even lower.The higher the embankment,the larger the low-wind-speed zone at the windward and leeward slope toes.As the longitudinal slope increases,the extent of the lowwind-speed zone at the same location along the route also increases.(2)Along the route direction,the wind speed at the windward and leeward slope toes decreases as embankment height increases.At the embankment toe,sand particles are transported from the top to the bottom of the longitudinal slope,and the greater the longitudinal slope,the stronger the transport effect.(3)Along the route direction,the sand accumulation around the embankment gradually gathers toward the bottom of the longitudinal slope as the slope increases.When the longitudinal slope is 3%and 4%,the trend of sand accumulation moving from the windward side at the end of the route to the leeward side at the start of the route is more significant.When the longitudinal slope is less than or equal to 3%,severe sand accumulation within the embankment range is reduced by 86.4%or more compared to when the slope is 4%.(4)Under the same longitudinal slope,the higher the embankment height,the smaller its transport rate.When the embankment height is the same,the greater the longitudinal slope,the greater the embankment transport rate.展开更多
The newly formulated non-Newtonian rivulet flows streaming down an inclined planar surface,with additional periodic perturbations arising from the application of the 2nd Stokes problem to the investigation of rivulet ...The newly formulated non-Newtonian rivulet flows streaming down an inclined planar surface,with additional periodic perturbations arising from the application of the 2nd Stokes problem to the investigation of rivulet dynamics,are demonstrated in the current research.Hereby,the 2nd Stokes problem assumes that the surface,with a thin shared layer of the fluid on it,oscillates in a harmonic manner along the x-axis of the rivulet flow,which coincides with the main flow direction streaming down the underlying surface.We obtain the exact extension of the rivulet flow family,clarifying the structure of the pressure field,which fully absorbs the arising perturbation.The profile of the velocity field is assumed to be Gaussian-type with a non-zero level of plasticity.Hence,the absolutely non-Newtonian case of the viscoplastic flow solution,which satisfies the motion and continuity equations,is considered(with particular cases of exact solutions for pressure).The perturbed governing equations of motion for rivulet flows then result in the Riccati-type ordinary differential equation(ODE),describing the dynamics of the coordinate x(t).The approximated schematic dynamics are presented in graphical plots.展开更多
Particle-and droplet-laden flows are central to many problems in mechanics and transport.They occur in sedimentladen boundary layers,gas-solid and gas-liquid dispersions,and surface-water films driven by external forc...Particle-and droplet-laden flows are central to many problems in mechanics and transport.They occur in sedimentladen boundary layers,gas-solid and gas-liquid dispersions,and surface-water films driven by external forcing.They also underpin practical applications ranging from environmental transport to high-speed and aerothermal systems.Despite decades of progress,prediction remains difficult.The physics spans a wide range of scales and often couples turbulence,interphase momentum exchange,collisions,and interfacial transport.Reliable computation therefore requires both robust numerical methodology and careful physical interpretation.展开更多
The Cooling-Storage-Ring External-target Experiment(CEE)at the Heavy Ion Research Facility in Lanzhou(HIRFL)is designed to study the properties of nuclear matter created in heavy-ion collisions at beam energies from a...The Cooling-Storage-Ring External-target Experiment(CEE)at the Heavy Ion Research Facility in Lanzhou(HIRFL)is designed to study the properties of nuclear matter created in heavy-ion collisions at beam energies from a few hundred MeV/u up to 1 GeV/u.It aims to facilitate research on the quantum chromodynamics(QCD)phase structure in the high-baryondensity region.Collective flow is a fundamental observable in heavy-ion collision experiments,providing information on the bulk properties of the produced matter.Although the standard event plane method has been widely used to measure collective flow,it is still important to validate and optimize this method for the CEE spectrometer.In this paper,we study the experimental procedures for measuring directed flow in^(238)U+^(238)U collisions at 500 MeV/u,using event planes reconstructed by Multi-Wire Drift Chamber and Zero Degree Calorimeter,respectively.Jet AA Microscopic(JAM)transport generator is used to generate events,and the detector response is simulated by the CEE Fast Simulation(CFS)package.Finally,the optimal kinematic region for proton directed flow measurements is discussed for the future CEE experiment.展开更多
The Reynolds-averaged Navier-Stokes(RANS)technique enables critical engineering predictions and is widely adopted.However,since this iterative computation relies on the fixed-point iteration,it may converge to unexpec...The Reynolds-averaged Navier-Stokes(RANS)technique enables critical engineering predictions and is widely adopted.However,since this iterative computation relies on the fixed-point iteration,it may converge to unexpected non-physical phase points in practice.We conduct an analysis on the phase-space characteristics and the fixed-point theory underlying the k-ε turbulence model,and employ the classical Kolmogorov flow as a framework,leveraging its direct numerical simulation(DNS)data to construct a one-dimensional(1D)system under periodic/fixed boundary conditions.The RANS results demonstrate that under periodic boundary conditions,the k-ε model exhibits only a unique trivial fixed point,with asymptotes capturing the phase portraits.The stability of this trivial fixed point is determined by a mathematically derived stability phase diagram,indicating the fact that the k-ε model will never converge to correct values under periodic conditions.In contrast,under fixed boundary conditions,the model can yield a stable non-trivial fixed point.The evolutionary mechanisms and their relationship with boundary condition settings systematically explain the inherent limitations of the k-ε model,i.e.,its deficiency in computing the flow field under periodic boundary conditions and sensitivity to boundary-value specifications under fixed boundary conditions.These conclusions are finally validated with the open-source code OpenFOAM.展开更多
This study investigates the enhancement of convective heat transfer in a serpentine pipe using ferrofluid flow influenced by dual non-uniform magnetic sources.The primary objective is to improve thermal performance in...This study investigates the enhancement of convective heat transfer in a serpentine pipe using ferrofluid flow influenced by dual non-uniform magnetic sources.The primary objective is to improve thermal performance in compact cooling systems,such as those used in heat exchangers.A two-dimensional,steady-state Computational Fluid Dynamic(CFD)model is developed in ANSYS Fluent to simulate the behavior of an incompressible ferrofluid under applied constant heat flux and magnetic fields.The magnetic force is modeled using the Kelvin force,which acts on magnetized nanoparticles in response to spatially varying electromagnetic fields generated by two strategically positioned current-carrying wires.The effects of magnetic field strength,quantified by the magnetic number(Mn),on flow behavior and temperature distribution are thoroughly analyzed.The results indicate that increasing Mn leads to higher Nusselt numbers,demonstrating enhanced convective heat transfer.Secondary vortices induced by magnetic forcing improve fluid mixing,particularly in curved regions of the pipe.A mesh-independence study and model validation with benchmark data support the reliability of the numerical framework.This work highlights the potential of magnetic-field-assisted thermal control in energy-efficient cooling applications and provides a foundation for the further development of advanced ferrofluid-based heat transfer systems.展开更多
We report an immobilized enzyme-catalyzed batch and continuous-flow synthesis of optically pure ethyl(R)-pantothenate((R)-PaOEt),the direct precursor of d-pantothenic acid.Firstly,a ketoreductase mutant designated as ...We report an immobilized enzyme-catalyzed batch and continuous-flow synthesis of optically pure ethyl(R)-pantothenate((R)-PaOEt),the direct precursor of d-pantothenic acid.Firstly,a ketoreductase mutant designated as M2,carrying two-point mutations of F97L and M242F relative to the wild-type SSCR,was constructed by site-directed mutagenesis,exhibited simultaneously improved activity toward ethyl 2′-ketopantothenate(K-PaOEt)and isopropanol,and could effectively catalyze the stereoselective reduction of K-PaOEt to(R)-PaOEt by using isopropanol as the sacrificial co-substrate to regenerate NADPH.After screening six commercially available carriers,an amino resin LXTE-700 was identified as the best solid support for the immobilization of M2 via the glutaraldehyde activation method.Upon optimization of the immobilization process and reaction conditions,the fabricated immobilized enzyme M2@amino resin demonstrated excellent recyclability and reusability,with the complete conversion of K-PaOEt to(R)-PaOEt being still realized after 12 cycles of reuse.Finally,M2@amino resin-catalyzed synthesis of(R)-PaOEt was successfully implemented in continuous-flow,accomplishing a 6.3 times higher space-time yield than that with the batch synthesis(529.2 versus 84 g L^(-1) d^(-1)).Our developed flow biocatalysis system also features an outstanding operational stability,as evidenced by the 100%conversion rate achieved after 15 consecutive days of operation.展开更多
The Rukwa Rift section of the East Africa Rift System presents a type setting for radiogenic helium accumulation in a petroleum free basin.As a prerequisite for accumulation,a considerable high heat flow anomaly is re...The Rukwa Rift section of the East Africa Rift System presents a type setting for radiogenic helium accumulation in a petroleum free basin.As a prerequisite for accumulation,a considerable high heat flow anomaly is required from tectonothermal events to drive the release and circulation of radiogenic helium in the continental crust.Here we apply statistical analysis on geochemical data observed in thermal springs and recorded heat flow to account for crustal helium mass balance for each tectonothermal event in the region.Our results demonstrate anomalously high heat flow~64-99 mW/m^(2) with a consistent trend of helium isotopic ratio and fluid chemistry in the Rukwa Rift.Mass balance calculation show that the whole crustal volume underlying the East Africa Helium Pool(EAHP)has a capability of producing radiogenic helium of about 9.9×10^(6) mol/yr(22×10^(-6) mol 4He/m^(2) yr)while the total radiogenic helium flux ranges between~2.39×10^(6) mol/yr and~2.68×10^(9) mol/yr.The Tanzania Craton contributes largely to radiogenic helium releasing up to 50% of the total capacity in the region.The total ^(4)He emission in the Rukwa Rift Basin is about 4.45×10^(5)-5.01×10^(8) mol/yr which is thus equivalent to 19%-21% of the total production capacity in the region.These results imply that the helium accumulation in the EAHP would have started as early as Paleoproterozoic(2.349 Ga).These results provide a qualitative and quantitative insight to assess both helium and geothermal potentiality in the region.展开更多
文摘The least-square gridless method was extended to simulate the compressible multi-material flows. The algorithm was accomplished to solve the Arbitrary Lagrange-Euler( ALE) formulation. The local least-square curve fits was adopted to approximate the spatial derivatives of a point on the base of the points in its circular support domain,and the basis function was linear. The HLLC( Harten-Lax-van Leer-Contact) scheme was used to calculate the inviscid flux. On the material interfaces,the gridless points were endued with a dual definition corresponding to different materials. The moving velocity of the interface points was updated by solving the Riemann problem. The interface boundary condition was built by using the Ghost Fluid Method( GFM).Computations were performed for several one and two dimensional typical examples. The numerical results show that the interface and the shock wave are well captured,which proves the effectiveness of gridless method in dealing with multi-material flow problems.
基金Project supported by the National Research Foundation of Korea(NRF)grant funded by the Korea government(MSIT)(No.RS-2025-02303676)。
文摘The stress minimization multi-material topology optimization(MMTO)approach has recently attracted significant attention because of its applications in aerospace and mechanical engineering.Nonetheless,the stress minimization MMTO approach may result in stress surpassing the material's tolerance limit,potentially culminating in failure.This research proposes a novel way for imposing stress constraints on each material to regulate their respective stress levels.The fundamental concept is that each material possesses its own interpolation function for the stress model.The maximum von Mises stress for each material can be established with the definition of an upper limit,ensuring that the materials will perform safely and effectively.This aids topological structures in resisting failure and augmenting strength.A multi-physics system including thermoelastic and self-weight loads is concurrently examined alongside stress limitations.The global stress constraint utilizes the p-norm function,and the adjoint method is used to derive sensitivity.This work employs a three-field strategy utilizing density filtering and Heaviside projection functions to mitigate the artificial stress in low density.The technique is assessed through two-dimensional(2D)and three-dimensional(3D)examples,illustrating the influence of stress limits on the compliance minimization under heat and self-weight loads.The optimized results indicate a substantial decrease in the stress levels accompanied by a minor gain in compliance,while maintaining the stress within the specified range for all materials.
基金supported by National Natural Science Foundation of China(Grant Nos.12071046,11571002,11772067).
文摘In this paper,a direct arbitrary Lagrangian-Eulerian(ALE)discontinuous Galerkin(DG)scheme is proposed for simulating compressible multi-material flows on the adaptive quadrilateral meshes.Our scheme couples a conservative equation related to the volume-fraction model with the Euler equations for describing the dynamics of the fluid mixture.The coupled system is discretized in the reference element and we use a kind of Taylor expansion basis functions to construct the interpolation polynomials of the variables.We show the property that the material derivatives of the basis functions in the DG discretization are equal to zero,with which the scheme is simplified.In addition,the mesh velocity in the ALE framework is obtained by using the adaptive mesh method from[H.Z.Tang and T.Tang,Adaptive mesh methods for one-and two-dimensional hyperbolic conservation laws,SIAMJ.NUMER.ANAL].This adaptivemesh method can automatically concentrate the mesh nodes near the regions with large gradient values and greatly reduces the numerical dissipation near the material interfaces in the simulations.With the help of this adaptivemesh method,the resolution of the solution near the target regions can be greatly improved and the computational efficiency of the simulation is increased.Our scheme can be applied in the simulations for the gas andwatermedia efficiently,and it ismore concise compared to some other methods such as the indirect ALE methods.Several examples including the gas-water flow problemare presented to demonstrate the efficiency of our scheme,and the results show that our scheme can capture the wave structures sharply with high robustness.
基金supported by National Natural Science Foundation of China(Grant Nos.62473277,62473275,62133004,52105072,and 62073230)Jiangsu Provincial Outstanding Youth Program(Grant No.BK20230072)+5 种基金National Key R&D Program of China(Grant Nos.2022YFC3802302 and 2023YFB4705600)Suzhou Industrial Foresight and Key Core Technology Project(Grant No.SYC2022044)Zhejiang Provincial Natural Science Foundation of China(Grant No.LZ24E050004)Shenzhen Polytechnic High-level Talent Start-up Project(Grant No.6023330006K)Shenzhen Science and Technology Program(Grant No.JCYJ20210324132810026)a Grant from Open Foundation of the State Key Laboratory of Fluid Power and Mechatronic Systems,Grants from Jiangsu QingLan Project and Jiangsu 333 high-level talents.
文摘Soft robots, inspired by the flexibility and versatility of biological organisms, have potential in a variety of applications. Recent advancements in magneto-soft robots have demonstrated their abilities to achieve precise remote control through magnetic fields, enabling multi-modal locomotion and complex manipulation tasks. Nonetheless, two main hurdles must be overcome to advance the field: developing a multi-component substrate with embedded magnetic particles to ensure the requisite flexibility and responsiveness, and devising a cost-effective,straightforward method to program three-dimensional distributed magnetic domains without complex processing and expensive machinery. Here, we introduce a cost-effective and simple heat-assisted in-situ integrated molding fabrication method for creating magnetically driven soft robots with three-dimensional programmable magnetic domains. By synthesizing a composite material with neodymium-iron-boron(NdFeB) particles embedded in a polydimethylsiloxane(PDMS) and Ecoflex matrix(PDMS:Ecoflex = 1:2 mass ratio, 50% magnetic particle concentration), we achieved an optimized balance of flexibility, strength, and magnetic responsiveness. The proposed heat-assisted in-situ magnetic domains programming technique,performed at an experimentally optimized temperature of 120℃, resulted in a 2 times magnetization strength(9.5 mT) compared to that at 20℃(4.8 m T), reaching a saturation level comparable to a commercial magnetizer. We demonstrated the versatility of our approach through the fabrication of six kinds of robots, including two kinds of two-dimensional patterned soft robots(2D-PSR), a circular six-pole domain distribution magnetic robot(2D-CSPDMR), a quadrupedal walking magnetic soft robot(QWMSR), an object manipulation robot(OMR), and a hollow thin-walled spherical magneto-soft robot(HTWSMSR). The proposed method provides a practical solution to create highly responsive and adaptable magneto-soft robots.
基金supported by the Jiangsu Provincial Basic Research Program(Natural Science Foundation)Youth Fund(Grant No.BK20230885)the International Joint Laboratory of Sustainable Manufacturing,Ministry of Education and Fundamental Research Funds for Central Universities(Grant No.NG2024012)Major Project on Fundamental Research of Aero-Engines and Gas Turbines,Ministry of Industry and Information Technology Special Project on High-Quality Development(Grant No.J2022-Ⅶ-0006-0048)。
文摘Multi-material 3D sand printing has gained significant attention;however,research has mainly focused on materials and mechanisms,with limited exploration of optimizing the sand-laying process through numerical simulations.In this study,we investigated the dynamic behavior of sand particles during a vibratory sand-laying process for multi-material additive manufacturing using discrete element simulations.The objective is to enable precise control over the amount and distribution of sand for multi-material printing.In this study,we combined experiments and simulations to calibrate the contact parameters of different sands and establish a relationship between the curing agent content and surface energy of sand particles.A model for the vibratory fall of multimaterial sand was developed to study the motion characteristics of sand particles.This allows for macro-control over the sand spreading flow and high-quality multi-material sand laying.The results show that the flow rate of falling sand increases with decreasing surface energy of the particles,wider spreader openings,and higher vibration frequencies.For silica and chromite sands,when their surface energy ranged from 0.15 to 25 J·m^(2)and0.01-0.03 J·m^(2),respectively,and the sand spreader opening was 6 mm with a vibration frequency of 500 Hz,the sand flow rates of both materials became nearly identical.However,a higher sand paving speed and height increased the scattering of sand particles outside the target area,thereby decreasing the paving quality.The results accomplished in this study enable precise and uniform sand particle deposition and offers guidelines for optimizing sand speed and height,thus expanding the application of multi-material sand 3D printing in complex and high-performance manufacturing.
基金supported by the National Natural Science Foundation of China(Grant No.52172356)the Hunan Provincial Natural Science Foundation of China(Grant No.2022JJ10012).
文摘The application of multi-material topology optimization affords greater design flexibility compared to traditional single-material methods.However,density-based topology optimization methods encounter three unique challenges when inertial loads become dominant:non-monotonous behavior of the objective function,possible unconstrained characterization of the optimal solution,and parasitic effects.Herein,an improved Guide-Weight approach is introduced,which effectively addresses the structural topology optimization problem when subjected to inertial loads.Smooth and fast convergence of the compliance is achieved by the approach,while also maintaining the effectiveness of the volume constraints.The rational approximation of material properties model and smooth design are utilized to guarantee clear boundaries of the final structure,facilitating its seamless integration into manufacturing processes.The framework provided by the alternating active-phase algorithm is employed to decompose the multi-material topological problem under inertial loading into a set of sub-problems.The optimization of multi-material under inertial loads is accomplished through the effective resolution of these sub-problems using the improved Guide-Weight method.The effectiveness of the proposed approach is demonstrated through numerical examples involving two-phase and multi-phase materials.
文摘To investigate the influence of Al-Zn-Mg-Cu alloy with as-homogenized and as-rolled initial microstructures on the tensile flow behavior,isothermal tensile tests were conducted on a GLEEBLE-3500 isothermal simulator at temperatures of 380-440℃and strain rates of 0.05-1 s^(−1).The Johnson-Cook model,Hensel-Spittel model,strain-compensated Arrhenius model,and critical fracture strain model were established.Results show that through the evaluation of the models using the correlation coefficient(R)and the average absolute relative error,the strain-compensated Arrhenius model can represent the flow behavior of the alloy more accurately.Shear bands are more pronounced in the as-homogenized specimens,whereas dynamic recrystallization is predominantly observed in as-rolled specimens.Fracture morphology analysis reveals that a mixed fracture mechanism is prevalent in the as-homogenized specimen,whereas a ductile fracture mechanism is predominant in the as-rolled specimen.The processing maps indicate that the unstable region is reduced in the as-rolled specimens compared with that in the as-homogenized specimens.The optimal hot working windows for the as-homogenized and as-rolled specimens are determined as 410-440℃/0.14-1 s^(−1)and 380-400℃/0.05-0.29 s^(−1),respectively.
基金supported by the University of Seoul’s 2025 Research Fund.
文摘Vanadium redox flow batteries(VRFBs)are a means of large-scale energy storage due to their excellent scalability,safety,long cycling life,and decoupled power and energy capacities.However,the slow redox kinetics of vanadium species on conventional carbon electrodes remains a major limitation to their performance.We investigated the deposition of carbon black,carbon nanotubes,and electrochemically exfoliated graphene(Exf-Gr)onto thermally-activated carbon paper(ACP)by spray coating to increase the electrode electrocatalytic activity.The modified electrodes were characterized using scanning electron microscopy,X-ray diffraction,Raman spectroscopy,X-ray photoelectron microscopy,and surface area analysis,while their electrochemical properties were evaluated by cyclic voltammetry,electrochemical impedance spectroscopy,and singlecell VRFB testing.Among the modified electrodes,Exf-Gr/ACP had the best performance,achieving a 2.9-fold reduction in charge transfer resistance compared to pristine ACP and delivering 2.5 times the discharge capacity in single-cell tests.This improvement is attributed to Exf-Gr’s high surface area,favorable catalytic activity,and excellent dispersion on the ACP substrate.Surface modification with electrochemically exfoliated graphene is a highly effective strategy for improving the electrode performance in VRFB systems,with significant implications for large-scale energy storage.
基金supported by the National Natural Science Foundation of China(NSFC)Basic Science Center Program for‘Multiscale Problems in Nonlinear Mechanics’(Grant No.11988102)NSFC project(Grant No.11972038)Chinese Academy of Sciences Project for Young Scientists in Basic Research(Grant No.YSBR-087).
文摘Large-eddy simulation(LES)is conducted to study the statistical properties of mixed-phase turbulence induced by the breaking of bow waves in flow past a partially submerged plate.The simulation is performed using a finite difference method,with the air-water interface captured by a coupled level-set and volume-of-fluid method.Four cases are conducted to investigate the effects of Froude number on turbulent statistics,including the mean velocity,turbulence kinetic energy,and turbulence mass flux(TMF),which is an additional unclosed term in the Reynolds-averaged momentum equation.The TMF,especially its vertical component,shows a complex behaviour with respect to the Froude number.This property of the TMF imposes high demands on the robustness of the closure model of TMF.The present LES data is further used to examine a closure model of the TMF production term,which shows a high correlation with the data obtained from LES.
基金supported by the National Science and Technology Major Project of China(Grant No.2024ZD 1004302)the Key Scientific and Technological Research project of SINOPEC(Grant No.P25186).
文摘To clarify fluid flow mechanisms and establish effective development conditions in continental shale oil reservoirs,a high-temperature,high-pressure steady-state flow system integrated with nuclear magnetic resonance(NMR)technology has been developed.The apparatus combines sample evacuation,rapid pressurization and saturation,and controlled displacement,enabling systematic investigation of single-phase shale oil flow under representative reservoir conditions.Related experiments allow proper quantification of the activation thresholds and relative contributions of different pore types to flow.A movable fluid index(MFI),defined using dual T_(2) cutoff values,is introduced accordingly and linked to key flow parameters.The results reveal distinct multi-scale characteristics of single-phase shale oil transport,namely micro-scale graded displacement and macro-scale segmented nonlinear behavior.As the injection-production pressure difference increases,flow pathways are activated progressively,beginning with fractures,followed by large and then smaller macropores,leading to a pronounced enhancement in apparent permeability.Although mesopores and micropores contribute little to direct flow,their indirect influence becomes increasingly important,and apparent permeability gradually approaches a stable limit at higher pressure difference.It is also shown that the MFI exhibits a strong negative correlation with the starting pressure gradient and a positive correlation with apparent permeability,providing a rapid and reliable indicator of shale oil flow capacity.Samples containing through-going fractures display consistently higher MFI values and superior flowability compared with those dominated by laminated fractures,highlighting the pivotal role of well-connected fracture networks generated by large-scale hydraulic fracturing in improving shale oil production.
基金financially supported by Youth Science“Research on Failure Mechanism and Evaluation Method of Sand Control Measures for Railway Machinery in Sandy Area”(12302511)Ningxia Transportation Department Science and Technology Project(20200173)The Central Guidance on Local Science and Technology Development Funds(22ZY1QA005)。
文摘Investigating the wind-sand flow response regularity in the longitudinal slope sections of desert highways provides a scientific basis for selecting the slope of desert roads.This study uses the Tengger Desert section of the Wuhai-Maqin Expressway as a case study,employing CFD numerical simulation methods to calculate and analyze the wind-sand flow field distribution characteristics in different longitudinal slope sections.The results show that:(1)Along with the direction of the incoming flow,the windward and leeward slope toes of the embankment are low-wind-speed zones,with the wind speed at the leeward slope toe being even lower.The higher the embankment,the larger the low-wind-speed zone at the windward and leeward slope toes.As the longitudinal slope increases,the extent of the lowwind-speed zone at the same location along the route also increases.(2)Along the route direction,the wind speed at the windward and leeward slope toes decreases as embankment height increases.At the embankment toe,sand particles are transported from the top to the bottom of the longitudinal slope,and the greater the longitudinal slope,the stronger the transport effect.(3)Along the route direction,the sand accumulation around the embankment gradually gathers toward the bottom of the longitudinal slope as the slope increases.When the longitudinal slope is 3%and 4%,the trend of sand accumulation moving from the windward side at the end of the route to the leeward side at the start of the route is more significant.When the longitudinal slope is less than or equal to 3%,severe sand accumulation within the embankment range is reduced by 86.4%or more compared to when the slope is 4%.(4)Under the same longitudinal slope,the higher the embankment height,the smaller its transport rate.When the embankment height is the same,the greater the longitudinal slope,the greater the embankment transport rate.
文摘The newly formulated non-Newtonian rivulet flows streaming down an inclined planar surface,with additional periodic perturbations arising from the application of the 2nd Stokes problem to the investigation of rivulet dynamics,are demonstrated in the current research.Hereby,the 2nd Stokes problem assumes that the surface,with a thin shared layer of the fluid on it,oscillates in a harmonic manner along the x-axis of the rivulet flow,which coincides with the main flow direction streaming down the underlying surface.We obtain the exact extension of the rivulet flow family,clarifying the structure of the pressure field,which fully absorbs the arising perturbation.The profile of the velocity field is assumed to be Gaussian-type with a non-zero level of plasticity.Hence,the absolutely non-Newtonian case of the viscoplastic flow solution,which satisfies the motion and continuity equations,is considered(with particular cases of exact solutions for pressure).The perturbed governing equations of motion for rivulet flows then result in the Riccati-type ordinary differential equation(ODE),describing the dynamics of the coordinate x(t).The approximated schematic dynamics are presented in graphical plots.
文摘Particle-and droplet-laden flows are central to many problems in mechanics and transport.They occur in sedimentladen boundary layers,gas-solid and gas-liquid dispersions,and surface-water films driven by external forcing.They also underpin practical applications ranging from environmental transport to high-speed and aerothermal systems.Despite decades of progress,prediction remains difficult.The physics spans a wide range of scales and often couples turbulence,interphase momentum exchange,collisions,and interfacial transport.Reliable computation therefore requires both robust numerical methodology and careful physical interpretation.
基金supported in part by the National Key R&D Program of China(No.2024YFA1610700)the National Natural Science Foundation of China(No.12475147)。
文摘The Cooling-Storage-Ring External-target Experiment(CEE)at the Heavy Ion Research Facility in Lanzhou(HIRFL)is designed to study the properties of nuclear matter created in heavy-ion collisions at beam energies from a few hundred MeV/u up to 1 GeV/u.It aims to facilitate research on the quantum chromodynamics(QCD)phase structure in the high-baryondensity region.Collective flow is a fundamental observable in heavy-ion collision experiments,providing information on the bulk properties of the produced matter.Although the standard event plane method has been widely used to measure collective flow,it is still important to validate and optimize this method for the CEE spectrometer.In this paper,we study the experimental procedures for measuring directed flow in^(238)U+^(238)U collisions at 500 MeV/u,using event planes reconstructed by Multi-Wire Drift Chamber and Zero Degree Calorimeter,respectively.Jet AA Microscopic(JAM)transport generator is used to generate events,and the detector response is simulated by the CEE Fast Simulation(CFS)package.Finally,the optimal kinematic region for proton directed flow measurements is discussed for the future CEE experiment.
基金Project supported by the National Natural Science Foundation of China(Nos.12372214 and U2341231)。
文摘The Reynolds-averaged Navier-Stokes(RANS)technique enables critical engineering predictions and is widely adopted.However,since this iterative computation relies on the fixed-point iteration,it may converge to unexpected non-physical phase points in practice.We conduct an analysis on the phase-space characteristics and the fixed-point theory underlying the k-ε turbulence model,and employ the classical Kolmogorov flow as a framework,leveraging its direct numerical simulation(DNS)data to construct a one-dimensional(1D)system under periodic/fixed boundary conditions.The RANS results demonstrate that under periodic boundary conditions,the k-ε model exhibits only a unique trivial fixed point,with asymptotes capturing the phase portraits.The stability of this trivial fixed point is determined by a mathematically derived stability phase diagram,indicating the fact that the k-ε model will never converge to correct values under periodic conditions.In contrast,under fixed boundary conditions,the model can yield a stable non-trivial fixed point.The evolutionary mechanisms and their relationship with boundary condition settings systematically explain the inherent limitations of the k-ε model,i.e.,its deficiency in computing the flow field under periodic boundary conditions and sensitivity to boundary-value specifications under fixed boundary conditions.These conclusions are finally validated with the open-source code OpenFOAM.
文摘This study investigates the enhancement of convective heat transfer in a serpentine pipe using ferrofluid flow influenced by dual non-uniform magnetic sources.The primary objective is to improve thermal performance in compact cooling systems,such as those used in heat exchangers.A two-dimensional,steady-state Computational Fluid Dynamic(CFD)model is developed in ANSYS Fluent to simulate the behavior of an incompressible ferrofluid under applied constant heat flux and magnetic fields.The magnetic force is modeled using the Kelvin force,which acts on magnetized nanoparticles in response to spatially varying electromagnetic fields generated by two strategically positioned current-carrying wires.The effects of magnetic field strength,quantified by the magnetic number(Mn),on flow behavior and temperature distribution are thoroughly analyzed.The results indicate that increasing Mn leads to higher Nusselt numbers,demonstrating enhanced convective heat transfer.Secondary vortices induced by magnetic forcing improve fluid mixing,particularly in curved regions of the pipe.A mesh-independence study and model validation with benchmark data support the reliability of the numerical framework.This work highlights the potential of magnetic-field-assisted thermal control in energy-efficient cooling applications and provides a foundation for the further development of advanced ferrofluid-based heat transfer systems.
基金the Science and Technology R&D Major Project of Jiangxi Province(No.20244AFI92001)the National Natural Science Foundation of China(Nos.22071033 and 21801047)for the financial supports.
文摘We report an immobilized enzyme-catalyzed batch and continuous-flow synthesis of optically pure ethyl(R)-pantothenate((R)-PaOEt),the direct precursor of d-pantothenic acid.Firstly,a ketoreductase mutant designated as M2,carrying two-point mutations of F97L and M242F relative to the wild-type SSCR,was constructed by site-directed mutagenesis,exhibited simultaneously improved activity toward ethyl 2′-ketopantothenate(K-PaOEt)and isopropanol,and could effectively catalyze the stereoselective reduction of K-PaOEt to(R)-PaOEt by using isopropanol as the sacrificial co-substrate to regenerate NADPH.After screening six commercially available carriers,an amino resin LXTE-700 was identified as the best solid support for the immobilization of M2 via the glutaraldehyde activation method.Upon optimization of the immobilization process and reaction conditions,the fabricated immobilized enzyme M2@amino resin demonstrated excellent recyclability and reusability,with the complete conversion of K-PaOEt to(R)-PaOEt being still realized after 12 cycles of reuse.Finally,M2@amino resin-catalyzed synthesis of(R)-PaOEt was successfully implemented in continuous-flow,accomplishing a 6.3 times higher space-time yield than that with the batch synthesis(529.2 versus 84 g L^(-1) d^(-1)).Our developed flow biocatalysis system also features an outstanding operational stability,as evidenced by the 100%conversion rate achieved after 15 consecutive days of operation.
基金funded by United Kingdom Commonwealth Scholarship Commission。
文摘The Rukwa Rift section of the East Africa Rift System presents a type setting for radiogenic helium accumulation in a petroleum free basin.As a prerequisite for accumulation,a considerable high heat flow anomaly is required from tectonothermal events to drive the release and circulation of radiogenic helium in the continental crust.Here we apply statistical analysis on geochemical data observed in thermal springs and recorded heat flow to account for crustal helium mass balance for each tectonothermal event in the region.Our results demonstrate anomalously high heat flow~64-99 mW/m^(2) with a consistent trend of helium isotopic ratio and fluid chemistry in the Rukwa Rift.Mass balance calculation show that the whole crustal volume underlying the East Africa Helium Pool(EAHP)has a capability of producing radiogenic helium of about 9.9×10^(6) mol/yr(22×10^(-6) mol 4He/m^(2) yr)while the total radiogenic helium flux ranges between~2.39×10^(6) mol/yr and~2.68×10^(9) mol/yr.The Tanzania Craton contributes largely to radiogenic helium releasing up to 50% of the total capacity in the region.The total ^(4)He emission in the Rukwa Rift Basin is about 4.45×10^(5)-5.01×10^(8) mol/yr which is thus equivalent to 19%-21% of the total production capacity in the region.These results imply that the helium accumulation in the EAHP would have started as early as Paleoproterozoic(2.349 Ga).These results provide a qualitative and quantitative insight to assess both helium and geothermal potentiality in the region.
