Multi-Physics Coupled Acoustic-Mechanics Analysis and Synergetic Optimization for a Twin-Fluid Atomization Nozzle

Fine particulate matter produced during the rapid industrialization over the past decades can cause significant harm to human health.Twin-fluid atomization technology is an effective means of controlling fine particulate matter pollution.In this paper,the influences of the main parameters on the droplet size,effective atomization range and sound pressure level(SPL)of a twin-fluid nozzle(TFN)are investigated,and in order to improve the atomization performance,a multi-objective synergetic optimization algorithm is presented.A multi-physics coupled acousticmechanics model based on the discrete phase model(DPM),large eddy simulation(LES)model,and Ffowcs Williams-Hawkings(FW-H)model is established,and the numerical simulation results of the multi-physics coupled acoustic-mechanics method are verified via experimental comparison.Based on the analysis of the multi-physics coupled acoustic-mechanics numerical simulation results,the effects of the water flow on the characteristics of the atomization flow distribution were obtained.A multi-physics coupled acoustic-mechanics numerical simulation result was employed to establish an orthogonal test database,and a multi-objective synergetic optimization algorithm was adopted to optimize the key parameters of the TFN.The optimal parameters are as follows:A gas flow of 0.94 m^(3)/h,water flow of 0.0237 m^(3)/h,orifice diameter of the self-excited vibrating cavity(SVC)of 1.19 mm,SVC orifice depth of 0.53 mm,distance between SVC and the outlet of nozzle of 5.11 mm,and a nozzle outlet diameter of 3.15 mm.The droplet particle size in the atomization flow field was significantly reduced,the spray distance improved by 71.56%,and the SPL data at each corresponding measurement point decreased by an average of 38.96%.The conclusions of this study offer a references for future TFN research.

Transient multi-physics behavior of an insert high temperature superconducting no-insulation coil in hybrid superconducting magnets with inductive coupling

A transient multi-physics model incorporated with an electromagneto-thermomechanical coupling is developed to capture the multi-field behavior of a single-pancake(SP)insert no-insulation(NI)coil in a hybrid magnet during the charging and discharging processes.The coupled problem is resolved by means of the finite element method(FEM)for the magneto-thermo-elastic behaviors and the Runge-Kutta method for the transient responses of the electrical circuits of the hybrid superconducting magnet system.The results reveal that the transient multi-physics responses of the insert NI coil primarily depend on the charging/discharging procedure of the hybrid magnet.Moreover,a reverse azimuthal current and a compressive hoop stress are induced in the insert NI coil during the charging process,while a forward azimuthal current and a tensile hoop stress are observed during the discharging process.The induced voltages in the insert NI coil can drive the currents flowing across the radial turns where the contact resistance exists.Therefore,it brings forth significant Joule heat,causing a temperature rise and a uniform distribution of this heat in the coil turns.Accordingly,a thermally/mechanically unstable or quenching event may be encountered when a high operating current is flowing in the insert NI coil.It is numerically predicted that a quick charging will induce a compressive hoop stress which may bring a risk of buckling instability in the coil,while a discharging will not.The simulations provide an insight of hybrid superconducting magnets under transient start-up or shutdown phases which are inevitably encountered in practical applications.

Subsurface multi-physical characterization of mountain excavation and city construction in loess plateau with a fiber-optic sensing system

Mountain excavation and city construction(MECC)projects being launched in the Loess Plateau in China involve the creation of large-scale artificial land.Understanding the subsurface evolution characteristics of the artificial land is essential,yet challenging.Here,we use an improved fiber-optic monitoring system for its subsurface multi-physical characterization.The system enables us to gather spatiotemporal distribution of various parameters,including strata deformation,temperature,and moisture.Yan’an New District was selected as a case study to conduct refined in-situ monitoring through a 77 m-deep borehole and a 30 m-long trench.Findings reveal that the ground settlement involves both the deformation of the filling loess and the underlying intact loess.Notably,the filling loess exhibits a stronger creep capability compared to underlying intact loess.The deformation along the profile is unevenly distributed,with a positive correlation with soil moisture.Water accumulation has been observed at the interface between the filling loess and the underlying intact loess,leading to a significant deformation.Moreover,the temperature and moisture in the filling loess have reached a new equilibrium state,with their depths influenced by atmospheric conditions measuring at 31 m and 26 m,respectively.The refined investigation allows us to identify critical layers that matter the sustainable development of newly created urban areas,and provide improved insights into the evolution mechanisms of land creation.

An internal ballistic model of electromagnetic railgun based on PFN coupled with multi-physical field and experimental validation

To accelerate the practicality of electromagnetic railguns,it is necessary to use a combination of threedimensional numerical simulation and experiments to study the mechanism of bore damage.In this paper,a three-dimensional numerical model of the augmented railgun with four parallel unconventional rails is introduced to simulate the internal ballistic process and realize the multi-physics field coupling calculation of the rail gun,and a test experiment of a medium-caliber electromagnetic launcher powered by pulse formation network(PFN)is carried out.Various test methods such as spectrometer,fiber grating and high-speed camera are used to test several parameters such as muzzle initial velocity,transient magnetic field strength and stress-strain of rail.Combining the simulation results and experimental data,the damage condition of the contact surface is analyzed.

Simultaneously realizing thermal andelectromagnetic cloaking by multi-physicalnull medium

Simultaneously manipulating multiple physical fields plays an important role in the increasingly complex integrated systems,aerospace equipment,biochemical productions,etc.For on-chip systems with high integration level,the precise and efficient control of the propagation of electromagnetic waves and heat fluxes simultaneously is particularly important.In this study,we propose a graphical designing method(i.e.,thermal-electromagnetic surface transformation)based on thermal-electromagnetic null medium to simultaneously control the propagation of electromagnetic waves and thermal fields according to the pre-designed paths.A thermal-electromagnetic cloak,which can create a cloaking effect on both electromagnetic waves and thermal fields simultaneously,is designed by thermal-electromagnetic surface transformation and verified by both numerical simulations and experimental measurements.The thermal-electromagnetic surface transformation proposed in this study provides a new methodology for simultaneous controlling on electromagnetic and temperature fields,and may have significant applications in improving thermal-electromagnetic compatibility problem,protecting of thermal-electromagnetic sensitive components,and improving efficiency of energy usage for complex onchip systems.

Verification of a self-developed CFD-based multi-physics coupled code MPC-LBE for LBE-cooled reactor

To perform an integral simulation of a pool-type reactor using CFD code,a multi-physics coupled code MPC-LBE for an LBE-cooled reactor was proposed by integrating a point kinetics model and a fuel pin heat transfer model into self-developed CFD code.For code verification,a code-to-code comparison was employed to validate the CFD code.Furthermore,a typical BT transient benchmark on the LBE-cooled XADS reactor was selected for verification in terms of the integral or system performance.Based on the verification results,it was demonstrated that the MPC-LBE coupled code can perform thermal-hydraulics or safety analyses for analysis for processes involved in LBE-cooled pool-type reactors.

Multi-physics analysis of permanent magnet tubular linear motors under severe volumetric and thermal constraints

Permanent magnet tubular linear motors(TLMs) arranged in multiple rows and multiple columns used for a radiotherapy machine were studied. Due to severe volumetric and thermal constraints, the TLMs were at high risk of overheating. To predict the performance of the TLMs accurately, a multi-physics analysis approach was proposed. Specifically, it considered the coupling effects amongst the electromagnetic and the thermal models of the TLMs, as well as the fluid model of the surrounding air. To reduce computation cost, both the electromagnetic and the thermal models were based on lumped-parameter methods. Only a minimum set of numerical computation(computational fluid dynamics, CFD) was performed to model the complex fluid behavior. With the proposed approach, both steady state and transient state temperature distributions, thermal rating and permissible load can be predicted. The validity of this approach is verified through the experiment.

Development and application of a multi-physics and multi-scale coupling program for lead-cooled fast reactor

In this study,a multi-physics and multi-scale coupling program,Fluent/KMC-sub/NDK,was developed based on the user-defined functions(UDF)of Fluent,in which the KMC-sub-code is a sub-channel thermal-hydraulic code and the NDK code is a neutron diffusion code.The coupling program framework adopts the"master-slave"mode,in which Fluent is the master program while NDK and KMC-sub are coupled internally and compiled into the dynamic link library(DLL)as slave codes.The domain decomposition method was adopted,in which the reactor core was simulated by NDK and KMC-sub,while the rest of the primary loop was simulated using Fluent.A simulation of the reactor shutdown process of M2LFR-1000 was carried out using the coupling program,and the code-to-code verification was performed with ATHLET,demonstrating a good agreement,with absolute deviation was smaller than 0.2%.The results show an obvious thermal stratification phenomenon during the shutdown process,which occurs 10 s after shutdown,and the change in thermal stratification phenomena is also captured by the coupling program.At the same time,the change in the neutron flux density distribution of the reactor was also obtained.

A robust multi-objective and multi-physics optimization of multi-physics behavior of microstructure

A new strategy is presented to solve robust multi-physics multi-objective optimization problem known as improved multi-objective collaborative optimization （IMOCO） and its extension improved multi-objective robust collaborative （IMORCO）. In this work, the proposed IMORCO approach combined the IMOCO method, the worst possible point （WPP） constraint cuts and the Genetic algorithm NSGA-II type as an optimizer in order to solve the robust optimization problem of multi-physics of microstructures with uncertainties. The optimization problem is hierarchically decomposed into two levels： a microstructure level, and a disciplines levels, For validation purposes, two examples were selected： a numerical example, and an engineering example of capacitive micro machined ultrasonic transducers （CMUT） type. The obtained results are compared with those obtained from robust non-distributed and distributed optimization approach, non-distributed multi-objective robust optimization （NDMORO） and multi-objective collaborative robust optimization （McRO）, respectively. Results obtained from the application of the IMOCO approach to an optimization problem of a CMUT cell have reduced the CPU time by 44% ensuring a Pareto front close to the reference non-distributed multi-objective optimization （NDMO） approach （mahalanobis distance, D2M =0.9503 and overall spread, So=0.2309）. In addition, the consideration of robustness in IMORCO approach applied to a CMUT cell of optimization problem under interval uncertainty has reduced the CPU time by 23% keeping a robust Pareto front overlaps with that obtained by the robust NDMORO approach （D2M =10.3869 and So=0.0537）.

Minimal Realization of Linear Graph Models for Multi-physics Systems

An engineering system may consist of several different types of components,belonging to such physical"domains"as mechanical,electrical,fluid,and thermal.It is termed a multi-domain(or multi-physics)system.The present paper concerns the use of linear graphs(LGs)to generate a minimal model for a multi-physics system.A state-space model has to be a minimal realization.Specifically,the number of state variables in the model should be the minimum number that can completely represent the dynamic state of the system.This choice is not straightforward.Initially,state variables are assigned to all the energy-storage elements of the system.However,some of the energy storage elements may not be independent,and then some of the chosen state variables will be redundant.An approach is presented in the paper,with illustrative examples in the mixed fluid-mechanical domains,to illustrate a way to recognize dependent energy storage elements and thereby obtain a minimal state-space model.System analysis in the frequency domain is known to be more convenient than in the time domain,mainly because the relevant operations are algebraic rather than differential.For achieving this objective,the state space model has to be converted into a transfer function.The direct way is to first convert the state-space model into the input-output differential equation,and then substitute the time derivative by the Laplace variable.This approach is shown in the paper.The same result can be obtained through the transfer function linear graph(TF LG)of the system.In a multi-physics system,first the physical domains have to be converted into an equivalent single domain(preferably,the output domain of the system),when using the method of TFLG.This procedure is illustrated as well,in the present paper.

Mechatronic Modeling and Domain Transformation of Multi-physics Systems

The enhanced definition of Mechatronics involves the four underlying characteristics of integrated,unified,unique,and systematic approaches.In this realm,Mechatronics is not limited to electro-mechanical systems,in the multi-physics sense,but involves other physical domains such as fluid and thermal.This paper summarizes the mechatronic approach to modeling.Linear graphs facilitate the development of state-space models of mechatronic systems,through this approach.The use of linear graphs in mechatronic modeling is outlined and an illustrative example of sound system modeling is given.Both time-domain and frequency-domain approaches are presented for the use of linear graphs.A mechatronic model of a multi-physics system may be simplified by converting all the physical domains into an equivalent single-domain system that is entirely in the output domain of the system.This approach of converting(transforming)physical domains is presented.An illustrative example of a pressure-controlled hydraulic actuator system that operates a mechanical load is given.

Multi-physics multi-scale simulation of unique equiaxed-to-columnar-to-equiaxed transition during the whole solidification process of Al-Li alloy laser welding

In this study,a novel multi-physics multi-scale model with the dilute multicomponent phase-field method in three-dimensional(3D)space was developed to investigate the complex microstructure evolu-tion in the molten pool during laser welding of Al-Li alloy.To accurately compute mass data within both two and three-dimensional computational domains,three efficient computing methods,including central processing unit parallel computing,adaptive mesh refinement,and moving-frame algorithm,were uti-lized.Emphasis was placed on the distinctive equiaxed-to-columnar-to-equiaxed transition phenomenon that occurs during the entire solidification process of Al-Li alloy laser welding.Simulation results indi-cated that the growth distance of columnar grains that epitaxially grew from the base metal(BM)de-creased as the nucleation rate increased.As the nucleation rate increased,the morphology of the newly formed grains near the fusion boundary(FB)changed from columnar to equiaxed,and newly formed equiaxed grains changed from having high-order dendrites to no obvious dendrite structure.When the nucleation rate was sufficiently high,non-dendritic equiaxed grains could directly form near the FB,and there was nearly no epitaxial growth from the BM.Additionally,simulation results illustrated the com-petition among multiple grains with varying orientations that grow in 3D space near the FB.Finally,how equiaxed grain bands develop was elucidated.The equiaxed band not only hindered the growth of early columnar grains but also some of its grains could grow epitaxially to form new columnar grains.These predicted results were in good agreement with experimental measurements and observations.

Microwave-induced thermoacoustic elastic imaging:A simulation study

Microwave-induced thermoacoustic imaging(MTI)has the advantages of high resolution,high contrast,non-ionization,and non-invasive.Recently,MTI was used in the¯eld of breast cancer screening.In this paper,based on the¯nite element method(FEM)and COMSOL Multiphysics software,a three-dimensional breast cancer model suitable for exploring the MTI process is proposed to investigate the in°uence of Young's modulus(YM)of breast cancer tissue on MTI.It is found that the process of electromagnetic heating and initial pressure generation of the entire breast tissue is earlier in time than the thermal expansion process.Besides,compared with normal breast tissue,tumor tissue has a greater temperature rise,displacement,and pressure rise.In particular,YM of the tumor is related to the speed of thermal expansion.In particular,the larger the YM of the tumor is,the higher the heating and contraction frequency is,and the greater the maximum pressure is.Di®erent Young's moduli correspond to di®erent thermoacoustic signal spectra.In MTI,this study can be used to judge di®erent degrees of breast cancer based on elastic imaging.In addition,this study is helpful in exploring the possibility of microwave-induced thermoacoustic elastic imaging(MTAE).

Thermal management by manipulating electromagnetic parameters

Electromagnetic absorbing materials may convert electromagnetic energy into heat energy and dissipate it.However,in a high-power electromagnetic radiation environment,the temperature of the absorbing material rises significantly and even burns.It becomes critical to ensure electromagnetic absorption performance while minimizing temperature rise.Here,we systematically study the coupling mechanism between the electromagnetic field and the temperature field when the absorbing material is irradiated by electromagnetic waves.We find out the influence of the constitutive parameters of the absorbing materials(including uniform and non-uniform)on the temperature distribution.Finally,through a smart design,we achieve better absorption and lower temperature simultaneously.The accuracy of the model is affirmed as simulation results aligned with theoretical analysis.This work provides a new avenue to control the temperature distribution of absorbing materials.

Understanding and improving Yangtze River Basin summer precipitation prediction using an optimal multi-Physics ensemble

