Recent Development in Partial Differential Equations and Their ApplicationsView this Special Issue
Numerical Modeling of the Photothermal Processing for Bubble Forming around Nanowire in a Liquid
An accurate computation of the temperature is an important factor in determining the shape of a bubble around a nanowire immersed in a liquid. The study of the physical phenomenon consists in solving a photothermic coupled problem between light and nanowire. The numerical multiphysic model is used to study the variations of the temperature and the shape of the created bubble by illumination of the nanowire. The optimization process, including an adaptive remeshing scheme, is used to solve the problem through a finite element method. The study of the shape evolution of the bubble is made taking into account the physical and geometrical parameters of the nanowire. The relation between the sizes and shapes of the bubble and nanowire is deduced.
In the last years, many researchers are interested in the use of nanomaterials. In the chemical industry and in the manufactures of nanotubes and nanowires, the usually used materials are and [1–3]. Such a use of these nanomaterials (natural or artificial) increases and these are dispersed in air or in water . Their impact on the environment and health must be evaluated (e.g., toxicity analysis) . Therefore, the detection of the presence of such nanomaterials in the environment becomes crucial. Two modes of detections of such a nanowire/nanotube can be achieved. The first one consists in a direct detection of the nanowire by optical microscopy through the measurement of the scattering of light emitted by the nanomaterial. Due to a weak signal/noise ratio, such a detection mode can be difficult. The second mode is an indirect method and consists in studying the bubble created by the photothermal response of the nanowire immersed in a liquid and illuminated by an electromagnetic wave. In such an approach, the nanowire absorbs the electromagnetic radiation (energy) for a range of wavelengths and heats and, for temperature exceeding the threshold of vaporization of the liquid, induces the creation of a nanobubble [6, 7]. The created bubble grows before being detected. The analysis of the shape and size of the bubble should permit studying the morphology of the nanowire. The studied problem consists in solving a photothermic coupled system (light, nanowire, and heat) taking into account the physical parameters of the system (i.e., permittivity of materials, material conductivity, laser wavelength, and laser power).
In that context, a numerical multiphysic model, allowing studying the behavior of the nanowires illuminated by an incident laser field, is presented. The formation of the bubble, associated with a nanowire of immersed in water and illuminated by a laser pulse, is studied. An optimization process, including adaptive remeshing scheme, is used to detect the variations of the temperature, the bubble shape evolution, ensuring the convergence of the solution to the physical solution [8, 9]. The paper is organized as follows. Section 2 describes the equations of the model and the numerical resolution method. The adaptive remeshing process and the optimization steps are presented in Section 3. In Section 4, the results of numerical simulations are presented before concluding.
2. Model and Numerical Methods
The section is devoted to presenting the equation systems modeling the photothermic process and the numerical method used to solve the system.
2.1. Electromagnetic Problem
In electromagnetic system, the partial differential equations are derived from Maxwell's equations. The problem can be reduced to Helmholtz equation for the harmonic electric and magnetic fields (i.e., in the form , where is the angular frequency of the harmonic wave) . In the 2D case of an infinity elliptical cylinder along the -axis, the unknown field is the magnetic component and, for a polarized illumination in the transverse magnetic mode TM, the magnetic field can be written as . Therefore, the electromagnetic problem is reduced to a scalar problem and the computation of , in a domain , allows deducing the electric field by using the Maxwell-Ampere equation : where is the rotational operator, the angular frequency, the permittivity of vacuum, and the relative complex permittivity of the considered materials which are functions of the spatial coordinates . The component of the magnetic field satisfies the scalar equation where is the wave number of the monochromatic incoming wave and the speed of light in vacuum. To compute the solution of the electromagnetic problem, a set of conditions on the boundary of the computational domain must be imposed. The natural boundary condition at the interface between materials is the continuity of the normal component of the electromagnetic excitation: where is the normal derivative operator, is the normal vector component along the -axis, and is the incident illumination field along the -axis with and being the permeability of vacuum. Such a boundary condition is used in problems of wave propagation [11–13].
2.2. Thermic Problem
Under illumination by an electromagnetic wave, the nanowire absorbs energy. That energy produces a heat source given by The resolution of the thermal problem requires solving the heat equation which is a partial differential parabolic equation describing the evolution of the temperature with a heat source . That equation is written as follows: with a Dirichlet boundary condition and is the thermal conductivity of the materials. The variation of the temperature depends on both imaginary part of the permittivity and the intensity of the electric field .
The resolution of the coupled electromagnetic and heat problems allows extracting the spatial distribution of the temperature in the computational domain. From the map of temperature and for a fixed threshold of vaporization , the identification of the shape and size of the bubble around the nanowire can be achieved. Such information on shape and size of the bubble would be used to construct a relation between the geometric characteristics of the bubble and the nanowire.
2.3. The Finite Element Method
The objective is to solve (2) and (5) for the coupled system in a domain whose geometry can be complex. The Finite Element Method (FEM) was applied since the 1940s in mechanics, thermodynamics, electromagnetics, and electrical engineering [14, 15]. The method is used to solve partial differential equation systems with boundary conditions in open or close domains. The resolution of problem necessitates a discrete domain, generally named mesh of the domain . The solutions of the problem are computed on the nodes of the mesh. In order to both control the error on the solution and to decrease the number of nodes, an improved method, including an iterative remeshing process, is developed and used. Such an improved FEM allows describing the complex structures with arbitrary shapes. Moreover, the stability of the FEM is also improved by using a weak formulation (or variational formulation) of (2) and (5). Therefore, the electromagnetic and thermic fields satisfy where and are test functions defined on (the linear space of the scalar functions and , being 2-integrable on ). The basis of polynomial functions provides an approximation of the solutions and in each node . The field (resp., ) is a linear combination of such basic polynomial functions (resp., ) and the problem consists in solving a linear system [15, 17]. The solution verifies exactly the partial differential equations on each node for the given boundary conditions. Ritz's formulation of the variational problem is used to satisfy the continuity of the tangential components of the electromagnetic field .
3. Optimization Process and Adaptive Remeshing
Partial differential equations (electromagnetic and thermic) are formulated and solved on the mesh of the computational domain through the FEM. But the accuracy of the computed solution depends on the quality of the mesh [9, 18, 19]. A remeshing process and adaptive loops have been developed in order to improve the quality of the solutions by adapting the size of the mesh elements to the physical solution [9, 20]. The mesh adaption is required to converge to a stable solution, in particular where strong variation of the electromagnetic or temperature fields occurred. For each step of the adaption process, the approximate solution of the Helmholtz equation, the electric field , the heat source , and the temperature are computed . The interpolation error, based on an estimation of the discrete Hessian of the solution, is used to limit the maximum deviation between the exact solution and the solution associated with the mesh [21, 22]. The a posteriori error estimator, based on the interpolation error, allows defining a physical size map such as where is the physical size defined at each node and is proportional to the inverse of the deviation of the Hessian. For a given maximum tolerance on the physical error , the size is given by where and are the minimum and maximum sizes of the elements and is an estimation of the maximum deviation obtained from the Hessian of the solution. The physical size map is used to govern the adaptive remeshing of the domain with the BL2D-V2 software (adaptive remeshing generating isotropic or anisotropic meshes) . The domain is then entirely remeshed and a new mesh is obtained. The resolution of the multiphysics problem is based on the computation of two physical size maps: the first one related to the heat source and the second one related to the temperature . The adaptive computational scheme consists in iterative and adaptive loops: initial mesh generated with triangular elements of the computational domain ,computation of the field (solution of (2)) on ,computation of the solutions and on ,physical error estimation: computation of the interpolation error of the physical solution ; definition of a physical size map connected to the field enabling relating the error to a given threshold ,remeshing of the domain conforming to the size map ,if the threshold is not satisfied loop to step , with , in order to obtain a new mesh , else , and ,computation of the solutions on ,physical error estimation: computation of the interpolation error of the physical solution ; definition of a physical size map connected to the field enabling relating the error to a given threshold ,remeshing of the domain conforming to the size map ,if the threshold is not satisfied loop to step , with , in order to compute the temperature on the new adapted mesh , else ,detection of the new domain (water vapor) on for a fixed threshold of vaporization in order to produce a mesh and ,computation of the physical solutions on ,physical error estimate: computation of the interpolation error of the physical solution ; definition of a physical size map connected to the field enabling relating the error to a given threshold ,remeshing of the domain conforming to the size map ,if the threshold is not satisfied loop to step , with , in order to compute the temperature on the last adapted mesh , else .
4. Numerical Results and Discussion
Here, we consider a elliptical nanowire of semiaxes ( nm and nm), with thermal conductivity immersed in water ( and ) at temperature . The nanowire is illuminated by a TM polarized laser pulse at wavelength of complex permittivity with a power density per area units [24, 25]. The materials of the system are considered isotropic and homogeneous.
The results of the adaptive process on mesh and on the temperature maps are illustrated in Figure 1. Figures 1(a) and 1(b) show the initial mesh and the associated temperature. The adaptive process on the temperature field (with , nm, nm, and ) produces the mesh and the temperature map. The mesh is adapted on the outline of the nanowire that presents strong variations of the temperature. For a water vaporization threshold ), the detection of the new material (water vapor) is obtained from the temperature map computed on the mesh . Figure 1(c) presents the areas of the three materials: (red), vapor (green), and water (blue). The computation of the temperature on the domain that contains the water vapor requires including the physical parameters of the vapor (permittivity and thermic conductivity ). The spatial distribution of the temperature field on the mesh after detection of the bubble produced around the nanowire is shown in Figure 1(d). The final mesh is obtained, after eight iterations, by applying the adaptive process on the field (with ) taking into account the bubble. That mesh is adapted in the bubble especially on its outline where variations in the temperature occur and relaxed inside the nanowire where the temperature is almost constant (Figure 1(e)). The remeshing process takes into account the shape and size of the bubble. Figure 1(f) shows the temperature map on the mesh after convergence to a stable solution. The level curves are smooth where a strong variation of the temperature is shown (in the vicinity of the boundary of the bubble and the nanowire). The map also shows an increase of the temperature in the nanowire due to the creation of the bubble. Such an increase is due to the diffusion of the temperature, produced by the nanowire after detection of the bubble (i.e., the water vapor has a smaller thermal conductivity than water).
In order to study the evolution of the shape and size of the bubble, we also consider the nanowire illuminated at three wavelengths nm, nm, and nm with physical parameters (, , and ). Figures 2(a), 2(b), 2(c), 2(d), 2(e), and 2(f) present the mesh after bubble detection and distribution of the temperature on the meshes for the three different wavelengths nm, nm, and nm, respectively. These show the evolution of the meshes (shape and size of the bubble) and the temperature as function of the wavelength. For increasing, the imaginary part of the complex permittivity of the decreases, leading to a decrease of the energy absorbed by the nanowire. Therefore, the temperature also decreases and the shape and size of bubble are changing. The created bubble follows the shape of the nanowire (elliptical) at the beginning and becomes circular with the increase of the temperature. Figure 3 shows the evolution of the mean temperature in the nanowire as function of the aspect ratio for three wavelengths , , and nm on the mesh .
Figure 4(a) shows the evolution of the aspect ratio of the bubble ( and being the semiaxes of the bubble) as function of the aspect ratio of the nanowire for the three wavelengths. From the computed data, a function , satisfying , can be obtained through a nonlinear least-squares fit method (LLS) by the Marquardt-Levenberg algorithm [26–29]. The method is used to find a set of the best parameters fitting the data. It is based on the sum of the squared differences or residuals (SSR) between the input data and the function evaluated at the data. The applied algorithm consists in minimizing the residual variance with NDF being the number of degrees of freedom (number of the data points (NDP) minus number of the estimated parameters) after a finite number of iterations (FNI). Therefore, the function can be written as follows: with , , and being set of parameters varying as function of the wavelength. The parameter concerns the asymptote value which is related to the maximum ratio for a circular bubble, is the inverse of the decay rate of the function which is related to the speed tending to the circular shape, and is the initial ratio from which the bubble begins to form. Table 1 shows the fit parameters for each wavelength. The function is continuous and strictly decreasing for in the interval ; therefore the inverse function exists. The measurement of the aspect ratio of the bubble allows predicting the aspect ratio of the nanowire through the relation . Figure 4(b) presents the evolution of the bubble volume (in 2D: ) as function of the volume of the nanowire (i.e., ) for each wavelength. From the computed data (vapor bubble for each nanowire and for each wavelength) and by using the same method and the same algorithm, a function can be obtained through a fit. That function allows obtaining the relation between the volumes . The function is continuous and strictly increasing; therefore the inverse function also exists. The measurement of the bubble volume can be used to determine the volume of the nanowire through the relation . With the two functions and , the size and shape of the nanowire can be obtained from the information on the bubble: consequently, Therefore, the measurement of the size and shape of the bubble can be used to obtain information on the geometry of the nanowire and to reconstruct the size and shape of the nanowire.
The paper focuses on the forming and the evolution of the shape and size of the bubble through a photothermal process between a nanowire of immersed in water and an electromagnetic wave. The increase of temperature is related to the geometry of the nanowire which leads to an increase in the shape and size of the bubble. That solution is computed by developing an adaptive remeshing method. That allows to compute with accuracy the temperature by adapting the mesh to the evolution of the bubble. The coupled problem (light, matter, heat) is solved through an adaptive loop process allowing converging to a stable solution and decreasing the number of nodes. The influence of the laser source and the geometrical parameters (wavelength, size, and shape of the nanowire) related to the size and shape of the bubble are presented and analyzed. The aspect ratio and the volume of the bubble can be expressed as function of the aspect ratio and the volume of the nanowire. By solving the inverse problem, two functions are obtained enabling finding the size and shape of the nanowire from the size and shape of bubble.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
The authors thank the ANR-2011-NANO-008 NANOMORPH for financial support.
M. Boroski, A. C. Rodrigues, J. C. Garcia, L. C. Sampaio, J. Nozaki, and N. Hioka, “Combined electrocoagulation and TiO2 photoassisted treatment applied to wastewater effluents from pharmaceutical and cosmetic industries,” Journal of Hazardous Materials, vol. 162, no. 1, pp. 448–454, 2009.View at: Publisher Site | Google Scholar
D. Barchiesi, T. Grosges, E. Kremer, and M. L. de la Chapelle, “Electromagnetic heat induced in meso-structures: computation of temperature in metallic dimers,” PIERS Online, vol. 7, no. 5, pp. 406–410, 2011.View at: Google Scholar
B. Born and E. Wolf, Principle of Optics, Pergamon Press, Oxford, UK, 1993.
T. Grosges, A. Vial, and D. Barchiesi, “Models of near-field spectroscopic studies: comparison between finite-element and finite-difference methods,” Optics Express, vol. 13, no. 21, pp. 8483–8497, 2005.View at: Google Scholar
J. Jin, The Finite Element Method in Electromagnetics, John Wiley & Sons, New York, NY, USA, 1993.
P. Silvester and G. Pelosi, Finite Elements for Wave Electromagnetics: Methods and Techniques, IEEE Press, New York, NY, USA, 1994.
I. Stakgold, Boundary Value Problems of Mathematical Physics, Volumes I and II, Macmillan, New York, NY, USA, 1969.
R. Radovitzky and M. Ortiz, “Error estimation and adaptive meshing in strongly nonlinear dynamic problems,” Computer Methods in Applied Mechanics and Engineering, vol. 172, no. 1–4, pp. 203–240, 1999.View at: Google Scholar
M. Ainsworth and J. T. Oden, “A posteriori error estimation in finite element analysis,” Computer Methods in Applied Mechanics and Engineering, vol. 142, no. 1-2, pp. 1–88, 1997.View at: Google Scholar
A. O'Connor, O. Morris, and E. Sokell, “Angular and energy distribution of Sn ion debris ejected from a laser-produced plasma source, for laser power densities in the range suitable for extreme ultraviolet lithography,” Journal of Applied Physics, vol. 109, no. 7, Article ID 073301, 2011.View at: Publisher Site | Google Scholar
K. Levenberg, “A method for the solution of certain problems in least squares,” Quarterly of Applied Mathematics, vol. 2, pp. 164–168, 1944.View at: Google Scholar
V. V. Vasin and G. Y. Perestoronina, “The Levenberg-Marquardt method and its modified versions for solving nonlinear equations with application to the inverse gravimetry problem,” Proceedings of the Steklov Institute of Mathematics, vol. 280, no. 1, supplement, pp. 174–182, 2013.View at: Publisher Site | Google Scholar