Simulation and Analysis of Cascading Faults in Integrated Heat and Electricity Systems Considering Degradation Characteristics

Cascading faults have been identified as the primary cause of multiple power outages in recent years.With the emergence of integrated energy systems(IES),the conventional approach to analyzing power grid cascading faults is no longer appropriate.A cascading fault analysis method considering multi-energy coupling characteristics is of vital importance.In this study,an innovative analysis method for cascading faults in integrated heat and electricity systems(IHES)is proposed.It considers the degradation characteristics of transmission and energy supply com-ponents in the system to address the impact of component aging on cascading faults.Firstly,degradation models for the current carrying capacity of transmission lines,the water carrying capacity and insulation performance of thermal pipelines,as well as the performance of energy supply equipment during aging,are developed.Secondly,a simulation process for cascading faults in the IHES is proposed.It utilizes an overload-dominated development model to predict the propagation path of cascading faults while also considering network islanding,electric-heating rescheduling,and load shedding.The propagation of cascading faults is reflected in the form of fault chains.Finally,the results of cascading faults under different aging levels are analyzed through numerical examples,thereby verifying the effectiveness and rationality of the proposed model and method.

A novel electron-phonon coupling thermoelasticity with Burgers electronic heat transfer

The electron-phonon interaction can reveal the microscopic mechanism of heat transfer in metals.The two-step heat conduction considering electron-phonon interaction has become an effective theoretical model for extreme environments,such as micro-scale and ultrafast processes.In this work,the two-step heat transfer model is further extended by considering the Burgers heat conduction model with the secondorder heat flux rate for electrons.Then,a novel generalized electron-phonon coupling thermoelasticity is proposed with the Burgers electronic heat transfer.Then,the problem of one-dimensional semi-infinite copper strip subject to a thermal shock at one side is studied by the Burgers two-step(BTS)model.The thermoelastic analytical solutions are systematically derived in the Laplace domain,and the numerical Laplace inversion method is adopted to obtain the transient responses.The new model is compared with the parabolic two-step(PTS)model and the hyperbolic two-step(HTS)model.The results show that in ultrafast heating,the BTS model has the same wave front jump as the HTS model.The present model has the faster wave speed,and predicts the bigger disturbed regions than the HTS model.More deeply,all two-step models also have the faster wave speeds than one-step models.This work may benefit the theoretical modeling of ultrafast heating of metals.

Coexisting fast–slow dendritic traveling waves in a 3D-array electric field coupled neuronal network

Coexistence of fast and slow traveling waves without synaptic transmission has been found in hhhippocampal tissues,which is closely related to both normal brain activity and abnormal neural activity such as epileptic discharge. However, the propagation mechanism behind this coexistence phenomenon remains unclear. In this paper, a three-dimensional electric field coupled hippocampal neural network is established to investigate generation of coexisting spontaneous fast and slow traveling waves. This model captures two types of dendritic traveling waves propagating in both transverse and longitude directions: the N-methyl-D-aspartate(NMDA)-dependent wave with a speed of about 0.1 m/s and the Ca-dependent wave with a speed of about 0.009 m/s. These traveling waves are synaptic-independent and could be conducted only by the electric fields generated by neighboring neurons, which are basically consistent with the in vitro data measured experiments. It is also found that the slow Ca wave could trigger generation of fast NMDA waves in the propagation path of slow waves whereas fast NMDA waves cannot affect the propagation of slow Ca waves. These results suggest that dendritic Ca waves could acted as the source of the coexistence fast and slow waves. Furthermore, we also confirm the impact of cellular spacing heterogeneity on the onset of coexisting fast and slow waves. The local region with decreasing distances among neighbor neurons is more liable to promote the onset of spontaneous slow waves which, as sources, excite propagation of fast waves. These modeling studies provide possible biophysical mechanisms underlying the neural dynamics of spontaneous traveling waves in brain tissues.

Analysis of influence of heat exchangerfouling on heat transfer performancebased on thermal fluid coupling

