A method for improving nanospraying precision based on molecular dynamics simulation

Abstract

The application discloses a kind of nanometer spraying precision promotion method based on molecular dynamics simulation, it is related to nanometer manufacturing and spraying technical field, to solve the technical problems of uncontrollable liquid film spreading form, serious sputtering, big parameter control blindness in existing nanometer spraying.The method is based on the dynamic mechanism of binary nanometer droplet impact liquid film revealed by molecular dynamics (MD) simulation, accurately identify three typical dynamic results of regular crown diffusion, columnar jet, central sputtering, by synergistically regulating droplet impact weber number (We) We )、Droplet dimensionless center distance (D L )And liquid film dimensionless thickness (H H ),Combination sub-region theory model (H H ≤0.5 and H >0.5) Accurately predict the maximum crown radius (R R c,max ), Realize the directional control of liquid film spreading form.The application makes full use of the dynamic law of droplet impact under nanometer scale, can significantly improve the uniformity and controllability of nanometer spraying, applicable to nanometer inkjet printing, nanometer spray cooling, nanometer 3D printing, nanometer precision coating preparation and other various scenes needing high-precision droplet control, material utilization rate and spraying precision are more than 25%-40% than traditional method.

Description

Technical Field

[0001] This invention relates to the field of nanofabrication and spraying technology, specifically to a method for improving the precision of nanospraying based on the dynamic control of binary nanodroplet impact. This method is applicable to various scenarios that rely on the precise manipulation of nanodroplets, including but not limited to nano-inkjet printing, nano-3D printing, nano-spray cooling, nano-precision coating preparation, and micro / nano spraying of pesticides. It is particularly suitable for high-precision manufacturing processes with strict requirements for liquid film uniformity, jet stability, or sputtering suppression. Background Technology

[0002] Nanospraying technology is a core supporting technology in the field of micro- and nano-manufacturing, and its precision directly determines the performance and reliability of nano-devices and functional coatings. In practical applications such as nano-inkjet printing, nano-spray cooling, and nano-3D printing, the process of nanodroplets impacting the substrate surface and covering it with a liquid film is a key step in forming a uniform coating, improving heat transfer efficiency, or constructing fine structures. This process requires directional control of the liquid film spreading morphology—for example, nano-inkjet printing requires regular crown-shaped uniform diffusion, nano-spray cooling requires a stable columnar jet to enhance heat transfer, and the preparation of precision nano-coatings requires strict suppression of central sputtering to reduce material waste. However, this process involves complex solid-liquid interface interactions, liquid-liquid extrusion coupling, and significant nanoscale effects, which makes the precise control of the liquid film spreading morphology face multiple technical bottlenecks. Existing technologies are insufficient to meet the demands of high-precision manufacturing. Specific problems are as follows:

[0003] 1. The dynamic mechanisms at the nanoscale differ fundamentally from those at the macroscopic scale, rendering traditional theories ineffective. Traditional droplet impact studies have largely focused on macroscopic scales from millimeters to centimeters, with their dynamics based on classical fluid dynamics. However, when droplet sizes shrink to the nanometer scale (e.g., ...), the impaction dynamics become more complex. D 0 = 8nm), the significant scale effect leads to a fundamental change in the dynamic mechanism: on the one hand, the dimensionless parameter Ohnesorge number ( Oh = μ / ( ρD 0 γ ) 1 / 2 From the macroscopic O (10 - ³) The increase in O (1) to the nanoscale represents an increase of three orders of magnitude, making viscous forces indispensable even in low-viscosity liquids (such as water), thus becoming a key factor in the dominant energy dissipation. Furthermore, the velocity gradient within the nanodroplet can penetrate the entire droplet body, unlike macroscopic droplets which are confined to the boundary layer, resulting in significantly different film spreading rates, jet formation, and sputtering thresholds compared to macroscopic patterns. In addition, Figure 5Data shows that viscous dissipation accounts for as much as 80%-90% of the total energy at the nanoscale, while this proportion is almost negligible at the macroscale. This further leads to the inability of traditional macroscopic fluid dynamics models to accurately describe the impact behavior of nanodroplets.

[0004] 2. There is a significant research gap in the study of nanodroplet impacts on liquid films, and existing data are insufficient to support this research. In both natural and industrial settings, the impact of binary or multi-component droplets on liquid films is far more common than that of single droplets—for example, continuous jetting of droplets in nano-inkjet printing and dense droplet swarms in nano-spray cooling both involve interactions between droplets. However, existing research has significant gaps: on the one hand, most nanoscale studies focus only on single droplet impacts on dry substrates, neglecting the influence of liquid films covering the substrate surface in real-world scenarios, while the presence of liquid films significantly alters energy transfer and morphological evolution after droplet impact; on the other hand, the few studies involving liquid films are either limited to macroscopic scales or only focus on the thickness of a single liquid film or the distance between droplets, failing to systematically explore the centroid distance of binary droplets at the nanoscale. L ), liquid film thickness ( H ) and impact velocity ( We The coupling effect of droplets and liquid films is significant. Due to the scale effect, research conclusions on macroscopic or single-droplet systems cannot be directly transferred to nanoscale droplet-liquid film systems, resulting in a lack of molecular-level dynamic mechanisms to support current technologies and extremely high degree of blindness in parameter control.

[0005] 3. The influence of liquid film thickness has not been fully revealed, and the control range lacks a basis. liquid film thickness ( H = H / D 0) is a key parameter affecting droplet impact morphology, but existing research has limitations in understanding its role: early studies only fixed H While focusing on ultrathin liquid film scenarios where the film thickness is comparable to the droplet size, this study does not cover nanoscale applications with thicknesses ≤ 0.1. H The entire range of values ​​≤1. System simulation revealed that... H It has a decisive influence on the propagation direction, spreading rate, and energy dissipation of the liquid cap: H When the value is ≤0.5, the liquid film crown spreads obliquely, with a fast early spreading rate but inhibited by solid-liquid interaction in the later stage; H When the value is >0.5, the liquid film propagates vertically, kinetic energy dissipation accelerates, and jetting and sputtering are significantly suppressed (e.g., Figure 5 (As shown in c). Furthermore... H It will also affect the thickness of the bottom of the liquid film depression—whenH When the size is small, the liquid film at the bottom of the depression is easily thinner than 1 nm (e.g. Figure 6 (As shown in b), this triggers direct solid-liquid interaction, further altering the kinetic behavior. Current technologies have not discovered this mechanism, resulting in a lack of theoretical guidance for selecting liquid film thickness and making targeted control difficult.

[0006] 4. Lack of accurate prediction models for nanoscale droplets Maximum crown radius ( R c,max Film spreading is a core indicator for measuring spraying accuracy, and accurate prediction of its properties is a prerequisite for controlling the film spreading range. Existing prediction models have two major shortcomings: first, they mostly target single droplets or macroscopic systems, failing to consider the squeezing effect between binary droplets; second, they do not distinguish between... H ≤0.5 and H The two propagation modes with a value greater than 0.5 result in insufficient prediction accuracy. Figure 9 The traditional model shows that... H When ≤0.5, as We As the temperature rises, the deviation between the predicted values ​​and the MD simulation results gradually increases; and even with the proposed regional model, in H ≤0.5 and high We Under these conditions, the accuracy still has room for improvement due to the insufficient quantification of the additional energy dissipation caused by solid-liquid interactions. The limitations of existing models make it impossible to predict the liquid film spreading range in advance during actual spraying, and parameters can only be repeatedly adjusted and optimized, which is not only inefficient but also results in a large amount of material waste.

[0007] 5. Poor controllability of sputtering and jetting makes it difficult to match the needs of different application scenarios. Different nano-spraying scenarios have vastly different requirements for droplet impact morphology: nano-inkjet printing needs to avoid jetting and sputtering, retaining only regular crown-shaped diffusion; nano-spray cooling requires enhanced columnar jetting to increase heat transfer area; and precision 3D printing requires complete suppression of center sputtering. However, current technologies cannot clearly define the control boundaries of jetting and sputtering: high We The center is prone to sputtering (e.g.) Figure 2 (As shown in c), this results in material waste; and when a jet is needed, there is a lack of parameter combinations that can stably form a columnar jet. Studies have shown that... We It is the key to triggering the form change — We When ∈ [77.34, 120.85], it is a regular crown-shaped diffusion. We When the velocity is ∈ [120.85, 174.02], a columnar jet is formed. We Center sputtering occurs at >174.02 (e.g.) Figure 2 (as shown); at the same time L and H The coupling will further modulate this boundary— L Increasing the size can suppress sputtering. H Increasing the size of the jet stream will simultaneously suppress both jetting and sputtering (e.g.) Figure 4 (As shown). Existing technologies have not grasped the coupling law of this parameter, making it impossible to control it in a targeted manner according to the target shape.

[0008] In summary, existing nano-spraying technologies suffer from uncontrollable film spreading morphology, low spraying accuracy, and poor material utilization due to a lack of in-depth understanding of the dynamics mechanism of binary droplet impacts on liquid films at the nanoscale and the absence of effective parameter optimization strategies and accurate prediction models. Therefore, developing a multi-parameter synergistic optimization and accurate prediction method based on the fundamental laws revealed by molecular dynamics simulations is crucial to overcoming existing technological bottlenecks and improving the accuracy of nano-spraying. Summary of the Invention

[0009] In view of this, the present invention discloses a method for improving the precision of nano-spraying based on molecular dynamics simulation.

[0010] It should be noted that this invention aims to solve the following core problems existing in current nano-spraying technology. All of these problems are based on the unique laws governing the impact of nanoscale droplets on liquid films revealed by research, and existing technologies cannot effectively address them: 1. The control boundaries of the three dynamic results of nanodroplet impact on liquid film (regular crown diffusion, columnar jet, and central sputtering) are unclear. Central sputtering is unique to binary nanodroplets (this phenomenon does not exist in single droplet systems), and its trigger threshold and parameter coupling relationship have not been quantified. 2. We , L , H The coupling mechanism of the three parameters has not been fully utilized—existing technologies have not discovered it. H The regulatory effect on the propagation direction of the liquid film crown (oblique / vertical, such as...) Figure 5 (as shown in c), also not clearly defined. L It only affects the jet and sputtering without affecting the crown diffusion rate (e.g.) Figure 4 The characteristics shown in e and 6d indicate that parameter optimization lacks theoretical support; 3. Different H The influence of the propagation direction of the lower liquid membrane crown on the spreading rate was not included in the prediction model. H When the value is ≤0.5, oblique propagation leads to rapid early spreading. H Vertical propagation at >0.5 results in slow spreading (e.g., Figure 5 As shown in a), traditional models do not distinguish this difference, resulting in large prediction errors; 4. High-viscosity dissipation at the nanoscale (accounting for 80%-90%, such as...) Figure 5 (As shown in b) and solid-liquid interaction (triggered when liquid film thickness < 1 nm, such as...) Figure 6 (as shown in b) leads to R c,max Low prediction accuracy – Existing models do not consider the combined effect of these two types of factors, especially in H ≤0.5 and high We Significant time deviation (e.g.) Figure 9 (as shown in b) 5. The evolution of the "half-moon" raised liquid film between nanodroplets is uncontrollable - This structure is formed by the extrusion between droplets, and its height and stability directly affect the uniformity of spraying. Existing technologies lack targeted control methods.

[0011] To achieve the above-mentioned technical objectives, the present invention adopts the following technical solution: A method for improving the precision of nano-spraying based on molecular dynamics simulation includes the following steps: (1) Determine the basic parameters of the target substrate, liquid film system and nanodroplets for nano-spraying, including the initial diameter of the droplets. D 0, liquid density ρ Viscosity μ and surface tension γ ; (2) Based on the results of molecular dynamics simulation, the core control parameter affecting the spraying effect is identified: Weber number. We (Characterizing the ratio of inertial force to capillary force,) We = ρD 0 V 0 2 / γ ), dimensionless centroid distance L ( L = L / D 0, L Let be the distance between the centroids of the droplets, 1 ≤ L ≤2) Dimensionless liquid film thickness H ( H = H / D 0, H The initial liquid film thickness is 0.1 ≤ H ≤1); (3) Select the optimal combination of core control parameters based on the target spraying morphology (regular crown diffusion, columnar jet, centerless sputtering): When the goal is to achieve uniform diffusion of regular crown-shaped structures, regulation We ∈[77.34, 120.85], L ∈[1.25, 1.5], H ∈[0.375, 0.5], at this point there is no jetting or sputtering of the liquid film (e.g. Figure 2 (as shown in a) When the goal is to stabilize the columnar jet, the control We ∈[120.85, 174.02], L =1.25, H ∈[0.25,0.375], at this time an edge columnar jet is formed (e.g. Figure 2 (as shown in b) When the goal is to suppress center sputtering, adjust We ≤174.02, L ≥1.5, H ≥0.75, at which point the liquid film crown diffusion is stable (e.g. Figure 4 (as shown in d and 4e). (4) Based on H The value of is used to predict the maximum crown radius using the corresponding regional theoretical model. R c,max : when H When >0.5, the liquid film crown propagates vertically (e.g. Figure 5 As shown in c), the formula is used. calculate; when H When ≤0.5, the liquid film crown propagates obliquely (e.g. Figure 5 As shown in c), the formula is used. calculate; (5) Start the nano-spraying equipment according to the optimized parameters, so that the nano-droplets hit the liquid film on the substrate surface. By comparing the real-time monitoring with the theoretical prediction value, fine-tune the parameters to complete the high-precision spraying.

[0012] Furthermore, the target substrate has a contact angle. θ For solid materials with a wettability of 0° to 40°, the substrate surface is reinforced with virtual springs to prevent impact deformation (e.g., Figure 1(As shown).

[0013] Furthermore, the nanodroplets were simulated using the mW single-atom water model, and the liquid film thickness was... H The value range is 0.8nm~8nm (corresponding to H The interaction between droplets and liquid films and substrates is described by the Leonard-Jones (LJ) potential (∈[0.1, 1]).

[0014] Furthermore, the theoretical model described in step (4) is based on the principle of energy conservation and takes into account the droplet kinetic energy. E d,k Surface energy E d,s Liquid film surface energy E f,s and viscous dissipation W The coupling effect.

[0015] Furthermore, the droplet ejection speed of the nano-coating device V 0 passed We The value was calibrated, and the simulation process was divided into two stages: NVT ensemble equilibration (2ns, 300K) and NVE ensemble impact (1ns).

[0016] Compared with existing technologies, the present invention has the following beneficial effects: The beneficial effects of this invention stem directly from the nanoscale laws revealed by the research and the precise control of this method. Furthermore, as verified by MD simulations and practical examples, it demonstrates significant innovation and practicality. High morphological orientation accuracy – achieving 100% orientation control of three dynamic results: thickness error ≤3.2% during uniform spraying (Example 1), jet stability ≥92% (Example 2), sputtering suppression rate ≥90% (Example 3), and the unique problem of center sputtering is completely solved; Prediction accuracy is significantly improved – regional models in H When the value is >0.5, the error is ≤0.45%. H When ≤0.5, the error is ≤10% (e.g. Figure 9 As shown in the figure, it is more than 30% better than the traditional model (error > 15%), and can predict the spreading range in advance, avoiding blind debugging; Energy utilization efficiency optimization—by combining parameters to stabilize the proportion of viscous dissipation at 80%-85% (e.g. Figure 5 (As shown in b), it utilizes high dissipation at the nanoscale to suppress sputtering while avoiding excessive dissipation that leads to insufficient spreading, thus improving energy utilization by 20%. Material utilization is greatly improved – after suppressing center sputtering, material waste is reduced by 25%-40% (up to 38% in Example 3), which is especially suitable for high-value systems such as precious metal nano-inks and bioactive materials; High versatility and compatibility – compatible with various substrates such as Pt and Si-based modified materials, and suitable for various droplet systems such as water-based and organic-based (as long as the LJ potential interaction is satisfied). No modification to existing nano-spraying equipment is required, only parameter calibration. Attached Figure Description

[0017] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on the provided drawings without creative effort.

[0018] Figure 1 This is an initial model diagram of a simulation system for binary nanodroplets impacting a liquid film.

[0019] Figure 2 The image shows the impact results of nanodroplets impacting the liquid film at different impact velocities.

[0020] Figure 3 This is a cloud map showing the velocity distribution inside the liquid during the droplet impact on the liquid film.

[0021] Figure 4 A comparison of the effects of droplet distance and liquid film thickness on the impact results.

[0022] Figure 5 This is a comparison diagram of nanodroplets and millimeter droplets impacting a liquid film.

[0023] Figure 6 The maximum propagation distance of a droplet impacting a liquid film R c Comparison chart showing the effects of liquid film thickness and droplet distance.

[0024] Figure 7 This is a schematic diagram of a newly predicted physical model for the maximum propagation distance of nanodroplets impacting a liquid film.

[0025] Figure 8 This is a comparison diagram of energy dissipation during droplet impact on the liquid film.

[0026] Figure 9 This figure compares the theoretical maximum propagation distance of nanodroplets impacting a liquid film with experimental results. Detailed Implementation

[0027] The technical solutions in the embodiments of the present invention will be clearly and completely described below. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0028] The term "embodiment" used herein, as an example, is not necessarily to be construed as superior to or better than other embodiments. Performance testing in the embodiments of this application, unless otherwise specified, employs conventional testing methods in the art. It should be understood that the terminology used in this application is merely for describing particular implementations and is not intended to limit the scope of this disclosure.

[0029] Unless otherwise stated, the technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application pertains; other experimental methods and technical means not specifically mentioned herein refer to experimental methods and technical means commonly used by one of ordinary skill in the art.

[0030] To better illustrate the content of this application, numerous specific details are provided in the following detailed embodiments. Those skilled in the art should understand that this application can be implemented even without certain specific details. In the embodiments, some methods, means, instruments, and devices well-known to those skilled in the art are not described in detail in order to highlight the main points of this application.

[0031] Without conflict, the technical features disclosed in the embodiments of this application can be combined arbitrarily, and the resulting technical solution belongs to the content disclosed in the embodiments of this application.

[0032] This invention discloses a method for improving the precision of nano-spraying based on molecular dynamics simulation.

[0033] Specifically, the accompanying diagram is explained below: Figure 1 The initial model diagram of the simulation system for binary nanodroplets impacting a liquid film is shown: It illustrates the binary nanodroplets ( ) in a 60×40×25 nm³ system. D The spatial relationship between the liquid film (0=8nm) and the Pt substrate, with the substrate atoms fixed by virtual springs, provides a geometric reference for parameter setting.

[0034] Figure 2 Figure 1 shows the impact results of nanodroplets impacting the liquid film at different impact velocities: Figure 2 a( We =120.85): Regular crown-shaped diffusion, double scallop structure + surrounding crown, no jet sputtering, for the target shape of uniform spraying; Figure 2 b ( We=174.02): Columnar jet, with a stable jet formed at the crown edge, which is the target morphology for nano-spray cooling; Figure 2 c ( We =236.86): Center sputtering, crown center ejection of sub-droplets, which is the form that needs to be suppressed for precision spraying.

[0035] Figure 3 A cloud map showing the velocity distribution inside the liquid during the droplet impact on the liquid film: Figure 3 a1, b1-d1 ( x (axial slice): + V h Migration towards the edge of the coronavirus explains the coronavirus spread mechanism; Figure 3 a2, b2-d2 ( y Axis slice): Jet region concentration + V v with + V h This reveals the principle behind the formation of columnar jets.

[0036] Figure 4 A comparison chart showing the effects of droplet distance and liquid film thickness on the impact results: Figure 4 a( We =120.85): Regular crown diffusion; Figure 4 b ( We =236.86): strong jet and splash, reflecting We The regulatory effect; Figure 4 c ( H =0.125): Vertical crown elevation is suppressed; Figure 4 d( H =0.75): jet sputtering is suppressed, reflecting H The regulatory effect; Figure 4 e ( L =1.75): Splash disappears, reflecting L The inhibitory effect.

[0037] Figure 5 A comparison image showing the impact of nanodroplets and millimeter droplets on a liquid film: Figure 5 a: Nanoscale R c The growth rate is lower than the macro level, and H When ≤0.5, early spreading is faster than H >0.5; Figure 5 b: Nanoscale viscous dissipation accounts for 80%-90%, providing an energy basis for model construction; Figure 5 c:H Oblique propagation when ≤0.5 H Vertical propagation is achieved when the value is >0.5, providing a basis for regional models.

[0038] Figure 6 The maximum propagation distance of a droplet impacting a liquid film R c Comparison chart showing the effects of liquid film thickness and droplet distance: Figure 6 a: We >100 hours R c The growth rate is divided into two stages, which are related to the thickness of the bottom of the liquid film depression being <1nm; Figure 6 b: different We The critical thickness (1 nm) at the bottom of the lower liquid film depression triggers solid-liquid interaction; Figure 6 c: H Increase inhibition R c increase; Figure 6 d: L right R c No significant effect.

[0039] Figure 7 A schematic diagram of the physical model for the newly predicted maximum propagation distance of nanodroplets impacting a liquid film: Figure 7 a: Initial state energy composition ( E d,k + E d,s + E f,s ); Figure 7 b, 7c: H ≤0.5 (oblique propagation) and H The maximum spread shape >0.5 (vertical propagation) provides a geometric basis for model derivation.

[0040] Figure 8 A comparison of energy dissipation during droplet impact on a liquid film: different L Down W dim Basically unchanged, proof L It has no impact on energy consumption.

[0041] Figure 9 A comparison of the theoretical maximum propagation distance of nanodroplets impacting a liquid film and experimental results: Figure 9 a( H ≥0.5): The model and MD results are highly consistent; Figure 9 b ( H ≤0.5): The prediction error of the corrected model is significantly reduced.

[0042] This invention is based on systematic results from molecular dynamics simulations (covering systems of 60×40×25 nm³, 1≤ L ≤2、0.1≤ H ≤1、77.34≤ We Using a full parameter space of ≤296.08, a four-in-one method for improving the precision of nano-spraying is constructed, which integrates "parameter control - morphology orientation - model prediction - process monitoring". Each step is closely combined with molecular-level mechanisms and data support. The specific technical solution is as follows: (I) Definition and control mechanism of core parameters (e.g.) Figure 2-6 (Quantitative regularity shown in Figure 8) This invention clarifies the physical meaning, control range, and mechanism of action of three core parameters, and all the laws have been verified by MD simulation: 1. Weber number We (77.34≤) We ≤296.08) Physical meaning: Characterizes the ratio of inertial force to capillary force in a droplet ( We = ρD 0 V 0 2 / γ ), which directly determines the intensity of kinetic energy input; Regulation mechanism: We It is the dominant parameter for morphological transformation, exhibiting a step-like regulatory effect (e.g. Figure 2 (as shown) We ∈[77.34, 120.85]: Inertial forces and capillary forces are balanced, and the liquid film forms a "double scallop + surrounding crown" structure (e.g. Figure 2 As shown in a), there is no jet or sputtering, corresponding to regular crown-shaped diffusion; We ∈[120.85, 174.02]: Increased inertial force causes the central "crescent-shaped" bulge liquid film to form a columnar protrusion, which then develops into a stable columnar jet (e.g., Figure 2 (as shown in b), jet region concentration + V v with + V h (like Figure 3(As shown in a2 and b2), the horizontal velocity decays more slowly than the vertical velocity, ensuring that jet extension takes precedence over lifting; We >174.02 (e.g., 236.86, 296.08): Inertial force far exceeds capillary force, the central liquid film is torn apart, and multiple sub-droplets are ejected (e.g. Figure 2 (as shown in c), i.e., center sputtering, which needs to be strictly suppressed in precision spraying; Key findings: We When >100, R c The growth rate is divided into two stages (e.g.) Figure 6 (As shown in a) — When the thickness of the liquid film at the bottom of the depression is >1nm, it grows rapidly; when it is <1nm, it grows rapidly due to solid-liquid interactions (such as... Figure 6 (As shown in b) It grows slowly, and this critical value (1nm) provides a basis for parameter calibration.

[0043] 2. Dimensionless parameter of centroid distance L (1≤ L ≤2) Physical meaning: L = L / D 0 ( L (where is the distance between the centroids of the nanodroplets), characterizing the mutual compression strength between droplets; Regulatory mechanisms (such as) Figure 4 e, 8 (as shown in Figure 8) L ∈[1, 1.5]: Strong inter-droplet compression facilitates the formation of a "crescent-shaped" raised liquid film, easily triggering a columnar jet ( L The jet is most stable when the value is 1.25, such as... Figure 2 (as shown in b) L ∈[1.5, 2]: The squeezing effect is weakened, the height of the "crescent-shaped" bulge liquid film is reduced, and the jet and sputtering are suppressed. L There is no splashing when =1.75, such as Figure 4 (as shown in e) Unique characteristics: L It has no significant effect on the diffusion rate of the liquid film crown (e.g.) Figure 6 As shown in d), and without changing the proportion of energy dissipation (different L Down W dim Deviation <1.3%, such as Figure 8 As shown in the figure, the jet / splash can be independently controlled without affecting the spreading range.

[0044] 3. Dimensionless parameter for liquid film thickness H (0.1≤ H ≤1) Physical meaning: H = H / D 0 ( H (where is the initial liquid film thickness), characterizing the liquid film's ability to buffer and transfer impact energy; Regulatory mechanisms (such as) Figure 4 (as shown in c, 4d, 5c, and 6c) Control of propagation direction: H When ≤0.5, the liquid film crown propagates obliquely (e.g. Figure 5 As shown in the left figure (c), the early spreading rate is fast ( τ <12.5 hours R c Larger, such as Figure 5 (as shown in the local evolution). H When the value is >0.5, it propagates vertically (e.g., Figure 5 (As shown in the right figure), the early spreading rate is slow; Morphological inhibition regulation: H ∈[0.1, 0.375]: The liquid film is thin, the solid-liquid interaction is strong, and the vertical lifting of the crown is suppressed, but jetting and sputtering are still likely to occur (e.g. Figure 4 (as shown in c) H ∈[0.5, 1]: Thick liquid film, rapid dissipation of kinetic energy ( H When the value is 0.75, the dissipation rate is 89.6%, such as... Figure 5 As shown in b), jetting and sputtering are significantly suppressed (e.g. Figure 4 (as shown in d) Critical value effect: H =0.5 is the dividing point. Below this value, the bottom of the liquid film depression is easily <1nm (triggering solid-liquid interaction), while above this value, it is difficult to contact the substrate (avoiding direct solid-liquid interaction), providing a basis for the partitioning of the theoretical model.

[0045] (II) Parameter Combination Optimization Strategy (Based on) Figure 2 , 4 (Simulation results of 6) For different spraying scenarios, precise parameter combination schemes are proposed. All combinations have been verified by MD simulation and can be directly applied. 1. Uniform spraying scenarios (nano inkjet printing, nano precision coating) Target morphology: Regular crown-shaped diffusion, no jet, no splash, "crescent-shaped" raised liquid film height ≤ 0.5 mm. D 0 (To avoid uneven coating thickness); Optimized combination: We =120.85±5、 L =1.25-1.5、 H =0.375-0.5; in accordance with: We =120.85 is the critical value between regular crown-shaped and columnar jets (e.g.) Figure 2 As shown in a), a fluctuation range of ±5 can balance spreading efficiency and stability. L =1.25-1.5 Ensure proper compression between droplets to avoid excessively high "crescent-shaped" bulges (e.g., Figure 4 (as shown in the comparison) H =0.375-0.5 can balance early spreading rate and later stability. H When the density is 0.5, the bottom of the liquid film depression is close to 1 nm, and the solid-liquid interaction is mild (e.g., ...). Figure 6 As shown in b), viscous dissipation accounts for 81.2% (e.g. Figure 5 (as shown in b), to avoid excessive dissipation leading to insufficient spreading.

[0046] 2. Jet enhancement scenarios (nano spray cooling, targeted nano spraying of pesticides) Target shape: stable columnar jet, jet duration ≥50ps, no center splash; Optimized combination: We =174.02±10、 L =1.25、 H =0.25-0.375; in accordance with: We =174.02 represents the stable stage of the columnar jet (e.g.) Figure 2 As shown in b), the jet height can be adjusted within a range of ±10 (as shown in b). We The larger the jet, the higher the jet. L At a concentration of 1.25, the interdroplet compression is strongest, and the "crescent-shaped" bulging liquid film is highest, providing a sufficient fluid source for jet formation (e.g., Figure 4 (as shown in b) H When the value is 0.25-0.375, the liquid film crown propagates obliquely (e.g. Figure 5 (as shown in c), + V h The migration rate towards the crown edge is fast (e.g.) Figure 3 As shown in d1), the jet extension effect is optimal, and the bottom of the liquid film depression is >1nm (to avoid solid-liquid interaction weakening the jet, such as...). Figure 6 (as shown in b).

[0047] 3. Sputtering suppression scenarios (nano 3D printing, nano-spraying of precious materials) Target morphology: Regular crown-shaped diffusion, no central sputtering, material utilization rate ≥90%; Optimized combination: We ≤120.85、 L ≥1.5、 H ≥0.75; in accordance with: We ≤120.85 can avoid center tearing caused by excessive inertial force (e.g.) Figure 2 (as shown in a) L ≥1.5 reduces droplet compression, lowers the height of the "crescent-shaped" raised liquid film, and makes the sputtering triggering conditions insufficient (e.g.) Figure 4 (as shown in e) H When ≥0.75, the liquid film crown propagates vertically (e.g. Figure 5 As shown in c), kinetic energy is dissipated quickly ( H When the value is 0.75, the dissipation rate is 89.6%, such as... Figure 5 As shown in b), the bottom of the liquid film depression makes it difficult to contact the substrate, thus failing to form the energy concentration required for sputtering.

[0048] 4. Extreme thin liquid film scenarios ( H ≤0.25, such as in the preparation of ultrathin nano-coatings) Target morphology: Regular crown-shaped diffusion, suppressing vertical lift and ensuring uniform coating thickness; Optimized combination: We =77.34-100、 L =1.25-1.5、 H =0.1-0.25; in accordance with: H When the value is ≤0.25, the hydrophilic substrate limits the horizontal velocity. V h Dissemination (e.g.) Figure 4 (as shown in c) We It needs to be lowered to 77.34-100 to avoid excessive vertical lifting; L =1.25-1.5 can partially offset the inhibition of spreading by solid-liquid interactions, ensuring R c To reach the target value (e.g.) Figure 5 (as shown in a).

[0049] (III) Construction of regional theoretical models (based on energy conservation and geometric models) This invention is based on the energy conservation framework, combined with H Influence on propagation direction and correction of solid-liquid interaction, establishing regional divisions R c,max The prediction model improves prediction accuracy by more than 30% compared to traditional models. 1. Basic Assumptions of the Model Initial total energy: E total = E d,k + E d,s + E f,s ( E d,k For the droplet kinetic energy, E d,s For the surface energy of the droplet, E f,s (for the surface energy of the liquid film). Energy at maximum spread: E total = E f,fs + W ( E f,fs The total surface energy of the system. W (for viscous dissipation); Viscous dissipation W The single-droplet impact formula is used as an approximation because L No impact W (like Figure 8 As shown), no additional corrections are needed.

[0050] 2.H >0.5 model (vertical propagation, such as Figure 7 (as shown in c) Applicable conditions: H >0.5, the liquid film crown rises vertically, and the solid-liquid interaction is negligible (bottom of the liquid film depression >1nm); Geometric relationship: Liquid film crown height H lc = H 2- ( H - H 1) H 1 represents the thickness of the fluid film at the bottom of the coronary cavity. H 2 represents the coronal cavity depth), during vertical propagation α =90° (the liquid film crown is perpendicular to the horizontal surface). Derivation process: Combining the formula ( E f,fs - E f,s = 2π R c,max H lc + π R c,max ²) and the energy balance equation, substituting into H lc The viscous dissipation characteristics of vertical propagation ultimately lead to: ; Accuracy verification: Compared with MD simulation results, the error is ≤0.45% (e.g. Figure 9 (as shown in a), such as We =120.85、 H When the value is 0.75, the model predicts a value of 2.21 and the MD value is 2.19.

[0051] 3. H ≤0.5 model (oblique propagation, such as) Figure 7 (as shown in b) Applicable conditions: H ≤0.5, the liquid film crown propagates obliquely, and the solid-liquid interaction is significant (when the bottom of the liquid film depression is <1nm). Geometric relationship: Angle between the liquid film cap and the horizontal surface α =60°±5° (e.g.) Figure 5 (as shown in the actual measurement) Effective spreading radius during oblique propagation R c,eff = R c,max ×cosα ; The correction introduces: solid-liquid interactions lead to additional viscous dissipation, according to... Figure 5 b data, H =0.25 dissipation accounts for a relatively large proportion H =0.75 higher than 8.4%, so a correction factor of 0.85 is introduced; Derivation process: Combining the geometric correction for oblique propagation and the energy correction for solid-liquid interaction, we finally obtain: ; Accuracy verification: Low We Time error ≤3%, high We ( We =174.02) When the error is ≤10% (e.g. Figure 9 As shown in b), it is a significant improvement over the uncorrected model (error > 15%).

[0052] 4. Prediction of liquid film height in crescent-shaped bulges Based on such Figure 4 Data b and e were fitted to obtain a height. h = h / D Relationship between 0 and parameters: h =0.12×( We / 120.85)×(1.5 / L )×(0.5 / H ); This model can predict the height of the bulge between droplets, preventing it from being too high and causing the coating to bulge, thus further improving uniformity.

[0053] (iv) Spraying Implementation Steps The implementation steps of this invention strictly follow the logic of thermodynamic equilibrium and impact processes in MD simulation, ensuring consistency between actual spraying and simulation results: 1. Substrate and liquid film pretreatment Substrate selection: Pt or Pt-like wettable materials ( θ 0=40°), the surface uses virtual springs to fix atoms (to avoid impact deformation, such as Figure 1 As shown in the figure, in practical applications, atomic-level smoothness can be achieved through surface modification; Liquid film preparation: An initial liquid film is formed using high-precision coating equipment, and the process is controlled. H Meet the goal H = H / D 0 (D 0=8nm), liquid film thickness uniformity error ≤5% (refer to the initial state of the liquid film simulated in the study). Interaction calibration: The interactions between droplets and liquid films, and between droplets and substrates, are calibrated using LJ potential parameters. ε w =0.26838eV, σ w =0.23925nm; ε s-l =0.022eV).

[0054] 2. Droplet parameter setting and calibration Droplet generation: Adjust the spraying equipment to generate droplets. D Binary nanodroplets with a density of 0=8nm were used, and the physical properties of the droplets were calibrated using the mW single-atom water model (density and surface tension are consistent with the simulation). Centroid distance control: The distance between the center of mass and the droplet jet spacing is adjusted via the device's droplet jet spacing adjustment function. L Meet the goal L = L / D 0, control accuracy ±0.1nm (ensuring) L Error ≤ 0.01); We Value calibration: based on We = ρD 0 V 0 2 / γ Calculation target V 0. The equipment's jet speed is calibrated using a laser velocimeter to ensure... We Value deviation ≤ ±2%.

[0055] 3. Thermodynamic Equilibrium and Impact Process Control Equilibrium Phase: The equipment simulates the NVT ensemble environment, controlling the temperature of the spraying area at 300K (Nose-Hoover thermostat), with an equilibrium time of 2ns (to ensure the energy stability of the liquid film and droplets). Impact Phase: Switch to NVE ensemble (without thermostat), initiate droplet ejection, impact process lasts 1 ns (sufficient to observe complete morphological evolution, such as...) Figure 2 The time series shown covers 14-34 ps. Data recording: Liquid film morphology and atomic positions were recorded every 2 ps using a high-speed visualization device (OVITO) (referencing the data recording frequency in the research simulation), and data was extracted. R c ( t )and h ( t(Height of the raised liquid film).

[0056] 4. Real-time monitoring and parameter fine-tuning Morphological monitoring: Comparing the real-time recorded liquid film morphology with the target morphology (e.g., Figure 2 (as shown in a, 2b, and 4d), determine whether a jet or splash occurs; Predicted value comparison: Extracting real-time values R c ( t ), compared with the regional model predictions R c,max When the deviation exceeds 5%, fine-tuning is initiated. Fine-tuning rules: like R c ( t Smaller than average: We Increase by 5-10, or H Reduce by 0.05-0.1 (not lower than 0.1); If splashing occurs: L Increase by 0.1-0.2, or We Reduce by 10-15; If the jet is unstable: L Adjust to 1.25, or H Adjusted to 0.375.

[0057] To better understand the present invention, the following embodiments are provided for further detailed description of the present invention, but they should not be construed as limiting the present invention. Any non-essential improvements and adjustments made by those skilled in the art based on the above-described invention are also considered to fall within the protection scope of the present invention.

[0058] Example 1: Regular coronal diffusion 1. Simulation Objective: To achieve regular crown-shaped diffusion of binary nanodroplets after impacting a liquid film, without jetting or sputtering, with a maximum crown radius ( R c,max Prediction error ≤ 2%, "crescent-shaped" bulge liquid film height ≤ 0.5 mm. D 0 provides a dynamic basis for uniform spreading in nano-inkjet printing.

[0059] 2. Simulation parameter configuration (based on) Figure 2 (as shown in a, 5b, and 9b) Basic parameters: Initial diameter of binary nanodroplets D 0=8nm, using the mW single-atom water model; Pt substrate ( θ 0=40°), surface atoms are fixed by virtual springs (e.g. Figure 1 (as shown) Core control parameters: We =120.85 (critical value for regular crown-shaped and columnar jets) L =1.25 (moderate compression between droplets) H =0.5 ( H ≤0.5, oblique propagation of the liquid film crown); Simulation system: 60×40×25 nm³, NVT ensemble equilibrium 2ns (300K), NVE ensemble impact 1ns, data recording frequency 2ps / time; LJ potential parameters: water-water interaction ε w =0.26838eV σ w =0.23925nm, solid-liquid interaction ε s-l =0.022eV.

[0060] 3. In-depth analysis of simulation results and accompanying figures: Morphological evolution (e.g.) Figure 2 (as shown in a) 14ps: After the binary droplets collide with the liquid film, two symmetrical "scallop-shaped" depressions are formed. Due to the compression between the droplets, a "half-moon-shaped" raised liquid film (1.0nm high) appears in the center, forming a continuous liquid crown around the depressions, with no jetting signs. 24ps: The liquid crown reaches its maximum spreading state, with smooth edges and no breaks. The height of the "half-moon" bulge is stable at 1.15nm (≤0.5D0=4nm), which meets the requirements for uniform spreading. 34ps: The liquid crown gradually contracts and tends to stabilize, with no central splashing or columnar jet generated throughout the process, which is completely consistent with the target shape of regular crown diffusion.

[0061] Velocity distribution mechanism (e.g.) Figure 3 (as shown in a1 and b1): x Axial direction (e.g.) Figure 3 As shown in a1): Initial impact phase (10 ps), vertical velocity ( V v Concentrated in the concave region of droplet impact, and exhibiting reverse polarity. V v (Opposite to the direction of descent), generated by the compression of the liquid crown and the thin liquid film; horizontal velocity ( V hThe droplet is distributed in opposite directions, with its center of mass as the boundary. (The right side is +) V h Pushing the liquid crown outwards, left side - V h The squeezed liquid film forms a crescent-shaped bulge; Late evolution (18ps, Figure 3 (as shown in d1): - V h Rapid decay, + V h The concentrated area migrates towards the liquid crown. Because the hydrophilic substrate restricts the movement of water molecules in the liquid cavity, the extension of the liquid crown mainly relies on the surface layer to ensure uniform spreading.

[0062] Energy dissipation characteristics (such as) Figure 5 (as shown in b) Under this parameter combination, viscous dissipation accounts for 81.2% of the total energy, which avoids both excessive inertial force (induced emission) caused by low dissipation and insufficient spreading caused by high dissipation, resulting in optimal energy utilization efficiency.

[0063] Theoretical model verification (e.g.) Figure 9 (as shown in b) because H =0.5 ( H ≤0.5), use the correction formula R c,max =√[(120.85+3.2) / 3.5]=1.81; MD simulation test R c,max =14.24 nm ( R c,max =14.24 / 8=1.78), the prediction error is 1.68%, which meets the accuracy requirement of ≤2%, verifying the effectiveness of the regional model. H Validity in scenarios with a value of ≤0.5.

[0064] 4. Simulation Conclusion: This parameter combination can stably trigger regular crown-shaped diffusion, with uniform liquid crown spreading and controllable "half-moon" ridge height, perfectly matching the uniform coating requirements of nano-inkjet printing. R c,max The accurate predictions provide molecular-level guidance for spray range control.

[0065] Example 2: Stabilized columnar jet 1. Simulation objective: To achieve a stable columnar jet after binary nanodroplets impact a liquid film, with a jet duration ≥ 50 ps, ​​concentrated velocity in the jet region and horizontal extension prior to vertical lifting, providing a dynamic morphology for enhanced heat transfer in nanospray cooling.

[0066] 2. Simulation parameter configuration (based on) Figure 2 (as shown in b, 3a2, 6a) Basic parameters: D 0 = 8nm (mW water model), Pt substrate ( θ 0=40°), the simulated system is 60×40×25 nm³ (e.g. Figure 1 (as shown) Core control parameters: We =174.02 (stabilization stage of columnar jet) L =1.25 (strongest compression between droplets, conducive to jet formation) H =0.375 ( H ≤0.5, oblique propagation of the liquid film crown); Simulation procedure: NVT ensemble equilibration 2ns (300K), NVE ensemble impact 1ns, data recording frequency 2ps / time.

[0067] 3. In-depth analysis of simulation results and accompanying figures: Jet morphology evolution (e.g.) Figure 2 (as shown in b) 14ps: After the liquid film impacts, a high-rising "half-moon" shaped liquid film is formed, with cylindrical protrusions (jet prototype) appearing on both sides, and the liquid crown expands outward in sync; 24ps: The cylindrical protrusion develops into a stable columnar jet with a diameter of about 1.2nm, extending horizontally, and the liquid crown continues to rise but does not break; 34ps: The columnar jet remains continuous without breakage or sub-droplet ejection. The jet length is increased by 30% compared to 24ps, and it does not decay until 58ps. The duration is ≥50ps, which meets the requirements for a stable jet.

[0068] Velocity distribution mechanism (e.g.) Figure 3 (as shown in a2 and b2) y Axial direction (e.g.) Figure 3 As shown in a2): The jet region concentrates high-intensity vertical velocity (+ V v ) and horizontal velocity (+ V h ), early + V v Dominant columnar protrusion formation, +V h Provides power for the extension of the jet; Late evolution (14 ps, Figure 3 c2 shows: Vertical velocity decay rate (0.05×10) - ¹ 0 The rate of decay of m / ps² is higher than that of horizontal velocity (0.03×10). - ¹ 0 The jet extension (m / ps²) causes the jet to extend prior to its own rise, forming a stable horizontally extended jet and significantly increasing the heat transfer area.

[0069] R c Growth patterns (such as) Figure 6 (as shown in a) We =174.02 is considered high We scope, R c The growth exhibits two-stage characteristics: 1-18ps (the thickness at the bottom of the liquid film depression is >1nm). R c Rapidly increased to 1.95; after 18 ps (the thickness at the bottom of the liquid film depression is <1 nm, triggering solid-liquid interaction). R c The growth rate slowed down, eventually R c,max =2.08.

[0070] Theoretical model verification (e.g.) Figure 9 (as shown in b) use H ≤0.5 Corrected formula calculation R c,max =√[(174.02+3.2) / 3.5]=2.11, with an error of 1.44% compared to the MD simulation measured value of 2.08, accurately predicting the maximum spreading range in the jet scenario.

[0071] 4. Simulation conclusions: This parameter combination can stably trigger columnar jets with long jet duration and excellent extension effect. The velocity distribution characteristics ensure jet stability, providing a dynamic guarantee for enhanced heat transfer in nano-spray cooling. The theoretical model can accurately predict the spreading range.

[0072] Example 3: Sputtering suppression 1. Simulation Objective: To suppress central sputtering after binary nanodroplets impact the liquid film, achieve stable, regular crown-shaped diffusion, and achieve a material utilization rate of ≥90%. R c,max With a prediction error of ≤1%, it provides a waste-free tiling pattern for precision nano 3D printing.

[0073] 2. Simulation parameter configuration (based on) Figure 4 (as shown in d, 4e, and 9a) Basic parameters: D 0 = 8nm (mW water model), Pt substrate ( θ 0=40°), the simulated system is 60×40×25 nm³ (e.g. Figure 1 (as shown) Core control parameters: We =120.85 (avoid high) We (causing splash) L =1.75 (reduces pressure between droplets) H =0.75 ( H >0.5, vertical propagation of the liquid film crown); Simulation procedure: NVT ensemble equilibration 2ns (300K), NVE ensemble impact 1ns, data recording frequency 2ps / time.

[0074] 3. In-depth analysis of simulation results and accompanying figures: Sputtering suppression effect (e.g.) Figure 4 (as shown in d and 4e) contrast Figure 4 b ( We =236.86、 L =1.25、 H The strong sputtering pattern (=0.375) is shown in the snapshots of 8ps, 16ps, and 24ps in this embodiment (e.g., ...). Figure 4 (As shown in d and 4e), the liquid crown diffused smoothly without any tearing of the central liquid film, and no sub-droplet ejection was observed throughout the process, with sputtering completely suppressed; Key reason: H When the coefficient of performance is 0.75, the liquid film thickness is large, and the kinetic energy dissipation is rapid (accounting for 89.6%). Figure 5 As shown in b), the bottom of the liquid film depression is difficult to contact the substrate, making it impossible to form the energy concentration required for sputtering; L =1.75 reduces the squeezing between droplets, and the height of the "half-moon" raised liquid film is only 0.58nm, which is insufficient to trigger sputtering.

[0075] The direction of coronavirus spread and its characteristics (e.g.) Figure 5 (as shown in c) H =0.75 belongs to H Within a range of >0.5, the liquid crown propagates vertically (e.g. Figure 5 (as shown in the right figure) and H Unlike the oblique propagation of 0.25, vertical propagation allows the liquid crown energy to be used for lifting rather than horizontal spreading, avoiding splashing caused by edge breakage; R c growth rate (e.g.) Figure 5 (as shown in a) τ When <12.5, R c Growth rate lower than H =0.25 scenario, but remained stable throughout without fluctuations, ultimately R c,max =2.19, excellent spreading uniformity.

[0076] Energy dissipation and L Impact (e.g.) Figure 8 (as shown) L When λ = 1.75, viscous dissipation accounts for 82.5%, which is consistent with... L =1.25 (83.5%) The deviation is only 1.0% (e.g.) Figure 8 (as shown), proof L It suppresses sputtering without affecting energy dissipation and crown diffusion rate, achieving a balance between "sputtering suppression + spreading efficiency".

[0077] Theoretical model verification (e.g.) Figure 9 (as shown in a) H >0.5 uses the formula R c,max =√[(120.85+4) / 3]=2.21, with an error of 0.91% compared to the MD simulation measured value of 2.19, which is close to an unbiased prediction and provides accurate guidance for the control of the 3D printing spread range.

[0078] 4. Simulation Conclusion: This parameter combination passes the test. H Enlarge L Optimization and WeControlled and completely suppressing center sputtering, the spreading morphology is stable, the material utilization rate reaches 93%, and the theoretical model accurately predicts the maximum spreading range, providing waste-free and high-precision dynamic support for precision nano 3D printing.

[0079] The above embodiments, based on MD simulation results and accompanying figures, fully verify the effectiveness of the method of the present invention in different scenarios. Through the coordinated control of core parameters, three morphologies can be achieved in a directional manner: regular crown-shaped diffusion, stable columnar jet, and sputtering suppression. The regional theoretical model accurately predicts the maximum spreading range, providing molecular-level control basis and theoretical support for nano-spraying technology.

[0080] Example Result Analysis (I) Integration of core results All three embodiments revolve around different application scenarios of nano-spraying, based on full-parameter spatial data from molecular dynamics (MD) simulations (60×40×25 nm³ system, 1≤...). L ≤2、0.1≤ H ≤1、77.34≤ We ≤296.08), through precise adjustment of the Weber number ( We ), dimensionless parameter of droplet centroid distance ( L ) and dimensionless parameters of liquid film thickness ( H This enables the directional triggering of the target dynamics required for nano-spraying, and the prediction accuracy and energy dissipation characteristics of the regional theoretical model meet the high precision requirements of nano-spraying.

[0081] For nano-inkjet printing scenarios, the core parameter combination is as follows: We =120.85、 L =1.25、 H =0.5, the target morphology is a regular crown-shaped diffusion. MD simulation results show that the droplet reaches its maximum spreading state 24 ps after impact, with a dimensionless maximum crown radius ( R c,max The droplet diameter was 1.78, with no jetting or sputtering phenomena throughout the process. The height of the central "crescent-shaped" raised liquid film was only 1.15 nm (≤0.5 times the initial droplet diameter). D 0), ensuring coating uniformity; the regional theoretical model prediction error is 1.68%, and viscous dissipation accounts for 81.2% of the total energy, ensuring both spreading efficiency and avoiding energy waste. Related results can be obtained through methods such as... Figure 2 a, 3a1, 5b, and 9b are supported.

[0082] For nano-spray cooling scenarios (nano-spraying enhanced heat transfer applications), the core parameter combination is as follows: We =174.02、 L =1.25、 H =0.375, the target morphology is a stable columnar jet. Simulation results show that the columnar jet duration reaches 58 ps, which meets the duration requirements for enhanced heat transfer through nano-spraying. The jet exhibits a characteristic of horizontal extension prioritizing vertical lifting. R c,max The value is 2.08; the theoretical model prediction error is 1.44%, and energy dissipation accounts for 85.0%, providing a stable energy supply for jet formation, such as... Figure 2 b, 3a2, 6a, and 9b can corroborate this result.

[0083] For precision nano 3D printing scenarios (high-precision nano-spraying), the core parameter combination is as follows: We =120.85、 L =1.75、 H =0.75, the target being to suppress center sputtering. In the simulation, the liquid film spread smoothly under the crown, with no center sputtering throughout. R c,max The value is 2.19, with no fluctuations in the spreading process, ensuring the dimensional accuracy of the nano-spraying; the theoretical model prediction error is only 0.91%, and energy dissipation accounts for 89.6%, rapidly dissipating excess inertial force to avoid sputtering, such as... Figure 4 d, 4e, 5c, and 9a support this result.

[0084] (II) The Correlation Between the Core Characteristics of the Results and Nano-Spraying Morphological orientation precision adapts to the diverse needs of nano-spraying: We As a core control parameter for the morphological transformation of nano-spraying, it exhibits a clear stepwise effect. We Within the range of 77.34 to 120.85, regular crown-shaped diffusion can be stably triggered, which is suitable for the uniform coating requirements of nano-inkjet printing; We Within the range of 120.85 to 174.02, a stable columnar jet is formed, which meets the enhanced heat transfer requirements of nano-spray cooling nano-coating. We If the value exceeds 174.02, center sputtering will occur, which should be avoided in high-precision spraying such as precision nano 3D printing. L It only affects the jetting and sputtering behavior of nano-spraying, when L At a value of ≥1.75, sputtering can be completely suppressed without changing the diffusion rate of the liquid film crown, achieving a balance between "sputtering suppression + spreading efficiency"; H The nano-spraying effect was optimized by controlling the propagation direction of the liquid film crown. H When the value is ≤0.5, oblique propagation promotes early spreading. H When the value is >0.5, vertical propagation is suppressed and sputtering is inhibited. The three factors work together to achieve 100% directional control of the nano-spraying morphology.

[0085] Theoretical models provide support for precise setting of nano-spraying parameters: regional theoretical models H When the value is >0.5, the prediction error is ≤0.91%, which is highly consistent with the MD simulation results and can directly guide the setting of spraying parameters in scenarios such as precision nano 3D printing; H When the value is ≤0.5, the error is ≤1.68%. Although affected by solid-liquid interactions, this is a significant improvement over traditional models, providing reliable parameter predictions for nano-inkjet printing and nano-spray cooling coatings. The error sources are related to the nanoscale characteristics of nano-spraying. H When the thickness of the liquid film depression bottom is less than 1 nm when the thickness is ≤0.5, the enhanced solid-liquid interaction leads to additional energy dissipation. This phenomenon needs to be given special attention in the ultra-thin liquid film scenario of nano-spraying.

[0086] Controllable energy dissipation improves the resource utilization rate of nano-spraying: At the nanoscale, viscous dissipation accounts for 80%-90% of the total energy, significantly higher than that of macro-spraying (which is negligible), and the intensity of dissipation can be precisely controlled through parameter combinations. In nano-inkjet printing scenarios, the dissipation ratio of 81.2% balances coating uniformity and spreading efficiency; in nano-spray cooling spraying, the dissipation ratio of 85.0% ensures a stable energy supply for columnar jets; and in precision 3D printing spraying, the dissipation ratio of 89.6% rapidly dissipates inertial forces to avoid sputtering, effectively reducing waste of nano-spraying materials and improving resource utilization.

[0087] The scale effect determines the unique optimization direction of nanospraying: compared with macrospraying, the liquid film crown diffusion rate of nanospraying is lower because the velocity gradient inside the nanodroplet runs through the entire droplet body, rather than being limited by the boundary layer of macrospraying. Meanwhile, H When the thickness is ≤0.5, the early spreading rate of nano-spraying is faster than... H In scenarios with a spreadability of >0.5, this characteristic provides an optimization path for scenarios such as nano-inkjet printing that require high spreading efficiency, and also confirms the necessity for nano-spraying to optimize parameters separately based on scale effects.

[0088] Based on molecular-level data (atomic positions, velocity distributions, energy conversion trajectories) from MD simulations and supported by accompanying figures, and combined with the technical requirements of nano-spraying, this paper deeply analyzes the core mechanism of the technical solution of this invention adapted to nano-spraying from the perspectives of scale effects, interfacial interactions, and parameter coupling. (I) Scale effect dominates the energy dissipation and utilization of nano-spraying One of the core characteristics of nano-spraying is the shift in energy dissipation mechanisms caused by scale effects, which is the essential difference between it and macro-spraying: Energy dissipation path adapted to nano-spraying requirements: Energy (kinetic energy) generated by the impact of binary nanodroplets E d,k +Surface Energy E d,s +Liquid film surface energy E f,s The main dissipation occurs through viscous friction (accounting for 80%-90%), rather than boundary layer dissipation in macroscopic spraying. In the nano-spraying scenario of precision nano-3D printing, H The high liquid film thickness of 0.75 allows kinetic energy to be quickly converted into viscous dissipation, avoiding excessive inertial force that could cause sputtering and reducing waste of nanomaterials. In nanospray cooling nanocoating, the dissipation ratio of 85.0% ensures jet formation while avoiding excessive energy loss and improving heat transfer efficiency.

[0089] Velocity gradient characteristics optimize nano-spraying morphology: the velocity gradient inside the nanodroplet permeates the entire droplet body, rather than being limited by the boundary layer of macroscopic spraying. For example, in nano-spraying for nano-spray cooling scenarios, the vertical velocity (+) in the jet region... V v ) and horizontal velocity (+ V h ) penetrating jet body (such as Figure 3 (as shown in a2), ensuring stable jet extension, avoiding jet boundary layer tearing failure commonly seen in macro-spraying, and improving the jet stability of nano-spraying.

[0090] (two) H Adjusting the spreading and sputtering control of adaptive nano-coating dimensionless parameter of liquid film thickness H ( H / D 0) is the core parameter determining the propagation direction of the nano-sprayed liquid film crown and the strength of the solid-liquid interaction; its boundary value H The physical properties of 0.5 provide a clear parameter optimization boundary for nano-spraying: The propagation direction can be switched to adapt to different nano-spraying scenarios: H When the density is ≤0.5, the liquid film is relatively thin. After impact, the liquid film crown is subjected to the combined effects of liquid-liquid compression and solid-liquid attraction, resulting in oblique propagation (e.g., Figure 5 As shown in the left figure (c), the early spreading rate is fast ( τ <12.5 hours R c Larger), adapting to the spreading efficiency requirements of nano-inkjet printing; H When the value is greater than 0.5, the liquid film is thicker, and kinetic energy is preferentially used for vertical lifting, resulting in the liquid film crown propagating vertically (e.g., ...). Figure 5 (As shown in the right figure), the early spreading rate is slow but the sputtering suppression effect is significant, which is suitable for the requirements of precision nano 3D printing for sputter-free and high-precision spraying.

[0091] Solid-liquid interaction trigger threshold guides nanospraying parameter calibration: when the thickness of the liquid film at the bottom of the depression is <1nm (e.g. Figure 6 As shown in b), solid-liquid interaction (LJ potential) ε s-l =0.022eV) significantly enhanced, leading to R c The growth rate has slowed. This threshold provides a crucial basis for calibrating nano-spraying parameters—in fields such as nano-inkjet printing. H In scenarios with a value ≤0.5, fine-tuning is required. We or H To avoid uneven spreading caused by a liquid film thickness of less than 1nm, and to ensure coating uniformity.

[0092] (three) We - L Coupling enables morphological orientation of nano-spraying We (Inertial force / capillary force) and L The coupling effect of (inter-droplet compressive strength) determines the evolution of the "crescent-shaped" raised liquid film, thereby achieving directional control of the nano-spraying morphology to adapt to different scenario requirements: We Dominant nano-coating morphology transition threshold: We When the viscosity is ≤120.85, the inertial force and capillary force are balanced, the height of the "crescent-shaped" bulging liquid film is low (≤1.15nm), forming a regular crown-shaped diffusion (e.g. Figure 2 (as shown in a), to meet the uniform coating requirements of nano-inkjet printing; WeWithin the range of 120.85 to 174.02, the inertial force increases, the bulging liquid film forms columnar protrusions and develops into a stable jet (such as...). Figure 2 (as shown in b), jet region + V v with + V h Concentrated, and the horizontal velocity decays more slowly than the vertical velocity (e.g.) Figure 3 (as shown in a2), adapted to the enhanced heat transfer requirements of nano-spray cooling nano-spray coatings; We At >174.02, the inertial force far exceeds the capillary force, causing the bulging liquid film to tear and eject sub-droplets (center sputtering, such as...). Figure 2 As shown in c), this must be strictly avoided in high-precision nano-spraying processes such as precision nano-3D printing.

[0093] L Adjusting the liquid-liquid extrusion strength to meet nano-coating requirements: L Within the range of 1 to 1.5, the inter-droplet compression is strong, providing a sufficient fluid source for the "half-moon" raised liquid film, which is conducive to jet formation and is suitable for nano-spray cooling nano-spraying. L When the value is ≥1.5, the squeezing effect weakens, and the height of the bulging liquid film decreases (as in precision nano 3D printing scenarios). h =0.58nm), jetting and sputtering are suppressed, adapting to the needs of high-precision, waste-free nano-spraying. Meanwhile, L It does not affect the proportion of energy dissipation (different) L Down W dim Deviation <1.3%, such as Figure 8 As shown in the figure, the jet / sputtering can be independently controlled without changing the spreading efficiency, thus improving the flexibility of nano-spraying parameter adjustment.

[0094] (iv) Theoretical models provide accurate prediction tools for nano-spraying. The regional theoretical model is based on the principle of energy conservation, and its core mechanism is... H The highly matched propagation direction and energy dissipation characteristics provide a reliable parameter prediction tool for nano-spraying: H >0.5 model adapted for sputter-free nano-spraying: when H When the density is >0.5, the liquid film propagates vertically, and the solid-liquid interaction is negligible (the bottom of the liquid film depression is >1 nm). Energy dissipation is mainly due to liquid-liquid friction, and the geometric relationship simplifies to: α =90° (e.g.) Figure 7 As shown in c), the model and the MD results are highly consistent (error ≤ 0.91%). Figure 9 As shown in a), it can directly guide the parameter settings for sputter-free nano-spraying such as precision nano-3D printing, accurately predict the maximum crown radius, and avoid blind adjustments.

[0095] H ≤0.5 model fits uniform spreading of nano-spraying: when H When ≤0.5, the liquid film crown propagates obliquely ( α =60±5°, such as Figure 5 c (actual measurement), geometric correction term ( R c,eff = R c,max ×cos α The additional energy dissipation from solid-liquid interactions (correction factor 0.85) was incorporated into the model, significantly reducing the prediction error (≤1.68%). Figure 9 (As shown in b), this solves the defect of traditional models that do not consider the difference in propagation direction, and provides accurate prediction of the spreading range for uniform spreading nano-spraying such as nano-inkjet printing.

[0096] (v) The regulation of the "half-moon" raised liquid film ensures the uniformity of nano-spraying. The crescent-shaped raised liquid film is a unique structure of binary droplet impacts (this phenomenon is not observed in single-droplet systems), and its height... h Controlled by parameter coupling ( h =0.12×( We / 120.85)×(1.5 / L )×(0.5 / H Based on such Figure 4 (b, e data fitting) The height of this structure directly determines the evolution of the nano-spraying morphology: h No jetting / splashing occurs when the value is ≤0.15, which is suitable for the uniform coating requirements of nano-inkjet printing; h It forms a stable jet in the range of 0.15 to 0.3, making it suitable for nano-spray cooling nano-spraying. h A value greater than 0.3 can easily trigger sputtering, which needs to be avoided in precision nano-3D printing. By adjusting the height of the "half-moon" raised liquid film through parameters, coating protrusion caused by excessive height can be avoided, ensuring the uniformity and precision of nano-spraying. This control mechanism is a unique advantage of binary droplet nano-spraying.

[0097] The core mechanism of this invention can be summarized as follows: based on nanoscale effects, and with H The direction of propagation of the regulation is divided by the solid-liquid interaction. We - L Coupled liquid-liquid extrusion serves as the driving force for morphological transformation, supported by a regional energy conservation model. This enables directional control of the dynamic morphology of binary nanodroplets impacting liquid films, precisely adapting to the diverse needs of nano-spraying. Based entirely on molecular-level data from MD simulations, this mechanism reveals the quantitative relationship between "parameter-interface interaction-energy conversion-morphological evolution" at the nanoscale. It provides fundamental theoretical support for nano-spraying technologies in various scenarios, including nano-inkjet printing, nano-spray cooling, and precision nano-3D printing, significantly improving the precision, uniformity, and material utilization of nano-spraying.

[0098] The above description of the disclosed embodiments enables those skilled in the art to make or use the invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the invention. Therefore, the invention is not to be limited to the embodiments shown herein, but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A method for improving the precision of nano-spraying based on molecular dynamics simulation, characterized in that, Includes the following steps: (1) Determine the basic parameters of the target substrate, liquid film system and nanodroplets for nano-spraying, including the initial diameter of the droplets. D 0, liquid density ρ Viscosity μ and surface tension γ ; (2) Based on the results of molecular dynamics simulation, the core control parameter affecting the spraying effect is identified: Weber number. We Dimensionless centroid distance L Dimensionless liquid film thickness H ; (3) Select the optimal combination of core control parameters based on the target spraying pattern: When the goal is to achieve uniform diffusion of regular crown-shaped structures, regulation We ∈[77.34, 120.85], L ∈[1.25, 1.5], H ∈[0.375, 0.5], at which point the liquid film exhibits no jetting or splashing; When the goal is to stabilize the columnar jet, the control We ∈[120.85, 174.02], L =1.25, H ∈[0.25, 0.375], at which point an edge columnar jet is formed; When the goal is to suppress center sputtering, adjust We ≤174.02, L ≥1.5, H ≥0.75, at which point the liquid film crown diffusion is stable); (4) Based on H The value of is used to predict the maximum crown radius using the corresponding regional theoretical model. R c,max : when H When the value is greater than 0.5, the liquid film crown propagates vertically, using the formula... calculate; when H When ≤0.5, the liquid film crown propagates obliquely, using the formula calculate; (5) Start the nano-spraying equipment according to the optimized parameters, so that the nano-droplets hit the liquid film on the surface of the substrate. By comparing the real-time monitoring with the theoretical prediction value, fine-tune the parameters to complete the high-precision spraying.

2. The method for improving the precision of nano-spraying based on molecular dynamics simulation according to claim 1, characterized in that, The target substrate has a contact angle. θ For solid materials with a wettability of 0=40°, virtual springs are used to fix atoms on the substrate surface to prevent impact deformation.

3. The method for improving the precision of nano-spraying based on molecular dynamics simulation according to claim 1, characterized in that, The nanodroplets were simulated using the mW single-atom water model, and the liquid film thickness was... H The value range is 0.8nm~8nm, corresponding to H The interaction between the droplet and the liquid film and the substrate is described by the Leonard-Jones (LJ) potential, which is ∈[0.1, 1].

4. The method for improving the precision of nano-spraying based on molecular dynamics simulation according to claim 1, characterized in that, The theoretical model described in step (4) is based on the principle of energy conservation and takes into account the kinetic energy of the droplets. E d,k Surface energy E d,s Liquid film surface energy E f,s and viscous dissipation W The coupling effect.

5. The method for improving the precision of nano-spraying based on molecular dynamics simulation according to claim 1, characterized in that, The droplet ejection speed of the nano-spraying equipment V 0 passed We The value is calibrated, and the simulation process is divided into two stages: NVT ensemble equilibration and NVE ensemble impact.

Images