Investigations of Energy Conversion and Surface Effect for Laser-Illuminated Gold Nanorod Platforms

: Achieving a quick temperature increase is a burning issue for biophysical applications, like germ inactivation and tumor ablation, and for energy performances, like solar collectors and steam generators. Based on the plasmon resonance phenomenon, noble metallic nanoparticles have emerged as promising weapons due to their very high biocompatibility, optical properties, and high surface-to-volume ratio, increasing energy conversion and allowing the maximum temperature to be reached faster. This work examines the energy conversion in sandwiched glassy platforms with gold nanorods. The platforms are kept vertically in the air and illuminated by a 0.5 W near-infrared laser (808 nm). To describe this aspect theoretically, the size and conversion efficiency of the electro-magnetic properties are compromised between the proposed model and the stability of the nano-rods. As a research approach, our model of cross-sections and polarizability for the surface effect is proposed, coupled with classical CFD numerical calculations. The results of the proposed model, validated by a thermal camera and spectroscopy measurements, indicate that as long as the energy conversion is visible with relatively low-power lasers (ΔT = 18.5 °C), the platforms do not offer fast heat dissipation. The results indicate that, despite the flow forcing by the air inflow, the entropy generation due to heat conduction is more than three orders higher than the dynamic entropy production. Flow forcing corresponds to the value of the velocity for classical convective motions. Therefore, the delivered heat flux must be distributed via convective transport or the associated high-conductive materials.