This study employs the regional Climate-Weather Research and Forecasting model(CWRF)to first investigate the primary physical mechanisms causing biases in simulating summer precipitation over the Yangtze River Basin(YRB),and then enhance its predictive ability through an optimal multi-physics ensemble approach.The CWRF 30-km simulations in China are compared among 28 combinations of varying physics parameterizations during 1980−2015.Long-term average summer biases in YRB precipitation are remotely correlated with those of large-scale circulations.These teleconnections of biases are highly consistent with the observed correlation patterns between interannual variations of precipitation and circulations,despite minor shifts in their primary action centers.Increased YRB precipitation aligns with a southward shifted East Asian westerly jet,an intensified low-level southerly flow south of YRB,and a south-eastward shifted South Asian high,alongside higher moisture availability over YRB.Conversely,decreased YRB precipitation corresponds to an opposite circulation pattern.The CWRF control configuration using the ensemble cumulus parameterization(ECP),compared to other cumulus schemes,best captures the observed YRB precipitation characteristics and associated circulation patterns.Coupling ECP with the Morrison or Morrison-aerosol microphysics and the CCCMA or CAML radiation schemes enhances the overall CWRF skills.Compared to the control CWRF,the ensemble average of these skill-enhanced physics configurations more accurately reproduces YRB summer precipitation’s spatial distributions,interannual anomalies,and associated circulation patterns.The Bayesian Joint Probability calibration to these configurations improves the ensemble’s spatial distributions but compromises its interannual anomalies and teleconnection patterns.Our findings highlight substantial potential for refining the representation of climate system physics to improve YRB precipitation prediction.This is notably achieved by realistically coupling cumulus,microphysics,and radiation processes to accurately capture circulation teleconnections.Further enhancements can be achieved by optimizing the multi-physics ensemble among skill-enhanced configurations.

Coupled Numerical Simulation of Electromagnetic and Flow Fields in a Magnetohydrodynamic Induction Pump

Magnetohydrodynamic(MHD)induction pumps are contactless pumps able to withstand harsh environments.The rate of fluid flow through the pump directly affects the efficiency and stability of the device.To explore the influence of induction pump settings on the related delivery speed,in this study,a numerical model for coupled electromagnetic and flow field effects is introduced and used to simulate liquid metal lithium flow in the induction pump.The effects of current intensity,frequency,coil turns and coil winding size on the velocity of the working fluid are analyzed.It is shown that the first three parameters have a significant impact,while changes in the coil turns have a negligible influence.The maximum increase in working fluid velocity within the pump for the parameter combination investigated in this paper is approximately 618%.As the frequency is increased from 20 to 60 Hz,the maximum increase in the mean flow rate of the working fluid is approximately 241%.These research findings are intended to support the design and optimization of these devices.

Research on heat dissipation characteristics of magnetic fluid bearings under multiple field coupling effects

This paper analyzes the sources of heat losses in magnetic fluid bearings,proposes various cou-pling relationships of physical fields,divides the coupled heat transfer surfaces while ensuring the continuity of heat flux density,and analyzes the overall heat dissipation pathways of the bearings.By changing parameters such as input current,rotor speed,and inlet oil flow rate,the study applies a multi-physics field coupling method to investigate the influence of different parameters on the temper-ature field and heat dissipation patterns of the bearings,which is then validated through experi-ments.This research provides a theoretical basis for the optimal design of magnetic fluid bearing sys-tems.

Integration of multi-physics and machine learning-based surrogate modelling approaches for multi-objective optimization of deformed GDL of PEM fuel cells