A study on heat transfer performance by thermal fluid coupling simulation for the fouling in a shell-tube heat exchanger used in engineering was presented. The coupling simulation was performed in a fluid and solid domains under three different fouling conditions: fouling inside the tube, fouling outside the tube, and fouling inside the shell. The flow field, temperature, and pressure distributions in the heat exchanger were solved numerically to analyze the heat transfer performance parameters, such as thermal resistance. It is found that the pressure drop of the heat exchanger and the thermal resistance of the tube wall increase by nearly 30% and 20%, respectively, when the relative fouling thickness reaches 10%. The fouling inside the tube has more impact on the heat transfer performance of the heat exchanger, and the fouling inside the shell has less impact.

Research on Operation Optimization of Heating System Based on Electric Storage Coupled Solar Energy and Air Source Heat Pump

For heating systems based on electricity storage coupled with solar energy and an air source heat pump(ECSA),choosing the appropriate combination of heat sources according to local conditions is the key to improving economic efficiency.In this paper,four cities in three climatic regions in China were selected,namely Nanjing in the hot summer and cold winter region,Tianjin in the cold region,Shenyang and Harbin in the severe cold winter region.The levelized cost of heat(LCOH)was used as the economic evaluation index,and the energy consumption and emissions of different pollutants were analyzed.TRNSYS software was used to simulate and analyze the system performance.The Hooke-Jeeves optimization algorithm and GenOpt software were used to optimize the system parameters.The results showed that ECSA systemhad an excellent operation effect in cold region and hot summer and cold winter region.Compared with ECS system,the systemenergy consumption,and the emission of different pollutants of ECSA system can be reduced by a maximum of 1.37 times.In cold region,the initial investment in an air source heat pump is higher due to the lower ambient temperature,resulting in an increase in the LOCH value of ECSA system.After the LOCH value of ECSA system in each region was optimized,the heating cost of the system was reduced,but also resulted in an increase in energy consumption and the emission of different pollutant gases.

Mixed Electric and Magnetic Coupling Design Based on Coupling Matrix Extraction

This paper proposes a design and fine-tuning method for mixed electric and magnetic coupling filters.It derives the quantitative relationship between the coupling coefficients(electric and magnetic coupling,i.e.,EC and MC)and the linear coefficients of frequencydependent coupling for the first time.Different from the parameter extraction technique using the bandpass circuit model,the proposed approach explicitly relatesEC and MC to the coupling matrix model.This paper provides a general theoretic framework for computer-aided design and tuning of a mixed electric and magnetic coupling filter based on coupling matrices.An example of a 7th-order coaxial combline filter design is given in the paper,verifying the practical value of the approach.

A Numerical Model Coupling Electromagnetism and Thermomechanics Application to Induction Heating Modeling

A numerical model coupling the various physical phenomena (electromagnetic, thermal and mechanical) taking place in the induction heating process has been developed. The mathematical model and the numerical methods are presented here, along with some results ( electric, thermal and mechanical fields in the workpiece)

Three Dimensional Coupling between Elastic and Thermal Fields in the Static Analysis of Multilayered Composite Shells

This new work aims to develop a full coupled thermomechanical method including both the temperature profile and displacements as primary unknowns of the model.This generic full coupled 3D exact shell model permits the thermal stress investigation of laminated isotropic,composite and sandwich structures.Cylindrical and spherical panels,cylinders and plates are analyzed in orthogonal mixed curved reference coordinates.The 3D equilibrium relations and the 3D Fourier heat conduction equation for spherical shells are coupled and they trivially can be simplified in those for plates and cylindrical panels.The exponential matrix methodology is used to find the solutions of a full coupled model based on coupled differential relations with respect to the thickness coordinate.The analytical solution is based on theories of simply supported edges and harmonic relations for displacement components and sovra-temperature.The sovra-temperature magnitudes are directly applied at the outer faces through static state hypotheses.As a consequence,the sovra-temperature description is assumed to be an unknown variable of themodel and it is calculated in the sameway as the three displacements.The final systemis based on a set of coupled homogeneous differential relations of second order in the thickness coordinate.This system is reduced in a first order differential relation system by redoubling the number of unknowns.Therefore,the exponential matrix methodology is applied to calculate the solution.The temperature field effects are evaluated in the static investigation of shells and plates in terms of displacement and stress components.After an appropriate preliminary validation,new benchmarks are discussed for several thickness ratios,geometrical data,lamination sequences,materials and sovra-temperature values imposed at the outer faces.Results make evident the accordance between the uncoupled thermo-mechanical model and this new full coupled thermo-mechanical model without the need to separately solve the Fourier heat conduction relation.Both effects connected with the thickness layer and the related embedded materials are included in the conducted thermal stress analysis.

Heat-fluid-solid coupling model for gas-bearing coal seam and numerical modeling on gas drainage promotion by heat injection

Improving the absorbed gas to active desorption and seepage and delaying gas drainage attenuation are considered as key methods for increasing drainage efficiency and gas output.According to the solid mechanics theory,the nonlinear Darcy seepage theory and thermodynamics,the heat-fluid-solid coupling model for gassy coal has been improved.The numerical model was founded from the improved multi-field coupling model by COMSOL Multiphysics and gas drainage by borehole down the coal seam enhanced by heat injection was modelled.The results show that the heatfluid-solid model with adsorption effects for gassy coal was well simulated by the improved multi-field model.The mechanism of coal seam gas desorption seepage under the combined action of temperature,stress and adsorption can be well described.Gas desorption and seepage can be enhanced by heat injection into coal seams.The gas drainage rate was directly proportional to the temperature of injected heat in the scope of 30-150 ℃ and increasing in the whole modelleddrainage process (0-1000 d).The increased level was maximum in the initial drainage time and decreasing gradually along with drainage time.The increasing ratio of drainage rate was maximum when the temperature raised from 30 to 60 ℃.Although the drainage rate would increase along with increasing temperature,when exceeding 60 ℃,the increasing ratio of drainage rate with rising temperature would decrease.Gas drainage promotion was more effective in coal seams with lower permeability than with higher permeability.The coal seam temperature in a 5 m distance surrounding the heat injection borehole would rise to around 60 ℃ in 3 months.That was much less than the time of gas drainage in the coal mines in sites with low permeability coal seams.Therefore,it is valuable and feasible to inject heat into coal seams to promote gas drainage,and this has strong feasibility for coal seams with low permeability which are widespread in China.

Optimal Operation of Integrated Heat and Electricity Systems:A Tightening McCormick Approach

Combined heat and electricity operation with variable mass flow rates promotes flexibility,economy,and sustainability through synergies between electric power systems(EPSs)and district heating systems(DHSs).Such combined operation presents a highly nonlinear and nonconvex optimization problem,mainly due to the bilinear terms in the heat flow model—that is,the product of the mass flow rate and the nodal temperature.Existing methods,such as nonlinear optimization,generalized Benders decomposition,and convex relaxation,still present challenges in achieving a satisfactory performance in terms of solution quality and computational efficiency.To resolve this problem,we herein first reformulate the district heating network model through an equivalent transformation and variable substitution.The reformulated model has only one set of nonconvex constraints with reduced bilinear terms,and the remaining constraints are linear.Such a reformulation not only ensures optimality,but also accelerates the solving process.To relax the remaining bilinear constraints,we then apply McCormick envelopes and obtain an objective lower bound of the reformulated model.To improve the quality of the McCormick relaxation,we employ a piecewise McCormick technique that partitions the domain of one of the variables of the bilinear terms into several disjoint regions in order to derive strengthened lower and upper bounds of the partitioned variables.We propose a heuristic tightening method to further constrict the strengthened bounds derived from the piecewise McCormick technique and recover a nearby feasible solution.Case studies show that,compared with the interior point method and the method implemented in a global bilinear solver,the proposed tightening McCormick method quickly solves the heat–electricity operation problem with an acceptable feasibility check and optimality.

Spin current and its heat effect in a multichannel quantum wire with Rashba spin-orbit coupling

Using the perturbation method, we theoretically study the spin current and its heat effect in a multichannel quantum wire with Rashba spin-orbit coupling. The heat generated by the spin current is calculated. With the increase of the width of the quantum wire, the spin current and the heat generated both exhibit period oscillations with equal amplitudes. When the quantum-channel number is doubled, the oscillation periods of the spin current and of the heat generated both decrease by a factor of 2. For the spin current js,xy, the amplitude increases with the decrease of the quantum channel; while the amplitude of the spin current js,yx remains the same. Therefore we conclude that the effect of the quantum-channel number on the spin current js,xy is greater than that on the spin current js,yx. The strength of the Rashba spin-orbit coupling is tunable by the gate voltage, and the gate voltage can be varied experimentally, which implies a new method of detecting the. spin current. In addition, we can control the amplitude and the oscillation period of the spin current by controlling the number of the quantum channels. All these characteristics of the spin current will be very important for detecting and controlling the spin current, and especially for designing new spintronic devices in the future.

