Next Article in Journal
Export Constraints Applicable to Renewable Generation to Enhance Grid Hosting Capacity
Previous Article in Journal
Experimental Study of Pollutant Emissions from Biomass Combustion and Modeling of PM Transportation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Piotr Radomski
1,
Federica Zaccagnini
2,
Paweł Ziółkowski
1,*,
Francesca Petronella
3,
Luciano De Sio
2,
Aimad Koulali
1 and
Dariusz Mikielewicz
1
1
Faculty of Mechanical Engineering and Ship Technology, Energy Institute, Gdańsk University of Technology, Narutowicza 11/12, 80-233 Gdańsk, Poland
2
Department of Medico-Surgical Sciences and Biotechnologies, Sapienza University of Rome, Corso della Repubblica 79, 04100 Latina, Italy
3
Department Chemical Sciences and Material Technologies, Institute of Crystallography of National Research Council, Area della Ricerca Roma 1, Strada Provinciale 35d, 9, 00010 Montelibretti, Italy
*
Author to whom correspondence should be addressed.
Energies 2024, 17(11), 2587; https://doi.org/10.3390/en17112587
Submission received: 8 April 2024 / Revised: 4 May 2024 / Accepted: 8 May 2024 / Published: 27 May 2024
(This article belongs to the Section A2: Solar Energy and Photovoltaic Systems)

Abstract

:
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 electromagnetic properties are compromised between the proposed model and the stability of the nanorods. 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.

1. 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 Section 4.4 and Section 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.

2. Theoretical Models

2.1. 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
Γ · 2 t 2 ρ · e + t ρ · e + d i v ρ · e v = d i v k e f f · g r a d T + τ p · v + S e f + S e n p
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
Γ = 3 · k e f f ρ · c p e f f · 1 v d 2
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 k - E computational fluid dynamics (CFD) model for flows and the Boussinesq approach, is represented in this work by the following equations [12]:
t ρ ρ v ρ e ρ k ρ E + d i v ρ v ρ v v ρ e v ρ k v ρ E v = d i v 0 τ p τ p · v + k e f f · g r a d T μ + μ t P r k · k μ + μ t P r E · E + 0 ρ ρ o · g S e f + S e n p G k ρ E + G Б C 1 ϵ · E k · G k + C 3 E · G Б C 1 ϵ · ρ · E 2 k
where
τ = g r a d X 2 3 μ d + 2 μ d
G k = 2 μ t · d · d T
G Б = Б · μ t P r · g · g r a d T
R = 2 3 ( ρ k + μ t d ) + 2 μ t d
Furthermore, the total entropy generation, S , is investigated as the sum dynamic and static entropy:
S = S s + S d  
Static entropy generation is defined in classical expression [28,29], as follows:
S s = ( k e f f + k t ) T · g r a d T 2
Using the Boussinesq and turbulent methods with Formula (7), it is possible to determine dynamic entropy generation [30,31], as follows:
S d   = μ T · 2 3 μ d + 2 d + 1 T · 2 3 ( ρ k + μ t d ) + 2 μ t d
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.

2.2. 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 S e f and S e n p , 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]
S e f = A a b s M · I a b s M = A a b s M · I o r · 1 R M · 1 exp A a b s M · δ
where the absorption coefficient is represented by
A a b s M = 4 π · i m n M λ λ
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:
S e n p = i W i w i · j W j w j · L W L w L · A a b s i , j , L · I a b s i , j , L = i W i w i · j W j w j · L W L w L · A a b s i , j , L · I o r · 1 R g l · 1 R i , j , L · 1 exp A a b s i , j , L · δ i
where w i , w j , and w L are polydispersive (distribution) parameters, respectively, for size, length, and distance, and δ i 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:
δ i = d s i + 2 · d g
Moreover, the optical coefficients, A a b s , 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:
η P T i , j , L = I a b s I e x t C a b s i , j , L C e x t i , j , L
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:
η T O T i , j , L = I a b s r I o r 1 R g l · 1 R i , j , L · 1 exp A a b s i , j , L · δ i
This is purposely emphasized to highlight what the standard laser beams imply, when assuming the normal distribution, according to the following formula:
I o r I o M a x · exp 2 · r R B 2
where R B is the beam size.

2.3. 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:
λ 2 π > d s 2 · d l ̿ 1 3
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]
A a b s i , j , L = η P T i , j , L · ξ i , j , L · C e x t i , j , L = ξ i , j , L · C e x t i , j , L C s c a i , j , L = ξ i , j , L · 4 π · 2 π λ · i m α i , j , L 8 π 3 · 2 π λ 4 · α e f f i , j , L 2
As can be seen, the cross-sections, C e x t i , j , L and C s c a i , j , L , as an average electric field distribution, can be described easily by the effective parameter, α e f f i , j , L . For capped ellipsoids, this is assumed as (for Cartesian coordinates)
α e f f i , j , L = 2 3 · α i , j , L , y + α i , j , L , z + 1 3 · α i , j , L , x , d l +
and for the considered axis, i ^ , it is
α i , j , L , i ^ = V c i , j + V g i , j 4 π · ε g ε h · ε s + ε c ε g · Л i , j , L , i ^ , c δ i · Л i , j , L , i ^ , g + V c i , j V c i , j + V g i , j · ε g · ε c ε g ε h + ε g ε h · Л i , j , i ^ , g · ε g + ε c ε g · Л i , j , L , i ^ , c δ i · Л i , j , L , i ^ , g + V c i , j V c i , j + V g i , j · Л i , j , L , i ^ , g · ε g · ε c ε g
where c , g , and h refer to core, agent, and host (working medium), respectively, and + indicates the oriented direction amongst the Cartesian coordinates, x , y , and z .
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 ( ϕ i , j ) or depolarization factor ( Л i , j ) 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:
α e f f i , j , L = 2 3 2 3 · α i , j , L , y + α i , j , L , z + 1 3 · α i , j , L , x , d l + + 1 3 · 2 3 · α i , j , L , y , d s + α i , j , L , z , d s + 1 3 · α i , j , L , x +
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.