The development of artificial intelligence(AI)greatly boosts scientific and engineering innovation.As one of the promising candidates for transiting the carbon intensive economy to zero emission future,proton exchange membrane(PEM)fuel cells has aroused extensive attentions.The gas diffusion layer(GDL)strongly affects the water and heat management during PEM fuel cells operation,therefore multi-variable optimization,including thickness,porosity,conductivity,channel/rib widths and compression ratio,is essential for the improved cell performance.However,traditional experiment-based optimization is time consuming and economically unaffordable.To break down the obstacles to rapidly optimize GDLs,physics-based simulation and machine-learning-based surrogate modelling are integrated to build a sophisticated M 5 model,in which multi-physics and multi-phase flow simulation,machine-learning-based surrogate modelling,multi-variable and multi-objects optimization are included.Two machine learning methodologies,namely response surface methodol-ogy(RSM)and artificial neural network(ANN)are compared.The M 5 model is proved to be effective and efficient for GDL optimization.After optimization,the current density and standard deviation of oxygen dis-tribution at 0.4 V are improved by 20.8%and 74.6%,respectively.Pareto front is obtained to trade off the cell performance and homogeneity of oxygen distribution,e.g.,20.5%higher current density is achieved when sacrificing the standard deviation of oxygen distribution by 26.0%.

Electro-chemo-mechanical design of polymer matrix in composited LiNi_(0.8)Co_(0.1)Mn_(0.1)O_(2) cathode endows solid-state batteries with superior performance

Nickel-rich LiNi_(0.8)Co_(0.1)Mn_(0.1)O_(2)(NCM811) cathode material has been widely concerned due to its high voltage,high specific capacity and excellent rate performance,which is considered as one of the most promising cathode materials for the next generation of high-energy-density solid-state lithium batteries.However,serious electro-chemo-mechanical degradation of Nickel-rich cathode during cycling,especially at a high voltage(over 4.5 V),constrains their large-scale application.Here,using the multiphysical simulation,highly-conductive polymer matrix with spontaneous stress-buffering effect was uncovered theoretically for reinforcing the electrochemical performance of composited NCM81 1 cathode through the visualization of uniform concentration distribution of Li-ion coupled with improved stress field inside NCM811 cathode.Thereupon,polyacrylonitrile(PAN) and soft polyvinylidene fluoride(PVDF) were selected as the polymer matrix to fabricate the composited NCM811 cathode(PVDFPAN@NCM811) for improving the electrochemical performance of the solid-state NMC811|Li full cells,which can maintain high capacity over 146.2 mA h g^(-1)after 200 cycles at a high voltage of 4.5 V.Suggestively,designing a multifunctional polymer matrix with high ionic conductivity and mechanical property can buffer the stress and maintain the integrity of the structure,which can be regarded as the door-opening avenue to realize the high electrochemical performance of Ni-rich cathode for solidstate batteries.

Research on Electromagnetic Loss Characteristics of Submarine Cables

The electromagnetic losses of submarine cables are mainly caused by the metal shielding layer to prevent the water tree effect and the armor layer that strengthens the strength of the submarine cables.While these losses cause the temperature of submarine cable to rise,and temperature variation will in turn change the conductivity of its metal layer material.In this paper,the electric-magnetic-thermal multi-physical field coupling of the electromagnetic loss variation of the submarine cable is realized by establishing a full coupling system containing Fourier’s law and Maxwell-Ampère’s Law for the photoelectric composite submarine cable.The multi-physical field coupling model is solved and analyzed by using the finite elementmethod.Firstly,the loss of each layer of the optoelectronic composite submarine cable is analyzed,and the lossof eachlayer of the submarine cable and themainfactors leading to the loss of the submarine cable are given.Secondly,the influence of environmental temperature,ampacity and armor layer on the electromagnetic loss of submarine cables is studied,and the main operating factors affecting the electromagnetic loss of submarine cables are summarized.The research shows that the influence of ambient temperature can be ignored,and the loss of shielding layer and armor layer increases with the increase of ampacity,but the impact of shielding layer loss is greater.Finally,this paper studies the electromagnetic loss of each metal layer of the submarine cable and the influence of the laying spacing on the electromagnetic loss.The research results show that the two ways of improving the conductivity of the armor layer and reducing the relative permeability of the armor layer can effectively reduce the loss of each metal layer in the cable structure and increase the current carrying capacity when the tensile strength of the armor layer meets the requirements for single-core and threecore photoelectric composite submarine cables laid horizontally.At the same time,increasing the laying spacing will increase the loss,but it can improve the overall current carrying capacity of the cable.The research in this paper provides a theoretical basis for the design of submarine cable carrying capacity,and also provides a reference for the optimization design of submarine cable structures.