Multi-dimensional Simulation of Phase Change by a 0D-2D Model Coupling via Stefan Condition

Considering phase changes associated with a high-temperature molten material cooled down from the outside,this work presents an improvement of the modelling and the numerical simulation of such processes for an application pertaining to the safety of light water nuclear reactors.Postulating a core meltdown accident,the behaviour of the core melt(aka corium)into a steel vessel is of tremendous importance when evaluating the vessel integrity.Evaluating correctly the heat fluxes requires the numerical simulation of the interaction between the liquid material and its solid counterpart which forms during the solidification process,but also may melt back.To simulate this configuration,encoun-tered in various industrial applications,one considers a bi-phase model constituted by a liquid phase in contact and interaction with its solid phase.The liquid phase may solidify in presence of low energetic source,while the solid phase may melt due to an intense heat flux from the high-temperature liquid.In the frame of the in-house legacy code,several simplifying assumptions(0D multi-layer discretization,instantaneous heat transfer via a quadratic temperature profile in solids)are made for the modelling of such phase changes.In the present work,these shortcomings are illustrated and further overcome by solving a 2D heat conduction model in the solid by a mixed Raviart-Thomas finite element method coupled to the liquid phase due to heat and mass exchanges through Stefan condition.The liquid phase is modeled with a 0D multi-layer approach.The 0D-liquid and 2D-solid mod-els are coupled by a Stefan like phase change interface model.Several sanity checks are performed to assess the validity of the approach on 1D and 2D academical configurations for which exact or reference solutions are available.Then more advanced situations(genu-ine multi-dimensional phase changes and an"industrial-like scenario")are simulated to verify the appropriate behavior of the obtained coupled simulation scheme.

City-Bus-Route Demand-based Efficient Coupling Driving Control for Parallel Plug-in Hybrid Electric Bus

Recently,plug?in hybrid electric bus has been one of the energy?e cient solutions for urban transportation. However,the current vehicle e ciency is far from optimum,because the unpredicted external driving conditions are di cult to be obtained in advance. How to further explore its fuel?saving potential under the complicated city bus driving cycles through an e cient control strategy is still a hot research issue in both academic and engineering area. To realize an e cient coupling driving operation of the hybrid powertrain,a novel coupling driving control strategy for plug?in hybrid electric bus is presented. Combined with the typical feature of a city?bus?route,the fuzzy logic inference is employed to quantify the driving intention,and then to determine the coupling driving mode and the gear?shifting strategy. Considering the response deviation problem in the execution layer,an adaptive robust controller for electric machine is designed to respond to the transient torque demand,and instantaneously compensate the response delay and the engine torque fluctuation. The simulations and hard?in?loop tests with the actual data of two typical driving conditions from the real?world city?bus?route are carried out,and the results demonstrate that the pro?posed method could guarantee the hybrid powertrain to track the actual torque demand with 10.4% fuel economy improvement. The optimal fuel economy can be obtained through the optimal combination of working modes. The fuel economy of plug?in hybrid electric bus can be significantly improved by the proposed control scheme without loss of drivability.

Review of ionospheric irregularities and ionospheric electrodynamic coupling in the middle latitude region

This paper briefly reviews ionospheric irregularities that occur in the E and F regions at mid-latitudes. Sporadic E(ES) is a common ionospheric irregularity phenomenon that is first noticed in the E layer. ES mainly appears during daytime in summer hemispheres, and is formed primarily from neutral wind shear in the mesosphere and lower thermosphere(MLT) region. Field-aligned irregularity(FAI) in the E region is also observed by Very High Frequency(VHF) radar in mid-latitude regions. FAI frequently occurs after sunset in summer hemispheres, and spectrum features of E region FAI echoes suggest that type-2 irregularity is dominant in the nighttime ionosphere. A close relationship between ES and E region FAI implies that ES may be a possible source of E region FAI in the nighttime ionosphere. Strong neutral wind shear, steep ES plasma density gradient, and a polarized electric field are the significant factors affecting the formation of E region FAI. At mid-latitudes, joint observational experiments including ionosonde, VHF radar, Global Positioning System(GPS) stations, and all-sky optical images have revealed strong connections across different scales of ionospheric irregularities in the nighttime F region, such as spread F(SF), medium-scale traveling ionospheric disturbances(MSTID), and F region FAI.Observations suggest that different scales of ionospheric irregularities are generally attributed to the Perkins instability and subsequently excited gradient drift instability. Nighttime MSTID can further evolve into small-scale structures through a nonlinear cascade process when a steep plasma density gradient exists at the bottom of the F region. In addition, the effect of ionospheric electrodynamic coupling processes, including ionospheric E-F coupling and inter-hemispheric coupling on the generation of ionospheric irregularities, becomes more prominent due to the significant dip angle and equipotentiality of magnetic field lines in the mid-latitude ionosphere. Polarized electric fields can map to different ionospheric regions and excite plasma instabilities which form ionospheric irregularities. Nevertheless,the mapping efficiency of a polarized electric field depends on the ionospheric background and spatial scale of the field.

