Adiabatic Boundary Conditions Research Articles

Thermodynamic characteristics of counter flow and cross flow plate fin heat exchanger based on distributed parameter model

The distributed parameter models of counter flow and cross flow plate fin heat exchanger (PFHE) are constructed to investigate their thermodynamic characteristics. In practical application, the cold end of PFHE may be directly connected to the cold equipment. Therefore, the effect of longitudinal heat conduction (LHC) under different thermal boundary conditions set on two ends of the separating plates (SPs) is investigated based on the model of counter flow PFHE. It is found that LHC effect transmits the effect of boundary condition on two ends of the SPs into the inner domain of PFHE. The LHC effect on the heat transfer performance of PFHE is positive when setting appropriate boundary condition on two ends of the SPs. When considering LHC, compared with the adiabatic boundary conditions on two ends of the SPs, the effectiveness of constant wall temperature boundary conditions increases by 31.7 ∼ 35.8 %, and the effectiveness of constant heat flux boundary conditions increases by 22.8 %. The model of cross flow PFHE is used to investigate the effect of configuration and operation parameters on the heat transfer amount and the pump power. The cross flow PFHE is optimized to arrive at the maximization of heat transfer amount per unit pump power (HTPUP). In the three-stream cross flow PFHE studied, the heat transfer amount almost doesn’t change when the xz projected area is maintained and one projected length (z side length) is changing. The pump power is the lowest and the HTPUP is the largest when z side length of PFHE (z side is the flow direction of the fluid with large mass flowrate) is 0.15 ∼ 0.2 m.

Effects of Uniform and Non-Uniform Temperature Gradients on Marangoni Convection in a Composite Layer with Variable Heat Sources

An investigation is carried out to determine the effect of uniform and non-uniform temperature profiles on single component Darcy-Benard Marangoni (DBM) convection in a composite layer system consisting of an incompressible fluid saturated porous layer on top of which is a layer of the same fluid with variable heat sources in both layers. The upper surface of the fluid layer is free with surface tension effects depending on temperature and the lower surface of the porous layer is rigid. The eigen value, thermal Marangoni number (TMN) is solved exactly for linear, parabolic and inverted parabolic temperature profiles for the adiabatic thermal boundary conditions at the horizontal boundaries of the composite layer. The influence of various dimensionless parameters on the eigen value against depth ratio is discussed in detail.

Turbulent kinetic energy evolution in turbulent boundary layers during head-on interaction of premixed flames with inert walls for different thermal boundary conditions

The statistical behaviour of turbulent kinetic energy, and its transport in turbulent boundary layers during premixed flame-wall interaction for isothermal and adiabatic chemically inert walls have been analysed using a Direct Numerical Simulation (DNS) database. The terms due to the mean velocity gradient (commonly known as the production term) and viscous dissipation remain leading order contributors to the turbulent kinetic energy transport when the flame is away from the wall, similar to that in the corresponding non-reacting turbulent boundary layer, but in addition to these terms, the pressure dilatation and pressure transport contributions play leading order roles when the flame interacts with the wall. It has been found that the gradient hypothesis-based Reynolds stress closures in the context of the k−ε model may yield an inaccurate prediction of the mean velocity gradient term due to the counter-gradient behaviour of Reynolds stresses in the premixed flame cases considered in this analysis. The thermal boundary condition at the wall does not have a major influence apart from a higher wall friction velocity for adiabatic boundary condition than in the case of isothermal boundary condition. The assumption of local equilibrium of production and dissipation of turbulent kinetic energy has been found to be rendered invalid when the flame is either near to the wall or interacting with the wall. This invalidates the conventional log-law for mean streamwise velocity variation, and the wall functions derived based on production-dissipation balance equipped by the usual gradient hypothesis in the context of the standard k−ε model for premixed flame-wall interaction in turbulent boundary layers. Therefore, improved wall functions are necessary for Reynolds averaged Navier-Stokes simulations of premixed flame-wall interaction in turbulent boundary layers.

An Advanced Model for Long-Term Localized Corrosion of Copper