2.4. 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, α e f f i , j , L , and the slab thickness, δ i . 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 x -axis, the substrate calculations can be presented as follows:
α i , j , L , x / z = α i , j , x / z s u r 1 + α i , x / z s u r 4 3 π · d s c , i 2 · d l c , j · β | | i , j , L
α i , j , L , y = α i , j , y s u r 1 + α i , j , y s u r 4 3 π · d s c , i 2 · d l c , j · β i , j , L
α i , j , x s u r = d s c , i 2 3 · 1 ϕ i , j 2 ϕ i , j 2 1 2 3 · q 11 i · 1 ϕ i , j 2 ϕ i , j 2 1 2 · ε c ε h ε c e 11 i · 1 ϕ i , j 2 ϕ i , j 2 1 2
α i , j , y s u r = 2 · d s c , i 2 3 · 1 ϕ i , j 3 · q 10 i · 1 ϕ i , j 2 ϕ i , j 2 1 2 · ε c ε h ε c e 10 i · 1 ϕ i , j 2 ϕ i , j 2 1 2
α i , j , z s u r = 2 · d s c , i 2 3 · 1 ϕ i , j 3 · q 11 i · 1 ϕ i , j 2 ϕ i , j 2 1 2 · ε c ε h ε c e 11 i · 1 ϕ i , j 2 ϕ i , j 2 1 2
e Λ Μ = p Λ Μ i · 1 ϕ i , j 2 ϕ i , j 2 1 2 · d q Λ Μ i · ψ d ψ ψ = 1 ϕ i , j 2 ϕ i , j 2 1 2 q Λ Μ i · 1 ϕ i , j 2 ϕ i , j 2 1 2 · d p Λ Μ i · ψ d ψ ψ = 1 ϕ i , j 2 ϕ i , j 2 1 2
and distance effect, according to the paper by Royer [21]:
β | | i , j , L = 4 3 π · d s c , i 2 · d l c , j ε s u r + ε h · ε s u r ε h d s c , i + d g 3 + 2 ε o · ε h · 1 3 · u L u L 2 + d s c , i + d g 2 2 u L 3
β i , j , L = 4 3 π · d s c , i 2 · d l c , j ε s u r + ε h · 2 · ε s u r ε h 4 π · ε o · ε w · d s c , i + d g 3 + 2 · ε s u r ε o · ε h · 1 u L 3
where the prolation parameter, ϕ i , j , specifies the selected projection, as highlighted in Figure 2. From Formula (22),   ϕ i , j , x , ϕ i , j , y , z , ϕ i , j , x ,   d l , ϕ i , j , y , z ,   d l , ϕ i , j , y , z ,   d s , and ϕ i , j , x ,   d s are determined as follows:
ϕ i , j , x = 1 d s i d l j 2
ϕ i , j , y , z = 1 0 d l j 2 = 0
ϕ i , j , x ,   d l = ϕ i , j , y , z ,   d s = 1 d l j d s i d l j d s i 2 + d s i 2 2
ϕ i , j , y , z ,   d l = ϕ i , j , x ,   d s = 1 d s i d l j d s i 2 + d s i 2 2
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.

3. Methodology

3.1. Considered System

Figure 4 and Figure 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.

3.2. 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 W i , W j , and W L are reduced and divided into 50 × 500 × 200 parts. The expected values are determined as d s i = 15 ± 3   n m , d l j = 55 ± 8   n m , and u L = 120 ± 30   n m . As a function which describes polydispersity, Gaussian and Lorentzian distributions are selected. The boundaries are fixed as 3 × d s i , 3 × d l j , and 3 × u L . Table 1 summarizes details of the considered nanoparticles, whereas Table 2 indicates electrical properties under 808 nm laser illumination.

3.3. 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 Figure 4 and Figure 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 k - E 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 second-order 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 f q u a 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:
f q u a = f H = 0 + C 1 · H + C 2 · H 3 +
The continuum value f H = 0 for the second-order scheme employed in this work yields:
f H = 0 = f 1 + f 2 f 1 H 2 H 1 θ 1
Similarly, for the three testing steps, Equation (25) is transformed into
f q u a = 2 n · f 1 f 2 2 n 1
n = log 2 f 2 f 3 f 1 f 2
which enables the investigated value to be extrapolated with a specified error, e.g., the relative ϑ and fractional error E 1 between two sets of data that follow
ϑ = f 2 f 1 f 1
E 1 = ϑ H 2 H 1 θ 1
Furthermore, the grid convergence index ( G C I ) 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 G C I parameter:
G C I = F S · E 1 = F S · ϑ H 2 H 1 θ 1
where
  • ϑ —Relative error between two sets/values of data;
  • f H = 0 —Continuum value as if the mesh or time size equals zero;
  • C 1,2 , 3 , —Constant or a function of the selected quantity that is other than H ;
  • F S —Factor of safety, here: 1.25.

3.4. 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:
X d σ d T · g r a d T
where the thermal coefficient of surface tension equals 5.263   · 10 5   N · m 1 · K 1 [31].
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   m W . 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.

4. Results

4.1. 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— t s 1 = 5.0 , t s 2 = 2.5 , and t s 3 = 1.0   m s . 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.

