Search

Background: During the dissolution of amorphous solid dispersion (ASD) formulations, the drug load (DL) often impacts the release mechanism and the occurrence of loss of release (LoR). The ASD/water interfacial gel layer and its specific phase behavior in connection with DL strongly dictate the release mechanism and LoR of ASDs, as reported in the literature. Thermodynamically driven liquid-liquid phase separation (LLPS) and/or drug crystallization at the interface are the key phase transformations that drive LoR. Methods: In this study, a combination of Perturbed-Chain Statistical Associating Fluid Theory (PC-SAFT) thermodynamic modeling and in silico molecular simulation was applied to investigate the release mechanism and the occurrence LoR of an ASD formulation consisting of ritonavir as the active pharmaceutical ingredient (API) and the polymer, polyvinylpyrrolidone-co-vinyl acetate (PVPVA64). A thermodynamically modeled ternary phase diagram of ritonavir (PVPVA64) and water was applied to predict DL-dependent LLPS in the ASD/water interfacial gel layer. Microscopic Erosion Time Testing (METT) was used to experimentally validate the phase diagram predictions. Additionally, in silico molecular simulation was applied to provide further insights into the phase separation, the release mechanism, and aggregation behavior on a molecular level. Results: Thermodynamic modeling, molecular simulation, and experimental results were consistent and complementary, providing evidence that ASD/water interactions and phase separation are essential factors driving the dissolution behavior and LoR at 40 wt% DL of the investigated ritonavir/PVPVA64 ASD system, consistent with previous studies. Conclusions: This study provides insights into the potential of blending thermodynamic modeling, molecular simulation, and experimental research to comprehensively understand ASD formulations. Such a combined approach can be leveraged as a computational framework to gain insights into the ASD dissolution mechanism, thereby facilitating in silico screening, designing, and optimization of formulations with the benefit of significantly reducing the number of experimental tests. [ABSTRACT FROM AUTHOR]

Dissipative particle dynamics (DPD) simulations is used to study the effect of Am/2BmAm/2 and H-shaped (Am/4)2Bm(Am/4)2 block copolymers on the interfacial properties of ternary blends. Our simulations show the following: (i) The capacity of block copolymers to diminish interfacial tension is closely linked to their compositions. With identical molecular weights and concentrations, H-shaped block copolymers outperform triblock copolymers in mitigating interfacial tension. (ii) The interfacial tension within the blends correlates positively with the escalation in H-shaped block copolymer molecular weight. This correlation suggests that H-shaped block copolymers featuring a low molecular weight demonstrate superior efficacy as compatibilizers when contrasted with those possessing a high molecular weight. (iii) Enhancing the concentration of H-shaped block copolymers fosters their accumulation at the interface, leading to a reduction in correlations between immiscible homopolymers and a consequent decrease in interfacial tension. (iv) As the length of the homopolymer chains increases, there is a concurrent elevation in interfacial tension, suggesting that H-shaped block copolymers perform more effectively as compatibilizers in blends characterized by shorter homopolymer chain lengths. These findings elucidate the associations between the efficacy of H-shaped block copolymer compatibilizers and their specific molecular characteristics. [ABSTRACT FROM AUTHOR]

Ultraviolet (UV) ink is a primary type of ink used in additive manufacturing with 3D inkjet printing. However, ink aggregation presents a challenge in nano-inkjet printing, affecting the stability and quality of the printing fluid and potentially leading to the clogging of nanometer-sized nozzles. This paper utilizes a Dissipative Particle Dynamics (DPD) simulation to investigate the aggregation behavior of alumina in a blend of 1,6-Hexanediol diacrylate (HDDA) and Trimethylolpropane triacrylate (TMPTA). By analyzing the effects of solid content, polymer component ratios, and dispersant concentration on alumina aggregation, the optimal ink formulation was identified. Compared to traditional experimental methods, DPD simulations not only reduce experimental costs and time but also reveal particle aggregation mechanisms that are difficult to explore through experimental methods, providing a crucial theoretical basis for optimizing ink formulations. This study demonstrates that alumina ceramic ink achieves optimal performance with a solid content of 20%, an HDDA-to-TMPTA ratio of 4:1, and 9% oleic acid as a dispersant. [ABSTRACT FROM AUTHOR]

Triblock copolymers such as styrene-b-(ethylene-co-butylene)-b-styrene (SEBS) have been widely used as an anion exchange membrane for fuel cells due to their phase separation properties. However, modifying the polymer architecture for optimized membrane properties is still challenging. This research develops a strategy to control the membrane morphology based on quaternized SEBS (SEBS-Q) by dual-tapering the interfacial block sequences. The structural and transport properties of SEBS-Q with various tapering styles at different hydration levels are systematically investigated by coarse-grained molecular simulations. The results show that the introduction of the tapered regions induces the formation of a bicontinuous water domain and promotes the diffusivity of the mobile components. The interplay between the solvation of the quaternary groups and the tapered fraction determines the conformation of polymer chains among the hydrophobic–hydrophilic subdomains. The strategy presented here provides a new path to fabricating fuel cell membranes with controlled microstructures. [ABSTRACT FROM AUTHOR]

An automated pipeline for comprehensive calculation of intermolecular interaction energies based on molecular force-fields using the Tinker molecular modelling package is presented. Starting with non-optimized chemically intuitive monomer structures, the pipeline allows the approximation of global minimum energy monomers and dimers, configuration sampling for various monomer–monomer distances, estimation of coordination numbers by molecular dynamics simulations, and the evaluation of differential pair interaction energies. The latter are used to derive Flory–Huggins parameters and isotropic particle–particle repulsions for Dissipative Particle Dynamics (DPD). The computational results for force fields MM3, MMFF94, OPLS-AA and AMOEBA09 are analyzed with Density Functional Theory (DFT) calculations and DPD simulations for a mixture of the non-ionic polyoxyethylene alkyl ether surfactant C10E4 with water to demonstrate the usefulness of the approach. Scientific Contribution To our knowledge, there is currently no open computational pipeline for differential pair interaction energies at all. This work aims to contribute an (at least academically available, open) approach based on molecular force fields that provides a robust and efficient computational scheme for their automated calculation for small to medium-sized (organic) molecular dimers. The usefulness of the proposed new calculation scheme is demonstrated for the generation of mesoscopic particles with their mutual repulsive interactions. [ABSTRACT FROM AUTHOR]