Study on Electromagnetic-Fluid-Temperature Multiphysics Field Coupling Model for Drum of Mine Cable Winding Truck

Aiming at solving problems of low efficiency,low cable capacity in current 300m open-pit mine cable winding truck,a 900 m cable winding plan was proposed.In this paper,the mechanism of the thermal effect of the cable was described,and a two-dimensional axisymmetric electromagnetic-fluid-temperature multiphysics coupling model of the cable reel was established regarding the 900m cable reel as independent system.Considering the structure of the drum,the number of cable winding layers,the factors of heat conduction,heat radiation and convective heat transfer in the actual working process,the steady state analysis of the multi-physical field coupling was carried out.The sum of the losses of each part of the cable was obtained through the calculation of electromagnetic field,which was used as a heat source to calculate and analyze the temperature distribution of different layers of cable winding,as well as the temperature distribution and heat dissipation characteristics of different structures of the drum.The results show that three layers of cable winding is the best design.The lowest temperature of closed cylindrical drum is 70℃after heat dissipation,which has obvious advantages compared with the lowest temperature of 85℃after heat dissipation of squirrel-cage cylindrical drum.The results provide a reliable theoretical basis for the research and development of a new type of mine cable winding truck with 900 m cable capacity.

Impact of Heat Treatments and Hole Density (p) on the Structural, Electrical, and Superconducting Properties of LnSrBaCu3O6+z (Ln = Eu, Sm, Nd) Compounds

In this study, we thoroughly examined the impact of heat treatments and hole count (p) on the properties of LnSrBaCu3O6+z (Ln = Eu, Sm, Nd) compounds. We focused on preparation, X-ray diffraction with Rietveld refinement, AC susceptibility, DC resistivity measurements, and heat treatment effects. Two heat treatment types were applied: oxygen annealing [O] and argon annealing followed by oxygen annealing [AO]. As the rare earth Ln’s ionic radius increased, certain parameters notably changed. Specifically, c parameter, surface area S, and volume V increased, while critical temperature Tc and holes (p) in the CuO2 plane decreased. The evolution of these parameters with rare earth Ln’s ionic radius in [AO] heat treatment is linear. Regardless of the treatment, the structure is orthorhombic for Ln = Eu, tetragonal for Ln = Nd, orthorhombic for Ln = Sm [AO], and pseudo-tetragonal for Sm [O]. The highest critical temperature is reached with Ln = Eu (Tc [AO] = 87.1 K). Notably, for each sample, Tc [AO] surpasses Tc [O]. Observed data stems from factors including rare earth ionic size, improved cationic and oxygen chain order, holes count p in Cu(2)O2 planes, and in-phase purity of [AO] samples. Our research strives to clearly demonstrate that the density of holes (p) within the copper plane stands as a determinant impacting the structural, electrical, and superconducting properties of these samples. Meanwhile, the other aforementioned parameters contribute to shaping this density (p).

Coupled heat-fluid-solid numerical study on heat extraction potential of hot dry rocks based on discrete fracture network model