4.2. 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.

4.3. 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 Figure 8 and Figure 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 Figure 4 and Figure 5 were selected to present the results.

4.3.1. Heating Stage

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 Figure 11 and Figure 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 Figure 8 and Figure 9.

4.3.2. Cooling Stage

Similarly, Figure 14, Figure 15, Figure 16 and Figure 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.

4.3.3. 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 represents the maximum velocity fluctuations over time. Results consider the cross-sections discussed in Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16 and Figure 17, located near the AuNRs’ platforms, at the bottom and top sides of the upper glass (4).

4.4. 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 Figure 11 and Figure 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 Figure 12 and Figure 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 Figure 12 and Figure 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 or 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.

5. 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 T 18.5 °C and T 35.5 °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, low-thermal 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.

Author Contributions

Conceptualization, P.Z. and L.D.S.; methodology, P.R.; software, P.R. and A.K.; validation, F.Z., F.P., P.R. and A.K.; formal analysis, P.R., P.Z. and A.K.; investigation, P.R., F.Z., F.P. and P.Z.; data curation, P.R.; writing—original draft preparation, P.R. and P.Z.; writing—review and editing, F.Z., F.P., A.K. and D.M.; visualization, P.R. and F.Z.; supervision, D.M., P.Z., F.P. and L.D.S.; project administration, P.Z. and L.D.S.; funding acquisition, P.Z. and L.D.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research of Gdańsk University of Technology was funded in whole by the National Science Center in Poland under the project “Shape and displacement optimization of gold nanorods in the killing chamber in order to phototermoablation processes”, no. UMO-2021/43/D/ST8/02504. For the purpose of open access, the author has applied a CC-BY public copyright license to any Author Accepted Manuscript (AAM) version arising from this submission. Italian work in this article has been supported: by the “NATO—Science for Peace and Security Programme (SPS-G5759, NANO-LC)”; by the Air Force Office of Scientific Research, Air Force Material Command, U.S. Air Force. “Digital optical network encryption with liquid-crystal grating metasurface perfect absorbers” FA8655-22-1-7007 (P. I. L. De Sio, EOARD 2022–2025); by the Italian PON project TITAN “Nanotechnology for cancer immunotherapy”, 2021–2023 ARS01_00906.

Data Availability Statement

The data presented in this study are available on request from the corresponding author ongoing doctoral dissertations in this area. However, some of the metadata are available under the DOI https://doi.org/10.34808/8rbr-1c71 (accessed on 7 May 2024).

Acknowledgments

The first author of the paper is grateful to the Doctoral School at Gdańsk University of Technology for a scholarship.

Conflicts of Interest

The authors declare no conflicts of interest.

Nomenclature

Roman letters
A a b s i , j , L absorption coefficient of nanoparticles, m 1
A a b s M absorption coefficient of the continuous material, m 1
C constant
C a b s i , j , L absorption cross-section of the i -particle, m 2
C e x t i , j , L extinction cross-section of the i -particle, m 2
C s c a i , j , L scattered cross-section of the i -particle, m 2
c p e f f effective heat capacity, J · K 1 · k g 1
c p specific heat capacity, J · K 1 · k g 1
d symmetric rate of deformation, s 1
d l c , j longer diameter of the i -nanoparticle, m
d g size of the nanoparticle’s coating, m
d s c , i shorter diameter of the i -nanoparticle, m
e specific energy, J · k g 1
F S factor of safety, here: 1.25
f selected quantity or testing parameter, here: temperature, K
G k turbulence kinetic energy generation due to velocity gradients, k g · m 1 · s 3
G Б turbulence kinetic energy generation due to buoyancy forces, k g · m 1 · s 3
g gravity, m · s 2
H numerical parameters to be compared (here: mesh size or timestep)
h heat transfer coefficient, W · m 2 · K 1
I o M a x laser intensity profile at its maximum
I o r initial laser intensity profile, W · m 2
I a b s i , j , L absorbed part of the i -particle, W · m 2
k e f f effective thermal conductivity coefficient, W · m 1 · K 1
n M refractive index of selected material, -
n numerical constant for 3- H correction, -
P r k turbulent Prandtl number for kinetic energy, -
P r ϵ turbulent Prandtl number for rate of dissipation, -
P r Prandtl number, -
p pressure, P a
R i , j , L reflection coefficient of nanoparticles, -
R M reflection coefficient of the continuous material, -
R B beam size of laser source, here: 0.0015 m
r radius in spherical coordinates, m
S e f source of energy for fluids, W · m 3
S e n p source of energy for the i -particle, W · m 3
T o initial temperature at the t = 0   s , K
T temperature, K
t time, s
u L distance between particles from their boundaries, m
v d speed of sound, m · s 1
v velocity of the fluid, m · s 1
w i size distribution (polydispersity), according to the Gaussian distribution, -
w j length distribution (polydispersity), according to the Gaussian distribution, -
w L distance distribution, according to the Lorentzian distribution, -
X diffusive momentum flux, k g · m · s 2
x , y , z Cartesian’s coordinates, m
Greek letters
α e f f i , j , L effective polarizability of the i -particle, m 3
α i , j , L , i ^ polarizability of the i -particle, m 3
β anisotropic factor for distance and surface effects,
δ i gold slab thickness, here: size of capped nanoparticles, m
δ thickness of the considered material, m
ε o permittivity in vacuum, -
ε A relative permittivity of the considered metal, -
ε c relative permittivity of the core, here: ε c = ε A ,
ε g relative permittivity of the capping agent,
ε h relative permittivity of the working (host) medium, -
ε s u r relative permittivity of the base where NRs are deposited, -
ϵ relative emissivity, -
θ rank of accuracy, here: θ = 2
ϑ relative error, -
λ incident wavelength, m
μ molecular viscosity, P a · s
ξ i , j , L nanoparticles concentration, m 3
ρ density of a material, k g · m 3
σ surface tension parameter, N · m 1
ϕ i , j elongation parameter, -
Slavic letters
Б thermal expansion coefficient, K 1
Γ relaxation time, s
Others
symbol of partial derivative
E turbulent energy dissipation, m 2 · s 3
E fractional error
e symmetrical rate of prolate rod–surface interaction
d first invariant of the strain rate, s 1
unit tensor,
i symbol of imaginary unit
i m imaginary part of a complex expression
k turbulent kinetic energy
N number of NP-NP pairs
p associated Legendre polynomials of a first kind
q associated Legendre polynomials of a second kind
S total specific entropy, J · k g 1 · K 1
r e real part of a complex expression
W i boundary for size distribution, -
W j boundary for length distribution, -
W L boundary for distance distribution
Subscripts and superscripts
A aurum, gold
a air
c core
f fluid
gagent (shell)
i size distribution
i ^ selected axis
j length distribution
L distance distribution between nanoparticles
l long
M material
o in reference to initial conditions
S safety
s short
s step
t turbulent
Λ , Μ integers of Legendre polynomials’ association
transposition
| | parallel
perpendicular
        ̿ average value