Crude oil desalination and dehydration are necessary for storage, transportation, and processing procedures. However, the behaviour of fine emulsion droplets under an electric field has always been questioned. This paper modified the dissipative particle dynamics method (DPD) to study the deformation process of fine emulsion droplets under a high-strength electric field. Compared with the literature data, the reliability of the DPD method is confirmed. The influence of the crude oil properties and the electric field characteristics on the behaviour of the emulsion droplet was analysed, and the effect factors included electric field intensity, electric field frequency, emulsion droplet size, centre distance ratio, conservative force intensity, dissipative strength, and crude oil density. The relationship between critical electric field intensity and emulsion droplet deformation was formulated based on the simulational dates. [ABSTRACT FROM AUTHOR]

To enhance the mechanical properties of the Nitrate Ester Plasticized Polyether solid propellant matrix, the uniaxial tension of multi-component systems is simulated and the factors influencing the mechanical properties of the propellant matrix are investigated. First, mesoscale models of five types of systems include poly alpha olefin (PAO(3)), polyethylene glycol (PEG200, PEG400, PEG600), and 1,4-butanediol (BDO) are established, followed by uniaxial tensile simulations. The results show PEG600, PEG400, PEG200, BDO, and PAO(3) in order of enhancing the mechanical performance of the matrix. Second, the diffusion behavior of nitroglycerin (NG) and butanetriol trinitrate (BTTN) in various systems is investigated. The results show that NG exhibits higher diffusion capacity than BTTN, and the diffusion coefficient increases with an increment in the molecular weight of PEG. Additionally, the influence of different plasticizer ratios (2.8–3.0), curing parameters (1.58–1.62), and chain extension parameters (0.08–0.10) on the mechanical properties of the PEG600 system are investigated. The results demonstrate that as the plasticizer ratio increases, there is a gradual decrease in the modulus of the matrix. Additionally, an increase in the curing parameter leads to a substantial enhancement in the tensile strength of the matrix, while increasing the chain extension parameter significantly expands the maximum tensile length of the matrix. Finally, employing the Slip-Spring model, the effects of the physical and chemical cross-linked network of the propellant are simulated. The result shows that increasing the content of a chemical cross-linked network significantly improves the tensile strength of the matrix.

Y-shaped polymer brushes represent a special class of binary mixed polymer brushes, in which a combination of different homopolymers leads to unique phase behavior. While most theoretical and simulation studies use monodisperse models, experimental systems are always polydisperse. This discrepancy hampers linking theoretical and experimental results. In this theoretical study, we employed dissipative particle dynamics to study the influence of polydispersity on the phase behavior of Y-shaped brushes grafted to flat surfaces under good solvent conditions. Polydispersity was kept within experimentally achievable values and was modeled via Schulz–Zimm distribution. In total, 10 systems were considered, thus covering the phase behavior of monodisperse, partially polydisperse and fully polydisperse systems. Using such generic representation of real polymers, we observed a rippled structure and aggregates in monodisperse systems. In addition, polydisperse brushes formed a stable perforated layer not observed previously in monodisperse studies, and influenced the stability of the remaining phases. Although the perforated layer was experimentally observed under good solvent conditions and in the melt state, further confirmation of its presence in systems under good solvent conditions required mapping real polymers onto mesoscale models that reflected, for example, different polymer rigidity, and excluded volume effects or direct influence of the surface, just to mention a few parameters. Finally, in this work, we show that mesoscale modeling successfully describes polydisperse models, which opens the way for rapid exploring of complex systems such as polydisperse Y-shaped brushes in selective or bad solvents or under non-equilibrium conditions. [ABSTRACT FROM AUTHOR]

Shape-memory polyurethanes (SMPUs) are promising materials that change shape in response to external heat. These polymers have a dual-segment structure: a hard segment for netpoint and a soft segment for molecular switch. Understanding the molecular behavior of each segment and microphase-separated morphology is crucial for comprehending the shape-memory mechanism. This study aimed to understand the shape-memory behavior by observing the phase separation of SMPU using mesoscale models based on dissipative particle dynamics (DPD) simulations. The SMPU copolymer was modeled using 4,4′-diphenylmethane diisocyanate (MDI, hard segment) and poly(ethylene oxide) (PEO, soft segment). By calculating segment solubility and repulsion parameters, we found that the hard-segment domain changes from isolated form to a lamellar and interconnected structure and eventually to a continuous form as its content increases. Combining these insights with shape-memory performance models can enhance our understanding of better SMPU design and contribute significantly to the optimization of smart stimuli-responsive materials.

This publication continues the cycle of our work aimed at improving the methodology for constructing mesoscale models of network polymers and characterizing their physical properties. As the object of study, the epoxy resin of bisphenol A, diglycidyl ether, and a tricarboxylic fatty acid hardener was chosen. Its structure is the result of three parallel reactions. For their correct reproduction, an algorithm was proposed, which allows to take into account the peculiarities of the relationship of all ongoing processes. The system model was constructed by mapping the chemical structure of the monomers onto an equivalent mesoscale representation. It was used to study the relationship between the structure and mechanical properties of feeding networks as a function of the ratio of volume fractions of comonomers in the initial reaction mixture. All calculations were performed within the reaction version of the dissipative particle dynamics method. The structure of polymer networks in the constructed samples was characterized by topological analysis. The study of mechanical properties was carried out by constructing the "stress-strain" dependencies. The results obtained show a good correlation between the density of the load-bearing chains and the mechanical properties of the resulting networks. It is shown that the material samples with the highest degree of transformation and the density of the number of load-bearing chains have the highest stiffness.

Dissipative particle dynamics (DPD) simulations is used to study the effect of Am/2BmAm/2 and H-shaped (Am/4)2Bm(Am/4)2 block copolymers on the interfacial properties of ternary blends. Our simulations show the following: (i) The capacity of block copolymers to diminish interfacial tension is closely linked to their compositions. With identical molecular weights and concentrations, H-shaped block copolymers outperform triblock copolymers in mitigating interfacial tension. (ii) The interfacial tension within the blends correlates positively with the escalation in H-shaped block copolymer molecular weight. This correlation suggests that H-shaped block copolymers featuring a low molecular weight demonstrate superior efficacy as compatibilizers when contrasted with those possessing a high molecular weight. (iii) Enhancing the concentration of H-shaped block copolymers fosters their accumulation at the interface, leading to a reduction in correlations between immiscible homopolymers and a consequent decrease in interfacial tension. (iv) As the length of the homopolymer chains increases, there is a concurrent elevation in interfacial tension, suggesting that H-shaped block copolymers perform more effectively as compatibilizers in blends characterized by shorter homopolymer chain lengths. These findings elucidate the associations between the efficacy of H-shaped block copolymer compatibilizers and their specific molecular characteristics.

Ultraviolet (UV) ink is a primary type of ink used in additive manufacturing with 3D inkjet printing. However, ink aggregation presents a challenge in nano-inkjet printing, affecting the stability and quality of the printing fluid and potentially leading to the clogging of nanometer-sized nozzles. This paper utilizes a Dissipative Particle Dynamics (DPD) simulation to investigate the aggregation behavior of alumina in a blend of 1,6-Hexanediol diacrylate (HDDA) and Trimethylolpropane triacrylate (TMPTA). By analyzing the effects of solid content, polymer component ratios, and dispersant concentration on alumina aggregation, the optimal ink formulation was identified. Compared to traditional experimental methods, DPD simulations not only reduce experimental costs and time but also reveal particle aggregation mechanisms that are difficult to explore through experimental methods, providing a crucial theoretical basis for optimizing ink formulations. This study demonstrates that alumina ceramic ink achieves optimal performance with a solid content of 20%, an HDDA-to-TMPTA ratio of 4:1, and 9% oleic acid as a dispersant.

Self-assembly of surfactants into complex structures is key to the performance of many formulated products, which form a significant fraction of the world’s manufactured goods. Here we adopt the dissipative particle dynamics simulation approach to explore the self-assembly process of surfactants, with the aim of understanding what information can be obtained that may correlate with an increased zero-shear viscosity of surfactant based products. To this end we experimentally measured the zero-shear viscosity of mixed micelle systems comprised of cocoamidopropyl betaine (CAPB) and sodium lauryl sarcosinate (SLSar), as a function of the CAPB/SLSar mass ratio and pH, and characterised the early stages of self-assembly of the same systems computationally. From simulation we identify three distinct behaviors in the micellar self-assembly process (logarithmic, linear and cubic growth) which we find show some degree of correlation with the experimental zero-shear viscosity. Owing to the relatively short simulation times required, this may provide formulation scientists with a practical route to identify regions of interest (i. e. those with a desired zero-shear viscosity) prior to synthesising de novo (potentially natural) surfactants. [ABSTRACT FROM AUTHOR]

The commercialized lithium‐ion battery separators are mass‐produced by mixing ultra‐high molecular weight polyethylene (UHMWPE) and paraffin oil (PO). The dissipative particle dynamics method is utilized to investigate the extrinsic factors (shear rate and cooling rate) and the intrinsic factors (the molecular chain length) on the microstructure of the UHMWPE‐PO mixture. For the mixture with UHMWPE possessing the same chain length, the high shear rate promoted a lower porosity (~28%) and smaller pores. In contrast, the slow shearing led to a high porosity (~40%) and larger pores. For the mixture with UHMWPE possessing short and long chains, the shear rate hardly affects the porosity and the pore size: the porosity was kept at ~30%, and the pore size was reduced by ~35% compared to the model with the same‐chain‐length UHMWPE. The cooling rate after shearing is the dominant factor in determining the porosity and pore size: the fast cooling raised the porosity by ~33% but hardly increased the pore size, while the slow cooling raised the porosity by ~74%, and the pore size by ~105%. The current study provided a deep understanding of the pore structure evolution in separator processing. Highlights: The effects of processing parameters on the pore structures are numerically illustrated.MD simulation and rheometer measurement assist DPD interaction parameters calibration for UHMWPE and PO.The low shear rate leads to a higher porosity and pore size.At the high shear rate, short UHMWPE chains reduce porosity but do not increase pore size.The fast‐cooling process slightly increases the porosity while keeping the pore size. [ABSTRACT FROM AUTHOR]

In this research, we have used the dissipative particle dynamics (DPD), a mesoscopic simulation technique, in order to investigate the dynamics of lipid domains in near critical temperature. Our specific focus has been on exploring the influence of lipid domain size on its lifetime, which mimics the behavior of lipid rafts within cellular membranes. The lipid membranes used in this study were composed of saturated and unsaturated lipids, which have been immersed in water. Through the simulation of these membranes close to their critical temperature, we have successfully generated fluctuating domains that mimic the lipid rafts observed in cellular systems. We have proposed a method to obtain the lifetime of the fluctuating domains by analyzing the sizes of the lipid domains at specific intervals of time. Our investigations have revealed a linear correlation between the initial size of the lipid domain and its lifetime. Our research finding give an insight into the underlying mechanisms that govern lipid rafts and their vital role in various cellular processes.

Triblock copolymers such as styrene-b-(ethylene-co-butylene)-b-styrene (SEBS) have been widely used as an anion exchange membrane for fuel cells due to their phase separation properties. However, modifying the polymer architecture for optimized membrane properties is still challenging. This research develops a strategy to control the membrane morphology based on quaternized SEBS (SEBS-Q) by dual-tapering the interfacial block sequences. The structural and transport properties of SEBS-Q with various tapering styles at different hydration levels are systematically investigated by coarse-grained molecular simulations. The results show that the introduction of the tapered regions induces the formation of a bicontinuous water domain and promotes the diffusivity of the mobile components. The interplay between the solvation of the quaternary groups and the tapered fraction determines the conformation of polymer chains among the hydrophobic–hydrophilic subdomains. The strategy presented here provides a new path to fabricating fuel cell membranes with controlled microstructures.

Self-assembly of surfactants into complex structures is key to the performance of many formulated products, which form a significant fraction of the world’s manufactured goods. Here we adopt the dissipative particle dynamics simulation approach to explore the self-assembly process of surfactants, with the aim of understanding what information can be obtained that may correlate with an increased zero-shear viscosity of surfactant based products. To this end we experimentally measured the zero-shear viscosity of mixed micelle systems comprised of cocoamidopropyl betaine (CAPB) and sodium lauryl sarcosinate (SLSar), as a function of the CAPB/SLSar mass ratio and pH, and characterised the early stages of self-assembly of the same systems computationally. From simulation we identify three distinct behaviors in the micellar self-assembly process (logarithmic, linear and cubic growth) which we find show some degree of correlation with the experimental zero-shear viscosity. Owing to the relatively short simulation times required, this may provide formulation scientists with a practical route to identify regions of interest (i. e. those with a desired zero-shear viscosity) prior to synthesising de novo (potentially natural) surfactants.

Although the behavior of the confined semi-dilute solutions of self-assembling copolymers represents an important topic of basic and applied research, it has eluded the interest of scientists. Extensive series of dissipative particle dynamics simulations have been performed on semi-dilute solutions of A5B5 chains in a selective solvent for A in slits using a DL-MESO simulation package. Simulations of corresponding bulk systems were performed for comparison. This study shows that the associates in the semi-dilute bulk solutions are partly structurally organized. Mild steric constraints in slits with non-attractive walls hardly affect the size of the associates, but they promote their structural arrangement in layers parallel to the slit walls. Attractive walls noticeably affect the association process. In slits with mildly attractive walls, the adsorption competes with the association process. At elevated concentrations, the associates start to form in wide slits when the walls are sparsely covered by separated associates, and the association process prevents the full coverage of the surface. In slits with strongly attractive walls, adsorption is the dominant behavior. The associates form in wide slits at elevated concentrations only after the walls are completely and continuously covered by the adsorbed chains. [ABSTRACT FROM AUTHOR]

A hexameric quaternary ammonium salt surfactant Tz-6C12QC featuring a rigid triazine spacer and six ammonium groups was synthesized. The molecular structure and aggregation behavior of Tz-6C12QC were characterized by nuclear magnetic resonance spectroscopy, surface tension, conductivity, dynamic light scattering, and transmission electron microscopy, etc. Dissipative particle dynamics (DPD) simulation was employed to investigate the self-assembly behavior of Tz-6C12QC at different concentrations. The rheological behavior of the polyacrylamide/Tz-6C12QC system was characterized by shear rheology. The results indicated that Tz-6C12QC exhibited superior surface activity and lower surface tension compared to conventional surfactants. Rheology analysis revealed that Tz-6C12QC had a significant viscosity reduction effect on polyacrylamide. DLS and TEM indicated that, as the concentration of Tz-6C12QC increased, monomer associations, spherical aggregations, vesicles, tubular micelles, and bilayer vesicles were sequentially formed in the solution. This study presents a synthetic approach for polysurfactants with a rigid spacer and sheds light on the self-assembly process of micelles. [ABSTRACT FROM AUTHOR]

The reaction of a biphasic catalysis with microgels bearing catalytic groups adsorbed at liquid interface was simulated for the first time using dissipative particle dynamics. It was shown that the rate of the catalytic process increases with the degree of deformation of the polymer network, which depends on the fraction of the crosslinker and the solubility of the polymer in both phases. It was found that the highest reaction rate was observed when the microgel was soluble in both phases due to an increase in its porosity (in comparison with amphiphilic microgels) and in the water–microgel–oil contact area with a simultaneous decrease in the time for the reagents to reach the catalytic groups due to the flattening of the microgel. The results obtained can be useful for increasing the efficiency of a wide range of catalytic reactions of the type considered through the use of network-like macromolecules. [ABSTRACT FROM AUTHOR]

Y-shaped polymer brushes represent a special class of binary mixed polymer brushes, in which a combination of different homopolymers leads to unique phase behavior. While most theoretical and simulation studies use monodisperse models, experimental systems are always polydisperse. This discrepancy hampers linking theoretical and experimental results. In this theoretical study, we employed dissipative particle dynamics to study the influence of polydispersity on the phase behavior of Y-shaped brushes grafted to flat surfaces under good solvent conditions. Polydispersity was kept within experimentally achievable values and was modeled via Schulz–Zimm distribution. In total, 10 systems were considered, thus covering the phase behavior of monodisperse, partially polydisperse and fully polydisperse systems. Using such generic representation of real polymers, we observed a rippled structure and aggregates in monodisperse systems. In addition, polydisperse brushes formed a stable perforated layer not observed previously in monodisperse studies, and influenced the stability of the remaining phases. Although the perforated layer was experimentally observed under good solvent conditions and in the melt state, further confirmation of its presence in systems under good solvent conditions required mapping real polymers onto mesoscale models that reflected, for example, different polymer rigidity, and excluded volume effects or direct influence of the surface, just to mention a few parameters. Finally, in this work, we show that mesoscale modeling successfully describes polydisperse models, which opens the way for rapid exploring of complex systems such as polydisperse Y-shaped brushes in selective or bad solvents or under non-equilibrium conditions.

The temporal and spatial evolution of shear banding during startup and steady-state shear flow was studied for solutions of entangled, linear, monodisperse polyethylene C 3000 H 6002 dissolved in hexadecane and benzene solvents. A high-fidelity coarse-grained dissipative particle dynamics method was developed and evaluated based on previous NEMD simulations of similar solutions. The polymeric contribution to shear stress exhibited a monotonically increasing flow curve with a broad stress plateau at intermediate shear rates. For startup shear flow, transient shear banding was observed at applied shear rates within the steady-state shear stress plateau. Shear bands were generated at strain values where the first normal stress difference exhibited a maximum, with lifetimes persisting for up to several hundred strain units. During the lifetime of the shear bands, an inhomogeneous concentration distribution was evident within the system, with higher polymer concentration in the slow bands at low effective shear rate; i.e., γ ˙ < τ R − 1 , and vice versa at high shear rate. At low values of applied shear rate, a reverse flow phenomenon was observed in the hexadecane solution, which resulted from elastic recoil of the molecules within the slow band. In all cases, the shear bands dissipated at high strains and the system attained steady-state behavior, with a uniform, linear velocity profile across the simulation cell and a homogeneous concentration. [ABSTRACT FROM AUTHOR]

Most of the existing hemolysis mechanism studies are carried out on the macro flow scale. They assume that the erythrocyte membranes with different loads will suffer the same damage, which obviously has limitations. Thus, exploring the hemolysis mechanism through the macroscopic flow field information is a tough challenge. In order to further understand the non-physiological shear hemolysis phenomenon at the cell scale, this study used the coarse-grained erythrocytes damage model at the mesoscopic scale based on the transport dissipative particle dynamics (tDPD) method. Combined with computational fluid dynamics the hemolysis of scalarized shear stress (1) in the clearance of "Impella 5.0" was evaluated under the Lagrange perspective and Euler perspective. The results from the Lagrange perspective showed that the change rate of scaled shear stress (...τ) was the most critical factor in damaging RBCs in the rotor region of "Impella 5.0 and other transvalvular micro-axial blood pumps. Then, we propose a dimensionless number Dk with time integration based on ...τ to evaluate hemolysis. The Dissipative particle dynamics simulation results are consistent with the Dk, evaluation results, so ...τ may be an important factor in the hemolysis of VADs. Finally, we tested the hemolysis of 30% hematocrit whole blood in the "Impella 5.0" shroud clearance from the Euler perspective. Relevant results indicate that because of the wall effect, the RBCs near the impeller side are more prone to damage, and most of the cytoplasm is also gathered at the rotor side. [ABSTRACT FROM AUTHOR]

