Multi-Material Topology Optimization for Spatial-Varying Porous Structures

This paper aims to propose a topology optimization method on generating porous structures comprising multiple materials.The mathematical optimization formulation is established under the constraints of individual volume fraction of constituent phase or total mass,as well as the local volume fraction of all phases.The original optimization problem with numerous constraints is converted into a box-constrained optimization problem by incorporating all constraints to the augmented Lagrangian function,avoiding the parameter dependence in the conventional aggregation process.Furthermore,the local volume percentage can be precisely satisfied.The effects including the globalmass bound,the influence radius and local volume percentage on final designs are exploited through numerical examples.The numerical results also reveal that porous structures keep a balance between the bulk design and periodic design in terms of the resulting compliance.All results,including those for irregular structures andmultiple volume fraction constraints,demonstrate that the proposedmethod can provide an efficient solution for multiple material infill structures.

Nano-Au-decorated hierarchical porous cobalt sulfide derived from ZIF-67 toward optimized oxygen evolution catalysis:Important roles of microstructures and electronic modulation

Enhancing both the number of active sites available and the intrinsic activity of Co-based electrocatalysts simultaneously is a desirable goal.Herein,a ZIF-67-derived hierarchical porous cobalt sulfide decorated by Au nanoparticles(NPs)(denoted as HP-Au@CoxSy@ZIF-67)hybrid is synthesized by low-temperature sulfuration treatment.The well-defined macroporous-mesoporous-microporous structure is obtained based on the combination of polystyrene spheres,as-formed CoxSy nanosheets,and ZIF-67 frameworks.This novel three-dimensional hierarchical structure significantly enlarges the three-phase interfaces,accelerating the mass transfer and exposing the active centers for oxygen evolution reaction.The electronic structure of Co is modulated by Au through charge transfer,and a series of experiments,together with theoretical analysis,is performed to ascertain the electronic modulation of Co by Au.Meanwhile,HP-Au@CoxSy@ZIF-67 catalysts with different amounts of Au were synthesized,wherein Au and NaBH4 reductant result in an interesting“competition effect”to regulate the relative ratio of Co^(2+)/Co^(3+),and moderate Au assists the electrochemical performance to reach the highest value.Consequently,the optimized HP-Au@CoxSy@ZIF-67 exhibits a low overpotential of 340 mV at 10 mA cm^(-2)and a Tafel slope of 42 mV dec-1 for OER in 0.1 M aqueous KOH,enabling efficient water splitting and Zn-air battery performance.The work here highlights the pivotal roles of both microstructural and electronic modulation in enhancing electrocatalytic activity and presents a feasible strategy for designing and optimizing advanced electrocatalysts.

A novel multi-channel porous structure facilitating mass transport towards highly efficient alkaline water electrolysis

An advantageous porous architecture of electrodes is pivotal in significantly enhancing alkaline water electrolysis(AWE)efficiency by optimizing the mass transport mechanisms.This effect becomes even more pronounced when aiming to achieve elevated current densities.Herein,we employed a rapid and scalable laser texturing process to craft novel multi-channel porous electrodes.Particularly,the obtained electrodes exhibit the lowest Tafel slope of 79 mV dec^(-1)(HER)and 49 mV dec^(-1)(OER).As anticipated,the alkaline electrolyzer(AEL)cell incorporating multi-channel porous electrodes(NP-LT30)exhibited a remarkable improvement in cell efficiency,with voltage drops(from 2.28 to 1.97 V)exceeding 300 mV under 1 A cm^(-1),compared to conventional perforated Ni plate electrodes.This enhancement mainly stemmed from the employed multi-channel porous structure,facilitating mass transport and bubble dynamics through an innovative convection mode,surpassing the traditional convection mode.Furthermore,the NP-LT30-based AEL cell demonstrated exceptional durability for 300 h under 1.0 A cm^(-2).This study underscores the capability of the novel multi-channel porous electrodes to expedite mass transport in practical AWE applications.

Modeling of multiphase flow in low permeability porous media:Effect of wettability and pore structure properties