+ oriented axis
unoriented axis
Abbreviations
a b s absorption
a v average
CFDcomputational fluid dynamics
CPUcentral processing unit
d i v divergence
e f f effective
e x p exponent
e x t extinction
FVMfinite volume method
G C I grid convergence index
g l glass
g r a d gradient
m a x maximum
m i n minimum
n p nanoparticle
NRsnanorods
PELpolyelectrolyte
p h phonon
q u a quantified
s c a scattering
s u r surface
Y a in reference to the Yamaguchi approach

References

  1. Coyle, E.D.; Simmons, R.A. Understanding the Global Energy Crisis; Purdue University Press: West Lafayette, IN, USA, 2014. [Google Scholar]
  2. Huang, Y.; Lee, C.K.C.; Yam, Y.-S.; Mok, W.-C.; Zhou, J.L.; Zhuang, Y.; Surawski, N.C.; Organ, B.; Chan, E.F.C. Rapid detection of high-emitting vehicles by on-road remote sensing technology improves urban air quality. Sci. Adv. 2022, 8, eabl7575. [Google Scholar] [CrossRef] [PubMed]
  3. Petronella, F.; De Biase, D.; Zaccagnini, F.; Verrina, V.; Lim, S.; Jeong, K.; Miglietta, S.; Petrozza, V.; Scognamiglio, V.; Godman, N.P.; et al. Label-free and reusable antibody-functionalized gold nanorod arrays for the rapid detection of Escherichia coli cells in a water dispersion. Environ. Sci. Nano 2022, 9, 3343–3360. [Google Scholar] [CrossRef]
  4. Zaccagnini, F.; Radomski, P.; Sforza, M.L.; Ziółkowski, P.; Lim, S.I.; Jeong, K.U.; Mikielewicz, D.; Godman, N.P.; Evans, D.R.; Slagle, J.E.; et al. White light termoplasmonic activated gold nanorod arrays enable the photo-thermal disinfection of medical tools from bacterial contamination. J. Mater. Chem. B 2023, 11, 6823–6836. [Google Scholar] [CrossRef] [PubMed]
  5. Fang, X.; Wang, Y.; Wang, Z.; Jiang, Z.; Dong, M. Microorganism Assisted Synthesized Nanoparticles for Catalytic Applications. Energies 2019, 12, 190. [Google Scholar] [CrossRef]
  6. Pontico, M.; Conte, M.; Petronella, F.; Frantellizzi, V.; De Feo, M.S.; Di Luzio, D.; Pani, R.; De Vincentis, G.; De Sio, L. 18F-fluorodeoxyglucose (18F-FDG) Functionalized Gold Nanoparticles (GNPs) for Plasmonic Photothermal Ablation of Cancer: A Review. Pharmaceutics 2023, 15, 319. [Google Scholar] [CrossRef] [PubMed]
  7. Radomski, P.; Ziolkowski, P.; Mikielewicz, D. Energy conversion in systems-contained laser irradiated metallic nanoparticles—Comparison of results from analytical solutions and numerical methods. Acta Mech. Autom. 2023, 17, 540–549. [Google Scholar] [CrossRef]
  8. Tofani, K.; Tiari, S. Nano-Enhanced Phase Change Materials in Latent Heat Thermal Energy Storage Systems: A Review. Energies 2021, 14, 3821. [Google Scholar] [CrossRef]
  9. Mebared-Oudina, F.; Chabani, I. Revie on Nano Enhanced PCMs: Insight on nePCM Application in Thermal Management/Storage Systems. Energies 2023, 16, 1066. [Google Scholar] [CrossRef]
  10. Radomski, P.; Ziolkowski, P.; De Sio, L.; Mikielewicz, D. Computational fluid dynamics simulation of heat transfer from densely packed gold nanoparticles to isotropic media. Arch. Thermodyn. 2021, 42, 87–114. [Google Scholar]
  11. Mikielewicz, D. Hydrodynamics and heat transfer in bubbly flow in the turbulent boundary layer. Int. J. Heat Mass Transf. 2002, 46, 207–220. [Google Scholar] [CrossRef]
  12. Ziolkowski, P.; Badur, J. A theoretical, numerical and experimental verification of the Reynolds thermal transpiration law. Int. J. Numer. Methods Heat Fluid Flow 2018, 28, 64–80. [Google Scholar] [CrossRef]
  13. De Sio, L.; Placido, T.; Comparelli, R.; Curri, M.L.; Striccoli, M.; Tabiryan, N.; Bunning, T.J. Next-generation thermo-plasmonic technologies and plasmonic nanoparticles in optoelectronics. Prog. Quantum Electron. 2015, 41, 23–70. [Google Scholar] [CrossRef]
  14. Sobhan, C.B.; Peterson, G.P. Microscale and Nanoscale Heat Transfer: Fundamentals and Engineering Applications; CRC Press: Boca Raton, FL, USA; Taylor & Francis Group: Abingdon, UK, 2008. [Google Scholar]
  15. Cattaneo, M.C. A form of heat conduction equation which eliminates the paradox of instantaneous propagation. Compte Rendus 1958, 247, 431–433. [Google Scholar]
  16. Radomski, P.; Zaccagnini, F.; Ziolkowski, P.; Petronella, F.; De Sio, L.; Mikielewicz, D. Heat Transfer of the Multicolor-Laser-Sources-Irradiated Nanoparticles in Reference to Thermal Processes. In Proceedings of the 36th International Conference on Efficiency, Cost, Optimization, Simulation and Environmental Impact on Energy Systems (ECOS), Las Palmas de Gran Canaria, Spain, 25–30 June 2023. [Google Scholar]
  17. Bohren, C.F.; Huffman, D.R. Absorption and Scattering of Light by Small Particles; A Wiley-Interscience Publication: Toronto, ON, Canada, 1998. [Google Scholar]
  18. Strutt, H.J.W. On the scattering of light by small particles. Lond. Edinb. Dublin Philos. Mag. J. Sci. 1871, 41, 447–454. [Google Scholar] [CrossRef]
  19. Smoluchowski, M. On conduction of heat by rarefied gases. Ann. Phys. Chem. 1898, 64, 101–130. (In German) [Google Scholar]
  20. Kreibig, U.; Vollmer, M. Optical Properties of Metal Clusters; Springer Series in Material Science; Springer: New York, NY, USA, 1995. [Google Scholar]
  21. Yamaguchi, T.; Yoshida, S.; Kinbara, A. Optical Effect of the Substrate on the Anomalous Absorption of Aggregated Silver Films. Thin Solid Film. 1974, 21, 173–187. [Google Scholar] [CrossRef]
  22. Fuchs, R. Theory of the optical properties of ionic crystal cubes. Phys. Rev. B 1975, 11, 1732–1740. [Google Scholar] [CrossRef]
  23. Roache, P.J. Quantification of Uncertainty in Computational Fluid Dynamics. Annu. Rev. Fluid Mech. 1997, 29, 123–160. [Google Scholar] [CrossRef]
  24. Richardson, L.F. The approximate arithmetical solution by finite differences of physical problems including differential equations, with an application to the stresses in a masonry dam. Phylosophical Trans. R. Soc. A 1911, 210, 307–357. [Google Scholar]
  25. Domański, R. Laser Radiation—Solid-State Interaction; Polish Scientific-Technical Journals: Warsaw, Poland, 1990; pp. 130–220. ISBN 83-204-1383-4. (In Polish) [Google Scholar]
  26. Carslaw, H.S.; Jaeger, J.C. Conduction of Heat in Solids, 2nd ed.; Oxford University Press: London, UK, 1959; pp. 50–132. [Google Scholar]
  27. Romanchuk, B.J. Computational Modeling of Buble Growth Dynamics in Nucleate Pool Boiling for Pure Water and Aqueous Surfactant Solutions. Master’s Thesis, Division of Research and Advanced of the University of Cincinnati, Cincinnati, OH, USA, 2014. [Google Scholar]
  28. Darbari, B.; Rashidi, S.; Esfahani, J.A. Sensitivity Analysis of Entropy Generation in Nanofluid Flow inside a Channel by Response Surface Methodology. Entropy 2016, 18, 52. [Google Scholar] [CrossRef]
  29. Nisar, K.S.; Khan, D.; Khan, A.; Khan, W.A.; Khan, I.; Aldawsari, A.M. Entropy Generation and Heat Transfer in Drilling Nanoliquids with Clay Nanoparticles. Entropy 2019, 21, 1226. [Google Scholar] [CrossRef]
  30. Zhang, L.; Bhatti, M.M.; Marin, M.; Mekheimer, K.S. Entropy Analysis on the Blood Flow through Anisotropically Tapered Arteries Filled with Magnetic Zinc-Oxide (ZnO) Nanoparticles. Entropy 2020, 22, 1070. [Google Scholar] [CrossRef] [PubMed]
  31. Ziolkowski, P.; Badur, J. On Navier slip and Reynolds transpiration numbers. Arch. Mech. 2018, 70, 269–300. [Google Scholar]
  32. Vial, A.; Laroche, T. Description of dispersion properties of metals by means of the critical points model and application to the study of resonant structures using the FDTD method. J. Phys. D Appl. Phys. 2007, 40, 7152–7158. [Google Scholar] [CrossRef]
  33. Bengte, D.L. The Refractive Index of Air. Metrologia 1966, 8, 71–80. [Google Scholar]
  34. Bansal, N.P.; Doremus, R.H. Handbook of Glass Properties; Academic Press: Cambridge, MA, USA; Materials Engineering Department Rensselaer Polytechnic Institute: Troy, NY, USA, 1986. [Google Scholar]
  35. Zaitlin MP Anderson, A.C. Thermal Conductivity of Borosilicate Glass. Phys. Rev. Lett. 1974, 33, 1158–1161. [Google Scholar] [CrossRef]
  36. Mark, J.E. Polymer Data Handbook; Oxford University Press: London, UK, 1999; Volume 131, p. 44. [Google Scholar]
  37. Reddy, H. Temperature-dependent optical properties of gold thin films. Opt. Mater. Express 2016, 6, 2776–2802. [Google Scholar] [CrossRef]
  38. Siegel, R.; Howell, J.R. Thermal Radiation Heat Transfer; McGraw-Hill Book Company: New York, NY, USA, 1972. [Google Scholar]
  39. Engineering ToolBox. Air—Density, Specific Weight and Thermal Expansion Coefficient vs. Temperature and Pressure. 2003. Available online: https://www.engineeringtoolbox.com/air-density-specific-weight-d_600.html (accessed on 17 February 2023).
  40. Engineering ToolBox. Air—Thermal Conductivity vs. Temperature and Pressure. 2009. Available online: https://www.engineeringtoolbox.com/air-properties-viscosity-conductivity-heat-capacity-d_1509.html (accessed on 17 February 2023).
  41. Engineering ToolBox. Air—Dynamic and Kinematic Viscosity. 2003. Available online: https://www.engineeringtoolbox.com/air-absolute-kinematic-viscosity-d_601.html (accessed on 17 February 2023).
  42. Fernández-Prini, R. Release on the Refractive Index of Ordinary Water Substance as a Function of Wavelength, Temperature and Pressure. J. Phys. Chem. Ref. Data 1998, 27, 761. [Google Scholar]
  43. Koulali, A.; Ziółkowski, P.; Radomski, P.; De Sio, L.; Zieliński, J.; Mikielewicz, D. Single-Phase CFD Approach for Investingating Bacterial Inactivation and Heat Transfer in a Microchamber; Book of extended papers from XXV Jubilee Congress of Thermodynamics in Gdańsk (Poland), 11–14 September 2023; Gdańsk University of Technology Redaction: Gdańsk, Poland, 2023; pp. 167–174. [Google Scholar]
  44. Koulali, A.; Abderrahmane, A.; Jamshed, W.; Hussain, S.M.; Nisar, K.S.; Abdel-Aty, A.H.; Yahia, I.S.; Eid, M.R. Comparative study on effects of thermal gradient direction on heat exchange between a pure fluid and a nanofluid: Employing finite volume method. Coatings 2021, 11, 1481. [Google Scholar] [CrossRef]
  45. Piotrowicz, M.; Flaszynski, P.; Doerffer, P. Effect of Hot Spot Location on Flow Structure in Nozzle Guide Vane. IOP Conf. Ser. J. Phys. Conf. Ser. 2018, 1101, 012025. [Google Scholar] [CrossRef]
  46. Patthew, P.; Garnett, B. Introduction to Metal-Nanoparticle Plasmonics; A Wiley–ScienceWise Publishing Co-Publication: Hoboken, NJ, USA, 2013. [Google Scholar]
  47. Zhengdao, W.; Wei, Y.; Qian, Y. Numerical study on entropy generation in thermal convection with differentially discrete heat boundary conditions. Entropy 2018, 20, 351. [Google Scholar] [CrossRef] [PubMed]
  48. Al-Kouz, W. MHD darcy-forchheimer nanofluid flow and entropy optimization in an odd-shaped enclosure filled with a (MWCNT-Fe3O4/water) using galerkin finite element analysis. Sci. Rep. 2021, 11, 2263. [Google Scholar] [CrossRef] [PubMed]
  49. Drożyński, Z. Entropy increase as a measure of energy degradation in heat transfer. Arch. Thermodyn. 2013, 34, 147–160. [Google Scholar] [CrossRef]
  50. Bejan, A. Fundamentals of Exergy Analysis, Entropy Generation Minimization, and the Generation of Flow Architecture. Int. J. Energy Res. 2002, 26, 1–43. [Google Scholar] [CrossRef]
  51. Sciubba, E. Use of Exergy Analysis To Compute The Resource Intensity of Biological Systems and Human Societies. In Proceedings of the 12th Joint European Thermodynamics Conference, Brescia, Italy, 1–5 July 2013; pp. 268–273. [Google Scholar]
  52. Jou, D. Entropy, entropy flux, temperature, and second law in extended irreversible thermodynamics. In Proceedings of the 12th Joint European Thermodynamics Conference, Brescia, Italy, 1–5 July 2013; pp. 211–216. [Google Scholar]
  53. Flaszynski, P.; Wasilczuk, F.; Piotrowicz, M.; Telega, J.; Mitraszewski, K.; Hansen, K.S. Numerical simulations for a parametric study of blockage effect on offshore wind farms. Wind Energy 2023, 27, 53–74. [Google Scholar] [CrossRef]
  54. Radomski, P.; De Sio, L.; Ziółkowski, P.; Mikielewicz, D. Heat Transfer in a Metallic-Nanoparticles-Contained System Using the Discrete-Ordinates and Thin-Film-Adjusted Models; Book of papers from XXV Jubilee Congress of Thermodynamics in Gdańsk (Poland), 11–14 September 2023; Gdańsk University of Technology Redaction: Gdańsk, Poland, 2023; p. 91. [Google Scholar]
  55. Zuo, W.; Li, F.; Li, Q.; Chen, Z.; Huang, Y.; Chu, H. Multi-objective optimization of micro planar combustor with tube outlet by RSM and NSGA-II for thermophotovoltaic applications. Energy 2024, 291, 130396. [Google Scholar] [CrossRef]
Figure 1. Scheme of energy conversion calculations for (a) thin films and (b) nanoparticles.
Figure 1. Scheme of energy conversion calculations for (a) thin films and (b) nanoparticles.
Energies 17 02587 g001
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 i ^ axis are marked by a violet color.
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 i ^ axis are marked by a violet color.
Energies 17 02587 g002
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 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’.
Energies 17 02587 g003
Figure 4. Scheme and dimensions of the gold nanorod platform from side view (scale has not been saved).
Figure 4. Scheme and dimensions of the gold nanorod platform from side view (scale has not been saved).
Energies 17 02587 g004
Figure 5. Scheme and dimensions of the gold nanorod platform from top view (scale has not been saved).
Figure 5. Scheme and dimensions of the gold nanorod platform from top view (scale has not been saved).
Energies 17 02587 g005
Figure 6. Three discretized meshes employed in computations on both sides. Mesh at the bottom has been assigned to demonstrate theoretical results.
Figure 6. Three discretized meshes employed in computations on both sides. Mesh at the bottom has been assigned to demonstrate theoretical results.
Energies 17 02587 g006
Figure 7. Experimental validation of the absorption spectra using a spectrophotometer for three different working media, ε h 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 7. Experimental validation of the absorption spectra using a spectrophotometer for three different working media, ε h 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.
Energies 17 02587 g007
Figure 8. Experimental validation of (a) the temperature distribution and (b) time response from a thermal camera for a 1 W laser approach.
Figure 8. Experimental validation of (a) the temperature distribution and (b) time response from a thermal camera for a 1 W laser approach.
Energies 17 02587 g008
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 9. Experimental validation from a thermal camera of (a) temperature distributions and (b) time response for a white-light approach based on [4].
Energies 17 02587 g009
Figure 10. Temperature contours, according to Figure 4 and Figure 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 10. Temperature contours, according to Figure 4 and Figure 5. A distance of 5 μm under (a), and 5 μm above (b) the AuNR platforms (5) after 70 s of the laser illumination.
Energies 17 02587 g010
Figure 11. Static entropy contours, according to Figure 4 and Figure 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.
Figure 11. Static entropy contours, according to Figure 4 and Figure 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.
Energies 17 02587 g011
Figure 12. Dynamic entropy contours, according to Figure 4 and Figure 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.
Figure 12. Dynamic entropy contours, according to Figure 4 and Figure 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.
Energies 17 02587 g012
Figure 13. Velocity vectors, according to Figure 4 and Figure 5. A distance of 5 μm under (a), and 5 μm above (b) the AuNR platforms (5) after 70 s of laser illumination.
Figure 13. Velocity vectors, according to Figure 4 and Figure 5. A distance of 5 μm under (a), and 5 μm above (b) the AuNR platforms (5) after 70 s of laser illumination.
Energies 17 02587 g013
Figure 14. Temperature contours, according to Figure 4 and Figure 5. A distance of 5 μm under (a), and 5 μm above (b) the AuNR platforms (5) after 10 seconds of cooling.
Figure 14. Temperature contours, according to Figure 4 and Figure 5. A distance of 5 μm under (a), and 5 μm above (b) the AuNR platforms (5) after 10 seconds of cooling.
Energies 17 02587 g014
Figure 15. Static entropy, according to Figure 4 and Figure 5. A distance of 5 μm under (a), and 5 μm above (b) the AuNR platforms (5) after 10 s of cooling. Gray color indicates places where glue has been dislocated.
Figure 15. Static entropy, according to Figure 4 and Figure 5. A distance of 5 μm under (a), and 5 μm above (b) the AuNR platforms (5) after 10 s of cooling. Gray color indicates places where glue has been dislocated.
Energies 17 02587 g015
Figure 16. Dynamic entropy, according to Figure 4 and Figure 5. A distance of 5 μm under (a), and 5 μm above (b) the AuNR platforms (5) after 10 s of cooling. Gray color indicates places where glue has been dislocated.
Figure 16. Dynamic entropy, according to Figure 4 and Figure 5. A distance of 5 μm under (a), and 5 μm above (b) the AuNR platforms (5) after 10 s of cooling. Gray color indicates places where glue has been dislocated.
Energies 17 02587 g016
Figure 17. Velocity vectors, according to Figure 4 and Figure 5. A distance of 5 μm under (a), and 5 μm above (b) the AuNR platforms (5) after 10 s of cooling.
Figure 17. Velocity vectors, according to Figure 4 and Figure 5. A distance of 5 μm under (a), and 5 μm above (b) the AuNR platforms (5) after 10 s of cooling.
Energies 17 02587 g017
Figure 18. Time responses of (a) maximum, (b) average temperature, and (c) magnitude of velocity vectors 5 μm above (top plane) (red color) and 5 μm under (bottom plane) (blue color) the glassy domain (4) of the AuNR platforms, respectively.
Figure 18. Time responses of (a) maximum, (b) average temperature, and (c) magnitude of velocity vectors 5 μm above (top plane) (red color) and 5 μm under (bottom plane) (blue color) the glassy domain (4) of the AuNR platforms, respectively.
Energies 17 02587 g018
Table 1. Parameters of gold nanorods.
Table 1. Parameters of gold nanorods.
General ParameterAssumed Value
Expected size of nanorods, d s 15 ± 3   n m
Expected length of nanorods, d l 55 ± 8   n m
Interdistance,   u L 120 ± 30   n m
Capping agent thickness, d g 4 nm
Nanoparticle concentration, ξ ( 1.63 ± 0.13 ) · 10 21 m 3
Table 2. Electrical properties for the utilized materials.
Table 2. Electrical properties for the utilized materials.
Permittivities λ = 808   n m References
Air, ε a 1.011639[29]
Glass, ε b g 2.281912[30,31]
Glue, ε g l 2.400680[32]
Polyelectrolyte, ε s u r 2.274245[4]
Gold nanoparticles r e −25.035329[33,34]
i m 1.564609
Table 3. Material properties in the considered simulations.
Table 3. Material properties in the considered simulations.
Material Density ,   ρ
k g · m 3
Specific   Heat   Capacity ,   c p e f f
J · k g 1 · K 1
Thermal
Conductivity
Coefficient, k e f f
W · m 1 · K 1
Dynamic
Viscosity, μ
P a · s
References
Glass 2124.9 779.74 0.9245251
+ 0.0004777689 · T
+ 9.178795 · 10 7 · T 2
-[33,34,35]
Glue 1231 1100 0.50 0.30 [36]
Gold nanorods19,320 129.81 317 -[37,38]
Air 352.965 · T 1 1006.43 0.004204762
+ 7.242857 · 10 5 · T
7.29 · 10 6
+ 4.0 · 10 8 · T
[39,40,41,42]
Table 4. Specified boundary conditions in the simulations carried out.
Table 4. Specified boundary conditions in the simulations carried out.
No.DomainBoundary Conditions
(1)Bottom glass T o h ϵ S e f 297.15   K 13.6   W · m 2 · K 1 ϵ = 0.957 0   W · m 3
(2)Air L a i r 0.01   m   ( 0.3937   i n ) (shell conduction)
S e f 0   W · m 3
v o 1   μ m · s 1 (at inlet)
(3)Glue T o 297.15   K
(4)Top glass S e f T o h ϵ 1.07516 · 10 5   W · m 3 297.15   K 13.6   W · m 2 · K 1 0.957
(5)Gold slabs δ ̿ i 26.0   n m
d σ d T 0.00005263   N · m 1 · K 1
S e n p U p 4.71214 · 10 11   W · m 3
S e n p D o w n 4.43412 · 10 11   W · m 3
System T o 297.15   K
p o 1013.25 hPa
Table 5. Details of the applied meshes.
Table 5. Details of the applied meshes.
ParameterMesh 1Mesh 2Mesh 3
Mean mesh size 0.1000   m m 0.0800   m m 0.0625   m m
Maximum aspect ratio 27.4834 30.7582 31.3168
Mean volume cell 6.2966 · 10 12   m 3 5.6025 · 10 12   m 3 3.4022 · 10 12   m 3
Orthogonal quality 1.0 1.0 1.0
Number of cells 4,985,344 8,324,230 12,473,784
Table 6. Parameters under Richardson and Roacher extrapolation proceedings after 30 s of illumination by a 1 W laser.
Table 6. Parameters under Richardson and Roacher extrapolation proceedings after 30 s of illumination by a 1 W laser.
Extrapolated ErrorMesh 1–Mesh 2Mesh 2–Mesh 3 t s 1 t s 2 t s 2 t s 3
ϑ T m a x 0.14526 % 0.05424 % 0.007190 % 0.002261 %
T a v 7.69209 · 10 5 % 3.536 · 10 5 % 0.001117 % 0.002921 %
| v m a x | 0.7104 % 0.2263 % 0.2265 % 0.1396 %
f H = 0 T m a x 51.5538 °C51.8365 °C52.2260 °C52.1884 °C
T a v 25.7018 °C25.7018 °C25.7018 °C25.7046 °C
| v m a x | 0.44305   c m · s 1 0.43170   c m · s 1 0.43267   c m · s 1 0.43268   c m · s 1
G C I ( f i n e ) T m a x 0.3721 % 0.1296 % 0.01027 % 0.003020 %
T a v 1.970 · 10 6 % 8.448 · 10 6 % 0.001596 % 0.003901 %
| v m a x | 0.01820 % 0.5407 % 0.3236 % 0.5290 %
f q u a T m a x 51.8277 °C51.7884 °C
T a v 25.7018 °C25.7018 °C
| v m a x | 0.43214   c m · s 1 0.43322   c m · s 1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Radomski, P.; Zaccagnini, F.; Ziółkowski, P.; Petronella, F.; De Sio, L.; Koulali, A.; Mikielewicz, D. Investigations of Energy Conversion and Surface Effect for Laser-Illuminated Gold Nanorod Platforms. Energies 2024, 17, 2587. https://doi.org/10.3390/en17112587

AMA Style

Radomski P, Zaccagnini F, Ziółkowski P, Petronella F, De Sio L, Koulali A, Mikielewicz D. Investigations of Energy Conversion and Surface Effect for Laser-Illuminated Gold Nanorod Platforms. Energies. 2024; 17(11):2587. https://doi.org/10.3390/en17112587

Chicago/Turabian Style

Radomski, Piotr, Federica Zaccagnini, Paweł Ziółkowski, Francesca Petronella, Luciano De Sio, Aimad Koulali, and Dariusz Mikielewicz. 2024. "Investigations of Energy Conversion and Surface Effect for Laser-Illuminated Gold Nanorod Platforms" Energies 17, no. 11: 2587. https://doi.org/10.3390/en17112587

APA Style

Radomski, P., Zaccagnini, F., Ziółkowski, P., Petronella, F., De Sio, L., Koulali, A., & Mikielewicz, D. (2024). Investigations of Energy Conversion and Surface Effect for Laser-Illuminated Gold Nanorod Platforms. Energies, 17(11), 2587. https://doi.org/10.3390/en17112587

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop