Search NASA⌕ Search

SEARCH · Search NASA

Results for “INTERACTION RANGE”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Application of the GBFP Method to Electron Partial-Wave Expansion Elastic Scattering Differential Cross Sections Within the Geant4 Toolkit - Paper 69

Analog Monte Carlo simulation of electrons in complex material geometries is generally impractical because the long range Coulomb interactions result in highly peaked differential cross sections (DCS) for scattering and energy-loss and extremely small collision mean free paths (mfp). The condensed history (CH) Monte Carlo method attempts to alleviate this limitation by advancing electrons in fixed steps and sampling angular deflections and energy losses from precomputed distributions obtained from approximate solutions to the underlying transport equation that are valid only for thin layers. Though CH is efficient and is widely employed in production codes (ITS, MCNP, EGS, Geant), step size restrictions and inconsistent handling of material and free surface boundaries limit the ultimate accuracy possible with this approach, especially in applications to highly heterogeneous media. Recently proposed alternatives to the CH method, such as the Generalized Boltzmann Fokker-Planck (GBFP) method, restore the correct transport mechanics by simulating single events based on a Boltzmann transport equation but with efficiency and accuracy achieved through the use of reduced-physics interaction models. Motivated by Lewis theory, various moment-preserving strategies are employed in the GBFP approach to construct DCSs that are less singular, and with mfps shorter than analog. Specifically, the approximate cross sections are constructed to rigorously preserve a finite number of angular and energy-loss moments of the analog DCSs while all higher moments are approximated in terms of the exact lower moments. By systematically increasing the number of moments retained, any desired accuracy can be obtained for integrated quantities, such as dose, reflection and transmission yields, as well as distributed quantities such as energy and angular distributions. For relativistic Rutherford cross section models, the GBFP method has been demonstrated in Monte Carlo simulation to yield accuracy comparable to analog but at a fraction of the cost. One advantage of the GBFP approach is that angular deflection and energy-loss interactions depend only on Legendre moments and energy-loss moments and not the detailed form of the analog DCSs. This feature allows moment-preservation to be enforced through DCSs expressed as a superposition of discrete (i.e., using delta-functions) angular deflections and discrete energy losses. This promotes extremely efficient computation when compared to analog and completely continuous or hybrid discrete-continuous moment-based DCS reconstructions, while yielding highly accurate results for integrated and distributed quantities. In this paper, we demonstrate that the moment-preserving algorithms central to our fully discrete GBFP method can be readily extended to apply to tabulated DCS data, specifically to the partial-wave expansion (PWE) elastic scattering DCSs. Since data libraries for the PWE DCSs for energy and material combinations of interest can be prohibitively large, the GBFP method provides a means to reduce data storage requirements. We implemented our algorithm in the Geant4 Monte Carlo code and all numerical results are obtained using this implementation. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

Design and simulation of a gyroklystron amplifier

In the present paper, a design methodology of the gyroklystron amplifier has been described and subsequently used for the design of a typically selected 200 kW, Ka-band, four-cavity gyroklystron amplifier. This conceptual device design has been validated through the 3D particle-in-cell (PIC) simulation and nonlinear analysis. Commercially available PIC simulation code “MAGIC” has been used for the electromagnetic study at the different location of the device RF interaction structure for the beam-absent case, i.e., eigenmode study as well as for the electron beam and RF wave interaction behaviour study in the beam present case of the gyroklystron. In addition, a practical problem of misalignment of the RF cavities with drift tubes within the tube has been also investigated and its effect on device performance studied. The analytical and simulation results confirmed the validity of the gyroklystron device design. The PIC simulation results of the present gyroklystron produced a stable RF output power of ∼218 kW for 0% velocity spread at 35 GHz, with ∼45 dB gain, 37% efficiency, and a bandwidth of 0.3% for a 70 kV, 8.2 A gyrating electron beam. The simulated values of RF output power have been found in agreement with the nonlinear analysis results within ∼5%. Further, the PIC simulation has been extended to study a practical problem of misalignment of the cavities axis and drift tube axis of the gyroklystron amplifier and found that the RF output power is more sensitive to misalignments in comparison to the device bandwidth. The present paper, gyroklystron device design, nonlinear analysis, and 3D PIC simulation using commercially available code had been systematically described would be of use to the high-power gyro-amplifier tube designers and research scientists.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

RADIO CONSTRAINTS ON LONG-LIVED MAGNETAR REMNANTS IN SHORT GAMMA-RAY BURSTS

The merger of a neutron star (NS) binary may result in the formation of a rapidly spinning magnetar. The magnetar can potentially survive for seconds or longer as a supramassive NS before collapsing to a black hole if, indeed, it collapses at all. During this process, a fraction of the magnetar’s rotational energy of ∼10{sup 53} erg is transferred via magnetic spin-down to the surrounding ejecta. The resulting interaction between the ejecta and the surrounding circumburst medium powers a year-long or greater synchrotron radio transient. We present a search for radio emission with the Very Large Array following nine short-duration gamma-ray bursts (GRBs) at rest-frame times of ≈1.3–7.6 yr after the bursts, focusing on those events that exhibit early-time excess X-ray emission that may signify the presence of magnetars. We place upper limits of ≲18–32 μ Jy on the 6.0 GHz radio emission, corresponding to spectral luminosities of ≲(0.05–8.3) × 10{sup 39} erg s{sup -1}. Comparing these limits to the predicted radio emission from a long-lived remnant and incorporating measurements of the circumburst densities from broadband modeling of short GRB afterglows, we rule out a stable magnetar with an energy of 10{sup 53} erg for half of the events in our sample. A supramassive remnant that injects a lower rotational energy of 10{sup 52} erg is ruled out for a single event, GRB 050724A. This study represents the deepest and most extensive search for long-term radio emission following short GRBs to date, and thus the most stringent limits placed on the physical properties of magnetars associated with short GRBs from radio observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

2D to 3D transition of polymeric carbon nitride nanosheets

The transition from a prevalent turbostratic arrangement with low planar interactions (2D) to an array of polymeric carbon nitride nanosheets with stronger interplanar interactions (3D), occurring for samples treated above 650 °C, was detected by terahertz-time domain spectroscopy (THz-TDS). The simulated 3D material made of stacks of shifted quasi planar sheets composed of zigzagged polymer ribbons, delivered a XRD simulated pattern in relatively good agreement with the experimental one. The 2D to 3D transition was also supported by the simulation of THz-TDS spectra obtained from quantum chemistry calculations, in which the same broad bands around 2 THz and 1.5 THz were found for 2D and 3D arrays, respectively. This transition was also in accordance with the tightening of the interplanar distance probably due to an interplanar π bond contribution, as evidenced also by a broad absorption around 2.6 eV in the UV–vis spectrum, which appeared in the sample treated at 650 °C, and increased in the sample treated at 700 °C. The band gap was calculated for 1D and 2D cases. The value of 3.374 eV for the 2D case is, within the model accuracy and precision, in a relative good agreement with the value of 3.055 eV obtained from the experimental results. - Graphical abstract: 2D lattice mode vibrations and structural changes correlated with the so called “2D to 3D transition”. - Highlights: • A 2D to 3D transition has been detected for polymeric carbon nitride. • THz-TDS allowed us to discover and detect the 2D to 3D transition of polymeric carbon nitride. • We propose a structure for polymeric carbon nitride confirming it with THz-TDS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Design and Beam Test Results for the sPHENIX Electromagnetic and Hadronic Calorimeter Prototypes

The super Pioneering High Energy Nuclear Interaction eXperiment (sPHENIX) at the Relativistic Heavy Ion Collider (RHIC) will perform high precision measurements of jets and heavy flavor observables for a wide selection of nuclear collision systems, elucidating the microscopic nature of strongly interacting matter ranging from nucleons to the strongly coupled quark-gluon plasma. A prototype of the sPHENIX calorimeter system was tested at the Fermilab Test Beam Facility as experiment T-1044 in the spring of 2016. The electromagnetic calorimeter (EMCal) prototype is composed of scintillating fibers embedded in a mixture of tungsten powder and epoxy. The hadronic calorimeter (HCal) prototype is composed of tilted steel plates alternating with plastic scintillator. Results of the test beam reveal the energy resolution for electrons in the EMCal is $2.8\%\oplus~15.5\%/\sqrt{E}$ and the energy resolution for hadrons in the combined EMCal plus HCal system is $13.5\%\oplus 64.9\%/\sqrt{E}$. These results demonstrate that the performance of the proposed calorimeter system satisfies the sPHENIX specifications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Collaborative Research Proposal: Time Resolved Optical Emissions Spectroscopy and Laser Induced Fluorescence Spectroscopy of Nanosecond Pulsed Discharges in a Gas-Liquid Water Film Reactor