Fracture networks within hot dry rock(HDR)geothermal reservoirs are complex,and heat extraction via water injection is thus a coupled process of heat-fluid-solid multifield.In this paper,utilizing the theory of normally distributed random functions,we develop a corresponding pre-processing subprogram to establish a discrete network model of complex fracture distribution in HDR reservoirs;then construct a heat-fluid-solid finite element model for heat extraction via water injection and compare the numerical solution with the analytical solution of the one-dimensional non-isothermal consolidation problem for verification.The numerical simulation results show that the main factors affecting the heat extraction efficiency of HDR reservoirs include fracture width,fracture density,fracture permeability,and matrix permeability.When a HDR reservoir is injected with water for heat extraction,there is a certain threshold value of these influential parameters,beyond which the outlet temperature drops significantly,resulting in an obvious thermal breakthrough.When injecting water for heat extraction,the values of these parameters should be controlled and kept at a reasonable level,otherwise,the HDR reservoir may enter a thermal breakthrough stage in advance,which is not conducive for long-period heat extraction.Influenced by the random distribution of complex fractures,the leading edge of the cold front may present an irregular distribution.During the process of heat extraction,the stress gradually changes from a compressional state to a tensile state,which induces further damage to the HDR reservoir.

A new type of conserved quantity of Mei symmetry for the motion of mechanico electrical coupling dynamical systems

We obtain a new type of conserved quantity of Mei symmetry for the motion of mechanico--electrical coupling dynamical systems under the infinitesimal transformations. A criterion of Mei symmetry for the mechanico-electrical coupling dynamical systems is given. Simultaneously, the condition of existence of the new conserved quantity of Mei symmetry for mechanico-electrical coupling dynamical systems is obtained. Finally, an example is given to illustrate the application of the results.

Effect of driving frequency on electron heating in capacitively coupled RF argon glow discharges at low pressure

A one-dimensional（1D） fluid model on capacitively coupled radio frequency（RF） argon glow discharge between parallel-plates electrodes at low pressure is established to test the effect of the driving frequency on electron heating. The model is solved numerically by a finite difference method. The numerical results show that the discharge process may be divided into three stages： the growing rapidly stage, the growing slowly stage, and the steady stage. In the steady stage,the maximal electron density increases as the driving frequency increases. The results show that the discharge region has three parts： the powered electrode sheath region, the bulk plasma region and the grounded electrode sheath region. In the growing rapidly stage（at 18 μs）, the results of the cycle-averaged electric field, electron temperature, electron density, and electric potentials for the driving frequencies of 3.39, 6.78, 13.56, and 27.12 MHz are compared, respectively. Furthermore,the results of cycle-averaged electron pressure cooling, electron ohmic heating, electron heating, and electron energy loss for the driving frequencies of 3.39, 6.78, 13.56, and 27.12 MHz are discussed, respectively. It is also found that the effect of the cycle-averaged electron pressure cooling on the electrons is to ＂cool＂ the electrons; the effect of the electron ohmic heating on the electrons is always to ＂heat＂ the electrons; the effect of the cycle-averaged electron ohmic heating on the electrons is stronger than the effect of the cycle-averaged electron pressure cooling on the electrons in the discharge region except in the regions near the electrodes. Therefore, the effect of the cycle-averaged electron heating on the electrons is to ＂heat＂ the electrons in the discharge region except in the regions near the electrodes. However, in the regions near the electrodes, the effect of the cycle-averaged electron heating on the electron is to ＂cool＂ the electrons. Finally, the space distributions of the electron pressure cooling the electron ohmic heating and the electron heating at 1/4 T, 2/4 T, 3/4 T, and 4/4 T in one RF-cycle are presented and compared.

Understanding the combustion behavior of electric bicycle batteries and unveiling its relationship with fire extinguishing

In this study,a detailed analysis of the combustion behaviors of the lithium iron phosphate(LFP)and lithium manganese oxide(LMO)batteries used in electric bicycles was conducted.This research included quantitative measurements of the combustion duration,flame height,combustion temperature,heat release rate,and total heat release.The results indicated that LMO batteries exhibited higher combustion temperatures of 600–700°C,flame heights of 70–75 cm,a significantly higher heat release rate of40.1 k W(12 Ah),and a total heat release of 1.04 MJ(12 Ah)compared to LFP batteries with the same capacity.Based on these experimental results,a normalized total heat release(NORTHR)parameter was proposed,demonstrating good universality for batteries with different capacities.Utilizing this parameter,quantitative calculations and optimization of the extinguishing agent dosage were conducted for fires involving these two types of batteries,and the method was validated by extinguishing fires for these two types of battery packs with water-based extinguishing fluids.