We consider a functional surface comprising thermoresponsive polymer chains, the material that has found numerous technological and biomedical applications. However, to achieve the required time and length scales for computer modelling of such applications, one is compelled to use coarse-grained mesoscopic modelling approaches. The model used here is based on the previous work [Soto-Figueroa et al., Soft Matter, 8, 1871 (2012)], and it mimics the principal feature of the poly(N-iso-propylacrylamide) (PNIPAM), namely, the rapid change of its hydrophilicity at the lower critical solution temperature (LCST). For the case of an isolated chain, we discuss scaling properties of the radius of gyration, end-to-end distance, various distribution functions, and the density profile of monomers below and above the LCST. For the case of the model thermoresposive brush, we search for the optimum grafting density at which the change in the brush height, upon crossing the LCST, reaches its maximum value. The interpretation of the thermoresponse, in terms of the Alexander-de Gennes blobs and the level of solvation of polymer chains in a brush, is provided. [ABSTRACT FROM AUTHOR]

In this paper, the two-dimensional Kelvin–Helmholtz (KH) instability occurring in the shear flow of polymer fluids is modeled by the dissipative particle dynamics (DPD) method at the coarse-grained molecular level. A revised FENE model is proposed to properly describe the polymer chains. In this revised model, the elastic repulsion and tension are both considered between the adjacent beads, the bond length of which is set as one segment's equilibrium length. The entanglements between polymer chains are described with a bead repulsive potential. The characteristics of such a KH instability in polymer fluid shear flow can be successfully captured in the simulations by the use of the modified FENE model. The numerical results show that the waves and vortexes grow more slowly in the shear flow of the polymer fluids than in the Newtonian fluid case, these vortexes become flat, and the polymer impedes the mixing of fluids and inhibits the generation of turbulence. The effects of the polymer concentration, chain length, and extensibility are also investigated regarding the evolution of KH instability. It is shown that the mixing of two polymer fluids reduces, and the KH instability becomes more suppressed as the polymer concentration increases. The vortexes become much longer with the evolution of the elongated interface as the chain length turns longer. As the extensibility increases, the vortexes become more flattened. Moreover, the roll-up process is significantly suppressed if the polymer has sufficiently high extensibility. These observations show that the polymer and its properties significantly influence the formation and evolution of the coherent structures such as the waves and vortexes in the KH instability progress. [ABSTRACT FROM AUTHOR]

Abstract The influence of molecular fragmentation and parameter settings on a mesoscopic dissipative particle dynamics (DPD) simulation of lamellar bilayer formation for a C10E4/water mixture is studied. A “bottom-up” decomposition of C10E4 into the smallest fragment molecules (particles) that satisfy chemical intuition leads to convincing simulation results which agree with experimental findings for bilayer formation and thickness. For integration of the equations of motion Shardlow’s S1 scheme proves to be a favorable choice with best overall performance. Increasing the integration time steps above the common setting of 0.04 DPD units leads to increasingly unphysical temperature drifts, but also to increasingly rapid formation of bilayer superstructures without significantly distorted particle distributions up to an integration time step of 0.12. A scaling of the mutual particle–particle repulsions that guide the dynamics has negligible influence within a considerable range of values but exhibits apparent lower thresholds beyond which a simulation fails. Repulsion parameter scaling and molecular particle decomposition show a mutual dependence. For mapping of concentrations to molecule numbers in the simulation box particle volume scaling should be taken into account. A repulsion parameter morphing investigation suggests to not overstretch repulsion parameter accuracy considerations.

A self-healing epoxy material is considered, based on bisphenol A diglycidyl ether and atricarboxylic fatty acid hardener, belonging to a new class of polymers called vitrimers. The res toration of the integrity of such systems in the case of a damage occurs due to the exchange reaction of covalent bonds between the comonomers forming a polymer network. In our previous work, we have developed a model of this material based on the method of reactive dissipative particle dynamics. In this work, we apply our model to study the welding process of vitrimer samples cut into two parts. The control of the integrity of the structure of the systems was carried out using a topological analysis by calculating the distributions over the lengths of simple cycles and the density of the number of load-bearing circuits. It has been shown that the rate of restoration of the integrity of the systems is determined by the concentration of the catalyst and the degree of crosslinking of the polymer. The results obtained also indicate that in the case of a high degree of crosslinking of the polymer, as well as a low catalyst concentration, the structure of the system is highly inhomogeneous.

The transient elongational viscosity η e (t) of the polymer melt is known to exhibit strain hardening, which depends on the strain rate ε ˙ . This phenomenon was elucidated by the difference of chain stretching in the entanglement network between extension and shear. However, to date, the microscopic evolution of polymer melt has not been fully statistically analyzed. In this work, the radial size distributions P( R g , t ) of linear polymers are explored by dissipative particle dynamics during the stretching processes. In uniaxial extensional flow, it is observed that the mean radius of gyration R ¯ g (t) and standard deviation σ (t) remain unchanged until the onset of strain hardening, corresponding to linear viscoelasticity. Both R ¯ g and σ rise rapidly in the non-linear regime, and bimodal size distribution can emerge. Moreover, the onset of strain hardening is found to be insensitive to the Hencky strain ( ε ˙ H t ) and chain length (N). [ABSTRACT FROM AUTHOR]

Polymers can increase the viscosity of water, reduce the relative permeability of the water phase, and enhance the flowability of the oil phase; surfactants can form molecular films at the oil–water interface boundaries, thereby reducing interfacial tension. Surfactant/polymer (S/P) flooding technology for enhancing oil recovery has become a major way to increase crude oil production. This study used dissipative particle dynamics (DPD) technology to simulate the emulsification process of a four-component composite system consisting of oil, water, sodium dodecylbenzene sulfonate (SDBS), and partially hydrolyzed polyacrylamide (HPAM). By changing the concentration of the S/P system, the effect on emulsification behavior was analyzed. Combined with particle distribution diagrams and interfacial tension parameters, the effect of the emulsification behavior on the performance of the S/P binary system was analyzed. On this basis, the effect of different emulsion performances on the recovery factor was evaluated using micro-experiments. The study found that the S/P system that produced stable emulsification had a lower interfacial tension and relatively good effect on improving the recovery factor. Increasing the concentration of the polymer and surfactant may cause changes in the interfacial film of the emulsion, thereby affecting the ability of the S/P system to reduce interfacial tension and may not improve the oil recovery factor. The research results help to better analyze and screen the S/P system used for oil extraction and improve crude oil recovery. [ABSTRACT FROM AUTHOR]

[Display omitted] Multicore block copolymer aggregates correspond to self-assembly such that the polymer system spontaneously phase separates to multiple, droplet-like cores differing in the composition from the polymer surroundings. Such multiple core aggregates are highly useful capsules for different applications, e.g., drug transport, catalysis, controlled solvation, and chemical reactions platforms. We postulate that polymer system composition provides a direct means for designing polymer systems that self-assemble to such morphologies and controlling the assembly response. Using dissipative particle dynamics (DPD) simulations, we examine the self-assembly of a mixture of highly and weakly solvophobic homopolymers and an amphiphilic block copolymer in the presence of solvent. We map the multicore vs single core (core–shell particles) assembly response and aggregate structure in terms of block copolymer concentration, polymer component ratios, and chain length of the weakly solvophobic homopolymer. For fixed components and polymer chemistries, the amount of block copolymer is the key to controlling single core vs multicore aggregation. We find a polymer system dependent critical copolymer concentration for the multicore aggregation and that a minimum level of incompatibility between the solvent and the weakly solvophobic component is required for multicore assembly. We discuss the implications for polymer system design for multicore assemblies. In summary, the study presents guidelines to produce multicore aggregates and to tune the assembly from multicore aggregation to single core core–shell particles. [ABSTRACT FROM AUTHOR]

Hierarchical self-assembly is efficient in transferring functionality from small to large scales. Polymer nanoparticles with regularity are highly desirable because they are ideal units of hierarchical self-assembly. Herein, we proposed a crosslinking strategy for generating nanoparticles from ABC linear triblock copolymers and conducted a computational study on the self-assembly of the nanoparticles. The effect of nanoparticle symmetry on aggregation behaviors was examined. We found that the symmetric patchy nanoparticles can form various low-curvature aggregates tuned by the crosslinking numbers and A/C block lengths. Increasing the asymmetry of patchy nanoparticles causes the increased curvature of assemblies, forming well-defined structures such as loops, vesicles, and wavy ribbons. Furthermore, the self-assembly kinetics reveals a step-growth nature, resembling the step polymerization for one-dimensional growth. The work demonstrates that ABC-type polymeric nanoparticles are versatile units for building hierarchical assemblies, considerably enriching the self-assembly strategies.

Polyolefins account for more than half of global primary polymer production, however only a small fraction of these polymers are currently being recycled. Fragmentation of polymer chains into shorter chains with a targeted molecular weight distribution with the goal of reusing these fragments in subsequent chemical synthesis can potentially introduce an alternative approach to polyolefins recycling. Herein we develop a mesoscale framework to model degradation of polyethylene melts at a range of high temperatures. We use the dissipative particle dynamics approach with modified segmental repulsive potential to model the process of random scission in melts of linear polymer chains. We characterize the fragmentation process by tracking the time evolution of the distribution of degrees of polymerization of chain fragments. Specifically, we track the weight average and the number average degrees of polymerization and dispersity of polymer fragments as a function of the fraction of bonds broken. Furthermore, we track the number fraction distribution and the weight fraction distribution of polymer fragments with various degrees of polymerization as functions of the fraction of bonds broken for a range of high temperatures. Our results allow one to quantify to what extent the distribution of polymer chain fragments during random scission can be captured by the respective analytical distributions for the range of conversions considered. Understanding the thermal degradation of polyolefins on the mesoscale can result in the development of alternative strategies for recycling a range of thermoplastics.

The influence of molecular fragmentation and parameter settings on a mesoscopic dissipative particle dynamics (DPD) simulation of lamellar bilayer formation for a C10E4/water mixture is studied. A "bottom-up" decomposition of C10E4 into the smallest fragment molecules (particles) that satisfy chemical intuition leads to convincing simulation results which agree with experimental findings for bilayer formation and thickness. For integration of the equations of motion Shardlow's S1 scheme proves to be a favorable choice with best overall performance. Increasing the integration time steps above the common setting of 0.04 DPD units leads to increasingly unphysical temperature drifts, but also to increasingly rapid formation of bilayer superstructures without significantly distorted particle distributions up to an integration time step of 0.12. A scaling of the mutual particle–particle repulsions that guide the dynamics has negligible influence within a considerable range of values but exhibits apparent lower thresholds beyond which a simulation fails. Repulsion parameter scaling and molecular particle decomposition show a mutual dependence. For mapping of concentrations to molecule numbers in the simulation box particle volume scaling should be taken into account. A repulsion parameter morphing investigation suggests to not overstretch repulsion parameter accuracy considerations. [ABSTRACT FROM AUTHOR]

Context: Lignin has an immense potential for the production of lignin-based functional materials. In this work, effect of 2-chloro-ethyltrimethyl ammonium chloride (AC)-grafted alkali lignin (AL) on the morphologies in water was investigated by dissipative particle dynamics (DPD) simulations. The results showed that AL molecules formed spherical micelles, but the corresponding phenylpropane units of AL were randomly distributed in spherical micelles. However, AC-grafted modification of phenolic hydroxyl groups in AL led to the formation of multilamellar spherical micelles. The formation of multilamellar spherical micelles of AL mainly went through four stages: small clusters, larger aggregates with a core-shell structure, trilaminar, and multilamellar spherical micelles. AL molecules resulted in dimethomorph molecules being randomly distributed in the spherical micelle, while the dimethomorph molecules were perfectly entrapped into the spherical micelles of AC-grafted AL. Various molecular weights of AL had no effect on the formation and size of multilamellar spherical micelles. With increasing the content of AC-grafted AL, small clusters, multilamellar spherical micelles, tube-like, and lamellar aggregates were observed successively. This work highlights the potential of lignin to prepare monodispersed lignin-based functional colloidal spheres. Methods: Coarse-grained beads were performed energy minimization, geometric optimization, NPT ensemble (298 K and 1.0 bar), and NVT ensemble (298 K) calculations. DPD simulations were carried out at 300,000 steps in a 30×30×30 Rc3 cubic box with Materials Studio 7.0 program. [ABSTRACT FROM AUTHOR]

[Display omitted] Polymer microgels have proven to be highly promising macromolecular objects for a wide variety of applications. In particular, the soft particles of an anisotropic (rod-like) shape are of special interest because of their potential use in tissue engineering or materials design. However, a little is known about the physical behavior of such microgels in solution, which inspired us to study them using mesoscopic computer simulations. For single networks, depending on the solvent quality, the dimensional characteristics were obtained for microgels of different molecular weight, crosslinking density and aspect ratio. In particular, the conditions for the rod-to-rod (preserving the nonspherical shape) and rod-to-sphere collapse were found. In addition, the effect of the liquid-crystalline (LC) ordering was demonstrated for the ensemble of rod-like microgels at different swelling ratios, and the influence of microgel aspect ratio on the volume fraction of the LC transition was shown. [ABSTRACT FROM AUTHOR]

In the current work, atom transfer radical polymerization-induced self-assembly (ATRP PISA) phase diagrams were obtained by the means of dissipative particle dynamics simulations. A fast algorithm for determining the equilibrium morphology of block copolymer aggregates was developed. Our goal was to assess how the chemical nature of ATRP affects the self-assembly of diblock copolymers in the course of PISA. We discovered that the chain growth termination via recombination played a key role in determining the ATRP PISA phase diagrams. In particular, ATRP with turned off recombination yielded a PISA phase diagram very similar to that obtained for a simple ideal living polymerization process. However, an increase in the recombination probability led to a significant change of the phase diagram: the transition between cylindrical micelles and vesicles was strongly shifted, and a dependence of the aggregate morphology on the concentration was observed. We speculate that this effect occurred due to the simultaneous action of two factors: the triblock copolymer architecture of the terminated chains and the dispersity of the solvophobic blocks. We showed that these two factors affected the phase diagram weakly if they acted separately; however, their combination, which naturally occurs during ATRP, affected the ATRP PISA phase diagram strongly. We suggest that the recombination reaction is a key factor leading to the complexity of experimental PISA phase diagrams. [ABSTRACT FROM AUTHOR]

The intricate process of bubble growth dynamics involves a broad spectrum of physical phenomena from microscale mechanics of bubble formation to macroscale interplay between bubbles and surrounding thermo-hydrodynamics. Traditional bubble dynamics models including atomistic approaches and continuum-based methods segment the bubble dynamics into distinct scale-specific models. To bridge the gap between microscale stochastic fluid models and continuum-based fluid models for bubble dynamics, we develop a composite neural operator model to unify the analysis of nonlinear bubble dynamics across microscale and macroscale regimes by integrating a many-body dissipative particle dynamics (mDPD) model with a continuum-based Rayleigh–Plesset (RP) model through a novel neural network architecture, which consists of a deep operator network for learning the mean behavior of bubble growth subject to pressure variations and a long short-term memory network for learning the statistical features of correlated fluctuations in microscale bubble dynamics. Training and testing data are generated by conducting mDPD and RP simulations for nonlinear bubble dynamics with initial bubble radii ranging from 0.1 to 1.5 micrometers. The results show that the trained composite neural operator model can accurately predict bubble dynamics across scales, with a 99% predictive accuracy for the time evolution of the bubble radius under varying external pressure while containing correct size-dependent stochastic fluctuations in microscale bubble growth dynamics. The composite neural operator is the first deep learning surrogate for multiscale bubble growth dynamics that can capture correct stochastic fluctuations in microscopic fluid phenomena, which sets a new direction for future research in multiscale fluid dynamics modeling. [Display omitted] • We developed a composite neural network model for multiscale bubble dynamics. • Our model offers seamless scale-bridging in multiscale bubble dynamics. • It enables learning of scale-dependent stochastic processes on top of mean-field bubble growth. [ABSTRACT FROM AUTHOR]

As a complex component of crude oil, asphaltenes determine the stability of emulsion. However, the mechanism of asphaltenes stabilize emulsion is still insufficient recognition due to the complexity of the molecular structure of asphaltenes. In this paper, the experiments and dissipative particle dynamics simulations were used to study the interfacial properties of asphaltenes. The adsorption of asphaltenes at the oil–water interface was studied by interfacial tension (IFT) experiments. A realistic molecular structure of asphaltene was constructed on the results of the characterization of the asphaltenes. The effects of asphaltenes concentration, temperature, and shear on the formation and stability of emulsion based on this molecular structure were investigated by simulations. The results show that the IFT decreases as asphaltenes concentration increases, while the droplet size of the emulsion decreases and becomes more uniformly distributed. Increasing the temperature enhances the interfacial activity of the asphaltenes, but also promotes its dissolution in the oil phase, which has a negative effect on reducing the IFT value and the stability of the emulsion. As the shear rate increases, the configuration of the emulsion changes from droplets to ellipsoids and rods, and the structure of the interfacial film is disrupted. This work provides useful insights into the mechanisms of asphaltenes concentration, temperature, and shear on the stability of the emulsion. [Display omitted] [ABSTRACT FROM AUTHOR]

In this study, submicron poly(methyl methacrylate) (PMMA) particles with surfaces modified using poly(vinyl alcohol) (PVA) (PMMA-PVA) were produced to formulate solid-stabilized oil-in-water (O/W) emulsions including high internal phase emulsions with a dispersed phase of over 74 %. The high-internal phase emulsions stabilized by PMMA–PVA exhibited gel-like behavior, and the oil droplets exhibited faceted shapes. PVA with a degree of saponification of 80 and 88 mol% was used as the surface modifier rather than PVA with > 96 mol%, assisting in the formation of a robust adsorbed layer on the PMMA particles. Coarse-grained molecular simulations revealed that the partially-hydrolyzed PVA exhibited a better adsorption behavior compared to the fully-hydrolyzed PVA, being thermodynamically adsorbed onto the PMMA surface in water owing to stable adsorption via its hydrophobic acetyl groups. Furthermore, PMMA–PVA reduced the interfacial tension as the degree of saponification of PVA decreased. Simulation and experimental results suggested that the acetyl groups of PVA not only adsorbed onto the PMMA surface but also promoted adsorption between PMMA and the oil–water interface. Overall, the findings of this study underscore the potency of PMMA–PVA as a particulate stabilizer for formulating O/W high internal phase emulsions. [Display omitted] • Oil-in-water emulsions were successfully formulated using PMMA–PVA particles. • The gel-like high-internal-phase emulsion contained polyhedral oil droplets. • Partially hydrolyzed PVA adsorbed more strongly on PMMA than fully hydrolyzed PVA. • Acetyl groups of PVA facilitated particle adsorption onto the oil–water interface. [ABSTRACT FROM AUTHOR]

[Display omitted] • Using dissipative particle dynamics, investigate how the alkanes in crude oil are extracted. • The effect of extraction has on the miscible and immiscible displacement of carbon dioxide and crude oil. • The images of the extraction are exhibited at the molecular level. Understanding the phase behavior of a CO 2 -oil–water surfactant system is critical for carbon dioxide flooding engineering. Using dissipative particle dynamics, we investigate how the alkanes in crude oil are extracted after the injection of supercritical carbon dioxide into the formation, and what effect this extraction has on the miscible and immiscible displacement of carbon dioxide and crude oil. The images of the extraction are exhibited at the molecular level, and the effects of temperature and pressure on the extraction are investigated. The phase distribution of carbon dioxide, oil, and water under immiscible phase is simulated using actual reservoir fluid data based on a comprehensive understanding of extraction. Oil tends to be spread at the interface of carbon dioxide and water when it is immiscible. The miscible pictures were generated under conditions of reduced CO 2 -oil–water miscible pressure by adding surfactants to the CO 2 -oil–water system, and the fluid distribution characteristics under immiscible and miscible phases were examined. dissipative particle dynamics (DPD) simulation demonstrates that when the molar proportion of carbon dioxide is less than 0.4, the solid phase surface can always create a water film of varied thickness, regardless of whether the surface charge of the solid phase is positive or negative. As the molar fraction of carbon dioxide increases, the water film's thickness decreases. [ABSTRACT FROM AUTHOR]

A hexameric quaternary ammonium salt surfactant Tz-6C12QC featuring a rigid triazine spacer and six ammonium groups was synthesized. The molecular structure and aggregation behavior of Tz-6C12QC were characterized by nuclear magnetic resonance spectroscopy, surface tension, conductivity, dynamic light scattering, and transmission electron microscopy, etc. Dissipative particle dynamics (DPD) simulation was employed to investigate the self-assembly behavior of Tz-6C12QC at different concentrations. The rheological behavior of the polyacrylamide/Tz-6C12QC system was characterized by shear rheology. The results indicated that Tz-6C12QC exhibited superior surface activity and lower surface tension compared to conventional surfactants. Rheology analysis revealed that Tz-6C12QC had a significant viscosity reduction effect on polyacrylamide. DLS and TEM indicated that, as the concentration of Tz-6C12QC increased, monomer associations, spherical aggregations, vesicles, tubular micelles, and bilayer vesicles were sequentially formed in the solution. This study presents a synthetic approach for polysurfactants with a rigid spacer and sheds light on the self-assembly process of micelles.

Although the behavior of the confined semi-dilute solutions of self-assembling copolymers represents an important topic of basic and applied research, it has eluded the interest of scientists. Extensive series of dissipative particle dynamics simulations have been performed on semi-dilute solutions of A5B5 chains in a selective solvent for A in slits using a DL-MESO simulation package. Simulations of corresponding bulk systems were performed for comparison. This study shows that the associates in the semi-dilute bulk solutions are partly structurally organized. Mild steric constraints in slits with non-attractive walls hardly affect the size of the associates, but they promote their structural arrangement in layers parallel to the slit walls. Attractive walls noticeably affect the association process. In slits with mildly attractive walls, the adsorption competes with the association process. At elevated concentrations, the associates start to form in wide slits when the walls are sparsely covered by separated associates, and the association process prevents the full coverage of the surface. In slits with strongly attractive walls, adsorption is the dominant behavior. The associates form in wide slits at elevated concentrations only after the walls are completely and continuously covered by the adsorbed chains.

Multigeometry nanoparticles with high complexity in composition and structure have attracted significant attention for enhanced functionality. We assess a simple but versatile strategy to construct hybrid nanoparticles with subdivided geometries through the cooperative assembly of diblock copolymer blends with asymmetric architectures. We report the formation of multicompartmental, vesicular, cylindrical, and spherical structures from pure AB systems. Then, we explore the assemblies of binary AB/AC blends, where the two incompatible, hydrophobic diblock copolymers subdivide into self-assembled local geometries, and the complexity of the obtained morphologies increases. We expand the strategy to ternary AB/AC/AD systems by tuning the effect of phase separation of different hydrophobic domains on the surface or internal region of the nanoparticle. The kinetic control of the coassembly in the initial stage is crucial for controlling the final morphology. The interactions of copolymers with different block lengths and chemistries enable the stabilization of interfaces, rims and ends of subdomains in the hybrid multigeometry nanoparticles. With further exploration of size and shape, the dependence of local geometry on the volume fraction is discussed. We show an efficient approach for controllable multigeometry nanoparticle construction that will be useful for multifunctional and hierarchical nanomaterials.

Rational design of thermotropic liquid crystal (LC)-based sensors utilizing different mesophases holds great promise to open up novel detection modalities for various chemical and biological applications. In this context, we present a dissipative particle dynamics study to explore the unique anchoring behavior of nematic and smectic LCs at amphiphile-laden aqueous-LC interface. By increasing the surface coverage of amphiphiles, two distinct anchoring sequences, a continuous planar-tilted-homeotropic transition and a discontinuous planar-to-homeotropic transition, can be observed for the nematic and smectic LCs, respectively. More importantly, the latter occurs at a much lower surface coverage of amphiphiles, demonstrating an outstanding sensitivity for the smectic-based sensors. The dynamics of reorientation further reveals that the formation of homeotropic smectic anchoring is mainly governed by the synchronous growth of smectic layers through the LCs, which is significantly different from the mechanism of interface-to-bulk ordering propagation in nematic anchoring. Furthermore, the smectic LCs have also been proven to possess a potential selectivity in response to a subtle change in the chain rigidity of amphiphiles. These simulation findings are promising and would be valuable for the development of novel smectic-based sensors. [ABSTRACT FROM AUTHOR]

Magnetically optimization of mini-MR damper with focused on MR fluid properties and damper construction is investigated using dissipative particle dynamics method as a molecular scale modeling technique. To select a suitable MR fluid, the effect of diameter and weight of magnetic particles on damping force of 10 N as set point is studied. The results show damping force increases by increasing diameter of magnetic particles and finally trends to a constant value while changes nonparabolic by enhancement of their weight. The results of studies on structural parameters of damper show that by increasing gap size of flow passage, damping force increases while enhancement in inner diameter of cylinder and piston length has reverse effect. Optimum design of MR damper completed by investigation on electrical coils in terms of their arrangement and step wise distribution of magnetic field strength as our major innovation. To optimum operating conditions of MR damper at minimum electrical energy consumption and hysteresis level, the scenario of three segment coils in length of 3, 5, and 7 mm at base relative magnetic strength of 40% with 15% difference in steps by utilizing 140CG MR fluid as agent fluid are selected. [ABSTRACT FROM AUTHOR]