Multiphase flow in low permeability porous media is involved in numerous energy and environmental applications.However,a complete description of this process is challenging due to the limited modeling scale and the effects of complex pore structures and wettability.To address this issue,based on the digital rock of low permeability sandstone,a direct numerical simulation is performed considering the interphase drag and boundary slip to clarify the microscopic water-oil displacement process.In addition,a dual-porosity pore network model(PNM)is constructed to obtain the water-oil relative permeability of the sample.The displacement efficiency as a recovery process is assessed under different wetting and pore structure properties.Results show that microscopic displacement mechanisms explain the corresponding macroscopic relative permeability.The injected water breaks through the outlet earlier with a large mass flow,while thick oil films exist in rough hydrophobic surfaces and poorly connected pores.The variation of water-oil relative permeability is significant,and residual oil saturation is high in the oil-wet system.The flooding is extensive,and the residual oil is trapped in complex pore networks for hydrophilic pore surfaces;thus,water relative permeability is lower in the water-wet system.While the displacement efficiency is the worst in mixed-wetting systems for poor water connectivity.Microporosity negatively correlates with invading oil volume fraction due to strong capillary resistance,and a large microporosity corresponds to low residual oil saturation.This work provides insights into the water-oil flow from different modeling perspectives and helps to optimize the development plan for enhanced recovery.

Dependence of lithium metal battery performances on inherent separator porous structure regulation

Boosting of rechargeable lithium metal batteries(LMBs) holds challenges because of lithium dendrites germination and high-reactive surface feature.Separators may experience structure-determined chemical deterioration and worsen Li plating-stripping behaviors when smoothly shifting from lithium-ion batteries(LIBs) to LMBs.This study precisely regulations the crystal structure of β-polypropylene and separator porous construction to investigate the intrinsic porous structure and mechanical properties determined electrochemical performances and cycling durability of LMBs.Crystal structure characterizations,porous structure analyses,and electrochemical cycling tests uncover appropriate annealing thermal stimulation concentrates β-lamellae thickness and enhances lamellae thermal stability by rearranging molecular chain in inferior β-lamellae,maximally homogenizing biaxial tensile deformation and resultant porous constructions.These even pores with high connectivity lower ion migration barriers,alleviate heterogeneous Li^(+) flux dispersion,stabilize reversible Li plating-stripping behaviors,and hinder coursing and branching of Li dendrites,endowing steady cell cycling durability,especially at higher currents due to the highlighted uncontrollable cumulation of dead Li,which offers new insights for the current pursuit of high-power density battery and fast charging technology.The suggested separator structure-chemical nature functions in ensuring cyclic cell stability and builds reliable relationships between separator structure design and practical LMBs applications.

High-Throughput Growth of Hexagonal Boron Nitride Film Using Porous-Structure Isolation Layer

Chemical vapor deposition is considered as the most hopeful method for the synthesis of large-area high-quality hexagonal boron nitride on the substrate of catalytic metal. However, the size the hexagonal boron nitride films are limited to the size of growth chamber, which indicates a lower production efficiency. In this paper, the utilization efficiency of growth chamber is highly improved by alternately stacking multiple pieces of Cu foils and carbon fiber surface felt with porous structure. Uniform and continuous hexagonal boron nitride films are prepared on Cu foils through chemical vapor deposition utilizing ammonia borane as the precursor. This work develops a simple and practicable method for high-throughput preparation of hexagonal boron nitride films, which could contribute to the industrial application of hexagonal boron nitride. .

Velocity-Free MS/AE Source Location Method for Three-Dimensional Hole-Containing Structures

Microseismic/acoustic emission(MS/AE)source localization method is crucial for predicting and controlling of potentially dangerous sources of complex structures.However,the locating errors induced by both the irregular structure and pre-measured velocity are poorly understood in existing methods.To meet the high-accuracy locating requirements in complex three-dimensional hole-containing structures,a velocity-free MS/AE source location method is developed in this paper.It avoids manual repetitive training by using equidistant grid points to search the path,which introduces A*search algorithm and uses grid points to accommodate complex structures with irregular holes.It also takes advantage of the velocity-free source location method.To verify the validity of the proposed method,lead-breaking tests were performed on a cubic concrete test specimen with a size of 10 cm10 cm10 cm.It was cut out into a cylindrical empty space with a size of/6cm10 cm.Based on the arrivals,the classical Geiger method and the proposed method are used to locate lead-breaking sources.Results show that the locating error of the proposed method is 1.20 cm,which is less than 2.02 cm of the Geiger method.Hence,the proposed method can effectively locate sources in the complex three-dimensional structure with holes and achieve higher precision requirements.

Observed spatiotemporal variation of three-dimensional structure and heat/salt transport of anticyclonic mesoscale eddy in Northwest Pacific

As in-situ observations are sparse,targeted observations of a specific mesoscale eddy are rare.Therefore,it is difficult to study the three-dimensional structure of moving mesoscale eddies.From April to September 2014,an anticyclonic eddy located at 135°E-155°E,26°N-42°N was observed using 17 rapidsampling Argo floats,and the spatiotemporal variations in the three-dimensional structure were studied.The results are as follows:(1)the eddy was identified and tracked using satellite altimeter data.It had a lifetime of 269 days and an average radius of 91.5 km.The lifetime of the eddy can be divided into three phases,i.e.,the initiation,maturity,and termination phases.The depth of its influence reached 1000 m;(2)the Argo profiles were divided into seven periods(approximately 20 days in each)for composite analysis,and the composite Argo profiles and CARS2009(CSIRO Atlas of Regional Seas)climatology data were merged following the data-interpolating variational analysis(DIVA)method to reconstruct the three-dimensional structure.The temperature and salinity anomaly cores of the anticyclonic mesoscale eddy are located from 400 to 600 m.From 800 to 900 m,there is an area of low salinity at the center of the eddy.A high concentration anomaly of dissolved oxygen was located at approximately 250 m;(3)to better understand the features of the eddy and its interaction with the surroundings,we calculated the anomalous velocity of the geostrophic flow and the heat,salt,dissolved oxygen transport anomaly,and discussed the eddy's origin and its adjustments to topography.The maximum heat,salt,and oxygen transport caused by eddy were 9.37×10^11 W,3.08×10^3 kg/s,and 2.70×10^2 kg/s,which all occurred during the termination phase.This study highlights the applicability of using Argo floats to understand the three-dimensional structure thermohaline features of eddies in the North Pacific.

Three-dimensional simulation of a Ka-band relativistic Cherenkov source with metal photonic-band-gap structures

This paper presents a three-dimensional particle-in-cell （PIC） simulation of a Ka-band relativistic Cherenkov source with a slow wave structure （SWS） consisting of metal photonic band gap （PBG） structures. In the simulation, a perfect match layer boundary is employed to absorb passing band modes supported by the PBG lattice with an artificial metal boundary. The simulated axial field distributions in the cross section and surface of the SWS demonstrate that the device operates in the vicinity of the π point of a TM01-1ike mode. The Fourier transformation spectra of the axial fields as functions of time and space show that only a single frequency appears at 36.27 GHz, which is in good agreement with that of the intersection of the dispersion curve with the slow space charge wave generated on the beam. The simulation results demonstrate that the SWS has good mode selectivity.

The three-dimensional structure and seasonal variation of the North Pacific meridional overturning circulation

The three-dimensional structure and the seasonal variation of the North Pacific meridional overturning circulation （NPMOC） are analyzed based on the Simple Ocean Data Assimilation data and Argo profiling float data. The NPMOC displays a multi-cell structure with four cells in the North Pacific altogether. The TC and the STC are a strong clockwise meridional cell in the low latitude ocean and a weaker clockwise meridional cell between 7°N and 18°N, respectively, while the DTC and the subpolar cell are a weaker anticlockwise meridional cell between 3°N and 15°N and a weakest anticlockwise meridional cell between 35°N and 50°N, respectively. The DTC, the TC and the STC are all of very strong seasonal variations. As to the DTC, the southward transport is strongest in fall and weakest in spring. For the TC, the northward transport is strongest in winter and weakest in spring, while the southward transport is strongest in fall and weakest in spring, which is associated with the strong southward fiow of the DTC in fall. As the STC, the northward transport is strongest in winter and weakest in summer, while the southward transport is strongest in summer and weakest in spring. This seasonal difference may be associated with the DTC. The zonal wind stress and the east-west slope of sea level play important roles in the seasonal variations of the TC, the STC and the DTC.

Three-dimensional crustal P-wave velocity structure in the Yangbi and Eryuan earthquake regions, Yunnan, China

A magnitude 5.5 earthquakes occurred in Eryuan County,Dali Bai Autonomous Prefecture,Yunnan Province,China,on March 3.And a magnitude 5.0 earthquake occurred in the same place on April 17,2013,i.e.,45 days later.Then,on May 21,2021,multiple earthquakes,one with magnitude 6.4 and several at 5.0 or above,occurred in Yangbi County,Dali Bai Autonomous Prefecture,Yunnan Province,China.All of these occurred in the Weixi-QiaohouWeishan fault zone.In this study,1,874 seismic events in Yangbi and Eryuan counties were identified by automatic micro-seismic identification technology and the first arrivals were picked up manually.Following this,a total of 11,968 direct P-wave absolute arrivals and 73,987 high-quality Pwave relative arrivals were collected for joint inversion via the double difference tomography method.This was done to obtain the regional three-dimensional fine crustal P-wave velocity structure.The results show that the travel time residuals before and after inversion decreased from the initial–0.1–0.1 s to–0.06–0.06 s.The upper crust in the study area,which exhibited a low-velocity anomaly,corresponded to the basin region;this indicated that the low-velocity anomaly in the shallow part of the study area was affected by the basin.Results also showed some correlation between the distribution of the earthquakes and velocity structure,as there was a lowvelocity body Lv1 with a wide distribution at depths ranging from 15–20 km in the Yangbi and Eryuan earthquake regions.In addition,earthquakes occurred predominantly in the highlow velocity abnormal transition zone.The low-velocity body in the middle and lower crust may be prone to concentrating upper crustal stress,thus leading to the occurrence of earthquakes.

Three-dimensional visualization of soil pore structure using computed tomography

The geometric and spatial characteristics of pore structures determine the permeability and water retention of soils, which have important effects on soil functional diversity and ecological restoration. Until recently, there have not been tools and methods to visually and quantitatively describe the characteristics of soil pores. To solve this problem, this research reconstructs the geometry and spatial distribution of soil pores by the marching cubes method, texture mapping method and the ray casting method widely used in literature. The objectives were to explore an optimal method for three-dimensional visualization of soil pore structure by comparing the robustness of the three methods on soil CT images with single pore structure and porosity ranging from low (2–5%) to high (12–18%), and to evaluate the reconstruction performance of the three methods with different geometric features. The results demonstrate that there are aliases (jagged edges) and deficiency at the boundaries of the model reconstructed by the marching cubes method and pore volumes are smaller than the ground truth, whereas the results of the texture mapping method lack the details of pore structures. For all the soil images, the ray casting method is preferable since it better preserves the pore characteristics of the ground truth. Furthermore, the ray casting method produced the best soil pore model with higher rendering speed and lower memory consumption. Therefore, the ray casting method provides a more advanced method for visualization of pore structures and provides an optional technique for the study of the transport of moisture and the exchange of air in soil.

The Shengli I Point Bar on the Yellow River Delta: Three-Dimensional Structures and Their Evolution

Point bars are well developed on the Yellow River delta, an~ which theShengli I point bar is the most typical. The point bar, being about 4 km in length and several tensto more than 100 meters in width, is located on the south side of the Shengli Bridge in KenliCounty, Dongying, Shandong. It is a typical fine-grained point bar with silt, which is predominant,some clay and minor plant debris and clay boulders. The Shengli I Point bar has complicated 3-Dstructures. Firstly, in a plane view, it comprises mainly eight sedimentary units, bar edge, baredge, bar platform, bar plain, bar channel, bar gully, bar pond and bar bay, developing side by sideand superimposed one by one m a complex way. Secondly, its vertical structures are very complex dueto the partial superimposition of the 8 sedimentary units. Besides hydatogenesis, very intensivewind erosion, eolian, ice and meltwater actions are also visible on the Shengli I point bar. Thecomplex form is made even more complicated because of the above co-actions.