Introduction
Fossil fuels have dominated the world's energy sources for ages due to their widespread availability.However, it has recently been discovered that their exploitation and widespread consumption are hastening the greenhouse effect.As a result, in recent years, the diligent search for new alternative energy sources and their diversification has become a burning issue [1,2].However, the majority of renewable energy sources do not provide an instant and powerful temperature increase that can last for hours, or they drastically raise the usage costs.Unfortunately, many of the available techniques are based on supplying heat flux to a fluid-solid boundary, implying the slow temperature increase throughout the system.
People realized the importance of hand and room disinfection during the COVID-19 pandemic.Unlike powerful chemical compounds, overheating can reduce the number of germs and viruses in a specific volume to nearly zero.Similarly, tumors and cancers may be totally eradicated; however, there exists the same challenge: how can we avoid destroying healthy cells whilst maintaining the noninvasive ones [3][4][5].
Metallic nanoparticles are emerging as a new strategy due to their own radiative properties.Moreover, their ability to travel through microscale pipes and their high heat transfer coefficient provides a quicker temperature rise than in standard metal sheet heat transfer.Therefore, the future vision generally is to apply these nanomaterials to solar collectors, steam production, and germ inactivation.Although several papers have proposed solutions for nanofluidic flows in the context of mass and heat transfer [6][7][8][9], heat generation in nanoparticles is still in its infancy, and the nanoparticle shape optimization necessitates the numerical computation of sophisticated formulae and equations, extending the time of a single simulation.
This paper extends the theoretical model presented in the authors' previous papers [4,[10][11][12][13] with continuum lasers and highlights the aspects that increase the energy conversion.The presented approach aims to use the validated theories, like the Cattaneo-Vernotte equation, the Rayleigh-Drude approximation, and the electromagnetic energy conversion into heat from the localized-surface plasmon resonance phenomenon.This allows for the study of the electrical field regardless of the nanoparticle's shape, although the model appears to produce dissatisfactory results in the blue and UV regions due to the Rayleigh limitation [14][15][16][17][18][19].Section 2 details these topics particularly for the computational fluid dynamics (CFD) (Section 2.1) equations, energy conversion (Section 2.2), and electric field distribution (Section 2.3) with the influence of the gold nanorods deposited on the charged surface (Section 2.4) [20][21][22].
Furthermore, Section 3.1 illustrates the numerical proceedings and discretized geometry, which were tested on mesh and with time independence to ensure minimal numerical errors [23,24].In Section 3.3, all the boundary conditions supporting the convection and radiation processes are presented, as well as the material properties that were used to solve all the necessary CFD equations [7,[25][26][27].
The mesh and time independence tests are described in Section 4.1.Section 4.2 presents the validation results confirming the validity of the adopted methodology.The results of further analysis (Section 4.3) take into account the movement of the fluid during the heating and cooling process.The discussion of the results and a summary are presented in Sections 4.4 and 5, respectively.
At this stage, it is worth emphasizing that despite the extraordinary importance of the experimental results carried out to obtain the temperature field from a thermal camera and spectroscopy measurements, they do not impart knowledge about the velocity fields and entropy generation in the area of the resulting system.Therefore, one of the goals of this work was the numerical analysis of the temperature field, velocity, and entropy generation in the area where the measurements were not made.
The main purpose of this study was to investigate how a quick temperature increase affects the heat and mass flow around gold nanorod platforms.Several investigations have already touched upon this topic [27][28][29][30], although they concern different nanomaterials under the forced convection process and not the natural one, which is stated to be dominant in this case.Moreover, in contrast to the well-known heat sources, laser beams encourage the use of Gaussian distribution in the space regime, which implies interesting and promising properties, that are worth studying.

CFD Governing Equations
The heat transfer of light-illuminated metallic nanoparticles is generally calculated using various methods, depending on the considered model and the working light source.
The majority of them require the use of the hyperbolic energy equation, also known as the Cattaneo-Vernotte equation [14,15], which is given by г This equation requires the establishment of the relaxation time, г, which defines the finite speed of heat propagation, and which is satisfied for crystalline solids by For noble metals, however, this value equals circa 10 −10 s, and, as may be realized, the Cattaneo-Vernotte equation can be easily reduced to the standard parabolic one if the CW-and millipulsed light sources are considered, as is performed in this paper.
Hence, the standard energy equation, associated with the continuity, momentum, turbulence, and dissipation equations, which are based on the standard -ℰ computational fluid dynamics (CFD) model for flows and the Boussinesq approach, is represented in this work by the following equations [12]: where Furthermore, the total entropy generation, , is investigated as the sum dynamic and static entropy: Static entropy generation is defined in classical expression [28,29], as follows: Using the Boussinesq and turbulent methods with Formula (7), it is possible to determine dynamic entropy generation [30,31], as follows: The total entropy generation parameter determines the irreversibility of each process and may also indicate, alongside the Richardson number, which type of convection, forced or natural, dominates in each simulation stage.

Energy Conversion Rate
The most relevant aspect in this work is to estimate this part of incident light energy that is directly responsible for energy conversion to heat.In the standard energy equation, this is described by a user-defined function, which is referred in this work as   f and   np , respectively, for fluids and nanoparticles.Some assumptions, however, should be considered in this case.Firstly, if a laser beam falls onto the surface at a low incident angle, the problem is easily reduced to the semi-infinite methods and Lambert-Beer-Bourger law, for which the heat generation formulae for thin films are satisfied by [4,7,10,13,16] where the absorption coefficient is represented by This approach is validated and commonly used, both for continuous media and porous structures (see Figure 1a).Hence, it has been assumed that nanoparticles can be treated as inclusions in a -thick slab, and each of them contributes to the total energy conversion (Figure 1b) by the following formula: where   ,   , and   are polydispersive (distribution) parameters, respectively, for size, length, and distance, and   is the thickness at which the electric dipoles can occur, declared as the slab thickness.Here, it is assumed as the size of capped nanoparticles, written by the following: Moreover, the optical coefficients,   , are calculated differently, depending on the materials.Classical Maxwell equations imply that the incident light may convert the electromagnetic energy if the absorbed light is determined.For nanoparticles, this is described by the photothermal efficiency equation: and it refers only to nanoparticles.Hence, this is not the parameter that reflects the total efficiency but only facilitates absorbance spectroscopy measurements.The total energy conversion efficiency indicates the absorbed part over the incoming light, which is presented as follows: This is purposely emphasized to highlight what the standard laser beams imply, when assuming the normal distribution, according to the following formula: where ℛ  is the beam size.

Cross-Sections and Polarizability
For a perfect crystal structure, the energy conversion is governed by the diffuse transport of electrons in lattice.Generally, two-temperature models are used to calculate the heat exchange between electrons and crystal lattice.For nanoparticles, however, boundary effects are extremely crucial to be included.In this work, the plasmon resonance phenomenon suggests that electrons cannot be treated separately, but as collective vibrations over the whole nanoparticle.Since the vibrations are confined by a capping agent or a surface, they convert and release the accumulated energy as heat.In this work, this phenomenon is determined by optical cross-sections, which were primarily studied in papers by Mie and Lorentz [17,18], but their theories were too complex for nanostructures other than spheres.Hence, the Rayleigh-Drude approximation is considered to be validated as long as the following condition is satisfied: Then, general assumptions of this theory refer particularly to nanoparticles illuminated by visible and infrared lasers.The representative optical cross-sections are described by [17,18]   ,, =   ,, As can be seen, the cross-sections,   ,, and   ,, , as an average electric field distribution, can be described easily by the effective parameter, α eff i,j, .For capped ellipsoids, this is assumed as (for Cartesian coordinates) where , , and ℎ refer to core, agent, and host (working medium), respectively, and + indicates the oriented direction amongst the Cartesian coordinates, , , and .
For rods, however, the polarizability is modified by including the correction for flat dimensions (Figure 1).It was introduced primarily in Fuchs' publication [19] that the elongation parameter (  , ) or depolarization factor ( Л , ) may indicate that the opposite charges prefer to be accumulated at edges, generating dipoles.Hence, as long as these parameters tend to 1.0 in ellipsoids, wires, or disks, the situation changes significantly in rods.Then, the flat, long dimension contributes to configuring new dipoles.This is clearly visible in Figure 2, where each axis projection is illustrated.Subsequently, an update to the effective polarizability yields the following: Although the effective parameter facilities including size, shape, distance, and surface effects, there is virtually zero examined reports studying the influence of these effects on the energy conversion.To calculate temperature distributions, however, the interaction between fluid and particles is primarily considered as an empirical boundary condition.For instance, the Smoluchowski approach investigates the heat transfer coefficient as the effective parameter for heat transfer [10,14,19,31].The model, however, is well described for two boundaries, and nanostructures deposited on a loaded surface necessitate the use of a surface effect where at least three different materials must be considered.

Surface Effect
The possibility of implementing the effective polarizability offers a way to adjust different effects to achieve the most satisfactory approach.Formulae (19) and ( 20) have considered the case where a single nanoparticle is located in a soluble medium or deposited on an unloaded surface, such as glass.In this work, however, a 3 nm polyelectrolyte material (PEL) whose preparation and chemical processing enable creating ions, which tend to dissociate at the boundaries, is used.The creation of a negative charge contributes to modifications in the electric field through quasi-dipole-induced dipole interactions, as seen in the PEL-AuNR case.This affects the energy conversion as well as increasing the effective polarizability,  eff ,, , and the slab thickness,   .The sandwiched structure and the scheme of electric field distribution are revealed in Figure 3.
The topic, known as a surface effect, has been studied particularly in the literature for different geometries and parameters.The most common approximation is to assume the arithmetic mean of either permittivities [17] or polarizabilities [20] between the working fluid and base material, which provides only qualitative results.Other models assume the electric field-based results, which become complex for shapes other than spheres [20].Here, the Royer-Yamaguchi approach [21] is claimed to provide satisfactory results due to the simplicity of combination surface and distance effect altogether.Although the model is dedicated to oblate structures, some corrections for rods based on Equations ( 20)-( 22) require transformations that are necessary to be accomplished.Having taken into account the nanoparticles' orientation and the projection on the axis, the substrate calculations can be presented as follows: ( ,,, ) − = ( ,, ) ) ) ) ) ) ) ) ) ) and distance effect, according to the paper by Royer [21]: where the prolation parameter,  , , specifies the selected projection, as highlighted in Figure 2. From Formula ( 22),  ,, ,  ,,, ,  ,,, ,  ,,,, ,  ,,,, , and  ,,, are determined as follows: ,,, =  ,,,, = √1 − ( In a similar way, having substituted the parameters ( 31)- (34) to the ( 23)- (30), and further to (22), the effective polarizability is possible to be calculated.

Considered System
Figures 4 and 5 depict a simple system consisting of two pieces of borosilicate glass surrounded by air and joined together by a transparent glue.This is to serve as a thermal transducer.On the top and bottom sides of the smaller piece, gold nanorods are dislocated.Glass is prepared here due to its high transparency, facilitating heat to be produced within the system.On the top of the glassy domain (4), where the gold nanorods are dislocated, parallel laser beams are directed.Due to the great distance between them, some light penetrates the top and reaches the second, bottom gold surface.

Calculations of Energy Conversion Rate
Calculations from a total number of 10 23 particles cannot be obviously accomplished due to the simulation time that would tend to infinity.This is why   ,   , and   are reduced and divided into 50 × 500 × 200 parts.The expected values are determined as    = 15 ± 3 nm ,    = 55 ± 8 nm , and   = 120 ± 30 nm .As a function which describes polydispersity, Gaussian and Lorentzian distributions are selected.The boundaries are fixed as 3 × ∆   , 3 × ∆   , and 3 × ∆  .Table 1 summarizes details of the considered nanoparticles, whereas

Numerical Model Implementations and Mesh Independences' Proceedings
After determining the heat production rate generated by nanoparticles under laser irradiation, as described in Equations ( 11) and ( 13), we carried out numerical simulations based on the problem presented in Figures 4 and 5, using the finite volume method (FVM).This computational technique discretizes the physical model into small control volumes and solves the governing equations at the centroids of these volumes.The calculations are run on the TrytonPlus supercomputer (CI TASK computing cluster, Gdansk, Poland, with the InfiniBand Fabric EDR 100 Gbps, fat tree topology network), which has 72 cores and supports the ANSYS Fluent software (versions 22.2 and 23.2).As outlined in the mathematical formulation section, we have treated Equations ( 2)-( 7) relating to turbulent flows in three dimensions, using the standard -ℰ turbulence model available in ANSYS Fluent software.For speed-pressure coupling, the SIMPLE algorithm is utilized, complemented by the PRESTO scheme for the spatial discretization of the pressure term, and a secondorder upwind scheme for discretizing convective terms in the momentum and energy equations ( 2)- (7).Under-relaxation factors of 0.3 for density and 0.5 for other parameters were used.Heat sources from Equations ( 11) and (13) were incorporated as boundary conditions on the considered surface, as indicated in Figure 4, using user-defined functions (UDFs) with the 'DEFINE_PROFILE' macro.UDFs dynamically distribute heat flux values over the surface under consideration, aligning them to the Gaussian distribution of the laser beam, as specified by Equation (17).These functions calculate the heat source term for each face of the control volumes during runtime, which enable us to model heat flux with spatial precision.The UDFs integrate the position vectors of the control volume faces and the laser beam intensity parameters to translate the physical beam profile into a precise heat source distribution during simulations.
Figure 6 reveals three meshes that were tested, providing minimal numerical error, using the Roacher and Richardson extrapolations [23,24].A similar procedure is utilized for a time regime.The similar procedure is used for a time regime.The Roacher technique determines the continuum values of a given quantity   and highlights the relative and fractional error between different mesh or timestep results.This method is also intended for lower-order discrete values and is based on the following equation: The continuum value  =0 for the second-order scheme employed in this work yields: Similarly, for the three testing steps, Equation ( 25) is transformed into which enables the investigated value to be extrapolated with a specified error, e.g., the relative  and fractional error  1 between two sets of data that follow Furthermore, the grid convergence index () methodology has been introduced to indicate how far the error may change as the mesh is finer, which, for Richardson extrapolation, is defined by the  parameter: where -Relative error between two sets/values of data;  =0 -Continuum value as if the mesh or time size equals zero;

Boundary Conditions
Boundary conditions allow the optimization of and reduction in the number of domains required to solve differential equations numerically.In this work, the system is assumed to be maintained under mixed conditions, incorporating both convection and radiation.To stabilize the computations and to omit the randomness and the initial zero problem, at starting point, a 1 μm/s velocity inlet and a pressure outlet along the gravity vector were assumed.Furthermore, due to the presence of an organic capping agent (citrate buffer), the application of the Marangoni stress to gold nanoparticles was necessary for managing the diffusive momentum flux: where the thermal coefficient of surface tension equals −5.263 As a heat source, a 808 nm laser was positioned on the top surface of the platform where gold nanorods are localized.This wavelength was selected due to the similar absorbance of gold nanorods for air, water, and glue as a working medium, as explained in further sections.In the spatial regime, the laser beam followed a Gaussian distribution, presented in Formula (15), and the laser output power was set to 500 mW.Higher powers were not considered due the power saturation, which would necessitate further modifications to the model.These would be subjects for future investigations.The laser-induced heat flux was delivered over a period of 70 s (irradiation time).After this time, the laser was off and the system cooled for 20 s.
Table 3 displays thermal properties of the materials used to prepare the system, whereas Table 4 summarizes the assumed boundary conditions.

No. Domain Boundary Conditions
(1) Bottom glass Top glass

Mesh and Time Independence Tests
Three meshes were tested to ensure the independence of the result to the mesh grid and time step size (Figure 6).Table 5 highlights details of all created grids.It should be added that the cells with maximum aspect ratio are located in the bottom glass (1) or in the glue zone (3) (see Figure 2).The maximum temperature in the entire system (max) and the average temperature of the top glass (4) were selected as quantity values for the Richardson extrapolation and GCI methodology, both in spatial and time regimes.The erroneous data are revealed in Table 6.
Given that the spatial errors were under 1%, the procedure was rerun for three timesteps- 1 = 5.0,  2 = 2.5, and  3 = 1.0 ms.The resulting comparison suggests that time errors of up to 0.1% are feasible.The 2.5 ms mesh results, being spatially and temporally independent, are now advanced to experimental validation.

Experimental Validation
To validate the theoretical model, AuNRs were deposited on a glass substrate using the immersive layer-by-layer assembly technique according to the protocol reported in reference [3].This technique allows for the immobilization of AuNRs on both sides of the glass substrate.After this step, a double glass cell was prepared by the scrupulously performing procedure reported in reference [4].Then, the refractive index surrounding the AuNRs was consecutively varied from the front to the back sides of the AuNR substrate.Additional details can be found in reference [3].
Figure 7 depicts the absorption spectra of the AuNRs' double cell, consecutively infiltrated with the glue (NOA-61) from the company Norland Products Inc. (Cranbury, NJ, USA), with a refractive index of 1.54 (from spectrophotometer PerkinElmer A Lambda 365 (Waltham, MA, USA)).The experimental optical behavior (continuous lines) agrees well with the theoretical model (dashed lines).
As may be noticed, the commendable agreement enables the ease of extension to almost the entire visible region and a part of the infrared region, up to 850 nm.
Figure 8 demonstrates the photothermal properties of the AuNR substrates upon laser light illumination, performed by using the optical setup and the experimental parameter reported in reference [3].
The results related to the temperature distribution, presented in Figure 8a, and the time-temperature profile (Figure 8b) demonstrate the full consistency between the experimental results and theoretical model.Validation of the model was conducted using a thermal imaging camera and the temperature field collected from the surface of the system under study [43].
As a consequence of this notification, this theoretical approach has been declared to be representative for the energy conversion and temperature distributions in this work.This decision is supported by the model's reliability and accuracy in capturing the intricate interactions of light absorption and heat transfer dynamics.Moreover,  = 808 nm wavelength has been declared to possess a similar energy conversion for all three cases and only the materials' properties and boundary conditions can affect heat transfer and the values of temperatures.Therefore, this laser is shifted to a further analysis.Furthermore, the photothermal behavior of the AuNRs' double cell was also investigated using a white light source, as outlined in the experimental procedure described in reference [4].Experimental data are presented in Figure 9.
The model's validity has also been raised and confirmed for the white-light lamp in [4] (Figure 9), notwithstanding the fact that the correlation is more pronounced for a single-laser approach.This can be attributed to the system's irregularities, which pose a challenge to achieving ideal geometry.Consequently, optimization processes are deemed necessary and beneficial.

CFD Results for 0.5 W Laser Illumination
It should be noted that the thermal imaging camera analyzed the temperature field on the surface of the system for both the systems in Figures 8 and 9. Consequently, the measurement tool employed does not permit a combined analysis of the temperature and velocity fields, and CFD calculations prove to be invaluable.The areas 5 μm above and 5 μm below the AuNRs' platform of Figures 4 and 5 were selected to present the results.As can be observed, the maximum temperature reaches 40.3 °C below and 42.5 °C above the top glass (4), which is a promising temperature increase provided by a low operating laser power.In the real situations, higher powers are anticipated to raise the system's maximum temperature significantly, as presented in Figures 8 and 9

Cooling Stage
Similarly, Figures 14-17 illustrate the abovementioned results after 10 s of cooling when the whole system was reasonably close to being thermally equilibrated.In contrast to the heating stage, some values decrease by an order of magnitude as the temperature reduces.The velocity vectors, on the other hand, continue being horizontally oriented and their values have not been modified, whereas the maximum value of dynamic entropy slightly increases.

Time Responses
Time responses are strongly dependent on the boundary conditions and the complicity of the assumed mathematical functions.Some nanoscale results are based on Green's or Gaussian functions, which are discussed in [7,[25][26][27].However, the presence of convection, which occurs in this work, encourages the use of numerical methods, which offer a similar effect.
Figure 18a depicts the maximum temperature, Figure 18b depicts the average temperature of the system, and Figure 18c

Discussion
The results in Figure 10 show a pronounced increase in temperature centered above and below the gold nanorod (AuNR) platforms, as expected due to the localized photothermal effect.The temperature contours suggest a strong radial distribution of heat from central hotspots corresponding to areas of direct laser irradiation.The maximum temperature recorded at these points corresponds to the known thermal properties of gold nanorods when subjected to laser radiation.It is also interesting to note that the temperature field decreases rapidly away from the center, indicating the efficiency of borosilicate glass's thermal insulation properties.The uniformity observed in the temperature distribution suggests a homogeneous thermal conductivity distribution in the glass medium.Moreover, the colder regions around the heated zones confirm the role of air as an insulator, which is supported by the distinct boundary between the high-temperature regions and the surrounding area.During the cooling phase (see Figure 14), the temperatures in the core where laser irradiation occurred demonstrate a clear reduction, indicating efficient heat dissipation in the borosilicate glass and surrounding air.The uniform temperature gradient reflects the glass's constant thermal conductivity and the air's insulating role.The corresponding static entropy fields for the observed temperature fields are displayed in Figures 11 and 15.As shown in Figure 11, the static entropy contours during the heating stage reflect the degree of disorder induced by laser heating for 70 s.Just above and below the AuNR platforms, there is a pronounced increase in entropy, indicating significant thermal activity.The central regions show the highest entropy values, revealing intense energy dissipation, which must be connected to the direct interaction of the laser with the AuNRs.As shown in Figure 15, during the transition to the cooling stage, the entropy levels decrease overall as the system loses thermal energy.However, the entropy distribution scheme remains similar, with the highest levels still concentrated around the central areas exposed to the laser.This region of constant high entropy suggests that AuNRs and associated regions retain heat longer or cool slower than surrounding materials.It should also be noted that the entropy does not immediately return to pre-heating values, suggesting a certain amount of thermal inertia or the presence of irreversible processes, even after the laser has been switched off.
As can be noticed in Figures 12 and 16, the dynamic entropy values are approximately three orders of magnitude lower than the static entropy, indicating that, although fixed inlet boundary conditions stabilize calculations, they do not play a dominant role in forced convection and do not contribute to the total entropy and the general case.Indeed, the low dynamic entropy demonstrates that the thermal gradients and localized photothermal effects induced by laser irradiation on the AuNR platforms are the primary drivers of entropy generation rather than the fluid motion or inlet flow conditions.This observation holds for both heating and cooling phases, reaffirming that heat transfer mechanisms primarily drive the entropy scenario in this study.Moreover, the flow is observed for the glassy domain (4) from above, in the immediate vicinity of gold nanorods and only from the top.Due to the highly viscous forces, this effect is not visible in the confined area (bottom) where external flow and gravity cannot disrupt the boundary, and only free convection prevails inside.Furthermore, in the entire air domain, velocity vectors are oriented along the gravity force (Figure 13).However, the magnitude of vectors is distributed similarly to the temperature gradient and the temperature wave.
Nevertheless, vectors are distributed by the Gaussian laser profile from Equation (17), and residual flow is elevated due to the buoyancy forces and air density dependence, as presented in Table 3.Hence, competition arises between fluid weight, buoyancy forces, and the assumed laser profile.This competition is not observed at the same rate if gravitational forces act in parallel to the heat source as they do when they act perpendicular to it.This is because gravity forces attempt to retain the fluid while buoyancy forces strive to circulate it.Since viscous forces prevent the flow from entering the confined volume, velocity changes are governed only by natural convection in this region (see Figure 18c), contrary to the plane (top) where the system is exposed directly to the oriented flow.It may be noticed by studying the differences in the dynamic entropy contours (see Figures 12 and 16).The topic is detailed in [44].
Furthermore, Gaussian (single-mode) and round (multi-mode) laser beams are distinguished by a similar radial temperature distribution, as seen in nozzles [45], where heat is delivered uniformly.Spatial distribution, characteristic of a Gaussian profile, however, offers an increase in the maximum temperature, although the energy conversion rate is on average comparable.It allows us to declare that a single-mode laser is sufficient to reach about 42.5 °C and 40.3 °C after 70 s of illumination (see Figure 18a) for a 0.5 W power and sparsely distributed nanoparticles if the surface is charged.Assuming the linear dependence between energy conversion rate and laser power, a two-watt near-infrared laser could potentially heat the system to reach 100 °C at this concentration.In reality, the energy conversion rate reaches a limit value that corresponds to the saturation temperature, as discussed in [46].
Additionally, produced heat appears to be distributed slowly due to low-diffusivity materials.This aspect is visible in Figure 18b, where both curves are similar.However, their tendencies are distinguishable due to the influence of natural convection and the magnitude of the velocity vectors.As long as the buoyancy effect intensifies heat transport, this influence is weak in the confined region where the region's size (10 μm) allows the flow to increase only during a few hundred milliseconds.This is why the time responses of velocities distinguish considerably from one another, whereas maximum and average temperatures remain comparable.In turn, while the irradiation is turned off, the system enters a cooling stage.The dissipation of heat exceeds the input during this phase, leading to a rapid decrease in temperature, and velocity.The diminishing thermal input allows the system to return to a more disordered state-the interaction of temperature, velocity, and entropy describes the system's dynamic response to a constant heat source (generated heat by AuNPs), demonstrating the complex balance between heat input, convective flow, and thermodynamic disorder.The complex interaction of several factors, including heat input, convective flow dynamics, and the system's response to spatial constraints, generates entropy, as discussed in [47][48][49].
The observed changes in the presented results can be attributed to the system's behavior during the simulation.The balance between heat input and output in the air medium influences the initial temperature rise.This phenomenon can be linked to the difference between the rate of heat absorption by the air and the dissipation rate.The velocity increases as the air near the wall heats up.This increase in velocity is accomplished by convective currents (free convection) induced by the temperature gradient.
The differences in parameter values observed between the air-exhibited area and the confined microscale volume can be attributed to the impact of spatial constraints on heat transfer and fluid dynamics.The restricted space in the confined microscale volume may result in enhanced heat conduction effects and a more localized influence of the heat source, resulting in different temperatures, velocities, and specific entropy patterns.Convective currents and heat dispersion, on the other hand, may operate more freely in the open area, where there are fewer spatial constraints, leading to higher values of temperature and velocity [50][51][52].
On the other hand, heat transfer in this case is not sufficiently fast to reach the steady state and a high temperature in less than one minute, which is indicated in tumors overheating.However, in this area, analyses of the material properties of biological structures and dedicated numerical calculations are necessary.Here, high-conductivity metallic sheets would reduce the effect of energy conversion, resulting in a low temperature increase.Therefore, glassy platforms with nanoparticles are declared only to convert energy and transduce it to the high-diffusivity materials, like metal sheets, that would distribute the heat throughout the entire system much faster than insulators.It is highly recommended studying this topic in the future.At this point, it is worth mentioning that in previous works, the authors have focused on comparing the CFD numerical model with simple geometry calculations [7,10].For classical geometries, large compatibilities are obtained when using CFD tools with analytical methods.However, due to the need to develop complex geometries with the full Navier-Stokes equation as in the case of wind farms [53], it is necessary to take into account the advantages of CFD and perform simulations using this approach.Also, in the presented geometry and in the given physical phenomenon in this paper, there are reasons to use the CFD approach.
In addition, the work [54] compared three different theories of energy conversion into heat due to unique properties of gold nanoparticles (AuNPs), which may boost the efficiency in selected wavelengths.White-light-lamp-illuminated nanoplatforms with AuNPs were investigated via radiative transport, which were solved numerically using the discrete ordinates method (DOM).Optical boundary conditions were calculated analytically from data in the literature using three different approaches to describe the absorption coefficient, which defines the efficiency of the conversion.The selected approaches supported either the classical Maxwell model, Rayleigh approximation for nanostructures the Mie-Lorenz theory for spheres.Likewise, in order to include substrate effect, a thin-film-adjusted model, developed by the authors of [4], was assumed.The results in [54] reveal a dozen Celsius-degrees increase for sparsely aligned gold nanoparticles, and vary significantly depending on the utilized model.Thus, the model presented in this paper is deemed the most appropriate compared to those mentioned above.

Conclusions and Perspectives
Although metallic nanoparticles can generally be modeled using standard transport equations described by exact differentials, they necessitate sophisticated mathematics, indicating the complexity of boundary conditions.According to this paper, metallic nanoparticles can be advantageous in applications where time is of the essence, such as cancer therapy and steam and bubble generation.The examined topic benefits from and encourages further research due to its promising temperature increase in energy or biomedical systems.The adopted model has been validated for a single-laser and a white-light source and could be further developed to accommodate more complex cases.Although NRs are not truly present, the interaction between NRs and light reflects the physical response, demonstrating that simple and well-known equations would unquestionably reduce simulation and CPU time.This work has also reported that a fixed inlet with a minimum velocity stabilizes calculations and improves results.Although it appears to be exaggerated and excessive, the gravitational forces overcome the inlet velocity after a few milliseconds, and the system follows the advisable conditions and results.According to the research, conduction is the most prevalent heat transfer mechanism.Despite the presence of viscous forces, their effect on overall entropy is minor.In addition, the low dynamic entropy observed compared to the static entropy indicates that the thermal effects of laser irradiation predominate over fluid dynamics in the entropy generation within the system.CFD models were compared with the analytical approach for solids in the works [7,10], confirming that CFD and analytical models give similar results for simple geometries.However, the CFD model used in this work is proven to be much more effective for complex geometries and fluid analysis [53], and therein lies the superiority of this tool over analytical solutions.
Additionally, the low dynamic entropy observed compared to the static entropy indicates that the thermal effects of laser irradiation predominate over fluid dynamics in the entropy generation within the system.
Our results reveal that the top plane of the gold nanorod platform exhibits a faster and more pronounced temperature and velocity response to irradiation than the bottom plane.In particular, the maximum temperature on the upper plane reaches around 43 °C after 70 s of irradiation, while the temperature on the lower plane increases gradually to around 38 °C.Both planes eventually exhibit similar behavior for the average temperature, reaching around 30 °C during the heating at 70 s.Velocity at the top plane quickly peaks at around 0.0008 m/s and remains at this level for 70 s, then rapidly decreases during the cooling phase.The theoretical results outline that even relatively low laser power values are sufficient to produce an impressive temperature increase.The maximum value in this work equals approximately ∆ ≈ 18.5 o C and ∆ ≈ 35.5 o C under a 808 nm illumination, respectively, for the 0.5 W and 1.0 W laser.However, challenges occur with heat transfer throughout the entire system.As long as convection and conduction can distribute heat, they cannot compete with changes in heat generation values.Consequently, lowthermal diffusivity and optically transparent materials, such as glass or epoxy, appear adequate for energy conversion, whereas other materials are preferable for intensifying heat transfer.The oriented flow, however, implies a heat wave whose front can affect and accelerate the achievement of high temperature on average using convective transportation methods.
Considering the present aspects of increased entropy generation with an increase in flow velocity, but an increase in heat transfer rates, there will be a need to perform an optimization task combining opposite objective functions.A multi-variant optimization approach similar to that described in the work of [55] is anticipated for further study.

Figure 1 .
Figure 1.Scheme of energy conversion calculations for (a) thin films and (b) nanoparticles.

Figure 2 .
Figure 2. Assumed electric field distribution (dipole configurations) for a 15-15-55 nm gold nanorod (orange color) capped by a 4-nanometer-thick citrate buffer (green color).Projections on the ̂ axis are marked by a violet color.

Figure 3 .
Figure 3. Scheme of the electrical field distribution with AuNRs deposited onto the PEL sandwiched functionalized substrate (side view).This reflects the mentioned name 'AuNRs platform'.

Figure 4 .
Figure 4. Scheme and dimensions of the gold nanorod platform from side view (scale has not been saved).

Figure 5 .
Figure 5. Scheme and dimensions of the gold nanorod platform from top view (scale has not been saved).

𝐶 1 ,
2,3,… -Constant or a function of the selected quantity that is other than ;   -Factor of safety, here: 1.25.

Figure 6 .
Figure 6.Three discretized meshes employed in computations on both sides.Mesh at the bottom has been assigned to demonstrate theoretical results.

Figure 7 .
Figure 7. Experimental validation of the absorption spectra using a spectrophotometer for three different working media,  ℎ deposited on the 3 nm PAH-PSS polyelectrolyte.The structure of the capping agent, as an organic material, is difficult to predict theoretically without stochastic methods, and is conceived to explain visible discrepancies.

Figure 8 .
Figure 8. Experimental validation of (a) the temperature distribution and (b) time response from a thermal camera for a 1 W laser approach.

Figure 9 .
Figure 9. Experimental validation from a thermal camera of (a) temperature distributions and (b) time response for a white-light approach based on [4].

Figure
Figure 10a,b exhibit the temperature after 70 s of light illumination, when the delivered heat begins to stabilize the temperature.Likewise, entropy contours are highlighted in Figures 11 and 12 to control if the forced convection and inlet boundary conditions significantly changes the temperature contours.The left side (a) presents the plane 5 μm under the glass (4) domain in the nearly confined (closed) volume, whereas on the right side (b) results are considered 5 μm above the glass (4) domain, as exhibited in the free space.Accordingly, Figure 13 reveals fluid flow via the velocity vectors.As can be observed, the maximum temperature reaches 40.3 °C below and 42.5 °C above the top glass (4), which is a promising temperature increase provided by a low operating laser power.In the real situations, higher powers are anticipated to raise the system's maximum temperature significantly, as presented in Figures8 and 9.
Figure 10a,b exhibit the temperature after 70 s of light illumination, when the delivered heat begins to stabilize the temperature.Likewise, entropy contours are highlighted in Figures 11 and 12 to control if the forced convection and inlet boundary conditions significantly changes the temperature contours.The left side (a) presents the plane 5 μm under the glass (4) domain in the nearly confined (closed) volume, whereas on the right side (b) results are considered 5 μm above the glass (4) domain, as exhibited in the free space.Accordingly, Figure 13 reveals fluid flow via the velocity vectors.As can be observed, the maximum temperature reaches 40.3 °C below and 42.5 °C above the top glass (4), which is a promising temperature increase provided by a low operating laser power.In the real situations, higher powers are anticipated to raise the system's maximum temperature significantly, as presented in Figures8 and 9.

Figure 10 .
Figure 10.Temperature contours, according to Figures 4 and 5.A distance of 5 μm under (a), and 5 μm above (b) the AuNR platforms (5) after 70 s of the laser illumination.

Figure 11 .Figure 12 .Figure 13 .
Figure 11.Static entropy contours, according to Figures 4 and 5.A distance of 5 μm under (a), and 5 μm above (b) the AuNR platforms (5) after 70 s of the laser illumination.Gray color indicates places where glue has been dislocated.
represents the maximum velocity fluctuations over time.Results consider the cross-sections discussed in Figures 10-17, located near the AuNRs' platforms, at the bottom and top sides of the upper glass (4).

Table 2
indicates electrical properties under 808 nm laser illumination.

Table 1 .
Parameters of gold nanorods.

Table 2 .
Electrical properties for the utilized materials.

Table 3 .
Material properties in the considered simulations.

Table 4 .
Specified boundary conditions in the simulations carried out.

Table 5 .
Details of the applied meshes.

Table 6 .
Parameters under Richardson and Roacher extrapolation proceedings after 30 s of illumination by a 1 W laser.
selected quantity or testing parameter, here: temperature,    turbulence kinetic energy generation due to velocity gradients,  •  − •  −  Б turbulence kinetic energy generation due to buoyancy forces,  •  − •  −  ⃗⃗⃗ gravity,  •  −  numerical parameters to be compared (here: mesh size or timestep)  heat transfer coefficient,  •  − •  −   laser intensity profile at its maximum   ( ⃗⃗) initial laser intensity profile,  •  −   ,, absorbed part of the -particle,  •  −   effective thermal conductivity coefficient,  •  − •  −   refractive index of selected material, - numerical constant for 3- correction, -  turbulent Prandtl number for kinetic energy, -  turbulent Prandtl number for rate of dissipation, - Prandtl number, - pressure,   ,, reflection coefficient of nanoparticles, -  reflection coefficient of the continuous material, -  beam size of laser source, here: 0.0015   ⃗⃗ radius in spherical coordinates,     source of energy for fluids,  •  −    source of energy for the -particle,  •  −   initial temperature at the  =  ,     ,, effective polarizability of the -particle,    ,,,̂ polarizability of the -particle,    anisotropic factor for distance and surface effects,   gold slab thickness, here: size of capped nanoparticles,   thickness of the considered material,    permittivity in vacuum, -  relative permittivity of the considered metal, -  relative permittivity of the core, here:   =   ,   relative permittivity of the capping agent,   relative permittivity of the working (host) medium, -  relative permittivity of the base where NRs are deposited, - relative emissivity, - rank of accuracy, here:  =   relative error, - incident wavelength,   molecular viscosity,  •   ,, nanoparticles concentration,  −  density of a material,  •  −  ⃡⃗ surface tension parameter,  •  −  , elongation parameter, -