The Nuclear Waste Management Organization (NWMO) is responsible for the implementation of the Adaptive Phased Management (APM) plan, the federally-approved plan for the safe long-term management of Canada’s used nuclear fuel.1 Under the APM plan, used nuclear fuel will ultimately be placed within a deep geological repository (DGR) in a suitable host rock formation. Part of evaluating the long-term performance and safety of the repository system is understanding the behavior of the copper-coated container with respect to localized corrosion. Even though localized corrosion is not expected in the DGR design system, there is still a need to develop a mathematical model for localized corrosion of copper.A 2-D axisymmetric time-dependent model for the localized corrosion of copper was developed using the concept of an Evans droplet.2 The model accounts for three reversible electrochemical reactions on the metal surface (copper dissolution/plating, cuprous ion oxidation (Cu+) /cupric ion (Cu2+) reduction, and hydrogen (H2) evolution/oxidation) and three irreversible electrochemical reactions (cupric ion reduction, oxygen (O2) reduction, and hydroperoxyl ion (HO2 -) reduction). Dissolved cuprous ion reacts in a second chemical step to form the cuprous chloride (CuCl) salt film treated as a diffusion barrier, or it reacts to form the cuprous oxide (Cu2O) film that partially blocked the active area on the metal surface. The model also accounts for fifteen homogeneous reactions involving copper chloride complex, copper hydroxide complex, copper carbonate complex, and superoxide species. The associated equilibrium constants were obtained from PHREEQC thermodynamic software.3 The mathematical model was developed using the finite-element method (COMSOL Multiphysics). The model includes coupled, nonlinear, diffusion equations for ionic species, which include the contribution of migration, local electroneutrality, and homogeneous reactions. The anodic and cathodic regions are not predefined but, instead, are determined by values of local concentration and potential from the simulation results. The growth of nm-scale films was calculated implicitly without using finite-element meshing and surface coverage was expressed in terms of the thickness in units of monolayers. The model also accounted for the effect of porous CuCl film thickness on the surface oxygen concentration and potential applied on electrochemical reactions. The influence of temperature was included on model parameters such as equilibrium rate constant, diffusion coefficient, Henry’s law constant, solubility product constant and kinetic parameters associated with electrochemical reactions. Transient droplet temperature was described by interpolating results from a one-million-year simulation for the temperature change on the used fuel container surface.4 A total of 28 dependent spatial-temporal variables including species’ concentrations, potentials and local corrosion rate were solved in this model.In the current version of model, two extreme conditions were chosen such that droplet boundary has a faster or slower rate of oxygen decay for an initially air-saturated droplet. The model shows a 10-year time-dependent radial distribution of anodic and cathodic current density, surface coverage of CuCl film and Cu2O film, and localized corrosion rates and depths. It also shows a distribution of pH, potential, and concentrations of dissolved gaseous and ionic species through the entire droplet. Temperature and oxygen concentration is shown to have a strong contribution to the simulation results. Simulations show that the depth of corrosion is almost uniform over the elapsed time simulated and the maximum deviation from the average corrosion depth was less than 0.0004% for slower oxygen decay and less than 0.002% for faster oxygen decay in 10 years.References NWMO, “Choosing a Way Forward. The Future Management of Canada’s Used Nuclear Fuel. Final Study,” Nuclear Waste Management Organization, Toronto, Ontario, 2005.R. Evans, The Corrosion of Metals, E. Arnold & Company, London, 1926.PHREEQC Version 3, United States Geological Survey, 202, https://www.usgs.gov/software/phreeqc-version-3Guo, “Thermal response of a Canadian conceptual deep geological repository in crystalline rock and a method to correct the influence of the near-field adiabatic boundary condition”, Engineering Geology, 218 (2017), p50-p62. AcknowledgementThis work was supported by the Nuclear Waste Management Organization, Canada, under project 2000904.

Prandtl number effects on the hydrodynamic stability of compressible boundary layers: flow–thermodynamics interactions

Hydrodynamic stability of compressible boundary layers is strongly influenced by the Mach number ( $M$ ), Prandtl number ( $Pr$ ) and thermal wall boundary condition. These effects manifest on the flow stability via the flow–thermodynamics interactions. Comprehensive understanding of stability flow physics is of fundamental interest and important for developing predictive tools and closure models for integrated transition-to-turbulence computations. The flow–thermodynamics interactions are examined using linear analysis and direct numerical simulations in the following parameter regime: $0.5 \leq M \leq 8$ ; and $0.5 \leq Pr \leq 1.3$ . For the adiabatic wall boundary condition, increasing Prandtl number has a destabilizing effect. In this work, we characterize the behaviour of production, pressure–strain correlation and pressure dilatation as functions of the Mach and Prandtl numbers. First and second instability modes exhibit similar stability trends but the underlying flow physics is shown to be diametrically opposite. The Prandtl number influence on instability is explicated in terms of the base flow profile with respect to the different perturbation mode shapes.

Water boiling time in a spent fuel pool with adiabatic and non-adiabatic boundary conditions

The Spent Fuel Pool (SFP) of Angra II, from Brazil, has received standard spent fuel (SF) assemblies of Uranium dioxide (UO2) discharged from Pressurized Water Reactors (PWR) of the Nuclear Power Plants (NPP) of Angra since the beginning of its operation. However, in case of using Mixed Oxide (MOX) or Thorium-based fuels, it would require further thermal studies of wet storage. It includes the determination of the water boiling time (Tb) of the SFP in case of breakdown of its external cooling system (ECS). This work presents studies of Tb of a simulated SFP storing mixed SF discharged from PWRs. The types of mixed SF studied include MOX plus UO2, oxide of thorium/uranium (U-Th)O2 plus UO2, and oxide of Thorium/transuranic (TRU-Th)O2 plus UO2. The simulations were implemented in CFX Ansys considering the top of the SFP as either adiabatic or non-adiabatic wall. Tb is considerably higher when the non-adiabatic boundary condition is used.

Entropy stability for the compressible Navier-Stokes equations with strong imposition of the no-slip boundary condition

We consider the compressible Navier-Stokes equations subject to no-slip adiabatic wall boundary conditions. The main goal is to investigate stability properties of schemes imposing the no-slip condition strongly (injection) and the temperature condition weakly by a simultaneous approximation term. To this end, we propose a low-order summation-by-parts scheme. By verifying the complete linearisation procedure, we prove linear stability for the scheme. In addition, and assuming that the interior scheme is entropy stable, we also prove entropy stability for the full scheme including the boundary treatment. Furthermore, we propose a linearly stable 3rd-order scheme with the same imposition of the wall conditions. However, the 3rd-order scheme is not provably non-linearly stable. A number of simulations show that the boundary procedure is robust for both schemes.

Modeling and analysis of heat transfer and fluid flow mechanisms in nanofluid filled enclosures irradiated from below

Radiation-driven transport mechanisms are ubiquitous in many natural flows and industrial processes. To mimic and to better understand these processes, recently, radiatively heated nanofluid filled enclosures have been extensively researched. The present work is essentially a determining step in quantifying and understanding the transport mechanisms involved in such enclosures. In particular, a two-dimensional square nanofluid filled enclosure irradiated from the bottom has been investigated in laminar flow situation. Effects of nanofluid optical depth, inclination angle of the enclosure, incident flux, and boundary conditions (adiabatic and isothermal) have been investigated. Moreover, the temperature and flow fields have been carefully analyzed in the situation ranging from “volumetric” to “mixed” to “surface” absorption modes. Under adiabatic boundary conditions, steady state is unconditionally achieved irrespective of the incident flux magnitude (varied between 5 W m-2 - 50 W m-2), enclosure inclination angle (varied between 0° to 60°) and mode of absorption (surface, mixed or volumetric). However, in case of isothermal boundaries; onset of natural convection and its transition into transient regime is significantly impacted by the mode of absorption and the enclosure inclination angle.

Forced Response of Thermoelectric Materials and Devices

The theory of the forced response of electric circuits applied to the study of the thermoelectricity here described allows the characterization of thermoelectric devices and materials determining the resistance of the thermal contacts, and the thermoelectric resistance. The successive computer analysis yields values for the Seebeck coefficient, electrical resistance, thermal conductance, and the figure of merit. The forced response of the thermoelectric materials and devices satisfies the Luttingers thermal transport coefficients theory, and the first-order electric circuitss behavior. Also, permits to find the thermoelectric time constant, and predicting the thermoelectric angular frequency which is necessary to determinate the complex impedance graphically through the Nyquist plot, due to that the thermoelectric time constant is inversely proportional to the thermoelectric angular frequency, as well as it makes accessible the prediction of the impedance spectroscopy measurements beyond the restrictive case of adiabatic boundary conditions regularly dicult to achieve experimentally, and therefore the characterization in situ. Like Harmans method, these parameters can be measured simultaneously on the same device or sample, and no requires neither reference nor standard material for comparison.

Turbulent melting characteristics of hybrid nano-enhanced molten salt phase change material in rectangular enclosure

The turbulent melting characteristics of hybrid nano-enhanced phase change material (HNEPCM) in a three-dimensional rectangular enclosure has been numerically investigated. The lower and upper walls are partially heated and cooled using constant surface temperature boundary condition and the other surfaces are specified with adiabatic boundary condition. The rectangular enclosure is filled with molten salt phase change material (PCM) of potassium nitrate (KNO3) with melting temperature of 334 °C. The transient melting rates of PCM are investigated by dispersing hybrid mono-dispersed spherical nanoparticles of Aluminum oxide (Al2O3) combined with Silicon dioxide (SiO2), Titanium dioxide (TiO2) and Multi-walled carbon nanotubes (MWCNT). The Rayleigh number (Ra) and particle volume fraction (ϕ) are varied in the range of 106≤Ra≤1010 and 2%≤ϕ≤5%. The turbulent phase change process of PCM is modeled using enthalpy–porosity approach with Lam–Bremhorst k−ϵ turbulence model. The three-dimensional governing equations are discretized using Finite difference technique and computations are carried out using Fortran 90 in-house code. The stronger thermal buoyancy force accelerates the melting rate and the melting pattern indicates the accumulation of large quantity of molten metal in the upper regions of the rectangular enclosure. The average melting velocity and turbulent kinetic energy of pure KNO3 PCM is increased by 51.17% and 77.84% by the addition of Al2O3−MWCNT hybrid nanoparticles. The impact of hybrid nanoparticles on the melting rate and thermal performance of PCM is significant in higher Rayleigh number flows than transitional flows. The present results are validated and are in close agreement with the benchmark experimental results.

Enstrophy evolution during head-on wall interaction of premixed flames within turbulent boundary layers

The statistical behaviors of mean enstrophy and its evolution during head-on interaction of premixed flames propagating toward a chemically inert flat wall across the turbulent boundary layer have been analyzed using direct numerical simulations for a friction velocity-based Reynolds number of Reτ=110. The enstrophy dynamics have been analyzed for both isothermal and adiabatic thermal wall boundary conditions. The contributions of vortex-stretching and viscous dissipation are found to be leading order source and sink, respectively, to the mean enstrophy transport in both non-reacting and reacting flows irrespective of the wall boundary condition. However, the contributions due to dilatation rate and baroclinic torque play important roles in addition to the leading order contributions of the vortex-stretching and viscous dissipation terms in the enstrophy transport in turbulent premixed flames. The thermal boundary condition has been demonstrated to affect the near-wall behavior of the enstrophy transport contribution due to dilatation rate, which also affects the near-wall distribution of the enstrophy. The magnitudes of the leading order contributors to the enstrophy transport decrease with the progress of head-on interaction for both wall boundary conditions. Moreover, the overall sink contributions to the enstrophy transport dominate over the source contributions, giving rise to a drop in the mean enstrophy with the progress of head-on interaction. The enstrophy distribution changes significantly during flame-wall interaction, which gives rise to a modification of the relative proportion of the coherent structures in the reacting flow turbulent boundary layer compared to the corresponding non-reacting flow features.

Simultaneous identification of the heat capacity and the anisotropic thermal conductivities of a Li-ion pouch cell by a non-destructive analytical approach

The specifications and aging of Li-ion batteries highly depends on temperature. Battery thermal models are therefore useful to improve the thermal management of battery packs. This study presents a method to determine the heat capacity and the anisotropic thermal conductivities of a Li-ion pouch cell. A single experimental setup is devised and involves local heating of the battery with a resistance heater. The thermal properties of the cell are then identified by calibrating an analytical thermal model with the experimental temperature data. A three-dimensional analytical model is developed for this purpose. It allows for simultaneous characterization of the cell anisotropy. The thermal properties are identified at 50% state of charge between 0 °C and 40 °C. Models with adiabatic and convection boundary conditions are compared and showed similar results. Overall, the values and tendencies agree with the literature. This non-destructive method can be directly applied to any pouch and prismatic cells, and adapted to cylindrical cell.

Numerical analysis of rapid solidification of NiTi alloy: Influence of boundary conditions

The objective of the present work is to analyze influence of temperature boundary conditions on the solidification process of NiTi alloy. Alloy is taken as an incompressible fluid where the heat is transferred by conduction and convection, including the thermal phase change phenomenon. The last one is modelled by the improvement procedure, so-called enthalpy-porosity formulation, where the liquid-solid mushy zone is treated as a porous zone with porosity equal to the liquid fraction. The numerical model is based on the finite volume method in body fitted coordinates with a SIMPLER scheme to couple the pressure and velocity fields. Simulation presents solidification for two cooling cases on the cylindrical part, where in the first case the adiabatic boundary condition is considered and in the second case the convective cooling over the wall is present. The results are presented for the velocity and temperature field as well as for the NiTi mass fraction during the solidification process. Results shows that velocity and temperature field is strongly affected by the different cooling condition on the cylindrical wall and therefore the solidification process of the alloy. Analysis shows that cooling on the cylindrical part is one of the major parameters for alloy solidification and therefore should not be neglected.

Super-suppression of long phonon mean-free-paths in nano-engineered Si due to heat current anticorrelations

The ability to minimize the thermal conductivity of dielectrics with minimal structural intervention that could affect electrical properties is an important capability for engineering thermoelectric efficiency in low-cost materials such as Si. We recently reported the discovery of special arrangements for nanoscale pores in Si that produce a particularly large reduction in thermal conductivity accompanied by strongly anticorrelated heat current fluctuations [1] – a phenomenon that is missed by the diffuse adiabatic boundary conditions conventionally used in Boltzmann transport models. This manuscript presents the results of molecular dynamics simulations and a Monte Carlo ray tracing model that teases apart this phenomenon to reveal that special pore layouts elastically backscatter long-wavelength heat-carrying phonons. This means that heat carriage by a phonon before scattering is undone by the scattered phonon, resulting in an effective mean-free-path that is significantly shorter than the geometric line-of-sight due to the pores. This effect is particularly noticeable for the long-wavelength, long mean-free-path phonons whose transport is impeded drastically more than is expected purely from the usual considerations of scattering defined by the distance between defects. This “super-suppression” of the mean-free-path below the characteristic length scale of the nanostructuring offers a route for minimizing thermal conductivity with minimal structural impact, while the stronger impact on long wavelengths offers possibilities for the design of band-pass phonon filtering. Moreover, the ray tracing model developed in this paper shows that different forms of correlated scattering imprint a unique signature in the heat current autocorrelation function that could be used as a diagnostic in other nanostructured systems.

Cold Start of PEMFC using Alternating Hydrogen Pump: Part I. Mechanistic Modeling

Cold start impedes the diffusion of fuel cell vehicles in regions with sub-zero temperatures. The alternating hydrogen pump (AHP) method has the advantages of fast, efficient and reliable startup without concerns for failure or degradation due to freezing of the product water. However, puzzling phenomena, including the variation of current in each half-cycle and the hook-shaped change of high frequency resistance (HFR), are observed during the AHP under constant voltage mode. These phenomena affect the efficiency and speed of cold starts, but are poorly understood. Herein, we develop a one-dimensional transient model, in which a mechanistic description of the water vapor transport enhanced by H2 convection is contemplated. The model is parameterized and validated using experiments under isothermal and adiabatic boundary conditions, respectively. The simulation results show that the puzzling phenomena are caused by the changes in the content and distribution of the water in the membrane, which in turn are induced by water transport in the membrane and water vapor removal from the membrane. The model has the potential to be applied in the optimization of the cold start process.

A novel thermal lattice Boltzmann model with heat source and its application in incompressible flow

In this paper, the source term of energy distribution in double distribution function approach is analysed in detail and a new discrete model is proposed to eliminate the error terms in previous models. The newly proposed model contains two coupling terms, i.e. a velocity-source coupling term and a temperature-force coupling term. The velocity-source coupling term is related to the transport of source which affects the accuracy of simulations of convection, chemical reaction and so on. The temperature-force coupling term is related to the inner product of force and temperature gradient which affects the accuracy of simulations of boundary layer and the edge of plumes in convection. We extended two new heuristic treatments to handle the adiabatic and isothermal boundary conditions. Their efficiency and accuracy are analyzed theoretically and proved by some numerical validations. Finally, we simulated natural convection flows with compression work and viscous dissipation over a wide range of Eckert number using the present model and obtained a scaling law between average Nusselt number and Eckert number.

Pressure-Tight and Non-stiff Volume Penalization for Compressible Flows

Embedding geometries in structured grids allows a simple treatment of complex objects in fluid simulations. Various methods for embedding geometries are available. The commonly used Brinkman-volume-penalization models geometries as porous media, and approximates a solid object in the limit of vanishing porosity. In its simplest form, the momentum equations are augmented by a term penalizing the fluid velocity, yielding good results in many applications. However, it induces numerical stiffness, especially if high-pressure gradients need to be balanced. Here, we focus on the effect of the reduced effective volume (commonly called porosity) of the porous medium. An approach is derived, which allows reducing the flux through objects to practically zero with little increase of numerical stiffness. Also, non-slip boundary conditions and adiabatic boundary conditions are easily constructed. The porosity terms allow keeping the skew symmetry of the underlying numerical scheme, by which the numerical stability is improved. Furthermore, very good conservation of mass and energy in the non-penalized domain can be achieved, for which the boundary smoothing introduces a small ambiguity in its definition. The scheme is tested for acoustic scenarios, for near incompressible and strongly compressible flows.