Three-dimensional structure of an observed cyclonic mesoscale eddy in the Northwest Pacific and its assimilation experiment

Mesoscale eddies play an important role in modulating the ocean circulation.Many previous studies on the threedimensional structure of mesoscale eddies were mainly based on composite analysis,and there are few targeted observations for individual eddies.A cyclonic eddy surveyed during an oceanographic cruise in the Northwest Pacific Ocean is investigated in this study.The three-dimensional structure of this cyclonic eddy is revealed by observations and simulated by the four-dimensional variational data assimilation(4 DVAR)system combined with the Regional Ocean Modeling System.The observation and assimilation results together present the characteristics of the cyclonic eddy.The cold eddy has an obvious dual-core structure of temperature anomaly.One core is at 50–150 m and another is at 300–550 m,which both have the average temperature anomaly of approximately-3.5℃.The salinity anomaly core is between 250 m and 500 m,which is approximately-0.3.The horizontal velocity structure is axis-asymmetric and it is enhanced on the eastern side of the cold eddy.In the assimilation experiment,sea level anomaly,sea surface temperature,and in situ measurements are assimilated into the system,and the results of assimilation are close to the observations.Based on the high-resolution assimilation output results,the study also diagnoses the vertical velocity in the mesoscale eddy,which reaches the maximum of approximately 10 m/d.The larger vertical velocity is found to be distributed in the range of 0.5 to 1 time of the normalized radius of the eddy.The validation of the simulation result shows that the 4 DVAR method is effective to reconstruct the three-dimensional structure of mesoscale eddy and the research is an application to study the mesoscale eddy in the Northwest Pacific by combining observation and assimilation methods.

Three-dimensional P velocity structure in Beijing area

A detail three-dimensional P wave velocity structure of Beijing, Tianjin and Tangshan area (BTT area) was deter-mined by inverting local earthquake data. In total 16 048 P wave first arrival times from 16048 shallow and mid-depth crustal earthquakes, which occurred in and around the BTT area from 1992 to 1999 were used. The first arrival times are recorded by Northern China United Telemetry Seismic Network and Yanqing-Huailai Digital Seismic Network. Hypocentral parameters of 1 132 earthquakes with magnitude ML=1.7~6.2 and the three-dimensional P wave velocity structure were obtained simultaneously. The inversion result reveals the com-plicated lateral heterogeneity of P wave velocity structure around BTT area. The tomographic images obtained are also found to explain other seismological observations well.

Three-dimensional discrete element numerical simulation of Paleogene salt structures in the western Kuqa foreland thrust belt

Taking the Paleogene salt strata in the west of Kuqa foreland thrust belt as study object, the deformation features of salt structure in the compression direction and perpendicular to the compression direction were examined to find out the control factors and formation mechanisms of the salt structures. By using the three-dimensional discrete element numerical simulation method, the formation mechanisms of typical salt structures of western Kuqa foreland thrust belt in Keshen and Dabei work areas were comprehensively analyzed. The simulation results show that the salt deformation in Keshen and Dabei work areas is of forward spread type, with deformation concentrated in the piedmont zone;the salt deformation is affected by the early uplift near the compression end, pre-existing basement faults, synsedimentary process and the initial salt depocenter;in the direction perpendicular to the compression direction, salt rocks near the compression end have strong lateral mobility with the velocity component moving towards the middle part, and the closer to the middle, the larger the velocity will be, so that salt rocks will aggregate towards the middle and deform intensely, forming complex folds and separation of salt structures from salt source, and local outcrop with thrust faults. Compared with 2 D simulation, 3 D simulation can analyze salt structures in the principal stress direction and direction perpendicular to the principal stress, give us a full view of the formation mechanisms of salt structures, and guide the exploration of oil and gas reservoirs related to salt structures.

Synthesis and Crystal Structure of a Three-dimensional Mn(Ⅱ) Coordination Polymer with 3-(Pyrazin-2-yloxy)-pyridine and Azide Anion as Mixed Bridge Ligand

A three-dimensional coordination polymer [Mn2（μ1.3-N3）4（μ-PP）2]n （PP = 3-（pyrazin-2-yloxy）-pyridine） has been synthesized with 3-（pyrazin-2-yloxy）-pyridine and azide anion as mixed bridge ligand, and its crystal structure was determined by X-ray crystallography. The crystal data： triclinic system, space group P1, with a = 6.794（4）, b = 9.885（6）, c = 9.947（6） A, α = 64.170（6）, β= 84.190（8）, γ= 85.319（8）°, V = 597.7（6）A^3, Z = 1, C18H14Mn2N18O2, Mr = 624.35, Dc = 1.735 g/cm^3, F（000） = 314 and μ = 1.117 mm^-1. In the crystal, the azide anion acts as a bridge ligand and makes adjacent Mn（Ⅱ） ions connect into a two-dimensional sheet on the ab plane, then 3-（pyrazin-2-yloxy）-pyridine serves as a bidentate bridge ligand to connect neighboring sheets along

Three-Dimensional Thermo-Elastic-Plastic Finite Element Method Modeling for Predicting Weld-Induced Residual Stresses and Distortions in Steel Stiffened-Plate Structures

The objective of the present paper is to develop nonlinear finite element method models for predicting the weld-induced initial deflection and residual stress of plating in steel stiffened-plate structures. For this purpose, three-dimensional thermo-elastic-plastic finite element method computations are performed with varying plate thickness and weld bead length (leg length) in welded plate panels, the latter being associated with weld heat input. The finite element models are verified by a comparison with experimental database which was obtained by the authors in separate studies with full scale measurements. It is concluded that the nonlinear finite element method models developed in the present paper are very accurate in terms of predicting the weld-induced initial imperfections of steel stiffened plate structures. Details of the numerical computations together with test database are documented.

Reconstruction of three-dimensional grain structure in polycrystalline iron via an interactive segmentation method

Using a total of 297 segmented sections,we reconstructed the three-dimensional(3D) structure of pure iron and obtained the largest dataset of 16254 3D complete grains reported to date.The mean values of equivalent sphere radius and face number of pure iron were observed to be consistent with those of Monte Carlo simulated grains,phase-field simulated grains,Ti-alloy grains,and Ni-based super alloy grains.In this work,by finding a balance between automatic methods and manual refinement,we developed an interactive segmentation method to segment serial sections accurately in the reconstruction of the 3D microstructure;this approach can save time as well as substantially eliminate errors.The segmentation process comprises four operations:image preprocessing,breakpoint detection based on mathematical morphology analysis,optimized automatic connection of the breakpoints,and manual refinement by artificial evaluation.

Three-dimensional thermohaline anomaly structures of rings in the Kuroshio Extension region

Using AVISO satellite altimeter observations during 1993–2015 and a manual eddy detection method, a total of 276 anticyclonic rings and 242 cyclonic rings shed from the Kuroshio Extension(KE) were identified, and their three-dimensional(3D) anomaly structures were further reconstructd based on the Argo float data and the Japan Agency for Marine-Earth Science and Technology(JAMSTEC) cruise and buoy data through an interpolation method. It is found that the cyclonic(anticyclonic) rings presented consistent negative(positive) anomalies of potential temperature;meanwhile the relevant maximum anomaly center became increasingly shallow for the cyclonic rings whereas it went deeper for the anticyclonic rings as the potential temperature anomaly decreased from the west to the east. The above deepening or shoaling trend is associated with the zonal change of the depth of the main thermocline. Moreover, the composite cold ring between 140° and 150°E was found to exhibit a double-core vertical structure due to the existence of mode water with low potential vorticity. Specifically, a relatively large negative(positive) salinity anomaly and a small positive(negative) one appeared for the composite cyclonic(anticyclonic) ring at the depth above and below 600 m, respectively. The underlying driving force for the temperature and salinity anomaly of the composite rings was also attempted, which varies depending on the intensity of the background current and the temperature and salinity fields in different areas of the KE region, and the rings’ influences on the temperature and salinity could reach deeper than 1 000 m on average.