Electrical discharge plasma formed in contact with liquid water is of interest for a wide range of applications in chemical, biomedical, agricultural, electrical, and materials science and engineering. Such plasma reactors are of very timely importance since they also have significant disinfection capability by inactivating bacteria, viruses and other pathogens. Many types of plasma sources including those driven by AC, DC, RF, microwave, and pulsed electrical power supplies coupled to a wide range of different electrode configurations and reactor designs have been developed and explored which contact the plasma with liquid water. Recent roadmaps have recommended that further work is needed to develop our understanding of the fundamental chemical and physical process which occur at the interface of non-thermal plasma with liquid water solutions in order to advance this large diversity of applications which ultimately depend upon efficient production of key reactive chemical species. There is a wide range of interacting factors that affect the chemical reactions that occur in the plasma, in the liquid phase, and at the interface. These factors ultimately govern the key reactive chemical species formed and used in the various applications and they include a) the reactor design and input parameters, b) the discharge and transport processes, c) the plasma properties, and d) the resulting chemical reactions. For pulsed discharges, the power supply design and output parameters control the applied voltage, frequency, rise time, and width (duration) of the applied pulses. The reactor design involves specification of the gas-liquid contacting methods, electrode gap distance, reactor volume and shape, and gas and liquid flow rates achievable. The gas and liquid compositions as well as the liquid properties such as pH and conductivity are also of key importance in determination of the resulting chemical reactions. In addition, the important plasma properties include plasma gas temperature, electron density, electron energy (distribution), and size of the plasma channels which are all affected by the discharge and transport properties. Many studies have focused on specific aspects of these various processes. Of particular important and relevance to the proposed work is the utilization of nanosecond pulses to generate plasma in gas-liquid systems, liquid bubbles, underwater, and in gases. Recent advances in nanosecond pulses provide significant advantages in utilization with liquid water. For example, fast rise time and short pulses are less sensitive to water conductivity, such short pulses may provide advantages in fast temporal quenching of the plasma, and fundamental analysis of pulse properties, including pulse shape and width, may be facilitated by investigation of single filamentary fast pulses. Many studies have also dealt with the role of the gas composition on formation of reactive oxygen species (ROS) (i.e., hydroxyl radicals – ·OH, hydrogen peroxide – H 2 O 2 , various atomic oxygen species, ozone-O 3 , hydroperoxyl radicals HO 2 ·) and reactive nitrogen species (RNS) (i.e., nitrogen oxides – NO, NO 2 , N 2 O - collectively termed NO X , nitrite-NO 2 - , nitrate-NO 3 - , peroxynitrite-ONOO - ). The hydroxyl radical is the critical species in many chemical oxidation reactions for chemical degradation of toxic compounds in water and gases and for synthesis of some compounds. The mixture of various nitrogen oxide species is important for many biochemical and biological processes involved in biomedical and agricultural applications including disinfection and fertilizer production. The present proposal focuses on determination of the effects of time resolved electron density and hydroxyl radicals on plasma chemical reactions through collaboration with the Princeton Collaborative Low Temperature Plasma Research Facility (PCRF) at the Princeton Plasma Physics Laboratory (PPPL). In order to further investigate the role of the plasma generated electrons and hydroxyl radicals on the overall formation of the key species including hydrogen peroxide, hydroxyl radicals, and nitrogen oxides, the proposed work, thus seeks to determine high resolution time resolved electron density by optical emissions spectroscopy and time resolved hydroxyl radicals using laser induced fluorescence in the nanosecond discharge reactor. The combination of data on electron density and hydroxyl radical concentration will be utilized in the present work to more fully characterize the chemical reaction processes in this system and to advance the design, development, and operation of such chemical reactors for a wide range of applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Non-Invasive Biophysical Techniques to Monitor the Structural Plasticity of the Photosynthetic Machinery of Live Diatom Cells

The photosynthetic performance of diatoms depends largely on the organization and structural flexibility of their thylakoid membranes, the densely packed, highly organized membrane vesicles in which light reactions of photosynthesis occur. Different regulatory mechanisms that fine tune the photosynthetic functions affect the organization of the photosynthetic machinery at different levels of structural complexity, from the level of individual protein complexes to the macroarray of membrane proteins and the remodeling of the entire thylakoid membrane system. To monitor these reorganizations, non-invasive techniques are of special value. In this chapter, we focus our attention on three of these techniques, which have been demonstrated to provide unique and useful information on the structure and structural and functional plasticity of live diatom cells: (i) circular dichroism (CD) spectroscopy, which has provided unique information on the chiral (macro-)organization of protein complexes and on their rapid, reversible reorganizations, fine-tuning the light-harvesting processes, as well as on variations in the short-range excitonic interactions in the antenna complexes; (ii) small-angle neutron scattering (SANS), which has been used to determine the periodic organization of the thylakoid membranes and to monitor reversible ultrastructural changes on the time-scale of minutes, induced by variations in the environmental conditions such as changes in temperature or light intensity; and (iii) electrochromic shift absorbance transients (ΔA ECS ), a spectroscopic tool which has been shown to be capable of identifying distinct functional groups of the light-harvesting carotenoid fucoxanthin in different diatoms and in cells exposed to different light intensities. Future use of these techniques will most certainly contribute to the deeper understanding of key regulatory mechanisms of photosynthesis in diatoms.

Szabo, Milan↗
DOE OSTI2017

Branched Aramid Nanofibers

Interconnectivity of components in three-dimensional networks (3DNs) is essential for stress transfer in hydrogels, aerogels, and composites. Entanglement of nanoscale components in the network relies on weak short-range intermolecular interactions. The intrinsic stiffness and rod-like geometry of nanoscale components limit the cohesive energy of the physical crosslinks in 3DN materials. Nature realizes networked gels differently using components with extensive branching. In this work, branched aramid nanofibers (BANFs) mimicking polymeric components of biological gels were synthesized to produce 3DNs with high efficiency stress transfer. Individual BANFs are flexible, with the number of branches controlled by base strength in the hydrolysis process. The extensive connectivity of the BANFs allows them to form hydro- and aerogel monoliths with an order of magnitude less solid content than rod-like nanocomponents. Branching of nanofibers also leads to improved mechanics of gels and nanocomposites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Time-Dependent Density Functional Theory for Fermionic Superfluids: From Cold Atomic Gases - To Nuclei and Neutron Stars Crust

We report in cold atoms and in the crust of neutron stars, the pairing gap can reach values comparable with the Fermi energy. While in nuclei the neutron gap is smaller, it is still of the order of a few percent of the Fermi energy. The pairing mechanism in these systems is due to short range attractive interactions between fermions and the size of the Cooper pair is either comparable to the inter-particle separation or it can be as big as a nucleus, which is still relatively small in size. Such a strong pairing gap is the result of the superposition of a very large number of particle–particle configurations, which contribute to the formation of the Cooper pairs. These systems have been shown to be the host of a large number of remarkable phenomena, in which the large magnitude of the pairing gap plays an essential role: quantum shock waves, quantum turbulence, Anderson–Bogoliubov–Higgs mode, vortex rings, domain walls, soliton vortices, vortex pinning in neutron star crust, unexpected dynamics of fragmented condensates and role of pairing correlations in collisions on heavy-ions, Larkin–Ovchinnikov phase as an example of a Fermi supersolid, role of pairing correlations in controlling the dynamics of fissioning nuclei.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Navigating collinear superspace

We introduce a new set of effective field theory rules for constructing Lagrangians with \( \mathcal{N} \) = 1 supersymmetry in collinear superspace. In the standard superspace treatment, superfields are functions of the coordinates \( \left({x}^{\mu },{\theta}^{\alpha },{\theta}^{\dagger \overset{\cdot }{\alpha }}\right) \) , and supersymmetry preservation is manifest at the Lagrangian level in part due to the inclusion of auxiliary F - and D -term components. By contrast, collinear superspace depends on a smaller set of coordinates ( x μ , η , η † ), where η is a complex Grassmann number without a spinor index. This provides a formulation of supersymmetric theories that depends exclusively on propagating degrees of freedom, at the expense of obscuring Lorentz invariance and introducing inverse momentum scales. After establishing the general framework, we construct collinear superspace Lagrangians for free chiral matter and non-Abelian gauge fields. For the latter construction, an important ingredient is a superfield representation that is simultaneously chiral, anti-chiral, and real; this novel object encodes residual gauge transformations on the light cone. Additionally, we discuss a fundamental obstruction to constructing inter- acting theories with chiral matter; overcoming these issues is the subject of our companion paper, where we introduce a larger set of superfields to realize the full range of interactions compatible with \( \mathcal{N} \) = 1. Along the way, we provide a novel framing of reparametrization invariance using a spinor decomposition, which provides insight into this important light-cone symmetry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Signatures and detection prospects for sub-GeV dark matter with superfluid helium

We explore the possibility of using superfluid helium for direct detection of sub-GeV dark matter (DM). We discuss the relevant phenomenology resulting from the scattering of an incident dark matter particle on a Helium nucleus. Rather than directly exciting quasi-particles, DM in this mass range will interact with a single He atom, triggering an atomic cascade which eventually also includes emission and thermalization of quasi-particles. We present in detail the analytical framework needed for modeling these processes and determining the resulting flux of quasi-particles. We propose a novel method for detecting this flux with modern force-sensitive devices, such as nanoelectro-mechanical system (NEMS) oscillators, and derive the sensitivity projections for a generic sub-GeV DM detection experiment using such sensors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Vibrational, Elastic Properties and Sound Velocities of MgTi{sub 2}O{sub 4} Spinel

Zone-centre phonon frequencies, elastic properties and sound velocities along high-symmetry directions and Debye temperature of MgTi{sub 2}O{sub 4} spinel is calculated using rigid ion model. The interatomic interactions up to third nearest neighbours are also calculated. The main results of the present paper reveals that the tetrahedral bonding is stronger than the octahedral bonding. Calculated results are compared with previously available experimental and theoretical data. It is observed that the obtained results are in good agreement with the available results in literature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Theory of multi-orbital topological superconductivity in transition metal dichalcogenides

We study possible superconducting states in transition metal dichalcogenide (TMD) monolayers, assuming an on-site pairing potential that includes both intra- and inter-orbital terms. We find that if the mirror symmetry with respect to the system’s plane is broken (e.g., by a substrate), this type of pairing can give rise to unconventional superconductivity, including time-reversal-invariant nodal and fully gapped topological phases. Using a multi-orbital renormalization group procedure, we show how these phases may result from the interplay between the local Coulomb repulsion, Hund’s rule coupling, and phonon-mediated attraction. In particular, for a range of interaction parameters, the system transitions from a trivial phase to a nodal phase and finally to a gapped topological phase upon increasing the strength of the mirror symmetry breaking term.

36 MATERIALS SCIENCE↗
DOE OSTI2017

OpenRBC: A Fast Simulator of Red Blood Cells at Protein Resolution

We present OpenRBC, a coarse-grained molecular dynamics code, which is capable of performing an unprecedented in silico experiment—simulating an entire mammal red blood cell lipid bilayer and cytoskeleton as modeled by multiple millions of mesoscopic particles—using a single shared memory commodity workstation. To achieve this, we invented an adaptive spatial-searching algorithm to accelerate the computation of short-range pairwise interactions in an extremely sparse three-dimensional space. The algorithm is based on a Voronoi partitioning of the point cloud of coarse-grained particles, and is continuously updated over the course of the simulation. The algorithm enables the construction of the key spatial searching data structure in our code, i.e., a lattice-free cell list, with a time and space cost linearly proportional to the number of particles in the system. The position and the shape of the cells also adapt automatically to the local density and curvature. The code implements OpenMP parallelization and scales to hundreds of hardware threads. It outperforms a legacy simulator by almost an order of magnitude in time-to-solution and >40 times in problem size, thus providing, to our knowledge, a new platform for probing the biomechanics of red blood cells.

Biophysics↗
DOE OSTI2024

From anti-Arrhenius to Arrhenius behavior in a dislocation-obstacle bypass: Atomistic simulations and theoretical investigation

Dislocations are the primary carriers of plasticity in metallic materials. Understanding the basic mechanisms for dislocation movement is paramount to predicting the material mechanical response. Relying on atomistic simulations, we observe a transition from non-Arrhenius to Arrhenius behavior in the rate of an edge dislocation overcoming the long-range elastic interaction with a prismatic loop in tungsten. Close to the critical resolved shear stress, the process shows a non-Arrhenius behavior at low temperatures. However, as the temperature increases, the activation entropy starts to dominate, leading to a traditional Arrhenius-like behavior. Here, we have computed the activation entropy analytically along the minimum energy path following Schoeck’s method, which captures the cross-over between anti-Arrhenius and Arrhenius domains. Also, the Projected Average Force Integrator (PAFI), another simulation method to compute free energies along an initial transition path, exhibits considerable concurrence with Schoeck’s formalism. We conclude that entropic effects need to be considered to understand processes involving dislocations bypassing elastic barriers close to the critical resolved shear stress. More work needs to be performed to fully understand the discrepancies between Schoeck’s and PAFI compared to molecular dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A step-wise ion hydration model of aqueous electrolyte solution: The 2:2, 2:1 and 1:2 punches

In prior work we showed that the current step-wise hydration framework shows promise for modeling the concentration-activity relationship in aqueous solutions of ions representing 1:1 electrolytes, alcohols, and sugars, and accurately predicting the partitioning of the water between the bound (hydrated) pool and an unbound (free) pool. Here, in the current work, we continue developing this framework by applying it to ions representing 2:2, 2:1 and 1:2 electrolyte in aqueous solutions. Properly treating ion hydration as an ideal solution along with employing mean spherical approximation of the long-range electrostatic interactions yields an accurate model for most electrolytes considered. This simple model of aqueous solutions is robust for electrolytes and organic solutes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Molecular-level understanding of gibbsite particle aggregation in water

Using molecular dynamics simulations, we investigate the molecular scale origin of crystal face selectivity when one gibbsite particle attaches to another in water. A comparison of the free energy per unit surface area of particle–particle attachment indicates that particle attachment through edge surfaces, where the edge surfaces are either (1 0 0) or (1 1 0) crystal faces, is more energetically favorable compared to attachment between two basal surfaces (i.e., (0 0 1) crystal faces) or between the basal surface of one particle and the edge surface of another. This result suggests that gibbsite crystals with low basal/edge surface area ratio will preferentially attach through edge surfaces, potentially helping the crystals grow laterally. However, for larger gibbsite particles (high basal/edge surface area ratio) the total free energy, not normalized by surface area, of particle attachment through the basal surfaces is lower (more negative) than attachment through the edge surfaces, indicating that larger gibbsite particles will preferentially aggregate through basal surface attachments. The short-range electrostatic interactions including the interparticle hydrogen bonds from surface –OH groups drive particle attachment, and the dominant contribution to the free energy minimum is enthalpic rather than entropic. However, the enthalpy of basal-edge attachment is significantly offset by the entropy leading to a higher free energy (less negative) compared to that of basal-basal attachment. Study of the free energy for a few imperfect attachments of two particles indicates a higher free energy (i.e., less negative, less stable), compared to a perfect attachment

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Resolved-particle simulation by the Physalis method: Enhancements and new capabilities

We present enhancements and new capabilities of the Physalis method for simulating disperse multiphase flows using particle-resolved simulation. The current work enhances the previous method by incorporating a new type of pressure-Poisson solver that couples with a new Physalis particle pressure boundary condition scheme and a new particle interior treatment to significantly improve overall numerical efficiency. Further, we implement a more efficient method of calculating the Physalis scalar products and incorporate short-range particle interaction models. We provide validation and benchmarking for the Physalis method against experiments of a sedimenting particle and of normal wall collisions. We conclude with an illustrative simulation of 2048 particles sedimenting in a duct. In the appendix, we present a complete and self-consistent description of the analytical development and numerical methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