Modeling and Analysis of Heat Transfer and Fluid Flow Mechanisms in Nanofluid Filled Enclosures Irradiated from Below

Radiation driven transport mechanisms are ubiquitous in many natural flows and industrial processes. To mimic and to better understand these processes, recently, radiatively heated nanofluid filled enclosures have been extensively researched. The present work is essentially a determining step in quantifying and understanding the transport mechanisms involved in such enclosures. In particular, a two dimensional square nanofluid filled enclosure irradiated from the bottom has been investigated in laminar flow situation. Effects of nanofluid optical depth, inclination angle of the enclosure, incident flux, and boundary conditions (adiabatic and isothermal) have been investigated. Moreover, the temperature and flow fields have been carefully analyzed in the situation ranging from volumetric to mixed to surface absorption modes. Under adiabatic boundary conditions, steady state is unconditionally achieved irrespective of the incident flux magnitude (varied between 5Wm-2 to 50Wm-2), enclosure inclination angle (varied between 0 to 60 degrees) and mode of absorption (surface, mixed or volumetric). However, in case of isothermal boundaries; onset of natural convection and its transition into transient regime is significantly impacted by the mode of absorption and the enclosure inclination angle.

Dynamical modeling of coupled heat and mass transfer process of coalbed methane desorption in porous coal matrix

Although it has been well known that heat transfer occurs with coalbed methane desorption in porous coal matrix, little is understood about the coupling interaction between the heat transfer and mass transfer in this process. In this work, considering the endothermic characteristics of gas desorption and the dependence of gas diffusion on temperature, a coupled gas diffusion and heat conduction model with isothermal or adiabatic boundary condition is developed to accurately describe the dynamical process of heat and mass transfer in coal matrix. The heat induced by gas desorption is treated as an unsteady volumetric heat sink and compared with the traditional uncoupled diffusion model. The governing equations are solved by an implicit finite-difference scheme with the method of Gauss-Seidel iterations. The variations of the fraction of gas desorbed, gas concentration and temperature distributions and heat flux as well as their parametric sensitivity analysis are systematically discussed. It is revealed that the heat transfer in coal matrix has significant effects on the behaviors of gas desorption. With the increasing desorption time, the excess temperature under isothermal boundary condition first drops rapidly and then rises slowly, while the excess temperature under adiabatic boundary condition monotonically drops and changes far more than the isothermal condition. Furthermore, the excess gas concentration decreases with the initial temperature, but increases with the activation energy of gas desorption. The differences in gas concentration under different thermal boundary conditions are increasing with time, but decrease with radial distance. Comparing with the coupled model, the uncoupled diffusion model overestimates the fraction of gas adsorbed. Finally, the heat flux through the coal matrix increases rapidly to the peak value and then decreases slowly to zero.

Optimal Design of Liquid-Cooled Plates for Lithium-Ion Batteries Using Multi-Objective Topology Optimization

This paper proposes a design method for liquid-cooled plates of lithium-ion batteries based on topology optimization. First, as a multi-objective optimization function, the Pareto algorithm is constructed as a weighted combination of viscous dissipation index and homogenous temperature index. Next, the optimal configuration of the channel with different import/export arrangements and different weighting factors is investigated using the fluid pressure drop and maximum temperature as evaluation indices, and the three-dimensional lithium-ion-battery liquid-cooled heat dissipation model is validated on this basis. The results show that under the adiabatic boundary conditions, the channel configuration with the import/export arranged in parallel diagonally on the long side has good overall performance in terms of fluid pressure drop, and so it has the lowest external pump power consumption, the maximum temperature of the channel configuration with the import/export arranged vertically diagonally is the lowest, and the cooling effect is best. Compared with the conventional fence-shaped channel, the optimized channel screened based on the Pareto optimal solution set reduces the maximum cell temperature by 0.21 K, the fluid pressure drop of the channel by 18.25%, and the cell temperature difference by 5.17%; compared with the serpentine channel, the fluid pressure drop is reduced by 85.55%.