Search NASA⌕ Search

SEARCH · Search NASA

Results for “DPD”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Anisotropic single-particle dissipative particle dynamics model

We have developed a new single-particle dissipative particle dynamics (DPD) model for anisotropic particles with different shapes, e.g., prolate or oblate spheroids. In particular, the conservative and dissipative interactions between anisotropic single DPD particles are formulated using a linear mapping from the isotropic model of spherical particles. The proper mapping operator is constructed between each interacting pair of particles at every time step. Correspondingly, the random forces are properly formulated to satisfy the fluctuation-dissipation theorem (FDT). Notably, the model exactly conserves both linear and angular momentum. We demonstrate the proposed model's accuracy and efficiency by applying it for modeling colloidal ellipsoids. Specifically, we show it efficiently captures the static properties of suspensions of colloidal ellipsoids. The isotropic-nematic transition in an ellipsoidal suspension is reproduced by increasing its volume fraction or the aspect ratio of ellipsoid particles. Moreover, the hydrodynamics and diffusion of a single colloidal ellipsoid (prolate or oblate with moderate aspect ratios) are accurately captured. The calculated drag force on the ellipsoid and its diffusion coefficients (both translational and rotational) agree quantitatively with the theoretical predictions in the Stokes limit.

Computer Science↗
DOE OSTI2024

Modeling hyperbranched polymer formation via ATRP using dissipative particle dynamics

Hyperbranched polymers (HBPs) offer distinguishing, advantageous properties that arise from their distinctive complex topology. One of the effective approaches to the synthesis of hyperbranched structures involves the use of a branching initiator (inibramer) that is activated only after incorporation into a polymer chain. There remain, however, challenges in determining and characterizing the structures of the synthesized HBPs. Dissipative particle dynamics (DPD) was used to probe the effects of inibramer concentration, solvent concentration, and inibramer reactivity on the kinetics, molecular weight, and dispersity of HBPs. Additionally, DPD allows for direct observation of branched structures, which was not possible in previously reported Monte Carlo type simulations. It was found that higher inibramer concentrations led to faster monomer consumption while forming more dendritic structures with fewer defects. Additionally, high dispersities characteristic of HBP systems were found to originate from asymmetric propagation rates between inibramer-inibramer and inibramer-monomer reactions.

Biswas, Santidan↗
DOE OSTI2017

Dislocation-pipe diffusion in nitride superlattices observed in direct atomic resolution

Device failure from diffusion short circuits in microelectronic components occurs via thermally induced migration of atoms along high-diffusivity paths: dislocations, grain boundaries, and free surfaces. Even well-annealed single-grain metallic films contain dislocation densities of about 10 14 m -2 ; hence dislocation-pipe diffusion (DPD) becomes a major contribution at working temperatures. While its theoretical concept was established already in the 1950s and its contribution is commonly measured using indirect tracer, spectroscopy, or electrical methods, no direct observation of DPD at the atomic level has been reported. We present atomically-resolved electron microscopy images of the onset and progression of diffusion along threading dislocations in sequentially annealed nitride metal/semiconductor superlattices, and show that this type of diffusion can be independent of concentration gradients in the system but governed by the reduction of strain fields in the lattice.

42 ENGINEERING↗
DOE OSTI2011

Comparative studies on charged particle flow in a double plasma device

An experiment is conducted to study the comparative effect of a transverse magnetic field (TMF) and stainless steel (ss) plates in a filament discharge in a double plasma device (DPD). A variable transverse magnetic field is setup between the source and the target regions of the DPD using permanent magnets. By changing the distance between magnetic pole separations, the magnetic field is varied in the range of 80 Gauss to 300 Gauss. Plasma diffuses from the source to the target region through this field. Changing the discharge voltage and the discharge current varies the discharge conditions. Langmuir probes are placed in the source and target regions to measure the plasma parameters. Ss plates of similar surface area then replace the TMF and these data are compared with the TMF data to show the effect of the TMF on plasma parameters in comparison to the ss plates. It is seen that the magnetic field causes the plasma parameters to vary in a consistent way when the discharge parameters are varied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Conservative and dissipative force field for simulation of coarse-grained alkane molecules: A bottom-up approach

We apply operational procedures available in the literature to the construction of coarse-grained conservative and friction forces for use in dissipative particle dynamics (DPD) simulations. The full procedure rely on a bottom-up approach: large molecular dynamics trajectories of n-pentane and n-decane modeled with an anisotropic united atom model serve as input for the force field generation. As a consequence, the coarse-grained model is expected to reproduce at least semi-quantitatively structural and dynamical properties of the underlying atomistic model. Two different coarse-graining levels are studied, corresponding to five and ten carbon atoms per DPD bead. The influence of the coarse-graining level on the generated force fields contributions, namely, the conservative and the friction part, is discussed. It is shown that the coarse-grained model of n-pentane correctly reproduces self-diffusion and viscosity coefficients of real n-pentane, while the fully coarse-grained model for n-decane at ambient temperature over-predicts diffusion by a factor of 2. However, when the n-pentane coarse-grained model is used as a building block for larger molecule (e.g., n-decane as a two blobs model), a much better agreement with experimental data is obtained, suggesting that the force field constructed is transferable to large macro-molecular systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Improvement of charged particles transport across a transverse magnetic filter field by electrostatic trapping of magnetized electrons

A study on the transport of charged particles across a magnetic filter field has been carried out in a double plasma device (DPD) and presented in this manuscript. The DPD is virtually divided into two parts viz. source and target regions by a transverse magnetic field (TMF) which is constructed by inserting strontium ferrite magnets into two stainless steel rectangular tubes. Plasma electrons are magnetized but ions are unmagnetized inside the TMF region. Negative voltages are applied to the TMF tubes in order to reduce the loss of electrons towards them. Plasma is produced in the source region by filament discharge method and allowed to flow towards the target region through this negatively biased TMF. It is observed that in the target region, plasma density can be increased and electron temperature decreased with the help of negatively biased TMF. This observation is beneficial for negative ion source development. Plasma diffusion across the negatively biased TMF follows Bohm or anomalous diffusion process when negative bias voltage is very less. At higher negative bias, diffusion coefficient starts deviating from the Bohm diffusion value, associated with enhanced plasma flow in the target region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Dissipative particle dynamics simulations of deformation and aggregation of healthy and diseased red blood cells in a tube flow

In this paper, we report simulation results assessing the deformation and aggregation of mixed healthy and malaria-infected red blood cells (RBCs) in a tube flow. A three dimensional particle model based on Dissipative Particle Dynamics (DPD) is developed to predict the tube flow containing interacting cells. The cells are also modelled by DPD, with a Morse potential to characterize the cell-cell interaction. As validation tests, a single RBC in a tube flow and two RBCs in a static flow are simulated to examine the cell deformation and intercellular interaction, respectively. The study of two cells, one healthy and the other malaria-infected RBCs in a tube flow demonstrates that the malaria-infected RBC (in the leading position along flow direction) has different effects on the healthy RBC (in the trailing position) at the different stage of parasite development or at the different capillary number. With parasitic development, the malaria-infected RBC gradually loses its deformability, and in turn the corresponding trailing healthy RBC also deforms less due to the intercellular interaction. With increasing capillary number, both the healthy and malaria-infected RBCs are likely to undergo an axisymmetric motion. The minimum intercellular distance becomes small enough so that rouleaux is easily formed, i.e., the healthy and malaria-infected RBCs are difficultly disaggregated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Shear-induced aggregation or disaggregation in edible oils: Models, computer simulation, and USAXS measurements

The effects of shear upon the aggregation of solid objects formed from solid triacylglycerols (TAGs) immersed in liquid TAG oils were modeled using Dissipative Particle Dynamics (DPD) and the predictions compared to experimental data using Ultra-Small Angle X-ray Scattering (USAXS). The solid components were represented by spheres interacting via attractive van der Waals forces and short range repulsive forces. A velocity was applied to the liquid particles nearest to the boundary, and Lees-Edwards boundary conditions were used to transmit this motion to non-boundary layers via dissipative interactions. The shear was created through the dissipative forces acting between liquid particles. Translational diffusion was simulated, and the Stokes-Einstein equation was used to relate DPD length and time scales to SI units for comparison with USAXS results. The SI values depended on how large the spherical particles were (250 nm vs. 25 nm). Aggregation was studied by (a) computing the Structure Function and (b) quantifying the number of pairs of solid spheres formed. Solid aggregation was found to be enhanced by low shear rates. As the shear rate was increased, a transition shear region was manifested in which aggregation was inhibited and shear banding was observed. Aggregation was inhibited, and eventually eliminated, by further increases in the shear rate. The magnitude of the transition region shear, $\dot{γ}_t$, depended on the size of the solid particles, which was confirmed experimentally.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Impact of cross-linking of polymers on transport of salt and water in polyelectrolyte membranes: A mesoscopic simulation study

Our recent atomistic simulation studies demonstrated that the transport properties of salt ions and water in non-crosslinked polymer electrolyte membrane exhibit an intriguing dependence on salt concentration that is opposite to that seen in electrolyte solutions. Here in this paper, we extend our study to probe the influence of the degree of cross-linking of the polymer on the transport properties of salt and water in polymer electrolyte membranes. Towards this objective, we use a coarse-grained model embedded within dissipative particle dynamics (DPD) mesoscale simulations, which allows us to access time scales necessary for studying crosslinked polymer systems. Our DPD simulations on non-crosslinked membranes reproduce results that are in qualitative agreement with our atomistic simulations. For the case of crosslinked membranes, our results demonstrate that the diffusion of salt ions and water is reduced significantly relative to crosslinked systems. However, the trends exhibited by the salt concentration dependence of diffusivities and the coordination of the cations with anions and with the polymer backbone remain qualitatively similar to those observed in non-crosslinked membranes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A modified many-body dissipative particle dynamics model for mesoscopic fluid simulation: methodology, calibration, and application for hydrocarbon and water

The many-body dissipative particle dynamics (mDPD) is a prominent mesoscopic multiphase model for fluid transport in mesoconfinement. However, it has been a long-standing challenge for mDPD (and other multiphase-enabled DPD models) to accurately predict real-fluid static and dynamic properties simultaneously. We have developed a modified mDPD model to overcome the issue and a rigorous calibration approach that uses reference data, including experimental and/or molecular dynamics (MD) simulations to parameterise the modified mDPD for real fluids. We choose heptane as a representative example of hydrocarbon in source rocks to demonstrate the model's capability to accurately predict the equation of state (EOS), free surface tension, diffusivity, and viscosity. Our timing test shows that the modified mDPD is 400–500 times faster than its MD counterpart for simulating bulk heptane in equivalent volumes. To further demonstrate the robustness of the model, we revisited the benchmark problem of mesoscopic modelling of water, in which all the previous works on DPD reported only a limited portion of the water properties. Here, we show that the modified mDPD can provide accurate modelling of water static and dynamic properties and an EOS that matches the experimental data to a large range of confinement pressure.

42 ENGINEERING↗
DOE OSTI2016

MD/DPD Multiscale Framework for Predicting Morphology and Stresses of Red Blood Cells in Health and Disease

Healthy red blood cells (RBCs) have remarkable deformability, squeezing through narrow capillaries as small as 3 microns in diameter without any damage. However, in many hematological disorders the spectrin network and lipid bilayer of diseased RBCs may be significantly altered, leading to impaired functionality including loss of deformability. We employ a two-component whole-cell multiscale model to quantify the biomechanical characteristics of the healthy and diseased RBCs, including Plasmodium falciparum-infected RBCs (Pf-RBCs) and defective RBCs in hereditary disorders, such as spherocytosis and elliptocytosis. In particular, we develop a two-step multiscale framework based on coarse-grained molecular dynamics (CGMD) and dissipative particle dynamics (DPD) to predict the static and dynamic responses of RBCs subject to tensile forcing, using experimental information only on the structural defects in the lipid bilayer, cytoskeleton, and their interaction. We first employ CGMD on a small RBC patch to compute the shear modulus, bending stiffness, and network parameters, which are subsequently used as input to a whole-cell DPD model to predict the RBC shape and corresponding stress field. For Pf-RBCs at trophozoite and schizont stages, the presence of cytoadherent knobs elevates the shear response in the lipid bilayer and stiffens the RBC membrane. For RBCs in spherocytosis and elliptocytosis, the bilayer-cytoskeleton interaction is weakened, resulting in substantial increase of the tensile stress in the lipid bilayer. Furthermore, we investigate the transient behavior of stretching deformation and shape relaxation of the normal and defective RBCs. Different from the normal RBCs possessing high elasticity, our simulations reveal that the defective RBCs respond irreversibly, i.e., they lose their ability to recover the normal biconcave shape in successive loading cycles of stretching and relaxation. Our findings provide fundamental insights into the microstructure and biomechanics of RBCs, and demonstrate that the two-step multiscale framework presented here can be used effectively for in silico studies of hematological disorders based on first principles and patient-specific experimental input at the protein level.

Biochemistry & Molecular Biology↗
DOE OSTI2021

Using Dissipative Particle Dynamics to Model Effects of Chemical Reactions Occurring within Hydrogels

Computational models that reveal the structural response of polymer gels to changing, dissolved reactive chemical species would provide useful information about dynamically evolving environments. However, it remains challenging to devise one computational approach that can capture all the interconnected chemical events and responsive structural changes involved in this multi-stage, multi-component process. Here, we augment the dissipative particle dynamics (DPD) method to simulate the reaction of a gel with diffusing, dissolved chemicals to form kinetically stable complexes, which in turn cause concentration-dependent deformation of the gel. Using this model, we also examine how the addition of new chemical stimuli and subsequent reactions cause the gel to exhibit additional concentration-dependent structural changes. Through these DPD simulations, we show that the gel forms multiple latent states (not just the “on/off”) that indicate changes in the chemical composition of the fluidic environment. Hence, the gel can actuate a range of motion within the system, not just movements corresponding to the equilibrated swollen or collapsed states. Moreover, the system can be used as a sensor, since the structure of the layer effectively indicates the presence of chemical stimuli.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Analytical and Computational Studies of Correlations of Hydrodynamic Fluctuations in Shear Flow

Correlations of hydrodynamic fluctuations in shear flow analytically and also by dissipative particle dynamics (DPD) simulations are studied. The hydrodynamic equations are linearized around the macroscopic velocity field and then solved by a perturbation method in Fourier-transformed space. The autocorrelation functions (ACFs) from the analytical method are compared with results obtained from DPD simulations under the same shear-flow conditions. Up to a moderate shear rate, various ACFs from the two approaches agree with each other well. At large shear rates, discrepancies between the two methods are observed, hence revealing strong additional coupling between different fluctuating variables, which is not considered in the analytical approach. In addition, the results at low and moderate shear rates can serve as benchmarks for developing multiscale algorithms for coupling of heterogeneous solvers, such as a hybrid simulation of molecular dynamics and fluctuating hydrodynamics solver, where thermal fluctuations are indispensable.

Physics↗
NASA NTRS2022

Un-Numbered Operational Areas PRL 229 Railroad Tie Disposal Area Pond Confirmation Sampling Report

This Confirmatory Sampling (CS) Report (CSR) was prepared for the National Aeronautics and Space Administration (NASA), Kennedy Space Center (KSC), Florida in accordance with the CS Work Plan (CSWP) (HGL, 2022a). The CSR summarizes investigation activities performed on the Railroad Tie Disposal Area (RTDA) of the Un-Numbered Operational Areas (UNOA), Potential Release Location (PRL) 229 located at KSC. HydroGeoLogic, Inc. (HGL) prepared this report under Contract Number 80KSC019D0012. The purpose of the CSR is to provide information to evaluate the level of potential contamination in surface water and sediment samples collected from the RTDA PRL 229 Pond and determine if contaminants pose unacceptable risks to ecological receptors. The CSR provides analytical data of 10 sets of surface water and sediment samples collected in 2021. Surface water samples were analyzed for volatile organic compounds (VOCs), semi-volatile organic compounds (SVOCs), polynuclear aromatic hydrocarbon (PAHs), metals, hardness, and salinity, while sediment samples were analyzed for metals, total organic carbon (TOC), SVOCs, PAHs, and total petroleum hydrocarbons (TPH). The CSR includes a Screening Level Ecological Risk Assessment (SLERA) conducted under KSC’s Resource Conservation and Recovery Act (RCRA) Corrective Action Program to determine if unacceptable risks to ecological receptors exist at PRL 229 Pond. Metals and PAHs were detected in surface water samples, including detections of antimony, arsenic, barium, calcium, copper, lead, magnesium, thallium, zinc, anthracene, benzo(a)anthracene, benzo(a)pyrene, benzo(b)fluoranthene, benzo(g,h,i)perylene, benzo(k)fluoranthene, chrysene, dibenzo(a,h)anthracene, fluoranthene, indeno(1,2,3-c,d)pyrene and pyrene. Based on measured salinity (1.03 to 1.21 parts per thousand), the pond was identified as freshwater. Sediment analytical results also contained detections of metals and PAHs, including antimony, arsenic, barium, beryllium, cadmium, chromium, copper, lead, mercury, nickel, selenium, thallium, zinc, anthracene benzo(a)anthracene, benzo(a)pyrene, benzo(b)fluoranthene, benzo(g,h,i)perylene, benzo(k)fluoranthene, chrysene, dibenzo(a,h)anthracene, fluoranthene, indeno(1,2,3-c,d)pyrene and pyrene. Potential risks to ecological receptors were evaluated in accordance with the SLERA Work Plan and the Decision Process Document for the Resource Conservation and Recovery Act (RCRA) Corrective Action Program at KSC (DPD) (Geosyntec, 2019), and using the 2021 sediment and surface water analytical results. Potential ecological receptors evaluated in the SLERA include the aquatic community, benthic invertebrates, aquatic birds, and aquatic mammals. The SLERA included Steps 2a through 2c of the Ecological Risk Assessment (ERA) process in accordance with the DPD. Using multiple lines of evidence, the SLERA concluded that many of the detected analytes are present at concentrations similar to those in background or reference locations at KSC, and/or at concentrations that would not pose a significant risk to ecological communities. There are no screening values available for the avian community, and PAHs are bioaccumulative compounds. For this reason, food web modeling was completed as a supplement to Step 2c to provide a quantitative evaluation of this exposure pathway. Food web modeling indicates that PAHs are unlikely to pose a risk to birds exposed to site sediment and surface water. Based on the results of the SLERA, the contaminants in sediment and surface water pose no to minimal risk to the ecological communities that could live or forage at the site. Based on the results of the SLERA, no further study is recommended for the RTDA Pond at PRL 229.

Confirmation Sampling↗
DOE OSTI2012

Multiscale Modeling of Hematologic Disorders

Parasitic infectious diseases and other hereditary hematologic disorders are often associated with major changes in the shape and viscoelastic properties of red blood cells (RBCs). Such changes can disrupt blood flow and even brain perfusion, as in the case of cerebral malaria. Modeling of these hematologic disorders requires a seamless multiscale approach, where blood cells and blood flow in the entire arterial tree are represented accurately using physiologically consistent parameters. In this chapter, we present a computational methodology based on dissipative particle dynamics (DPD) which models RBCs as well as whole blood in health and disease. DPD is a Lagrangian method that can be derived from systematic coarse-graining of molecular dynamics but can scale efficiently up to small arteries and can also be used to model RBCs down to spectrin level. To this end, we present two complementary mathematical models for RBCs and describe a systematic procedure on extracting the relevant input parameters from optical tweezers and microfluidic experiments for single RBCs. We then use these validated RBC models to predict the behavior of whole healthy blood and compare with experimental results. The same procedure is applied to modeling malaria, and results for infected single RBCs and whole blood are presented.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2011

Time-dependent and outflow boundary conditions for Dissipative Particle Dynamics

We propose a simple method to impose both no-slip boundary conditions at fluid-wall interfaces and at outflow boundaries in fully developed regions for Dissipative Particle Dynamics (DPD) fluid systems. The procedure to enforce the no-slip condition is based on a velocity-dependent shear force, which is a generalized force to represent the presence of the solid-wall particles and to maintain locally thermodynamic consistency. We show that this method can be implemented in both steady and time-dependent fluid systems and compare the DPD results with the continuum limit (Navier-Stokes) results. We also develop a force-adaptive method to impose the outflow boundary conditions for fully developed flow with unspecified outflow velocity profile or pressure value. We study flows over the backward-facing step and in idealized arterial bifurcations using a combination of the two new boundary methods with different flow rates. Finally, we explore the applicability of the outflow method in time-dependent flow systems. The outflow boundary method works well for systems with Womersley number of O(1), i.e. when the pressure and flowrate at the outflow are approximately in-phase.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Energy-conserving dissipative particle dynamics with temperature-dependent properties

The dynamic properties of fluid, including diffusivity and viscosity, are temperature-dependent and can significantly influence the flow dynamics of mesoscopic non-isothermal systems. To capture the correct temperature-dependence of a fluid, an energy-conserving dissipative particle dynamics (eDPD) model is developed by expressing the weighting terms of the dissipative force and the random force as functions of temperature. The diffusivity and viscosity of liquid water at various temperatures ranging from 273 K to 373 K are used as examples for verifying the proposed model. Simulations of a Poiseuille flow and a steady case of heat conduction for reproducing the Fourier law are carried out to validate the present eDPD formulation and the thermal boundary conditions. Results show that the present eDPD model recovers the standard DPD model when isothermal fluid systems are considered. For non-isothermal fluid systems, the present model can predict the diffusivity and viscosity consistent with available experimental data of liquid water at various temperatures. Moreover, an analytical formula for determining the mesoscopic heat friction is proposed. The validity of the formula is confirmed by reproducing the experimental data for Prandtl number of liquid water at various temperatures. The proposed method is demonstrated in water but it can be readily extended to other liquids. (C) 2014 Elsevier Inc. All rights reserved.

CONSERVATION↗