Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dissipative particle dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

Correction to “Frank–Kasper Phases of Diblock Copolymer Melts Studied with the DPD Model: SCF Results”

In a recent study, (1) we performed the self-consistent field (SCF) calculations of a dissipative particle dynamics (DPD) model of conformationally asymmetric diblock copolymer A-B melts to study the stability of various Frank–Kasper (FK) phases formed by such A-B melts. While all the equations in ref (1) are correct, we did not calculate the dimensionless (mean-field) Helmholtz free energies per chain βf c = βu c,χ + βu c,κ – s c /k B and the dimensionless stresses dβf c /dθ of ordered phases accurately. For example, the dimensionless (nonbonded) internal energy per chain due to the system compressibility βu c,κ = (N/2κ)∑ q βû 0 (|q|)(ϕ̂ A (q) + ϕ̂ B (q) – δ q,0 )(ϕ̂ A (−q) + ϕ̂ B (−q) – δ– q,0 ) should be calculated in the reciprocal space as given (note that βû 0 (|q|) is obtained analytically, and that since ϕ A (r) is real 𝜙̂ A (−𝐪)= $\overline{𝜙̂ A}$(𝐪), i.e., the complex conjugate of ϕ̂ A (q)); this corresponds to numerically evaluating βu c,κ = (N/2κ)∫dr(ϕ A (r) + ϕ B (r) – 1)ψ(r) with ψ(r) ≡ (∫dr′βu 0 (|r–r′|)(ϕ A (r′) + ϕ B (r′) – 1)/V in the real space via the (composite) trapezoidal rule. While the trapezoidal rule is less accurate than the Romberg integration (2) in general, the opposite occurs when the integration domain is periodic; this was pointed out by Matsen in ref (3) and explained in detail in refs (4and5). With our ignorance of this exceptional convergence property of the trapezoidal rule for periodic functions, we calculated all convolutions having the form of ψ(r) in the reciprocal space via the fast Fourier transforms (6) but those having the form of βu c,κ in the real space via Romberg integration; the latter leads to inaccurate results and unfortunately incorrect phase diagrams reported in ref (1), which are corrected here.

He, Juntong↗

Mesoscale Modeling of Hydrogels Under Frictional Shear Stress

Hydrogels are three-dimensional networks of hydrophilic polymers often used as a simplified model of hydrated biological materials, from cartilaginous joints to the ocular tear film. However, the lubrication mechanisms of hydrogels remain poorly understood, partly due to their complex polymeric structure, which creates blurred interfaces during sliding that are challenging to study experimentally. In this study, we employ dissipative particle dynamics (DPD) to investigate the frictional behavior of a polymeric hydrogel network sliding against a solid wall in an explicit viscous solvent. This computational approach enables us to model hydrodynamic interactions and mesoscale polymer dynamics, capturing key aspects of hydrogel friction. Our simulations reveal that hydrogel friction is governed by the interplay between polymer relaxation and viscous shear, characterized by the Weissenberg number (Wi). At low Wi, friction coefficient remain nearly constant, dominated by polymer relaxation. However, at higher Wi, friction is dominated by viscous drag within a near-wall solvent layer, leading to a linear increase in friction coefficient with Wi. Furthermore, our results demonstrate an inverse relationship between the friction coefficient and the applied normal load, consistent with experimental observations. This work provides new insights into the fundamental tribological properties of hydrogels, shedding light on the micromechanics of hydrogel friction. Improving our understanding of hydrogel structure and dynamics under friction advances our knowledge of the mechanisms regulating biological lubrication in health and disease.

36 MATERIALS SCIENCE↗

Kinetic Pathways to Diblock Copolymer Micelles via Thin Film Dissolution and Cosolvent Variation

We investigate how processing pathways control micelle formation during dissolution of lamellar diblock copolymers using dissipative particle dynamics (DPD) simulations and DPD-selfconsistent field theory (DPD-SCFT). Comparing thin film (TF) and cosolvent-assisted (CS) dissolution reveals fundamentally different micellization mechanisms and final states. In the TF pathway, dissolution proceeds through a metastable cylindrical intermediate that dictates final micelle dimensions. The lamellar melt first transforms into cylindrical micelles, which then fragment into spherical micelles. The final spherical radius is intrinsically linked to the cylinder radius by a constant ratio, Rsphere ≈ 1.4Rcylinder. DPD-SCFT confirms that both cylinders and spheres are near their equilibrium dimensions. In contrast, CS dissolution produces kinetically trapped micelles with sizes significantly below equilibrium. Micellization is hindered by high free energy barriers to single-chain exchange and fusion. Growth occurs primarily via fusion, but increasing corona repulsion suppresses further coarsening, preventing equilibration. These results demonstrate that, when starting from a bulk copolymer phase, the metastable cylindrical intermediate plays a central role in determining final micelle size.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Controlled Organization of Inorganic Materials Using Biological Molecules for Activating Therapeutic Functionalities

Precise control over the assembly of biocompatible three-dimensional (3D) nanostructures would allow for programmed interactions within the cellular environment. Nucleic acids can be used as programmable crosslinkers to direct the assembly of quantum dots (QDs) and tuned to demonstrate different interparticle binding strategies. In this work, morphologies of self-assembled QDs are evaluated via gel electrophoresis, transmission electron microscopy, small-angle X-ray scattering, and dissipative particle dynamics simulations, with all results being in good agreement. The controlled assembly of 3D QD organizations is demonstrated in cells via the colocalized emission of multiple assembled QDs, and their immunorecognition is assessed via enzyme-linked immunosorbent assays. RNA interference inducers are also embedded into the interparticle binding strategy to be released in human cells only upon QD assembly, which is demonstrated by specific gene silencing. The programmability and intracellular activity of QD assemblies offer a strategy for nucleic acids to imbue the structure and therapeutic function into the formation of complex networks of nanostructures, while the photoluminescent properties of the material allow for optical tracking in cells in vitro.

36 MATERIALS SCIENCE↗

Multi-scale computer-aided design and photo-controlled macromolecular synthesis boosting uranium harvesting from seawater

By integrating multi-scale computational simulation with photo-regulated macromolecular synthesis, this study presents a new paradigm for smart design while customizing polymeric adsorbents for uranium harvesting from seawater. A dissipative particle dynamics (DPD) approach, combined with a molecular dynamics (MD) study, is performed to simulate the conformational dynamics and adsorption process of a model uranium grabber, i.e., PAO m -b-PPEGMA n , suggesting that the maximum adsorption capacity with atomic economy can be achieved with a preferred block ratio of 0.18. The designed polymers are synthesized using the PET-RAFT polymerization in a microfluidic platform, exhibiting a record high adsorption capacity of uranium (11.4 ± 1.2 mg/g) in real seawater within 28 days. This study offers an integrated perspective to quantitatively assess adsorption phenomena of polymers, bridging metal-ligand interactions at the molecular level with their spatial conformations at the mesoscopic level. The established protocol is generally adaptable for target-oriented development of more advanced polymers for broadened applications.

97 MATHEMATICS AND COMPUTING↗

Onset of cavitation and vapor bubble development over hydrophilic and hydrophobic surfaces

Cavitation, the formation of vapor bubbles as the liquid pressure is reduced below the saturated vapor pressure, often requires a substantial negative relative pressure in a pure liquid. Classical nucleation theory (CNT) provides an estimate for the rate of cavitation but there is often a disconnect between the predictions at the molecular scale compared to observations at the macroscale. We report on mesoscale simulations of cavitation based on many-body dissipative particle dynamics (mDPD), a coarse-grained molecular dynamics (MD), which bridges the two scales. A liquid layer is confined between smooth planar walls at a constant temperature, while the pressure is reduced slowly by expanding the wall-bounded domain. The wetting properties of the liquid are determined by the parameters of the interaction potentials. With hydrophilic walls, homogeneous nucleation is observed in the liquid bulk. As a bubble forms and grows, it creates a strong pressure pulse and oscillations that cause other bubbles that may have formed slightly later to collapse. For a nearly neutral wall with a contact angle close to 90 ° , heterogeneous nucleation occurs at the walls at a smaller negative pressure and generates weaker pressure oscillations. With hydrophobic walls or seed particles, heterogeneous nucleation readily occurs, where fluctuations and the merger of transient surface bubbles are significant.

Science & Technology - Other Topics↗

Droplet bioprinting of acellular and cell-laden structures at high-resolutions

Advances in digital light projection(DLP) based (bio) printers have made printing of intricate structures at high resolution possible using a wide range of photosensitive bioinks. A typical setup of a DLP bioprinter includes a vat or reservoir filled with liquid bioink, which presents challenges in terms of cost associated with bioink synthesis, high waste, and gravity-induced cell settling, contaminations, or variation in bioink viscosity during the printing process. Here, we report a vat-free, low-volume, waste-free droplet bioprinting method capable of rapidly printing 3D soft structures at high resolution using model bioinks and model cells. A multiphase many-body dissipative particle dynamics model was developed to simulate the dynamic process of droplet-based DLP printing and elucidate the roles of surface wettability and bioink viscosity. Process variables such as light intensity, photo-initiator concentration, and bioink formulations were optimized to print 3D soft structures (∼0.4–3 kPa) with a typical layer thickness of 50 µm, an XY resolution of 38 ± 1.5 μm and Z resolution of 237 ± 5.4 µm. To demonstrate its versatility, droplet bioprinting was used to print a range of acellular 3D structures such as a lattice cube, a Mayan pyramid, a heart-shaped structure, and a microfluidic chip with endothelialized channels. Droplet bioprinting, performed using model C3H/10T1/2 cells, exhibited high viability (90%) and cell spreading. Additionally, microfluidic devices with internal channel networks lined with endothelial cells showed robust monolayer formation while osteoblast-laden constructs showed mineral deposition upon osteogenic induction. Overall, droplet bioprinting could be a low-cost, no-waste, easy-to-use, method to make customized bioprinted constructs for a range of biomedical applications.

DLP↗

Producing entangled photon pairs and quantum squeezed states in plasmas

Plasma is capable of mediating the conversion of two pump photons into two different photons through a relativistic four-wave mixing nonlinearity. Spontaneously created photon pairs are emitted at symmetric angles with respect to the colinear pump direction, and the emission rate is largest if they have identical frequency. Thus, two orthogonally polarized pumps can produce polarization-entangled photon pairs through a millimeter-long homogeneous plasma. Here, the noise from Raman scattering can be avoided if the pump detuning differs from twice the plasma frequency. However, pump detuning exactly equal to twice the plasma frequency can significantly enhance the interaction rate, which allows for the production of strong two-mode squeezed states. Remarkably, the amplified noise from Raman scattering are correlated and hence can be suppressed in one of the output quadratures, thereby maintaining the squeezing magnitude.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A versatile approach to stabilize liquid-liquid interfaces using surfactant self-assembly

Stabilizing liquid–liquid interfaces, whether between miscible or immiscible liquids, is crucial for a wide range of applications, including energy storage, microreactors, and biomimetic structures. In this study, a versatile approach for stabilizing the water-oil interface is presented using the morphological transitions that occur during the self-assembly of anionic, cationic, and nonionic surfactants mixed with fatty acid oils. The morphological transitions underlying this approach are characterized and extensively studied through small-angle X-ray scattering (SAXS), rheometry, and microscopy techniques. Dissipative particle dynamics (DPD) as a simulation tool is adopted to investigate these morphological transitions both in the equilibrium ternary system as well as in the dynamic condition of the water-oil interface. Such a versatile strategy holds promise for enhancing applications such as liquid-in-liquid 3D printing. Moreover, it has the potential to revolutionize a wide range of fields where stabilizing liquid–liquid interfaces not only offers unprecedented opportunities for fine-tuning nanostructural morphologies but also imparts interesting practical features to the resulting liquid shapes. These features include perfusion capabilities, self-healing, and porosity, which could have significant implications for various industries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Frank–Kasper Phases of Diblock Copolymer Melts: Self-Consistent Field Results of Two Commonly Used Models

We constructed phase diagrams of conformationally asymmetric diblock copolymer A-B melts using the polymer self-consistent field (SCF) calculations of both the dissipative particle dynamics chain (DPDC) model (i.e., compressible melts of discrete Gaussian chains with the DPD non-bonded potential) and the “standard” model (i.e., incompressible melts of continuous Gaussian chains with the Dirac δ-function non-bonded potential) in the χN-ε plane, where χN and ε characterize, respectively, the repulsion and conformational asymmetry between the A and B blocks, at the A-block volume fraction f = 0.2 and 0.3. Consistent with previous SCF calculations of the “standard” model, σ and A15 are the only stable Frank–Kasper (FK) phases among the five FK (i.e., σ, A15, C14, C15 and Z) phases considered. The stability of σ and A15 is due to their delicate balance between the energetic and entropic contributions to the Helmholtz free energy per chain of the system, which, within our parameter range, increases in the order of σ/A15, Z, and C14/C15. While in general the SCF phase diagrams of these two models are qualitatively consistent, A15 is not stable for the DPDC model at the copolymer chain length N = 10 and f = 0.3; any differences in the SCF phase diagrams are solely due to the differences between these two models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

“Transition from dissipative fluid dynamics to particles in heavy ion collisions”(Final Technical Report)

This report summarizes research activities and achievements made under Award No. DE-SC00016524. We advanced the following areas: i) self-consistent bulk viscous and baryon number diffusion corrections from kinetic theory, ii) hadron gas bulk viscosity determination, and iii) discovery of novel mechanisms of elliptic flow generation and suppression in heavy-ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quantum system dynamics with a weakly nonlinear Josephson junction bath

Here, we investigate the influence of a weakly nonlinear Josephson bath consisting of a chain of Josephson junctions on the dynamics of a small quantum system (LC oscillator). Focusing on the regime where the charging energy is the largest energy scale, we perturbatively calculate the correlation function of the Josephson bath to the leading order in the Josephson energy divided by the charging energy while keeping the cosine potential exactly. When the variation of the charging energy along the chain ensures fast decay of the bath correlation function, the dynamics of the LC oscillator that is weakly and capacitively coupled to the Josephson bath can be solved through the Markovian master equation. We establish a duality relation for the Josephson bath between the regimes of large charging and Josephson energies, respectively. The results can be applied to cases where the charging energy either is nonuniformly engineered or disordered in the chain. Furthermore, we find that the Josephson bath may become non-Markovian when the temperature is increased beyond the zero-temperature limit in that the bath correlation function gets shifted by a constant and does not decay with time.

1-dimensional systems↗

Collisional dynamics of perturbed particle disks in the solar system

Theoretical and computational studies were carried out on galactic and planetary disks. With the goal of addressing important open questions centered on galactic structure, the cloudy interstellar medium, giant molecular clouds, and star formation in galactic disks and the collisional dynamics of perturbed particulate matter in planetary disks, focus was largely on the self-gravitational effects, dissipative effects, and collisional dynamics of cloud-particle disks. N-body, 'cloud-particle' computational algorithms were developed for the purpose of isolating the role of gaseous self gravity from the roles of other dominant physical mechanisms and dynamical processes, e.g. the collisional dynamics and dissipative processes. The efforts focused largely on galactic disks show that observational constraints provide stringent tests of the numerical simulation techniques developed. Self gravitational effects of the galactic interstellar medium's gas clouds are included by means of Fourier Transform technique.

Roberts, William W., Jr.↗

Collaborative Research: Explosive reconnection in relativistic magnetically-dominated plasmas (Final Report)

Magnetic fields can be the dominant component of astrophysical plasmas, so that the magnetic energy density might exceed even the rest-mass energy density of matter. In this extreme (and largely unexplored) regime the magnetic field controls the overall plasma evolution, dissipation, and acceleration of non-thermal particles. This plasma regime, applicable to a variety of astrophysical sources - magnetars, pulsars and pulsar wind nebulae (PWNe), jets of Active Galactic Nuclei (AGNs) and Gamma-Ray Bursters (GRBs) - is dramatically different from laboratory plasmas, the magnetospheres of planets, and the interplanetary plasma. Relativistic astrophysical sources then provide an unique opportunity to study the fundamental plasma physics of magnetically-dominated plasmas; a novel and fast-evolving field of theoretical research which, by investigating energy conversion and particle energization processes in plasmas, is of vital importance to the Fusion Energy Sciences DoE program. Data coming from astrophysical high-energy missions, especially the Crab Nebula flares recently observed by the Fermi and AGILE satellites, suggest that the acceleration of non-thermal particles to the highest energies occurs in magnetic reconnection events - a major change of paradigm in high-energy astrophysics. Most importantly, observations demand that particle acceleration should proceed extremely fast (with accelerating electric field of the order of the magnetic field) and on macroscopic scales (much larger, e.g., than the microscopic plasma skin depth). We are conducting studies of the microphysics of magnetically-dominated plasmas focussing in particular on the highly dynamic regime of explosive reconnection and associated particle acceleration in relativistic plasmas. We are studying the stability and explosive plasma dynamics, particle acceleration and radiation production in a number of idealized plasma configurations that approximate relevant astrophysical sources (like the magnetic ABC structures and interacting flux tubes, as well as generalizations of analytical models of X-point collapse to relativistic plasmas). The well-studied case of the Crab Nebula is taken as a prototypical example for the application of the model. We are combining analytical studies of explosive magnetic dynamics and dissipation in relativistic plasmas with particle-in-cell (PIC) simulations and fluid simulations. The theoretical model, fluid and particle-in-cell simulations are cross-checked for agreement and convergence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