Loading [MathJax]/jax/output/SVG/jax.js
Research article Special Issues

Biochemical and biophysical mechanisms underlying the heart and the brain dialog

  • In this paper, within the so called “neuro-visceral integration model”, it is reviewed the ability of the heart to secrete numerous endocrine mediators, neurotransmitters and substances that regulate the immune function with repercussions on the central nervous system. The heart would also seem to be able to process various information independently, influencing the brain work through the “intrinsic cardiac nervous system” and baroreceptor pathways. In reviewing this matter, further physical mechanisms are also described, including mechanical contractions and deformations, that are involved in the “heart-brain symphony” based on intra-cardiac formation and propagation of blood vortices coupled to electrical signals. The relevance of the role of vorticity of the blood flow in the molecular dynamics and physiological activity is stressed. By resorting to some conceptual and formal aspects of the dissipative quantum model of brain, mechanisms such as the spontaneous breakdown of symmetry in many-body physics, the dynamical formation of long range correlations and their associated Nambu-Goldstone quanta, coherent states and fractal self-similarity are discussed with reference to the heart-brain dialog. Our discussion supports the view that the heart role is more than the one of a muscle responsible of the blood flow. Further related topics such as the formation of aneurysms and vein varices which in our modeling seem to be related to the weakening or loss of vorticity of the blood flow, the role of the recently discovered fluid-filled interstitial structure and the complex network of thick collagen bundles are finally briefly mentioned in the concluding remarks.

    Citation: C. Dal Lin, M. Falanga, E. De Lauro, S. De Martino, G. Vitiello. Biochemical and biophysical mechanisms underlying the heart and the brain dialog[J]. AIMS Biophysics, 2021, 8(1): 1-33. doi: 10.3934/biophy.2021001

    Related Papers:

    [1] Ping Zhu . Dynamics of the positive almost periodic solution to a class of recruitment delayed model on time scales. AIMS Mathematics, 2023, 8(3): 7292-7309. doi: 10.3934/math.2023367
    [2] Yan Yan . Multiplicity of positive periodic solutions for a discrete impulsive blood cell production model. AIMS Mathematics, 2023, 8(11): 26515-26531. doi: 10.3934/math.20231354
    [3] Guilin Tang, Ning Li . Chaotic behavior and controlling chaos in a fast-slow plankton-fish model. AIMS Mathematics, 2024, 9(6): 14376-14404. doi: 10.3934/math.2024699
    [4] Xiaofang Meng, Yongkun Li . Pseudo almost periodic solutions for quaternion-valued high-order Hopfield neural networks with time-varying delays and leakage delays on time scales. AIMS Mathematics, 2021, 6(9): 10070-10091. doi: 10.3934/math.2021585
    [5] Yi Zhu, Minghui Jiang . Prescribed-time synchronization of inertial memristive neural networks with time-varying delays. AIMS Mathematics, 2025, 10(4): 9900-9916. doi: 10.3934/math.2025453
    [6] Yun Xin, Hao Wang . Positive periodic solution for third-order singular neutral differential equation with time-dependent delay. AIMS Mathematics, 2020, 5(6): 7234-7251. doi: 10.3934/math.2020462
    [7] Călin-Adrian Popa . Synchronization of Clifford-valued neural networks with leakage, time-varying, and infinite distributed delays on time scales. AIMS Mathematics, 2024, 9(7): 18796-18823. doi: 10.3934/math.2024915
    [8] Jing Ge, Xiaoliang Li, Bo Du, Famei Zheng . Almost periodic solutions of neutral-type differential system on time scales and applications to population models. AIMS Mathematics, 2025, 10(2): 3866-3883. doi: 10.3934/math.2025180
    [9] Yanshou Dong, Junfang Zhao, Xu Miao, Ming Kang . Piecewise pseudo almost periodic solutions of interval general BAM neural networks with mixed time-varying delays and impulsive perturbations. AIMS Mathematics, 2023, 8(9): 21828-21855. doi: 10.3934/math.20231113
    [10] Liye Wang, Wenlong Wang, Ruizhi Yang . Stability switch and Hopf bifurcations for a diffusive plankton system with nonlocal competition and toxic effect. AIMS Mathematics, 2023, 8(4): 9716-9739. doi: 10.3934/math.2023490
  • In this paper, within the so called “neuro-visceral integration model”, it is reviewed the ability of the heart to secrete numerous endocrine mediators, neurotransmitters and substances that regulate the immune function with repercussions on the central nervous system. The heart would also seem to be able to process various information independently, influencing the brain work through the “intrinsic cardiac nervous system” and baroreceptor pathways. In reviewing this matter, further physical mechanisms are also described, including mechanical contractions and deformations, that are involved in the “heart-brain symphony” based on intra-cardiac formation and propagation of blood vortices coupled to electrical signals. The relevance of the role of vorticity of the blood flow in the molecular dynamics and physiological activity is stressed. By resorting to some conceptual and formal aspects of the dissipative quantum model of brain, mechanisms such as the spontaneous breakdown of symmetry in many-body physics, the dynamical formation of long range correlations and their associated Nambu-Goldstone quanta, coherent states and fractal self-similarity are discussed with reference to the heart-brain dialog. Our discussion supports the view that the heart role is more than the one of a muscle responsible of the blood flow. Further related topics such as the formation of aneurysms and vein varices which in our modeling seem to be related to the weakening or loss of vorticity of the blood flow, the role of the recently discovered fluid-filled interstitial structure and the complex network of thick collagen bundles are finally briefly mentioned in the concluding remarks.


    Euler-Lagrange fluid-particle simulation has been extensively used as a surrogate for experiments in respiratory aerosol dynamics and inhalation drug delivery. This approach's popularity is rooted in its high accuracy in resolving fine-scale flow details and tracking individual droplets/particles in a time-efficient manner and thus providing an estimation of the dosimetry of released pharmaceuticals [1]. However, an important issue often neglected in numerical simulations of pulmonary drug delivery was the momentum exchange from the discrete phase to the carrier fluid or the reciprocal interactions between aerosols and flow (i.e., two-way coupling, or 2WC). One-way coupling (1WC) was adopted in most previous inhalation dosimetry simulations based on the assumption of dilute aerosol concentrations, where the presence of the aerosols did not affect the carrier flow [2]. This assumption was valid only when aerosols had a sufficiently low concentration and low exiting velocities, such as those from a mesh nebulizer. One-way coupling was also valid for aerosols from a valved holding chamber (VHC) or spacer connected to a metered-dose inhaler (MDI) due to the mixing of the spray droplets within it, which greatly reduced droplets' speed and concentration [3].

    However, the local concentration of MDI spray aerosols can be high upon discharge close to the actuator orifice, which has a small diameter of 0.2–0.5 mm [4,5,6]. The pressurized propellant generates a high-speed jet flow from the orifice and atomizes the drug solution into micro-scale droplets, which develop into a spray plume by advancing along the streamwise direction and dispersing in transverse directions. Because of the large density ratio between the droplet and fluid (i.e., three orders of magnitude), the inertia and response time of the droplets differ notably from those of the carrier flow, which will induce large droplet-flow slip velocities. As a result, the momentum exchange from the spray droplets to the co-flow can be significant downstream of the actuator orifice and modify the flow behaviors there, which in turn will modify aerosol motions. These reciprocal interactions kick in from the onset of the actuation and occur throughout the actuation process (usually 0.2 s). Spatially, the 2WC effect is more significant near the orifice and persists till either the local mass loading or aerosol velocities become sufficiently low. Due to their higher inertia, aerosols often exhibit an integrative feature, i.e., the current aerosol distribution (or droplet location) depends on not only the instantaneous velocity distribution but also the prior histories of the droplets. By contrast, airflow can quickly respond to external changes and is predominantly determined by the instantaneous pressure field. As a result, the motions of individual droplets will be constantly modified by the carrier flow, and these modifications will accumulate throughout the time course, yielding even more different spray plume morphology and velocity than in the 1WC case.

    Many experimental and numerical studies have considered the 2WC effects on particle-laden flows in general. Early observations by Torobin and Gauvin demonstrated that particles modified the drag force and heat transfer rate on the pipe, which couldn't be explained unless the flow turbulence had been modulated by the entrained particles [7]. Questions also arose about whether, at what conditions, and how the turbulence would be reduced or enhanced by the dispersed particles. It was known that when the aerosol droplets are comparable in size to the smallest fluid length scale or larger, unsteady flows can be important in affecting the aerosol motions [8,9]. Gore and Crowe reviewed a wide spectrum of experimental studies on particle-modulated turbulence in jets and pipes at varying orientations, with a Reynolds number of 8000–105, volume fraction of 10-6–0.2, and density ratio of 10-3–2500 [10]. They first observed that particles smaller than 200 µm suppressed turbulence while particles larger than 200 µm enhanced turbulence. A non-dimensional parameter dp/lt (particle diameter to turbulence length scale) was then proposed, with turbulence attenuation at dp/lt < 0.1 and enhancement at dp/lt > 0.1 [10,11]. The percentage change in turbulence intensity varied from -50% to 400%, and exhibited a nonlinear relation to dp/lt [10,11]. Later, Elghobashi and Truesdell updated the criterion to dp/η, with η being the smallest turbulence length scale or Kolmogorov length scale [12].

    Note that the particle-induced turbulence modulation can reversely modify particle motions at varying degrees, depending on the volume fraction, mass fraction, and/or diameter of the particles [13]. Using direct numerical simulations (DNS), Monchaux & Dejoan showed that 2WC had a weak effect on the global flow statistics when the particle loading is low (φ = 1.5×10-5). They also reported that 2WC enabled the particles to drag the fluid down with them, which further enhanced the particle settling [14]. In a classic review, Elghobashi categorized the particle-laden turbulent flows into three regimes based on the particle volume fraction φ: one-way coupling (1WC) regime when φ = 10-8–10-6, two-way coupling (2WC) regime when φ = 10-6–10-3, and particle-particle collision regime when φ > 10-3 [15]. The 1WC and 2WC regimes had dilute suspensions while the particle-particle collision regime had dense suspension. In the 2WC regime (φ = 10-6–10-3), a sub-division was made based on the time scale ratio, τpη where τp was the particle residence time, and τη was the Kolmogorov time scale. For particles with τpη > 100, 2WC enhanced turbulence production, and for particles with τpη < 100, 2WC enhanced dissipation [15].

    Till now, studies of the 2WC effects on MDI spray flow and aerosol dynamics have not been reported. Considering the high droplet concentration close to the actuator orifice and the high-speed jet flow that generates small turbulence length scales, 2WC effects are expected to be significant in the orifice vicinity and decay from the orifice. It is unknown yet to what extent the spray plume will be affected and at what level it is affected in different directions (i.e., spray centerline vs. transverse directions) and under various conditions (i.e., mass loading, jet flow speed, and aerosol diameter). This study seeks to answer the above questions related to three popular MDIs, i.e., Ventolin, ProAir, and Qvar, which have different droplet size distributions and co-flow velocities. Specific aims include:

    1) Develop an LES-Lagrangian computational model for MDI spray plume development in open space and validate it with complimentary high-speed imaging of Ventolin sprays.

    2) Compare the spatial and temporal evolutions of the Ventolin spray plumes between 2WC and 1WC predictions.

    3) Evaluate the 2WC effects with mass loadings, droplet sizes, and inhaler types.

    4) Understand diameter-dependent spray plume behaviors by correlating to flow Kolmogorov length/time scales and respective droplet-flow length/time ratios.

    The results of this study will shed light on whether 2WC effects should be considered in MDI spray simulations and provide quantitative guidance in estimating the differences between 2WC and 1WC cases.

    The 2WC effects on MDI airflow and aerosol dynamics were numerically evaluated in an MDI-open-space model, as shown in Figure 1(a). The MDI was actuated, and the generated spray droplets were released into the open space. A cone-shaped space was adopted to account for the expansion of the spray plume and, at the same time, minimize the computational mesh size. An airspace geometry had a length of 0.5 m based on the observations that the maximal spray plume lengths from all three inhalers considered hereof were equal to or smaller than 0.5 m.

    Figure 1.  Study design and computational model: (a) computational domain consisting of an MDI and an open space with a length of 0.5 m, (b) model for the Ventolin inhaler geometry and polydisperse spray aerosols, (c) aerosol size distribution for ProAir and Qvar, and (d) multidomain computational mesh, with a zoomed view of refined meshes in the mouthpiece and close to the actuator orifice.

    Figure 1(b) shows the model for Ventolin, including both the geometrical model and the model for the initial airflow and spray aerosols at actuation. The inhaler geometry was developed in SolidWorks (Dassault Systèmes, Concord, MA) based on the Ventolin inhaler with the most essential details included. Examples included the exact size of the mouthpiece and the actuator orifice indented in the cylindrical reservoir with an orifice diameter of 0.3 mm (Figure 1(b)). The orifice diameter of 0.3 mm followed Gabrio et al. [16], who evaluated the plume impact forces from several prototype MDIs (including Ventolin) with two different orifice diameters of 0.32 and 0.48 mm. They also measured the throat deposition from model actuators with orifice diameters ranging from 0.29 to 0.56 mm and reported that the throat deposition was sensitive to orifice diameters up to approximately 0.4 mm [16]. In this study, only the prototype actuator diameter of 0.32 mm (rounded to 0.3 mm) was considered. The exiting velocity of the airflow and aerosol was 40 m/s, which was based on previous measurements and computations [17,18,19]. A polydisperse aerosol size distribution was adopted (Figure 1(b)) based on previous measurements, with the mean droplet size being 9.1 µm and the standard geometrical derivation (GSD) being 1.32 [17,19]. The inhaler actuation started from 0.1 s and lasted for 0.2 s. After that, no further airflow and aerosols were released from the orifice, while the aerosols that had been released between 0.1–0.3 s would continue to be simulated till their complete decay.

    An MDI is supplied with a canister that contains 8.5 g of formulations intended for 200 puffs [20,21]. Each actuation delivers 90–120 μg medications (i.e., 0.2–0.3% w/w), such as albuterol from Ventolin, albuterol sulfate from ProAir, and beclomethasone dipropionate from Qvar [22]. In addition to the 0.1–1.0% w/w medication, a formulation also contains 3–13% w/w ethanol to increase the drug solubility and propellent to atomize the drug formulation into micron-scaled droplets [23,24]. In this study, it was assumed that the propellant evaporated immediately and completely at the orifice, leaving the liquid droplets that contained only the medication and ethanol. As a result, the aerosol droplet mass per actuation was approximately 3.0 mg (8500 mg, 200 puffs, 7% medication plus ethanol). The computational model was first applied to simulate the Ventolin spray plume development and was validated against complimentary high-speed images. This validated case with 2WC was then used as the baseline case for later comparisons. Various spray properties were considered regarding their relative contributions to the 2WC effects. These included the presence of the propellant jet flow, the spray mass loading, and individual droplet sizes. All cases were simulated with both 2WC and 1WC to characterize the 2WC effects under different scenarios.

    To evaluate the influences of the propellant flow, the spray plume behaviors were simulated with a 40-m/s velocity for aerosols but a zero velocity for the airflow from the orifice during actuation (0.1–0.3 s) under both two-way and one-way coupling. To evaluate the mass loading effects, three other drug masses were considered, which were one-tenth, ten times, and one hundred times (i.e., 0.30, 30, and 300 mg) that of the baseline (3.0 mg). To study the droplet diameter effects, seven monodisperse aerosols with a dose of 3.0 mg were considered, with the aerosol sizes being 2, 3.5, 5, 7.5, 10, 15, and 20 µm. Two other inhalers, ProAir and Qvar, were considered in addition to Ventolin, with their aerosol size distribution shown in Figure 1(c). Based on the measurements of Liu et al., the exiting spray velocity was 26 m/s for ProAir and 23 m/s for Qvar [17]. A summary of the test cases is listed in Table 1, with the corresponding figure numbers of the results.

    Table 1.  Different categories in comparing the predicted spray plumes using two-way coupling (2WC) vs. one-way coupling (1WC).
    Baseline (Ventolin) Mass loading Droplet size Inhaler type
    Aim Ventolin spray plume simulation Mass loading effect Size-dependent spray plumes Inhaler-specific spray plumes
    Inputs Ventolin
    Polydisperse
    dp: 9.1 µm, GSD: 1.32 V0: 40, Up: 40 m/s
    Mose mass: 6 mg
    Dose mass:
    0.6, 6, 60,600 mg
    Monodisperse
    Droplet size: 2, 3.5, 5, 7.5, 10, 15, 20 µm
    ProAir
    dp: 5.3 µm
    V0 = Up: 26 m/s
    Qvar, dp: 2.9 µm
    V0 = Up: 23 m/s
    Outputs Spray evolution;
    Vortex;
    Velocity contour;
    Streamlines;
    Droplet dynamics.
    Vortex;
    Vel contour;
    Streamlines;
    Droplet dynamics.
    Vel contour;
    Droplet dynamics
    Kolmogorov scales: η, τη;
    Ratio: dp/η, τp/τη;
    (V-Up);
    Spray velocity.
    Vel contour;
    Vortex;
    Droplet dynamics;
    Streamlines.
    Figures 3–5 6 7–12 13, 14

     | Show Table
    DownLoad: CSV

    The high-speed camera used was Phantom VEO, with an acquisition rate of up to 11,000 frames per second (fps). In this study, the temporal and spatial evolutions of the Ventolin spray plume after actuation were recorded at 4000 fps in order to achieve the optimal quality of acquired images. A laser sheet (OXlaser, 488 nm, 100 mW) was implemented to enhance the spray-environment contrast by aligning the sheet with the spray plume centerline. To understand the spray plume development at the very beginning of the actuation, the mouthpiece was partially removed to reveal the orifice. Images of the spray from the orifice were thus directly captured without the blockage by the mouthpiece.

    The LES-WALE model was used to simulate the multi-regime flows because of its ability to accurately capture turbulent-laminar transitions [25]. Droplet dynamics were simulated with a Lagrangian approach. An in-house MATLAB code was used to generate MDI droplets at the orifice, which specified the droplet diameter, size distribution, exit velocities, and plume angle [26,27]. The 2WC between the airflow and aerosol droplets were considered through the force below,

    F=(18μCDReρpd2p24(upu)+Fother)˙mpΔt (1)

    where µ is the fluid viscosity, ρp is the particle density, dp is the particle diameter, Re is the relative Reynolds number (Re = ρ(upu)l/µ), CD is the drag coefficient, up is the particle velocity, u is the fluid velocity, p is the aerosol mass flow rate, Δt is the time step, and Fother is a collection of other interaction forces, which are zero in this study. The mutual momentum transfer between the fluid and aerosols was numerically considered in each control volume for a given instant.

    ANSYS Fluent 19.1 (Canonsburg, PA) was utilized to solve flow-particle governing equations. ANSYS ICEM CFD was utilized to generate the computational mesh (Figure 1(d)). The characteristic length scales in the MDI-airspace system differed by three orders of magnitude, i.e., 0.3-mm diameter for the actuator orifice, 2.54-cm width for the mouthpiece opening, and a 0.5-m length for the open space. To adequately resolve the flow details in all regions, multidomain meshes were generated by dividing the system into many zones, with each zone being filled with varying-sized cells according to the anticipated flow complexities. For instance, ultrafine meshes were generated close to the orifice and within the shear layers around the spray centerline, with the mesh density gradually decreasing from the centerline, which eventually became coarse in the far open space (Figure 1(d)). Grid independent studies were conducted by testing six meshes from 3.2 million to 12.0 million (Figure 1(d)). The parameter of interest was the airflow 3 cm from the mouthpiece at 0.5 s. The grid-independent mass-based DF was reached at 10.1 million. A total 2×105 seed droplets were released for each actuation during 0.1–0.3 s. To approximate the continuous releasing of sprays during actuation, a group of 5000 droplets was released every 5 ms for 40 times, as opposed to a single injection in most previous numerical studies [28,29,30]. For an applied mass of 3.0 mg, the mass flow rate for each droplet parcel in the injection file was calculated to be 3.0×10-9 kg/s based on a droplet parcel number of 5000 and a time step of 0.005 s per release for 40 releases. The mass loading of 0.30, 30, and 300 mg was obtained by scaling the baseline particle mass flow rate per parcel by a factor of 1/10, 10, and 100, respectively. A time step size of 5 ms was used for a flow time of 1.5 s, or a complete spray decay, which came first. One test case took 60 hours or so in a Ryzen 9 3960x workstation with 3.79 GHz frequency and 256G RAM.

    A deeper understanding of the fluid-aerosol interactions involves those at small scales, especially when the length and time scales are equivalent between the fluid and individual droplets. The Kolmogorov scales are the smallest in turbulent flows and describe the microscopic airflow behaviors being dominated by viscosity, such as how much energy is contained in the smallest eddies, how energy is transferred between small and large eddies, and how energy is dissipated by eddies [31]. The Kolmogorov length scale, η and time scale, τη of the airflow are modeled as [32],

    η=(v3ε)1/4;τη=(vε)1/2 (2)

    where ν is the molecular viscosity and ɛ is the dissipation rate per unit mass, as shown below.

    ε=2<SijSij>=v{uixj+ujxiuixj+ujxi} (3)

    The length scale of a particle or droplet is obviously its diameter, and the time scale is described by the residence time, i.e., the duration required to react to a local flow change,

    τp=ρpd2pρf18v=ρpd2p18μ;Stk=ρpd2pu018μl0=τpt0 (4)

    The particle/droplet Stokes number was also presented above, which is the ratio of the particle/droplet residence time over the characteristic macroscopic time scale of the flow (l0/u0). Here l0 and u0 are the macroscopic characteristic length and velocity scale in a local region. Considering the multiscale dimensions (three orders of magnitude differences) of the system in this study, the local characteristic scales also varied drastically, which further led to drastic changes in the particle/droplet Stokes number.

    The aerosol-fluid length scale ratio (dp/η) and time scale ratio (τp/τη) have been suggested as indexes to evaluate the effects of aerosol-to-fluid forces on the flow turbulence [15]. It has been demonstrated that the presence of aerosols suppressed turbulence when dp/η < 0.1 and enhanced turbulence when dp/η > 0.1 [10,11,12]. Both ratios were used to evaluate the spray plume morphology variations for different aerosol sizes.

    Figure 2(a) shows the recorded images of the spray plume development from Ventolin at two instants using the Phantom VEO high-speed camera at an acquisition rate of 4000 fps. Intensive turbulence and eddies were observed at the front of the spray plume and, to a lesser degree, at the transverse interface with the ambient air. Considering the plume front, an aerosol stripe developed and moved downward (red arrow, second panel, Figure 2(a)), which subsequently separated from the main plume and started to settle downward from the gravity (yellow arrow, third panel, Figure 2(a)). Simultaneously, the eddies at the transverse interface continued to grow by mixing with the ambient air due to large shear forces. Because of the front and shear forces, as well as the entrained mass of the air, the overall velocity of the spray plume decreased quickly, from 40 m/s at the orifice to nearly zero at 0.5 m. The morphology of the spray plume exhibited a highly dynamic nature even after the spray plume reached its maximal length, driven both by the decaying vortices and gravity.

    Figure 2.  High-speed imaging: (a) temporal development of the spray plume from Ventolin after actuation, and (b) spray actuation from the orifice viewed from a cut-open mouthpiece.

    Considering that the orifice was blocked by the mouthpiece, parts of the mouthpiece were removed, and the emanation process of the spray from the orifice was captured using the high-speed camera at 4000 fps (Figure 2(b)). It can be clearly seen that the spray plume exited the orifice as a very thin jet with a cross-section diameter comparable to that of the orifice (first recorded image, Figure 2(b)). The plume grew quickly in the axial direction, and at a slower speed, in the transverse direction as well (second and third recorded images, Figure 2(b)). However, no apparent recirculation or vortices were noticed close to the orifice and within the mouthpiece. These observations were used to validate the complimentary computational results in the later section.

    Figure 3 shows the temporal and spatial evolution of the spray plume at different instants under 2WC. There were three stages during the plume development. The first stage occurs during 0.1–0.3 s (actuation, upper row of Figure 3), with different dynamics between the carrier flow and spray droplets. Due to the large density and high inertia, droplets traveled a much longer distance (~ 13 mm) in the first 5 ms than the airflow, with a vortex ring forming close to the mouthpiece exit (red arrow in Figure 2(a)). However, due to the drag force from the ambient air, the spray plume slowed down quickly, with further advancement of 7.5 mm at the next 15 ms (from 0.105 to 0.120 s) and 3 mm at the next 80 ms (from 0.120 to 0.200 s), as shown in Figure 2(a)(c). By contrast, the coherent flow structures developed steadily. The single vortex ring generated from the mouthpiece at 0.105 s (red arrow in Figure 2(a)) moved forward, while new vortex rings were constantly generated at the mouthpiece (black and green arrows in Figure 2(b)). The green ring caught up with the red ring at 0.200 s and merged together, while younger rings (pink and black arrows, Figure 2(b)) were continuously generated and moved forward. At 0.3 s, the vortex rings caught up with the droplet plume and disturbed the shear-shaped spray plume into a more irregular morphology (Figure 3(d)).

    Figure 3.  Temporal evolution of the spray plume from Ventolin predicted with 2WC: (a) 0.105 s, (b) 0.120 s, (c) 0.200 s, (d) 0.300 s, (e) 0.400 s, (f) 0.600 s (zoomed by 2), (g) 0.900 s, (h) 1.200 s, and (i) 1.500 s. The MDI was actuated at 0.100 s and the applied dose was 3.0 mg.

    In the second stage (0.3–0.8 s), vortex rings formed and moved forward, carrying the spray droplets with them. These vertex rings interacted with each other, as well as with the ambient air, leading to increasingly irregular patterns. Recirculating flows were induced from the interactions between the vortex and ambient air, which slowed down the advancement of the vortices and distributed the entrained droplets in transverse directions. For the aerosol droplets that fell below the spray centerline, the gravitational force became more dominant and started to settle in the negative y-direction (filled hollow arrow).

    The third stage occurred after 0.8 s when no apparent vortices were generated at the mouthpiece. The vortices progressively decayed while all aerosol droplets started settling down under gravity, with larger droplets falling faster, as illustrated by the hollow arrow in Figure 3(h), (i).

    The airflow development under the 2WC effects is further visualized in Figure 4 in terms of the velocity contour lines and streamlines. The jet flow was apparent in Figure 4(a) at the beginning of the actuation, with recirculating flows forming at the mouthpiece. The recirculating flow increased its size from 1.05 to 1.20 s by mixing with the surrounding air (red arrow, Figure 4(b)). It was also noted that jet flow pushed the ambient air at the front (green arrow) and sucked in the surrounding air due to low pressures induced by the high-velocity jet (pink arrow, Figure 4(b)). The recirculating flows above and below the jet centerline were approximately symmetrical, indicating an insignificant effect of the current mass loading (3.0 mg) on the airflow dynamics along the jet centerline (Figure 4(c)(e)). The recirculating flow zones advanced much slowly than the initial discharge velocity (i.e., ~0.4 vs. 40 m/s), indicating the quick decay of the spray plume energy due to the strong mixing with the ambient air. Also, note that the recirculating flow below the jet centerline slightly lagged behind the one above the centerline; this was because there were more large droplets in the lower region, which slowed down the recirculating flow (red arrow, Figure 4(f)). Moreover, the streamlines in the left lower region (filled arrow) bent downward because of the predominate gravitational settling in this region.

    Figure 4.  Airflow development in terms of velocity contour and streamlines under the effects of 2WC: (a) 0.105 s, (b) 0.200 s, (c) 0.400 s, (d) 0.600 s, (e) 0.900 s, and (f) 1.500 s. The MDI was actuated at 0.1 s and the applied dose was 3.0 mg.

    Remarkable differences in spray dynamics were observed with and without the consideration of momentum exchange from aerosols to airflow during the first 0.3 s after actuation, as shown in Figure 5. The MDI was actuated at 0.1 s, and the applied dose was 3.0 mg. One obvious difference was the backflow into the mouthpiece when neglecting the aerosol-to-flow effects (i.e., 1WC from flow to aerosol only). This was reasonable because the spray droplets moved faster than the airflow and thus imparted a forward momentum onto the airflow, which effectively prevented the occurrence of backflows. By contrast, the airflow experienced a large resistance from the upfront pressure, which would push a fraction of airflow backward when neglecting the aerosol-to-flow momentum transfer. These backflows would entrain some droplets into the mouthpiece, as demonstrated in Figure 5(b). The upfront pressure resistance also gave rise to the mushroom-shaped vortex rings at the mouthpiece exit (blue arrow in Figure 5(b)), a phenomenon often reported in jet flows [33,34,35].

    Figure 5.  Comparison of the predicted spray plume development between 2WC and 1WC during the 0.1–0.4 s. The MDI was actuated at 0.1 s with an applied dose of 3.0 mg.

    The second major difference in spray plume evolution between two-way and one-way coupling is the vortex generation and transportation along the centerline. With 2WC, the airflow had a higher momentum and, therefore, exhibited a higher rigidity as a bulk. A more coherent vortex formed (hollow arrow, 0.150 s, Figure 5(a)) than its counterpart with 1WC (hollow arrow, 0.15 s, Figure 5(b)). The vortex with 2WC moved faster, which caught up with or even surpassed the large droplets (red spheres) at 0.3 s (or 0.2 s after actuation), as denoted by the filled arrow in Figure 5(a). By comparison, the vortex under 1WC was still 26 mm short of the large droplets (0.3 s, Figure 5(b)). As alluded to in Figure 3, the vortex would persist for more than one second, which moved forward at a much-reduced speed (~0.4 m/s) than the discharge speed (40 m/s); it transported the polydisperse aerosol droplets along the axial direction, mixed with the surrounding air in transverse directions, and led to complex spray plume morphologies, as illustrated at 0.4 s in Figure 5(a), (b).

    Effects of the mass loading on the airflow and aerosol dynamics are shown in Figure 6 by varying the applied dose mass by three orders of magnitude (0.25–250 mg). Both the streamlines, vortex structures, and aerosol droplets are shown at 0.4 and 0.9 s (i.e., 0.3 and 0.8 s after the actuation), with the vortices being 60% transparent to avoid blocking the aerosol droplets. At a very low mass loading (i.e., 0.25 mg, or one-tenth of the typical MDI dose), the airflow pattern and aerosol distribution resembled those under 1WC. The vortices above and below the axial centerline kept relatively symmetric at both 0.4 and 0.9 s, reflecting the negligible effects from the discrete droplets (Figure 6(a)). By contrast, a high mass loading (i.e., 25 mg, or ten times of typical MDI dose) led to a considerable delay of the vortex below the centerline because of the larger mass that, in turn, resulted from increased gravitational settling at this higher mass loading. Particularly, due to the faster motion of the upper vortex, the core flows bent downward, further escalating the spray plume asymmetry (Figure 6(b)). In Figure 6(c), increasing the applied spray dose to 250 mg resulted in even a higher level of asymmetry, with the majority of airflow streamlines bending downward at 0.9 s. As a result, the mass loading could have a significant impact on the spray plume evolution for a spray dose typical of MDI inhalers or higher; 2WC should be implemented in simulations of MDI drug deliveries.

    Figure 6.  Comparison of the airflow and aerosol dynamics among varying applied doses: (a) 0.30 mg, (b) 30 mg, and (c) 300 mg.

    We further studied the airflow dynamics for monodisperse aerosols with varying droplet diameters, as shown in Figure 7. Significant differences were observed in the airflow pattern among different droplets (Figure 7(a)(d)). For a given dose of 3.0 mg, the airflow with the 2-µm aerosol advanced a much short distance in comparison to the 1WC case (Figure 7(a) vs. 7(e)). The advancement distance increased with the droplet diameter of the monodisperse aerosol from 2 to 20 µm (Figure 7(a)(d)). Compared to 1WC, where the airflow was not affected by the presence of aerosols, the spray plume length was shorter for 5-µm aerosols but longer for 10-µm aerosols. Moreover, the spray plume persistently expanded in transverse directions with the increasing aerosol droplet size, with the 20-µm spray aerosols exhibiting a much wider plume angle than other aerosols (Figure 7(d) vs. 7(a)(c)).

    Figure 7.  Droplet diameter effects on airflow dynamics for monodisperse aerosols under 2WC: (a) 2 µm, (b) 5 µm, (c) 10 µm, (d) 20 µm, vs. (e) 1WC. The MDI was actuated at 0.1 s.

    The corresponding monodisperse aerosol dynamics under 2WC is shown in Figure 8(a)(d) at 0.4 s after actuation for droplet diameter of 2, 5, 10 and 20 µm, all with a dose of 3.0 mg. For comparison purposes, the polydisperse aerosol dynamics under 1WC were also plotted, as shown in Figure 8(e). Very different patterns of aerosol distributions were noted for different droplet sizes, which, however, were generally consistent with the velocity contours in Figure 7. As expected, the largest differences occurred for the largest droplets hereof (i.e., 20 µm, Figure 8(d)).

    Figure 8.  Droplet diameter effects on aerosol dynamics at t = 0.5 s for monodisperse aerosols under 2WC: (a) 2 µm, (b) 5 µm, (c) 10 µm, (d) 20 µm, vs. (e) 1WC (for polydisperse aerosols).

    To understand the strikingly distinct behaviors of the airflow and aerosols for different droplet sizes in Figures 7 and 8, the theory of Kolmogorov scales and ratios were extracted from the spray simulations at 0.5 s under 1WC (Figure 9). Examining 1WC simulations had the advantage of evaluating the tendency of flow/aerosol behaviors without complications caused by the fluid-aerosol interweaving under 2WC. For the spray flow at 0.5 s, the smallest Kolmogorov length scale, η, was 2.72 µm, which led to the maximal magnitude of the aerosol-flow length ratio, dp/η, being 0.74 for 2-µm droplets, 3.68 for 10-µm droplets, and 7.35 for 20-µm droplets (Figure 9(a)). Note that both scales, at their minimal magnitudes, occurred in a very small region immediately downstream of the actuator orifice (first panel, Figure 9(a)). Both scales increased quickly in magnitude at other locations away from the orifice.

    Figure 9.  Kolmogorov scales and ratios for aerosol droplets of varying diameters: (a) Kolmogorov length scale, η and aerosol-flow length scale, dp/η for droplets of 2, 10, and 20 µm, and (b) Kolmogorov time scale, τη and aerosol-flow time scale, tpη for droplets of 2, 10, and 20 µm.

    The distribution of the length ratio dp/η is shown in the second to fourth panels in Figure 9(a) based on a color code range of 0–0.1. For 2-µm droplets, the dp/η magnitude was less than 0.1 for almost all regions except for a very slim region close to the orifice (second panel, Figure 9(a)). Accordingly, the turbulence development was noticeable suppressed, leading to weaker turbulence and vortex intensities, which in turn yielded a shorter penetration distance than when aerosol-to-flow forces were neglected, as previously described in Figure 8(a) vs. 8(e). By contrast, the 10-µm droplets had a much larger region with dp/η > 0.1 (red zone, third panel, Figure 9(a)). In this region, the forces from the aerosols would enhance the turbulent flows and increase the flow strength, as aforementioned in Figure 8(c) vs. 8(e). The dp/η distribution largely resembled the morphology of the orifice jet flow, with the flow-enhancing zones along the axial centerline, which decayed quickly in the transverse directions (third panel, Figure 9(a)). Considering 20-µm droplets, the red zones enlarged in both axial and transverse directions, reflecting a much larger impact from the spray droplets. This observation was consistent with the highly expanded spray cloud for 20-µm droplets, as exhibited in Figures 7(d) and 8(d).

    Figure 9(b) shows the Kolmogorov time scale and the droplet-specific time-scale-ratio, τp/τη for 2, 10, and 20-µm aerosols. For the spray flow at 0.5 s, the smallest Kolmogorov time scale, τη, was 5.1e-7 s, which was smaller than the residence time, τp, for all droplet sizes considered (2–20 µm), as listed in Table 2. It is also noted that the smallest Kolmogorov time scale (τη = 5.1e-7) existed in a very small region only; there was an appreciable region along the spray centerline with τη ranging from 1.0e-5 to 1.0e-3, which was equivalent to the residence time τp for droplets ranging from 2 to 20 µm. As a result, the airflow dynamics would be noticeably modified by the entrained droplets downstream of the actuator orifice, which was also demonstrated to modify the droplet behaviors and fates, as discussed in detail in Figures 7 and 8.

    Table 2.  Droplet Residence time, τp and droplet Stokes number, Stk for different droplet sizes and with varying characteristic velocities and lengths (i.e., at different jet flow locations).
    Droplet size (µm) Residence time, τp (s) Stokes number, Stk
    Up: 40 m/s, L: 0.3 mm Up: 10 m/s, L = 0.3 mm Up: 10 m/s, L = 3 cm
    2 1.24E-05 1.66 0.41 0.004
    3.5 3.80E-05 5.07 1.27 0.01
    5 7.76E-05 10.35 2.59 0.03
    7.5 1.75E-04 23.28 5.82 0.06
    10 3.10E-04 41.39 10.35 0.10
    15 6.98E-04 93.12 23.28 0.23
    20 1.24E-03 165.55 41.39 0.41

     | Show Table
    DownLoad: CSV

    The Kolmogorov length scale, η, and ratio, dp/η, were further quantified in the axial and transverse directions downstream of the actuator orifice at 0.5 s (Figure 10(a)). As alluded to in Figure 9, the η magnitude increased with the distance from the orifice because of the decaying eddies or vortices away from the orifice (Figure 10(a)). Figure 10(b) shows the length scale ratio dp/η vs. the axial distance (i.e., Y-direction). As expected, in most regions, the dp/η magnitude for 2-µm droplets was lower than 0.1 (pink dashed line). By contrast, for 10-µm and 20-µm aerosols, the dp/η magnitude was larger than 0.1 for Y < 12 cm. A bump was observed at Y = 1–3 cm for 10-µm and 20-µm aerosols (Figure 10(b)) because the spray actuation had been completed at 0.5 s and the core jet flow started to decay thenceforth.

    Figure 10.  Kolmogorov length scale (η) and aerosol-flow ratio (dp/η): (a) η variation in the axial (Y) direction, (b) dp/η variation in the axial direction, (c) η variation in the transverse directions (Z, X) at Y = 3 and 6 cm, and (b) η variation in the transverse directions (Z, X) at Y = 9 and 12 cm.

    The distributions of the Kolmogorov length scale in the circumneutral directions (Z- and X- directions) are shown in Figure 10(c), (d) at different distances from the orifice. For all cross-sections considered, the Kolmogorov length scale η magnitude (the eddy scale) increased quickly from the axial centerline and became larger than 200 µm (10 times of the 20-µm droplets) when the radius from the centerline was larger than 0.5 cm. For a given radius, the η magnitude also decreased with the axial distance. At the cross-sections of Y = 3 and 6 cm, the η distributions were approximately symmetric upper-lower (Z) or left-right (X). This symmetry was also observed in the X-direction (left-right) at Y = 9 and 12 cm (two solid lines, Figure 10(d)). However, apparent asymmetry was observed in the gravitational Z-direction, with a much higher value above the axial centerline than below.

    The particle Reynolds number has been suggested as an index to differentiate the turbulence suppression and enhancement. In Figure 11, the flow and aerosol velocities were examined for the 2-µm aerosols under 2WC at t = 0.5 s. Figure 11(a) shows the airflow velocities at the locations that are coincident with the instantaneous particles. A wide range of velocities was observed, with low-speed flows in most of the regions downstream of the mouthpiece in contrast to a high-speed jet flow along the centerline. To better understand the velocity distributions, a slice view was presented in the lower panel of Figure 11(a), which only showed the aerosols within a 2-mm-thickness slice in the midplane. Without obscuring from surrounding droplets, the jet flow and vortex (pink square) were clearly displayed. The vortex was further visualized in a zoomed view, where flow recirculation was apparent.

    Figure 11.  Comparison of the airflow and aerosol droplet velocities for 2-µm aerosols with 2WC at t = 0.5 s: (a) airflow, (b) droplet velocities and relative velocities.

    The instantaneous velocities of the spray droplets within the midplane slice (2 mm thickness) are presented in Figure 11(b), which were much higher than the airflow velocities. Unlike the various directions of the airflow, most droplets exhibited a forward direction. The differences in the velocity and direction between airflow and aerosol droplets are shown in the right lower panel of Figure 11(a), where both the instantaneous velocity vectors of the flow and droplets are presented. Note that the flow and droplet velocity vectors were not to the scale in order to plot them in the same frame. It was observed that both the magnitude and direction could differ significantly between the airflow and an individual droplet, and that these differences varied among droplets. This made it challenging to use the particle Reynolds number as a practical criterion for turbulence-depression-enhancement.

    The effects of the monodisperse aerosol diameter on the spray plume evolution were evaluated by comparing the aerosol velocities at two points (i.e., 3 and 6 cm downstream of the mouthpiece). Both the maximal and average velocities were quantified and presented in Figure 12(a) as a function of the droplet diameter (2, 3.5, 5, 7.5, 10, 15, 20 µm). The aerosol velocity peaked around 5–7 µm for both the maximal and average values at both sampling points considered. The velocity decreased thereafter and then increased again from 15 µm (Figure 12(a)). The velocity decay from 3 to 6 cm had a similar magnitude for both the maximal and average velocities (Figure 12(a)).

    Figure 12.  Comparison of the airflow and aerosol velocities predicted using 2WC: (a) the maximal and average aerosol velocities vs. droplet size at two sampling points (3 and 6 cm downstream of the mouthpiece), and (b) airflow velocities vs. the distance from the mouthpiece for four monodisperse aerosols (i.e., 2, 5, 10, 20 µm).

    Figure 12(b) shows the velocities of the airflow along the axial direction (1–8 cm from the mouthpiece) for different aerosol diameters. Overall, the aerosol speeds are much higher than those of the carrier flow (Figure 12(a) vs. 12(b)). The airflow slowed down quickly after exiting the mouthpiece. However, the rate of decrease differed among aerosols, with 2-µm aerosols decaying most drastically and 5-µm aerosols decaying the slowest, which was consistent with the trend observed in Figure 12(a).

    The effects of the 2WC on different MDI types were evaluated by considering the other two inhalers: ProAir with a mean droplet size of 5.2 μm emitted at 26 m/s, and Qvar inhaler with a mean droplet size of 2.9 μm emitted at 23 m/s [17]. As expected, the predicted airflow dynamics without the aerosol impact were very similar regardless of the inhaler types (Figure 13(a) vs. 13(b)). The differences in the one-way-predicted spray morphologies at 0.5 s between ProAir and Qvar were noticeable but not significant. These slight differences were presumably attributed to the differences in the droplet sizes (5.2 vs. 2.9 µm) and initial velocities (26 vs. 23 m/s). Remarkable differences were predicted with 2WC in both the airflow and aerosol dynamics between Proair and Qvar (Figure 13(a) vs. 13(b)). Even though the airflow was suppressed by the entrained spray aerosols from both inhalers, that from Qvar was significantly more suppressed due to its smaller droplet sizes and lower dp/η ratio. As a result, the Qvar spray plume was much less dispersed and advanced a shorter distance (Figure 13(b)).

    Figure 13.  Comparison of the airflow and aerosol dynamics at t = 0.5 s between two-way and 1WC: (a) ProAir inhaler with a mean droplet size of 5.2 µm and an initial velocity of 26 m/s, and (b) Qvar inhaler with a mean droplet size of 2.9 µm and an initial velocity of 23 m/s.

    The airflow and aerosol dynamics under 2WC after 0.5 s are presented in Figure 14(a), (b) for ProAir and Qvar, respectively. For ProAir (with a mean aerosol diameter of 5.2 µm), a smaller vortex formed below the axial centerline, as denoted by the hollow red arrow in Figure 14(a). The airflow posterior to the vortex bent downward and became increasingly more apparent with time (hollow blue arrow in Figure 14(a)), indicating an increasing gravitational effect from large spray droplets. The bending-downward flow and the small vortex were not observed in Qvar. This was reasonable considering the much smaller aerosol mean diameter (2.9 µm) than that of ProAir (5.2 µm). Moreover, the spray plumes from Qvar were less dispersed than those from ProAir and advanced a shorter distance than their counterparts in ProAir. The dipole vortices below and above the axial centerline were slightly off asymmetric at t = 1.0 s, but overall, the bending-downward motions of airflow or aerosols were not observed (third panel, Figure 14(b)). The differences shown in this figure, together with Figure 7, indicated that the 2WC effects are sensitive to the spray droplet size, and thus it will be essential in numerical studies of MDI pulmonary drug delivery to consider the 2WC effects for different types of inhalers, which often had very different spray droplets and exit velocities.

    Figure 14.  Comparison of the airflow and aerosol dynamics under 2WC between: (a) ProAir and (b) Qvar at t = 0.6, 0.8, and 1.0 s.

    The magnitude of the 2WC force (Eq (1)) was proportional to both the slip velocity and local mass loading. It was demonstrated in this study (Figures 4 and 5) that a moderate mass loading typical of an MDI spray dose per puff (e.g., 3.0 mg) could alter the fluid flow and aerosol behaviors. For an applied dose ten times smaller (i.e., 0.30 mg), no significant differences were predicted between two-way and one-way coupling (Figure 7). However, for an applied dose ten times larger (i.e., 30 mg), both the flow and aerosols were remarkably modified under two-way-coupling. The jet flow became apparently asymmetric above and below the spray centerline while the front of the spray plume started to bend downward. In other words, the gravitational effect became more pronounced at higher mass loadings, which drove the droplet-laden flow toward the ground (Figure 7). These observations were in agreement with Tong and Wang, who reported that moderate loadings of particles could induce noticeable flow modulation [36]. Similarly, Monchaux and Dejoan reported that 2WC strongly enhanced particle settling than the 1WC case and that these particles dragged the fluid downward [14]. These also conformed to the three regimes of particle-laden flows proposed by Elghobashi et al. [15]. At very low volume fractions (φ < 10-6), the discrete phase did not affect the carrier flow and was regarded as the 1WC regime. In the second regime, the discrete phase had a low volume fraction but a moderate mass loading due to the high droplet-flow density ratio. In this case, droplets could perceivably modulate the co-flow (i.e., 2WC regime), but droplet-droplet interactions could be neglected. Note that the MDIs fell within this regime. In the third regime (i.e., dense flow regime), both the volume fraction and mass concentration were high, and collisions among droplets had to be included.

    To evaluate the droplet size effect on 2WC, monodisperse aerosols of varying sizes (2–20 μm) were simulated and compared for a given applied mass of 3.0 mg. Strikingly different spray plume morphologies were observed for different droplet sizes, with 2-μm aerosols significantly suppressed and 20-μm aerosols dispersed in comparison to their counterparts with 1WC simulations (Figures 8 and 9). To find out the underlying reasons for the sticky behavior of small particles and dispersing behavior of large parties, Kolmogorov scales for the airflow and scale ratios between droplets and airflow were quantified at t = 0.5 s. Turbulence modification in particle-laden flows has been acknowledged for several decades [7,12,37]. Early experimental observations by Torobin and Gauvin showed that the presence of particles changed the wall drag in pipes, as well as the heat transfer and chemical reaction rates, which could not be explained unless the fluid turbulence had been modified by particles [7]. In a comprehensive review, Gore and Crowe listed a spectrum of potential effects of the carrier flow depending on the particle size and volumetric loading. The ratio of the particle size to the turbulence integral length was proposed as an indicator of whether the turbulence was increased or decreased [10]. In this study, the ratio of the droplet diameter over the Kolmogorov length scale (dp/η) was quantified (Figures 10 and 11), which was found to be closely associated with the spray plume morphologies (i.e., advancement in the axial direction and dispersion in transverse directions) (Figures 8 and 9). For instance, the seemingly sticky aerosols of 2 μm had a dp/η < 0.1 in near all regions of the spray, while the highly dispersed and advanced aerosols of 20 μm had a dp/η < 0.1 only in the regions further than 15 cm (Figure 11(b)). This finding was consistent with Elghobashi and Truesdell, who reported that turbulence was suppressed when dp/η < 0.1 and enhanced when dp/η > 0.1 [12].

    It was also observed that the transverse dispersion of 20-μm aerosols in Figure 9 was much larger than the regions with dp/η > 0.1 (red zone) in Figure 10, indicating that the Kolmogorov scale ratio theory might not fully explain the spray suppression/expansion in the transverse directions. One salient feature of MDI sprays was mixing layers at the interface between the jet and ambient air, where intensive vortex formations occurred. Many experimental and numerical investigations have demonstrated that aerosol dynamics in the mixing layers differed significantly from that of the passive tracers if the particle/droplet response time was comparable to the characteristic time scale of the flow. Aerosols would accumulate in the outskirts of vortices, especially in the proximity of braid stagnation points where two counter-rotating vortices met [38]. Therefore, the strong mixing and vortex formations at the jet interfaces could have assisted the transverse dispersion in addition to the Kolmogorov-associated turbulence enhancement from the core flow of 20-μm aerosols.

    Turbulence modification by particles was also correlated to the particle Reynolds number Re,p based on the interphase slip velocity, with a low Re,p decreasing the fluid turbulence, and a high Re,p increasing the turbulence [39,40,41]. In this study, the Re,p was also sought by characterizing the slip velocity between an individual droplet and the local flow (Figure 12). However, both the airflow and droplets in the MDI sprays were highly transient in time and heterogeneous in space (Figure 12), making it impractical to quantify the Re,p for individual droplets based on their instantaneous slip velocities relative to the local flow. By contrast, the Kolmogorov length scale of the spray flow was more practical to quantify, and its ratio to droplet diameters has been demonstrated to be useful in predicting suppression or enhancement of the carrier flow and aerosol dispersion (Figures 8 and 9 vs. Figures 10 and 11).

    The high sensitivity of the MDI spray plume to its droplet sizes rested on several MDI-specific factors, such as the high-speed jet flow from the orifice that led to Kolmogorov length scales comparable to the droplet sizes (2–20 μm), and the applied mass loading (3.0 mg in this study) that was large enough to elicit perceivable aerosol-to-flow momentum exchanges. In this study, we demonstrated that reducing the mass loading by ten times (0.30 mg) led to negligible two-coupling effects while increasing the mass loading by ten times (30 mg) led to large variations in spray plume morphologies (Figure 7). It was reminded that the droplet-size-dependent spray behaviors were shown for monodisperse aerosols at a 3.0-mg mass loading (Figures 8 and 9). For polydisperse aerosols like an MDI spray, the transient morphology of the aerosol plume was an instantaneous collective manifestation of all droplets with varying diameters and probability densities. For Ventolin, 2WC provided a better approximation to the experimental aerosol behaviors close to the orifice (Figure 5 vs. Figure 2); it also predicted a more detached aerosol bolus from the mouthpiece and a slightly faster advancement than 1WC (Figure 5). By contrast, for ProAir and Qvar, 2WC predicted a spray plume advancing at lower speeds than 1WC (Figures 14 and 15). However, the Qvar spray development appeared more retarded in both axial advancement and transverse dispersion due to its smaller mean diameter (2.9 μm) and smaller length ratio, dp/η, than those of ProAir (5.2 μm).

    Even though only three types of inhalers were considered hereof, the spray behaviors of other MDIs can be roughly estimated based on the droplet-Kolmogorov length scale, which can be readily calculated from the aerosol size and discharge flow speed from the actuator orifice. Another implication was that previous one-way-coupled predictions of the delivery efficiencies from MDI delivery might have errors of varying degrees, especially for Qvar with small aerosol sizes (2.9 μm) and suppressed plume developments. Future studies are warranted to include the 2WC effects in predictive inhalation dosimetry in the respiratory tract for different inhalers.

    MDI actuation and the subsequent spray plume development are highly dynamic and involve multiscale and multi-physics. Multiple assumptions were made in this study to fill the gap of missing data in current literature or mitigate the prohibitive numerical requirements. These included the same discharge velocity for aerosol droplets and co-flow, a constant aerosol output during the 0.2-s actuation, a constant temperature, no evaporation, no electrostatic charge, and no droplet interactions. The MDI formulation contains approximately 80–90% w/w propellant, which vaporizes quickly to atomize the solution into micro-scale droplets; however, measurements of propellent flow speed at the actuator orifice are not yet available [24]. Liu et al. measured the spray (aerosol and co-flow) velocities at 3 and 6 cm from the mouthpiece of different inhalers, including Ventolin, ProAir, and Qvar [17]. These measurements were used to reversely estimate the discharge velocity at the orifice, as explained in Talaat et al. [18,30]. The spray output during actuation can also be time-dependent during the 0.2-s actuation, adding more flow fluctuations. Drastic temperature variation can occur during MDI actuation, leading to thermophoretic effects on droplets [42,43,44]. The ethanol droplets can evaporate and decrease in size, further altering the thermal environment and droplet behavior [45,46,47]. Particle-particle interactions can also be important at high mass loading, e.g., 300 mg, but should be insignificant due to the low volume fraction ratio of the liquid droplets in the open space. Electrostatic charges on the droplets can be another factor that complicates the trajectories of individual droplets and interactions among neighboring droplets [48,49,50,51,52]. An orifice diameter of 0.3 mm was used all three inhalers, which was larger than Qvar's actual diameter (0.25 mm) and smaller than Ventolin's diameter (0.48 mm). Its effects on aerosol dynamics were expected to be insignificant considering that the experimentally measured aerosol sizes and speeds were implemented for Qvar. For Ventolin, we acknowledged that our LES simulations could not be regarded as replicating the experimental observations, given the differences between the model and prototype orifice diameters. Rather, the general features of 2WC simulations seemed to be more representative of the Ventolin plume than 1WC simulations. Furthermore, this study considered only spray evolutions in the open space, which can be different from those in a closed space like the respiratory tract [53,54,55,56]. However, neglecting these compounding factors allowed us to isolate the 2WC effect and examine it vs. the 1WC counterpart systemically under controlled, parametrized conditions, i.e., with different mass loadings, droplet sizes, and inhaler types.

    In summary, spray plume evolutions in open space were simulated using the LES-Lagrangian approach under both 2WC and 1WC conditions. Differences in the temporal and spatial morphologies of the spray plumes were characterized, and the reasons underlying the differences were explored using Kolmogorov scale ratios relative to the spray droplets. Specific findings are listed below:

    1) Numerical simulations with 2WC seemed to provide better agreement with complimentary high-speed imaging of the spray from the orifice than 1WC.

    2) Increasing the mass loading altered the spray morphology, with both the airflow and aerosol droplets bending downward at the spray plume front.

    3) The droplet size had a significant impact on the spray evolution in both axial (advancement) and transverse (dispersion) directions, with suppressed sprays for 2-μm aerosols and more dispersed sprays for 20-μm aerosols.

    4) The Kolmogorov length scale, η of the spray jet flow were comparable to the MDI droplet sizes. The length scale ratio dp/η correlated well with spray suppression when less than 0.1 and with spray enhancement when larger than 0.1.

    5) Different levels of impact from 2WC were predicted on the spray plume evolution of Ventolin, ProAir, and Qvar in comparison to 1WC.

    Amr Seifelnasr was gratefully acknowledged for reviewing this manuscript.

    The authors declare there is no conflict of interest.


    Acknowledgments



    Thanks to the Pneumomeditazione© teachers and the Entolé staff for their support. Useful discussions with Professor Sergio Pagano are also acknowledged.

    Conflict of interest



    The authors declare no conflict of interest.

    [1] Kocica MJ, Corno AF, Carreras-Costa F, et al. (2006) The helical ventricular myocardial band: global, three-dimensional, functional architecture of the ventricular myocardium. Eur J Cardiothorac Surg 29: S21-S40. doi: 10.1016/j.ejcts.2006.03.011
    [2] Dal Lin C, Tona F, Osto E (2018) The heart as a psychoneuroendocrine and immunoregulatory organ. Sex-Specific Analysis of Cardiovascular Function 1065: 225-239. doi: 10.1007/978-3-319-77932-4_15
    [3] Marinelli RA, Penney DG, Marinelli W, et al. (1991) Rotary motion in the heart and blood vessels: a review. J Appl Cardiol 6: 421-431.
    [4] Sengupta PP, Narula J, Chandrashekhar Y (2014) The dynamic vortex of a beating heart: Wring out the old and ring in the new!. J Am Coll Cardiol 64: 1722-1724. doi: 10.1016/j.jacc.2014.07.975
    [5] Bottio T, Buratto E, Dal Lin C, et al. (2012) Aortic valve hydrodynamics: considerations on the absence of sinuses of Valsalva. J Heart Valve Dis 21: 718-723.
    [6]  Calgary University Medicine: Blood Moving Through the Heart - 4D Flow Available from: https://www.youtube.com/watch?v=sMeaD3Jh64E.
    [7] Dal Lin C, Marinova M, Rubino G, et al. (2018) Thoughts modulate the expression of inflammatory genes and may improve the coronary blood flow in patients after a myocardial infarction. J Tradit Complement Med 8: 150-163. doi: 10.1016/j.jtcme.2017.04.011
    [8] Armour JA (2007) The little brain on the heart. Cleve Clin J Med 74: S48-S51. doi: 10.3949/ccjm.74.Suppl_1.S48
    [9] Lane RD, Reiman EM, Ahern GL, et al. (1982) Activity in medial prefrontal cortex correlates with vagal component of heart rate variability during emotion. Brain Cognition 47: 97-100.
    [10] Jennings JR, Sheu LK, Kuan DCH, et al. (2016) Resting state connectivity of the medial prefrontal cortex covaries with individual differences in high-frequency heart rate variability. Psychophysiology 53: 444-454. doi: 10.1111/psyp.12586
    [11] Schandry R, Montoya P (1996) Event-related brain potentials and the processing of cardiac activity. Biol Psychol 42: 75-85. doi: 10.1016/0301-0511(95)05147-3
    [12] Garfinkel SN, Barrett AB, Minati L, et al. (2013) What the heart forgets: Cardiac timing influences memory for words and is modulated by metacognition and interoceptive sensitivity. Psychophysiology 50: 505-512. doi: 10.1111/psyp.12039
    [13] Azevedo RT, Garfinkel SN, Critchley HD, et al. (2017) Cardiac afferent activity modulates the expression of racial stereotypes. Nat Commun 8: 13854. doi: 10.1038/ncomms13854
    [14] Garfinkel SN, Minati L, Gray MA, et al. (2014) Fear from the heart: Sensitivity to fear stimuli depends on individual heartbeats. J Neurosci 34: 6573-6582. doi: 10.1523/JNEUROSCI.3507-13.2014
    [15] Montoya P, Schandry R, Müller A (1993) Heartbeat evoked potentials (HEP): topography and influence of cardiac awareness and focus of attention. Electroencephalogr Clin Neurophysiol Evoked Potentials 88: 163-172. doi: 10.1016/0168-5597(93)90001-6
    [16] Thayer JF, Lane RD (2000) A model of neurovisceral integration in emotion regulation and dysregulation. J Affect Disord 61: 201-216. doi: 10.1016/S0165-0327(00)00338-4
    [17] Park G, Thayer JF (2014) From the heart to the mind: cardiac vagal tone modulates top-down and bottom-up visual perception and attention to emotional stimuli. Front Psychol 5: 278. doi: 10.3389/fpsyg.2014.00278
    [18] Thayer JF, Hansen AL, Saus-Rose E, et al. (2009) Heart rate variability, prefrontal neural function, and cognitive performance: the neurovisceral integration perspective on self-regulation, adaptation, and health. Ann Behav Med 37: 141-153. doi: 10.1007/s12160-009-9101-z
    [19] Thayer JF, Sternberg E (2006) Beyond heart rate variability: vagal regulation of allostatic systems. Ann N Y Acad Sci 1088: 361-372. doi: 10.1196/annals.1366.014
    [20] Lin PF, Lo MT, Tsao J, et al. (2014) Correlations between the signal complexity of cerebral and cardiac electrical activity: a multiscale entropy analysis. PLoS One 9: e87798. doi: 10.1371/journal.pone.0087798
    [21] Aftanas LI, Brak IV, Reva NV, et al. (2013) Brain oscillations and individual variability of cardiac defense in human. Ross Fiziol Zh Im IM Sechenova 99: 1342-1356.
    [22] McCraty R, Atkinson M, Bradley RT (2004) Electrophysiological evidence of intuition: Part 2. A system-wide process? J Altern Complement Med 10: 325-336. doi: 10.1089/107555304323062310
    [23] Gray MA, Beacher FD, Minati L, et al. (2012) Emotional appraisal is influenced by cardiac afferent information. Emotion 12: 180-191. doi: 10.1037/a0025083
    [24] Craig ADB (2009) How do you feel—now? The anterior insula and human awareness. Nat Rev Neurosci 10: 59-70. doi: 10.1038/nrn2555
    [25] Craig ADB (2014)  How do you feel? An interoceptive moment with your neurobiological self Princeton: Princeton University Press. doi: 10.23943/princeton/9780691156767.001.0001
    [26] Grossmann I, Sahdra BK, Ciarrochi J, et al. (2016) Heart and a mind: Self-distancing facilitates the association between heart rate variability, and wise reasoning. Front Behav Neurosci 10: 68. doi: 10.3389/fnbeh.2016.00068
    [27] Rahman SU, Hassan M (2013) Heart's role in the human body: A literature review. ICCSS 2: 1-6.
    [28] McCraty R, Trevor Bradley R, Tomasino D (2004) The resonant heart. Front Counsciousness 5: 15-19.
    [29] McCraty R, Atkinson M, Tomasino D, et al. (2009) The coherent heart: Heart-brain interactions, psychophysiological coherence, and the emergence of system-wide order. Integr Rev 5: 10-115.
    [30] Goldstein DS (2012) Neurocardiology: therapeutic implications for cardiovascular disease. Cardiovasc Ther 30: 89-106. doi: 10.1111/j.1755-5922.2010.00244.x
    [31] Dal Lin C, Poretto A, Scodro M, et al. (2015) Coronary microvascular and endothelial function regulation: Crossroads of psychoneuroendocrine immunitary signals and quantum physics. J Integr Cardiol 1: 132-163.
    [32] Dal Lin C, Tona F, Osto E (2019) The crosstalk between the cardiovascular and the immune system. Vasc Biol 1: H83-H88. doi: 10.1530/VB-19-0023
    [33] Dal Lin C, Tona F, Osto E (2015) Coronary microvascular function and beyond: The crosstalk between hormones, cytokines, and neurotransmitters. Int J Endocrinol 2015: 1-17. doi: 10.1155/2015/312848
    [34] Lashley KS (1942) The problem of cerebral organization in vision. Visual Mechanisms. Biological Symposia Lancaster: Jaques Cattell Press, 301-322.
    [35] Pribram KH (1991)  Brain and Perception Hillsdale: Lawrence Erlbaum.
    [36] Freeman WJ (1975)  Mass Action in the Nervous System New York: Academic Press.
    [37] Freeman WJ (2000)  Neurodynamics: An Exploration of Mesoscopic Brain Dynamics Berlin: Springer. doi: 10.1007/978-1-4471-0371-4
    [38] Ricciardi LM, Umezawa H (1967) Brain and physics of many-body problems. Kybernetik 4: 44-48. doi: 10.1007/BF00292170
    [39] Goldstone J, Salam A, Weinberg S (1962) Broken Symmetries. Phys Rev 127: 965-970. doi: 10.1103/PhysRev.127.965
    [40] Umezawa H (1995)  Advanced field theory: Micro, macro, and thermal physics New York: AIP.
    [41] Blasone M, Jizba P, Vitiello G (2011)  Quantum Field Theory and its macroscopic manifestations: Boson condensation, ordered patterns, and topological defects London: Imperial College Press. doi: 10.1142/p592
    [42] Jibu M, Yasue K (1995)  Quantum brain dynamics and consciousness Amsterdam: John Benjamins Publ. doi: 10.1075/aicr.3
    [43] Umezawa H (1995) Development in concepts in quantum field theory in half century. Math Jpn 41: 109-124.
    [44] Vitiello G (1995) Dissipation and memory capacity in the quantum brain model. Int J Mod Phys B 9: 973-989. doi: 10.1142/S0217979295000380
    [45] Vitiello G (2001)  My double unveiled Amsterdam: John Benjamins Publ. doi: 10.1075/aicr.32
    [46] Freeman WJ, Vitiello G (2006) Nonlinear brain dynamics as macroscopic manifestation of underlying many-body dynamics. Phys Life Rev 3: 93-118. doi: 10.1016/j.plrev.2006.02.001
    [47] Freeman WJ, Livi R, Obinata M (2012) Cortical phase transitions, non-equilibrium thermodynamics and the time dependent Ginzburg-Landau equation. Int J Mod Phys B 26: 1250035. doi: 10.1142/S021797921250035X
    [48] Alfinito E, Vitiello G (2000) Formation and lifetime of memory domains in the dissipative quantum model of brain. Int J Mod Phys B 14: 853-868.
    [49] Vitiello G (2012) Fractals, coherent states and self-similarity induced noncommutative geometry. Phys Lett A 376: 2527-2532. doi: 10.1016/j.physleta.2012.06.035
    [50] Freeman W, Vitiello G (2016) Matter and mind are entangled in two streams of images guiding behavior and informing the subject through awareness. Mind Matter 14: 7-24.
    [51] Del Giudice E, Doglia S, Milani M, et al. (1985) A quantum field theoretical approach to the collective behavior of biological systems. Nucl Phys B 251: 375-400. doi: 10.1016/0550-3213(85)90267-6
    [52] Del Giudice E, Doglia S, Milani M, et al. (1986) Electromagnetic field and spontaneous symmetry breakdown in biological matter. Nucl Phys B 275: 185-199. doi: 10.1016/0550-3213(86)90595-X
    [53] Del Giudice E, Preparata G, Vitiello G (1988) Water as a free electric dipole laser. Phys Rev Lett 61: 1085-1088. doi: 10.1103/PhysRevLett.61.1085
    [54] Engel GS, Calhoun TR, Read EL, et al. (2007) Evidence for wavelike energy transfer through quantum coherence in photosynthetic systems. Nature 446: 782-786. doi: 10.1038/nature05678
    [55] Peacock JA (1990) An in vitro study of the onset of turbulence in the sinus of valsalva. Circ Res 67: 448-460. doi: 10.1161/01.RES.67.2.448
    [56] Mettauer B, Levy F, Richard R, et al. (2005) Exercising with a denervated heart after cardiac transplantation. Ann Transplant 10: 35-42.
    [57] Armour JA, Ardell JL (2004)  Basic and Clinical Neurocardiology Oxford: Oxford University Press.
    [58] Biasetti J, Hussain F, Gasser TC (2011) Blood flow and coherent vortices in the normal and aneurysmatic aortas: a fluid dynamical approach to intra-luminal thrombus formation. J R Soc Interface 8: 1449-1461. doi: 10.1098/rsif.2011.0041
    [59] Matsumoto H, Papastamatiou NJ, Umezawa H, et al. (1975) Dynamical rearrangement in the Anderson-Higgs-Kibble mechanism. Nucl Phys B 97: 61-89. doi: 10.1016/0550-3213(75)90215-1
    [60] Matsumoto H, Papastamatiou NJ, Umezawa H (1975) The boson transformation and the vortex solution. Nucl Phys B 97: 90-124. doi: 10.1016/0550-3213(75)90216-3
    [61] Manka R, Vitiello G (1990) Topological solitons and temperature effects in gauge field theory. Ann Phys 199: 61-83. doi: 10.1016/0003-4916(90)90368-X
    [62] Vitiello G (2000) Defect formation through boson condensation in quantum field theory. Topological Defects and the Non-Equilibrium Dynamics of Symmetry Breaking Phase Transitions Dordrecht: Springer, 171-191. doi: 10.1007/978-94-011-4106-2_9
    [63] Meyer G, Vitiello G (2018) On the molecular dynamics in the hurricane interactions with its environment. Phys Lett A 382: 1441-1448. doi: 10.1016/j.physleta.2018.03.044
    [64] Da Silva AF, Carpenter T, How TV, et al. (1997) Stable vortices within vein cuffs inhibit anastomotic myointimal hyperplasia? Eur J Vasc Endovasc Surg 14: 157-163. doi: 10.1016/S1078-5884(97)80185-2
    [65] Kefayati S, Amans M, Faraji F, et al. (2017) The manifestation of vortical and secondary flow in the cerebral venous outflow tract: An in vivo MR velocimetry study. J Biomech 50: 80-187. doi: 10.1016/j.jbiomech.2016.11.041
    [66] Lurie F, Kistner RL, Eklof B, et al. (2003) Mechanism of venous valve closure and role of the valve in circulation: a new concept. J Vasc Surg 38: 955-961. doi: 10.1016/S0741-5214(03)00711-0
    [67] Boisseau MR (1997) Venous valves in the legs: hemodynamic and biological problems and relationship to physiopathology. J Mal Vasc 22: 122-127.
    [68] Machi J, Sigel B, Ramos JR, et al. (1986) Sonographic evaluation of platelet aggregate retention in a vortex within a simulated venous sinus. J Ultrasound Med 5: 685-689. doi: 10.7863/jum.1986.5.12.685
    [69] Meyer G, Vitiello G (2019) On the hurricane collective molecular dynamics. J Phys Conf Ser 1275: 012017. doi: 10.1088/1742-6596/1275/1/012017
    [70] Higgs PW (1966) Spontaneous symmetry breakdown without massless bosons. Phys Rev 145: 1156-1163. doi: 10.1103/PhysRev.145.1156
    [71] Freeman W, Vitiello G (2010) Vortices in brain waves. Int J Mod Phys B 24: 3269-3295. doi: 10.1142/S0217979210056025
    [72] Del Giudice E, Vitiello G (2006) The role of the electromagnetic field in the formation of domains in the process of symmetry breaking phase transitions. Phys Rev A 74: 02210.
    [73] De Lauro E, De Martino S (2019) On the heart vibrations: Some insights from ECG and laser doppler vibrometry.p. atticon11705 V-C-5in press.
    [74] Scalise L (2012) Non contact heart monitoring. Advances in Electrocardiograms-Methods and Analysis Rijeka, Croatia: InTech, 81-106.
    [75] Tomasini E, Pinotti M, Paone N (1998) Carotid artery pulse wave measured by a laser vibrometer.3411: 611-616.
    [76] Morbiducci U, Scalise L, De Melis M (2006) Optical vibrocardiography: a novel tool for optical monitoring of cardiac activity. Ann Biomed Eng 35: 45-58. doi: 10.1007/s10439-006-9202-9
    [77] Buccheri G, De Lauro E, De Martino S, et al. (2016) Experimental study of self-oscillations of the trachea–larynx tract by laser doppler vibrometry. Biomed Phys Eng Express 2: 055009. doi: 10.1088/2057-1976/2/5/055009
    [78] Tavakolian K, Dumont GA, Blaber AP (2012) Analysis of seismocardiogram capability for trending stroke volume changes: A lower body negative pressure study. Computing in Cardiology 2012: 733-736.
    [79] Desjardins CL, Antonelli LT (2007) A remote and non-contact method for obtaining the blood pulse waveform with a laser Doppler vibrometer. Advanced Biomedical and Clinical Diagnostic Systems V, Proc. SPIE .
    [80] Hyvärinen A, Karhunen J, Oja E (2001)  Independent Component Analysis Hoboken, USA: John Wiley & Sons. doi: 10.1002/0471221317
    [81] Capuano P, De Lauro E, De Martino S, et al. (2017) Convolutive independent component analysis for processing massive datasets: a case study at Campi Flegrei (Italy). Nat Hazards 86: 417-429. doi: 10.1007/s11069-016-2545-0
    [82] Capuano P, De Lauro E, De Martino S (2016) Detailed investigation of long-period activity at Campi Flegrei by convolutive independent component analysis. Phys Earth Planet Inter 253: 48-57. doi: 10.1016/j.pepi.2016.02.003
    [83] Capuano P, De Lauro E, De Martino S (2011) Water-level oscillations in the Adriatic Sea as coherent self-oscillations inferred by independent component analysis. Prog Oceanogr 91: 447-460. doi: 10.1016/j.pocean.2011.06.001
    [84] De Lauro E, De Martino S, Falanga M, et al. (2006) Statistical analysis of Stromboli VLP tremor in the band [0.1–0.5] Hz: some consequences for vibrating structures. Nonlinear Process Geophys 13: 393-400. doi: 10.5194/npg-13-393-2006
    [85] Rüssel IK, Götte MJW, Bronzwaer JG, et al. (2009) Left ventricular torsion: an expanding role in the analysis of myocardial dysfunction. JACC Cardiovasc Imaging 2: 648-655. doi: 10.1016/j.jcmg.2009.03.001
    [86] Loppini A, Capolupo A, Cherubini C (2012) On the coherent behavior of pancreatic beta cell clusters. Phys Lett A 378: 3210-3217. doi: 10.1016/j.physleta.2014.09.041
    [87] Dal Lin C, Brugnolo L, Marinova M, et al. (2020) Toward a unified view of cognitive and biochemical activity: Meditation and linguistic self-reconstructing may lead to inflammation and oxidative stress improvement. Entropy 22: 818. doi: 10.3390/e22080818
    [88] Dal Lin C, Radu CM, Vitiello G, et al. (2020) In vitro effects on cellular shaping, contratility, cytoskeletal organization and mitochondrial activity in HL1 cells after different sounds stimulation. A qualitative pilot study and a theoretical physical model, 2020.
    [89] McFadden J, Al-Khalili J (2014)  Life on the edge: the coming of age of quantum biology London: Bantam Press.
    [90] McFadden J, Al-Khalili J (1999) A quantum mechanical model of adaptive mutation. Biosystems 50: 203-211. doi: 10.1016/S0303-2647(99)00004-0
    [91] Misra B, Sudarshan ECG (1977) The Zeno's paradox in quantum theory. J Math Phys 18: 756-763. doi: 10.1063/1.523304
    [92] Kraus K (1981) Measuring processes in quantum mechanics I. Continuous observation and the watchdog effect. Found Phys 11: 547-576. doi: 10.1007/BF00726936
    [93] Meloni M (2014) The social brain meets the reactive genome: neuroscience, epigenetics and the new social biology. Front Hum Neurosci 8: 309. doi: 10.3389/fnhum.2014.00309
    [94] Vitiello G (2014) On the isomorphism between dissipative systems, fractal self-similarity and electrodynamics. Toward an integrated vision of nature. Systems 2: 203-216. doi: 10.3390/systems2020203
    [95] Novack DH, Cameron O, Epel E, et al. (2007) Psychosomatic medicine: The scientific foundation of the biopsychosocial model. Acad Psychiatry 31: 388-401. doi: 10.1176/appi.ap.31.5.388
    [96] Dal Lin C, Gola E, Brocca A, et al. (2018) miRNAs may change rapidly with thoughts: The relaxation response after myocardial infarction. Eur J Integr Med 20: 63-72. doi: 10.1016/j.eujim.2018.03.009
    [97] Chrousos GP (2009) Stress and disorders of the stress system. Nat Rev Endocrinol 5: 374-381. doi: 10.1038/nrendo.2009.106
    [98] Muehsam D, Ventura C (2014) Life rhythm as a symphony of oscillatory patterns: electromagnetic energy and sound vibration modulates gene expression for biological signaling and healing. Glob Adv Health Med 3: 40-55. doi: 10.7453/gahmj.2014.008
    [99] (2015)  Heart Beat Made Visible on CymaScope Avaibale from: https://www.youtube.com/watch?v=2kuY98F7o_0.
    [100] Ingber DE, Wang N, Stamenović D (2014) Tensegrity, cellular biophysics, and the mechanics of living systems. Rep Prog Phys 77: 046603. doi: 10.1088/0034-4885/77/4/046603
    [101] Wang N, Tytell JD, Ingber DE (2009) Mechanotransduction at a distance: Mechanically coupling the extracellular matrix with the nucleus. Nat Rev Mol Cell Biol 10: 75-82. doi: 10.1038/nrm2594
    [102] Martino F, Perestrelo AR, Vinarsky V, et al. (2018) Cellular mechanotransduction: from tension to function. Frony Physiol 9: 824. doi: 10.3389/fphys.2018.00824
    [103] Buxbaum O (2016)  Key Insights into Basic Mechanisms of Mental Activity Switzerland: Springer International Publishing. doi: 10.1007/978-3-319-29467-4
    [104] Jamieson JP, Crum AJ, Goyer JP, et al. (2018) Optimizing stress responses with reappraisal and mindset interventions: an integrated model. Anxiety, Stress, Coping 31: 245-261. doi: 10.1080/10615806.2018.1442615
    [105] Pulvermüller F (2013) How neurons make meaning: Brain mechanisms for embodied and abstract-symbolic semantics. Trends Cogn Sci 17: 458-470. doi: 10.1016/j.tics.2013.06.004
    [106] Segall JM, Allen EA, Jung RE, et al. (2012) Correspondence between structure and function in the human brain at rest. Front Neuroinform 6: 10. doi: 10.3389/fninf.2012.00010
    [107] Alexander-Bloch A, Shou H, Liu S, et al. (2018) On testing for spatial correspondence between maps of human brain structure and function. Neuroimage 178: 540-551. doi: 10.1016/j.neuroimage.2018.05.070
    [108] Rebollo I, Devauchelle AD, Béranger B, et al. (2018) Stomach-brain synchrony reveals a novel, delayed-connectivity resting-state network in humans. Elife 7: e33321. doi: 10.7554/eLife.33321
    [109] Ventura C (2017) Seeing cell biology with the eyes of physics. NanoWorld J 3: S1-S8.
    [110] Fredericks S, Saylor JR (2013)  Shape oscillation of a levitated drop in an acoustic field, 2013 Available from: arXiv:1310.2967.
    [111] Zhang CY, Wang Y, Schubert R, et al. (2016) Effect of audible sound on protein crystallization. Cryst Growth Des 16: 705-713. doi: 10.1021/acs.cgd.5b01268
    [112] Guo F, Li P, French JB, et al. (2015) Controlling cell–cell interactions using surface acoustic waves. Proc Natl Acad Sci 112: 43-48. doi: 10.1073/pnas.1422068112
    [113] Vogel V, Sheetz M (2006) Local force and geometry sensing regulate cell functions. Nat Rev Mol Cell Biol 7: 265-275. doi: 10.1038/nrm1890
    [114] Shaobin G, Wu Y, Li K, et al. (2010) A pilot study of the effect of audible sound on the growth of Escherichia coli. Colloid Surface B 78: 367-371. doi: 10.1016/j.colsurfb.2010.02.028
    [115] Gu SB, Yang B, Wu Y, et al. (2013) Growth and physiological characteristics of E. coli in response to the exposure of sound field. Pakistan J Biol Sci 16: 969-975. doi: 10.3923/pjbs.2013.969.975
    [116] Sahu S, Ghosh S, Fujita D, et al. (2014) Live visualizations of single isolated tubulin protein self-assembly via tunneling current: effect of electromagnetic pumping during spontaneous growth of microtubule. Sci Rep 4: 7303. doi: 10.1038/srep07303
    [117] Acbas G, Niessen KA, Snell EH, et al. (2014) Optical measurements of long-range protein vibrations. Nat Commun 5: 3076. doi: 10.1038/ncomms4076
    [118] Christians ES, Benjamin IJ (2012) Proteostasis and REDOX state in the heart. Am J Physiol Heart Circ Physiol 302: H24-H37. doi: 10.1152/ajpheart.00903.2011
    [119] Christians ES, Mustafi SB, Benjamin IJ (2014) Chaperones and cardiac misfolding protein diseases. Curr Protein Pept Sci 15: 189-204. doi: 10.2174/1389203715666140331111518
    [120] Naviaux RK (2014) Metabolic features of the cell danger response. Mitochondrion 16: 7-17. doi: 10.1016/j.mito.2013.08.006
    [121] Crum A, Zuckerman B (2017) Changing mindsets to enhance treatment effectiveness. J Am Med Assoc 317: 2063-2064. doi: 10.1001/jama.2017.4545
    [122] Maas C, Belgardt D, Han KL, et al. (2009) Synaptic activation modifies microtubules underlying transport of postsynaptic cargo. Proc Natl Acad Sci 106: 8731-8736. doi: 10.1073/pnas.0812391106
    [123] Lo LP, Liu SH, Chang YC (2007) Assembling microtubules disintegrate the postsynaptic density in vitro. Cell Motil Cytoskeleton 64: 6-18. doi: 10.1002/cm.20163
    [124] Arimura N, Kaibuchi K (2007) Neuronal polarity: From extracellular signals to intracellular mechanisms. Nat Rev Neurosci 8: 194-205. doi: 10.1038/nrn2056
    [125] Macario AJL, Conway de Macario E (2000) Stress and molecular chaperones in disease. Int J Clin Lab Res 30: 49-66. doi: 10.1007/s005990070016
    [126] Dal Lin C, Marinova M, Brugnolo L, et al. Rapid senectome and alternative splicing miRNAs changes with the relaxation response: A one year follow-up study, 2020 Available from: Preprints doi:10.20944/preprints202007.0268.v1.
    [127] Picard M, McManus MJ, Gray JD, et al. (2015) Mitochondrial functions modulate neuroendocrine, metabolic, inflammatory, and transcriptional responses to acute psychological stress. Proc Natl Acad Sci 112: E6614-E6623. doi: 10.1073/pnas.1515733112
    [128] Piattelli-Palmarini M, Vitiello G (2015) Linguistics and some aspects of its underlying dynamics. Biolinguistics 9: 96-115.
    [129] Mańka R, Ogrodnik B (1991) A model of soliton transport along microtubules. J Biol Phys 18: 85-189. doi: 10.1007/BF00417807
    [130] Kučera O, Havelka D (2012) Mechano-electrical vibrations of microtubules-Link to subcellular morphology. BioSystems 109: 346-355. doi: 10.1016/j.biosystems.2012.04.009
    [131] Benias PC, Wells RG, Sackey-Aboagye B, et al. (2018) Structure and distribution of an unrecognized interstitium in human tissues. Sci Rep 8: 4947. doi: 10.1038/s41598-018-23062-6
    [132] Brizhik L, Chiappini E, Stefanini P, et al. (2019) Modeling meridians within the quantum field theory. J Acupunct Meridian Stud 12: 29-36. doi: 10.1016/j.jams.2018.06.009
    [133] Bentov I (1977)  Stalking the wild pendulum Glasgow: William Collins Sons & Co. Ltd.
    [134] Pavanello S, Campisi m, Tona F, et al. (2019) Exploring epigenetic age in response to intensive relaxing training: A pilot study to slow down biological age. Int J Env Res Pub He 16: 3074. doi: 10.3390/ijerph16173074
    [135] Dal Lin C, Grasso R, Scordino A, et al. (2020)  Ph, electric conductivity and delayed luminescence changes in human sera of subjects undergoing the relaxation response: A pilot study Available from: doi:10.20944/PREPRINTS202004.0202.V1.
    [136] Cifra M, Brouder C, Nerudová M, et al. (2015) Biophotons, coherence and photocount statistics: A critical review. J Lumin 164: 38-51. doi: 10.1016/j.jlumin.2015.03.020
    [137] Boveris A, Cadenas E, Reiter R (1980) Organ chemiluminescence: noninvasive assay for oxidative radical reactions. Proc Natl Acad Sci 77: 347-351. doi: 10.1073/pnas.77.1.347
    [138] Krasovitski B, Frenkel V, Shoham S (2011) Intramembrane cavitation as a unifying mechanism for ultrasound-induced bioeffects. Proc Natl Acad Sci 108: 3258-3263. doi: 10.1073/pnas.1015771108
    [139] Brujan EA (2000) Collapse of cavitation bubbles in blood. Europhys Lett 50: 175. doi: 10.1209/epl/i2000-00251-7
    [140] Brennen CE (2015) Cavitation in medicine. Interface Focus 5: 20150022. doi: 10.1098/rsfs.2015.0022
    [141] Didenko YT, Suslick KS (2002) The energy efficiency of formation of photons radicals and ions during single-bubble cavitation. Nature 418: 394-397. doi: 10.1038/nature00895
    [142] Sabbadini SA, Vitiello G (2019) Entanglement and phase-mediated correlations in quantum field theory. Application to brain-mind states. Appl Sci 9: 3203. doi: 10.3390/app9153203
    [143] Shaffer F, McCraty R, Zerr CL (2014) A healthy heart is not a metronome: an integrative review of the heart's anatomy and heart rate variability. Front Psychol 5: 1040. doi: 10.3389/fpsyg.2014.01040
    [144] Grippo A (2011) The utility of animal models in understanding links between psychosocial processes and cardiovascular health. Soc Pers Psychol Compass 5: 164-179. doi: 10.1111/j.1751-9004.2011.00342.x
    [145] Mensah G, Collins P (2015) Understanding mental health for the prevention and control of cardiovascular diseases. Glob Heart 10: 221. doi: 10.1016/j.gheart.2015.08.003
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(8017) PDF downloads(442) Cited by(0)

Figures and Tables

Figures(7)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog