Search NASA⌕ Search

SEARCH · Search NASA

Results for “solids flow”

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 361 records · Page 20

Universal behavior in fragmenting brittle, isotropic solids across material properties

A bonded particle model is used to explore how variations in the material properties of brittle, isotropic solids affect critical behavior in fragmentation. To control material properties, a model is proposed which includes breakable two- and three-body particle interactions to calibrate elastic moduli and mode I and mode II fracture toughnesses. In the quasistatic limit, fragmentation leads to a power-law distribution of grain sizes which is truncated at a maximum grain mass that grows as a nontrivial power of system size. In the high-rate limit, truncation occurs at a mass that decreases as a power of increasing rate. A scaling description is used to characterize this behavior by collapsing the mean-square grain mass across rates and system sizes. Consistent scaling persists across all material properties studied, although there are differences in the evolution of grain size distributions with strain as the initial number of grains at fracture and their subsequent rate of production depend on Poisson's ratio. Importantly, this evolving granular structure is found to induce a unique rheology where the ratio of the shear stress to pressure, an internal friction coefficient, decays approximately as the logarithm of increasing strain rate. The stress ratio also decreases at all rates with increasing strain as fragmentation progresses and depends on elastic properties of the solid.

36 MATERIALS SCIENCE↗

Model Validation and Uncertainty Quantification on the KRUSTY Microreactor Design Using GRIFFIN Neutron Transport Code [Poster]

Argonne National Laboratory (ANL) and INL have developed a GRIFFIN steady state neutronics model for the multiphysics simulations of the Kilopower Reactor Using Sterling TechnologY (KRUSTY) microreactor in the Multiphysics Object Oriented Simulation Environment (MOOSE). The reliability of such deterministic neutronics models can be validated by comparing with computations from Monte Carlo codes (e.g. MCNP, SERPENT, OpenMC, Shift, etc). Furthermore, potential modeling/design improvements can be identified by incorporating uncertainty quantification (UQ), which can be performed by MOOSE’s Stochastic Tools Module (STM). KRUSTY is a prototype for a 5-kW thermal nuclear-powered space reactor. Its primary components consist of nuclear fuel, heat pipes, a control rod, a reflector, and the shielding. The fuel consists of 3 stacked U-7.65Mo cylinders with a hole in the center for the control rod. 8 liquid sodium heat pipes transfer fission energy from the solid fuel block to the Sterling power conversion system where the energy is extracted, and the cooled sodium flows back to the core via capillary action . The movable Boron Carbide control rod regulates the neutron population during startup or when a reactor temperature boost is needed . The beryllium oxide reflector is in 3 places in the reactor; it surrounds the core axially, it lies beneath the core on a platen, and it is present in the shim. The axial and lower reflectors rest on an adjustable stainless-steel platen that moves upward to cover the fuel and help the reactor reach criticality. Lastly, radial stainless steel surrounds the core offering protection from radiation exposure .

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Method for solid-liquid separations

Devices, systems, and methods are disclosed to separate solids from liquids. A vessel includes a screw, an inlet, a first outlet, and a second outlet. The inlet receives a first slurry having an incoming solids concentration. The first slurry consists of solids and a contact liquid. The first outlet is disposed toward the first end of the vessel. The second outlet is disposed toward a second end of the vessel. The screw preferentially conveys the solids over the contact liquid towards the second end of the vessel, the solids displacing at least a portion of the contact liquid, causing that portion of the contact liquid to flow toward the first outlet. A melting device is included.

Baxter, Larry↗

Pushing the Envelope-Moving Dynamic Building Envelope Thermal Energy Storage Systems Mainstream: Preprint

Buildings contribute to nearly 40% of the U.S. national energy consumption and a significant portion of CO2 emissions. More importantly, disadvantaged communities are disproportionately affected by energy burden and thermal discomfort in their homes. This paper will discuss two novel DOE's BTO supported thermal energy storage (TES) integrated dynamic building envelope technologies, their ability to harvest ambient energy, and improve energy efficiency by reducing HVAC loads and peak electricity demand while enhancing energy and thermal resilience in buildings. The paper will also discuss the recent advancements that have made these systems more affordable and easier to integrate into new and existing buildings. The first solution is a thermally anisotropic building envelope (TABE) system that can redirect ambient thermal energy (heat or coolness) from diurnal outdoor conditions, solar irradiance, and night sky radiation from the envelope to a hydronic loop. The redirected thermal energy can be stored in a TABE-integrated thermal energy storage system and use the stored energy to offset HVAC energy use and peak demand. The second solution is an innovative plug-and-play thermal switch in the form of insertable plugs integrated with a phase change material (PCM). The plug can vary its thermal resistance based on the indoor and outdoor conditions, thus allowing preferential directional heat flow, and enhancing utilization of free ambient cooling and heating to charge/discharge the PCM, much like a solid-state economizer. While the fist solution can be actively controlled, the second solution is passive, requiring no external power, and work solely based on the ambient temperature.

anisotropic envelope↗

Pushing the Envelope—Moving Dynamic Building Envelope Thermal Energy Storage Systems Mainstream

Buildings contribute to nearly 40% of the U.S. national energy consumption and a significant portion of CO2 emissions. More importantly, disadvantaged communities are disproportionately affected by energy burden and thermal discomfort in their homes. This paper will discuss two novel DOE's BTO supported thermal energy storage (TES) integrated dynamic building envelope technologies, their ability to harvest ambient energy, and improve energy efficiency by reducing HVAC loads and peak electricity demand while enhancing energy and thermal resilience in buildings. The paper will also discuss the recent advancements that have made these systems more affordable and easier to integrate into new and existing buildings. The first solution is a thermally anisotropic building envelope (TABE) system that can redirect ambient thermal energy (heat or coolness) from diurnal outdoor conditions, solar irradiance, and night sky radiation from the envelope to a hydronic loop. The redirected thermal energy can be stored in a TABE-integrated thermal energy storage system and use the stored energy to offset HVAC energy use and peak demand. The second solution is an innovative plug-and-play thermal switch in the form of insertable plugs integrated with a phase change material (PCM). The plug can vary its thermal resistance based on the indoor and outdoor conditions, thus allowing preferential directional heat flow, and enhancing utilization of free ambient cooling and heating to charge/discharge the PCM, much like a solid-state economizer. While the fist solution can be actively controlled, the second solution is passive, requiring no external power, and work solely based on the ambient temperature.

Shrestha, Som↗

Viscometric flow of dense granular materials under controlled pressure and shear stress

This study examines the flow of dense granular materials under external shear stress and pressure using discrete element method simulations. In this method, the material is allowed to strain along all periodic directions and adapt its solid volume fraction in response to an imbalance between the internal state of stress and the external applied stress. By systematically varying the external shear stress and pressure, the steady rheological response is simulated for: (1) rate-independent quasi-static flow; and (2) rate-dependent inertial flow. The simulated flow is viscometric with non-negligible first and second normal stress differences. While both normal stress differences are negative in inertial flows, the first normal stress difference switches from negative to slightly positive, and second normal stress difference tends to zero in quasi-static flows. The first normal stress difference emerges from a lack of coaxiality between a second-rank contact fabric tensor and strain rate tensor in the flow plane, while the second normal stress difference is linked to an excess of contacts in the shear plane compared with the vorticity direction. A general rheological model of second order (in terms of strain rate tensor) is proposed to describe the two types of flow, and the model is calibrated for various values of interparticle friction from simulations on nearly monodisperse spheres. The model incorporates normal stress differences in both regimes of flow and provides a complete viscometric description of steady dense granular flows.

42 ENGINEERING↗

Interfacial Reactivity of Silicon Electrodes: Impact of Electrolyte Solvent and Presence of Conductive Carbon

Silicon (Si) is a promising high-capacity material for lithium-ion batteries; however, its limited reversibility hinders commercial adoption. Approaches such as particle and crystallite size reduction, introduction of conductive carbon, and use of different electrolyte solvents have been explored to overcome these electrochemical limitations. Herein, operando isothermal microcalorimetry (IMC) is used to probe the influence of silicon particle size, electrode composition, and electrolyte additives fluoroethylene carbonate and vinylene carbonate on the heat flow during silicon lithiation. In this work, the IMC data are complemented by X-ray photoelectron and Raman spectroscopies to elucidate differences in solid electrolyte interphase (SEI) composition. Nanosized (~50 nm, n-Si) and micrometer-sized (~4 μm, μ-Si) silicon electrodes are formulated with and without amorphous carbon and electrochemically lithiated in ethylene carbonate (EC), fluoroethylene carbonate (FEC), or vinylene carbonate (VC) based electrolytes. Notably, n-Si electrodes generate 53–61% more normalized heat relative to their μ-Si counterparts, consistent with increased surface area and electrode/electrolyte reactivity. Introduction of amorphous carbon significantly alters the heat flow profile where multiple exothermic peaks and increased normalized heat dissipation are observed for all electrolyte types. Notably, the VC-containing electrolyte demonstrates the greatest normalized heat dissipation of the electrode compositions tested showing as much as a 50% increase compared to the EC or FEC counterparts. The results are relevant to the understanding of silicon negative electrode function in the presence of electrolyte additives and provide insight relative to silicon containing cell reactivity and safety.

25 ENERGY STORAGE↗

Understanding effects of chemical complexity on helium bubble formation in Ni-based concentrated solid solution alloys based on elemental segregation measurements

Here, helium bubble formation and swelling were systematically studied in Ni-based concentrated solid solution alloys containing different numbers and types of elements. Our microscopy analysis showed that although increasing the alloy chemical complexity helps suppress bubble formation in general, there is no monotonic relationship between the bubble growth rate and the number of alloying elements. Certain elements (e.g., Fe and Pd) are more effective in suppressing bubble growth than others (e.g., Cr and Mn). Atom probe tomography was applied to accurately measure elemental segregation around bubbles, revealing unique effects of certain alloying elements on vacancy migration towards bubbles. More specifically, the high vacancy mobility via Cr sites leads to a large vacancy flux and an increased bubble size, while the high degree of atomic size mismatch introduced by Pd helps deflect vacancy flow away from bubbles and decrease the amount of swelling. The effects identified in this study provide new strategies to design concentrated solid solutions with superior resistance to swelling.

36 MATERIALS SCIENCE↗

Application of NEAMS Codes to Capture MSR Phenomena

This report documents work completed in FY21 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program’s Molten Salt Reactor (MSR) Application Drivers activity at Argonne. The common focus was on identifying the modeling and simulation functional requirements for designing, licensing, and operating MSRs and applying those capabilities already developed in NEAMS tools to example problems of interest. The four main parts of this report each focuses on a specific area of simulation physics as it relates to MSR phenomena: fuel evolution, chemistry, computational fluid dynamics (CFD), and systems analysis. In terms of fuel evolution, which includes depletion, decay, on-line separations, and transmutation, the current state of computational capabilities for modeling this behavior in liquid-fueled molten salt reactor is discussed. Some of the functional requirements to accomplish the various applications of MSR fuel depletion modeling are highlighted, followed by a summary of recent approaches and code development activities. The chemistry functional requirements were discussed in the context of several applications of high importance for MSRs, such as corrosion, salt chemistry, and species transport. Each of these types of chemistry modeling have considerable impact on various aspects of reactor applications, including informing on reactor designs, improving operational efficiencies, analyzing safety and reactivity concerns, and estimating the mechanistic source term of the reactor. A brief overview is also provided on code development activities ongoing under NEAMS relevant to chemistry modeling of MSRs. In terms of CFD applications, the state-of-the-art spectral element code, Nek5000 was used to model the fluid dynamics within a full core of the Molten Salt Fast Reactor (MSFR) concept designed as part of the Euratom EVOL project. This concept was selected as the challenge problem because of its similar features to several U.S. industry concepts. The goal was to model some of the fast MSR design challenges, including potential large internal re-circulations, the need of accurate tracking of delayed neutron precursors (DNP), and the lack of relevant thermal-hydraulics models/correlations, etc. Therefore, a series of CFD models were created for the MSFR core cavity using a k – τ model two-equation model for the turbulence. These first full core results demonstrated that any potential recirculation zones could be properly identified with the current NEAMS CFD capabilities. The development of these models will also set the stage for future testing of Nek5000’s functionalities to model other MSR thermal-fluid phenomena. Lastly, the validation of SAM against experimental data from the Molten Salt Reactor Experiment, which started in FY19, continues with the inclusion of modeling the reactivity insertion tests. This involved using SAM and its point kinetics model for flowing fuel salt to recreate the time dependent power changes and response after positive reactivity insertions at the 1, 5, and 8 MWt power levels. Through these exercises, several code modifications were suggested to the SAM development team and accommodated to enable closer agreement with solid technical and physical justifications. These include adding a moderator reactivity feedback coefficient as an available input and modifying the solution approach for the point kinetics model. Additional SAM development suggestions for flowing fuel MSRs include adding the capability to allow the moderator power change proportionally with the reactor power and enabling specification of the power and DNP distributions separately.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Confinement-Driven Heterogeneous Benzene Crystallization in Silica Nanopores

Nanoconfinement alters the thermodynamics, dynamics, and kinetics of fluids hosted in nanoscale solid nanopores to an extent that depends on the characteristics of the confining space and the chemistry of the confined fluids. Confinement-induced alterations in the phase behavior of confined energetic fluids under high pressure or low temperature are highly relevant to subsurface and subsea phenomena such as fluid flow in porous media, hydrate formation and dissociation, and gas storage capacity. Although extensive efforts have been directed toward understanding the phase behavior of confined fluids, the role of solid-liquid interfaces in the phase transitions of organic liquids has not been resolved yet. Here, we explore the onset and growth of benzene crystallization confined in 6 nm sized SBA-15 silica nanopores in the temperature range 300-200 K using in situ extended range small-angle and wide-angle X-ray scattering (SAXS/WAXS) measurements and atomistic classical molecular dynamics (MD) simulations. The crystallization onset of confined benzene depresses to 265 K compared to the freezing point of bulk benzene ( ~278 K), followed by the continuous growth of the emerged crystals in the pore space with complete crystallization at 200 K. The orientation of the emerged benzene crystals is dominated by parallel (π-π stacking) and perpendicular (T-shape stacking) orientations along the cylindrical pore radius and pore length, respectively. The onset of benzene crystals occurs heterogeneously on the pore surface and grows continuously toward the pore center. Further, confined benzene undergoes a dynamical crossover from fragile to strong dynamics behavior, inferred from the rotational and translational diffusion. The insights provided by this study have significant implications for the phase transitions of confined organic liquids that are relevant to a wide range of applications in the biological, geological, environmental, and chemical fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermal and hydraulic performance of volumetrically heated triply periodic minimal surface heaters

Advanced heat transfer surfaces, mainly Triply Periodic Minimal Surfaces (TPMSs) have great potential to increase heat transfer performance over traditional heat transfer surfaces due to their tortuous flow path and higher surface area-volume ratio. The emergence of additive manufacturing of nuclear fuel motivates experimental investigation of such advanced geometries in the context of a nuclear reactor core. Here, this study employed a novel method to explore the performance of gyroid and diamond TPMS geometries. Volumetrically heated heaters were additively manufactured from conductive polymer filament, and airflow was used to simulate the convective heat transfer from the core within a nuclear reactor. The mass flow and volumetric generation rates were varied, local Nusselt number correlations were obtained using embedded thermocouples to measure the solid and fluid temperature, and friction factor correlations were developed using differential pressure measurements. The gyroid TPMS presented a friction factor 1.5 times higher than the diamond TPMS and 90 times that of the rod bundle. TPMS geometries showed Nusselt numbers 8–10 times higher than that of the rod bundle, also maintaining colder internal temperatures, which suggested that they could sustain 250–275 % higher power density than the rod bundle for a given centerline temperature. Hence, having the potential to greatly increase the safety margins and power output of a nuclear reactor core with similar volume.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Interfacial and Activation Energies of Environmentally Abundant Heterogeneously Nucleated Iron(III) (Hydr)oxide on Quartz

Poorly crystalline iron(III) (hydr)oxide nanoparticles are ubiquitous in environmental systems and play a crucial role in controlling the fate and transport of contaminants. Yet, the thermodynamic and kinetic parameters, e.g., the effective interfacial (α^') and apparent activation (E a ) energies, of iron(III) (hydr)oxide nucleation on earth-abundant mineral surfaces have not been determined, which hinders an accurate prediction of iron(III) (hydr)oxide formation and its interactions with other toxic or reactive ions. Here, for the first time, we report experimentally obtained α^' and E a for iron(III) (hydr)oxide nucleation on quartz mineral surfaces, by employing a flow-through, time-resolved grazing incidence small angle X-ray scattering (GISAXS). GISAXS enabled the in situ detection of iron(III) (hydr)oxide nucleation rates under different supersaturations (σ, achieved by varying pH 3.3-3.6) and temperatures (12-35ºC). By quantitative analyses based on classical nucleation theory, α^' was obtained to be 34.6 mJ/m 2 and E a was quantified as 32.8 kJ/mol. The fundamental thermodynamic and kinetic parameters obtained here will advance our fundamental understanding of the surface chemistry and nucleation behavior of iron(III) (hydr)oxides in subsurface and water treatment systems, as well as their effects on the fate and transport of pollutants in natural and engineered water systems. The in situ flow-through GISAXS method can also be adapted to quantify thermodynamic and kinetic parameters at interfaces for many important solid–liquid environmental systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Energy Storage Technologies and U.S. Department of War Requirements

This report offers an overarching primer the energy storage market and assessment of each technology’s suitability for U.S. Department of War applications. We find that battery energy storage is the most promising technology for energy storage applications in terms of energy density and cost. While commercial and advanced Li-ion can satisfy certain electric mobility (e.g., cars, midsized vehicles) and electric flight application (e.g. drones) needs of DoW, next generation technologies like Li metal solid-state batteries, and conversion chemistries (Li-air, Li-CFx) are needed for heavy duty applications like armored vehicles, tanks, airplanes, jets. Batteries should also be supplemented with technologies such as supercapacitors and next-generation flywheels to provide short bursts of high energy for high power needs. Thermal energy storage solutions (sensible and latent heat) and fuel cells (PEMFC and SOFC) are good alternatives for supporting energy storage applications in terms of energy density and cost. Gravity-based energy storage promises the least energy density with high cost and should only be used for niche DoW applications. Of the different gravity energy storage technologies reviewed, pumped-hydro, flywheel, and compressed air energy storage are promising, while solid gravity energy storage appears the least attractive. Military facilities requiring 14 days of energy storage can benefit from flow batteries, especially commercial VFRBs and next generation iron-air systems. While VFRBs are expensive, iron-air batteries have the coupled benefits of low cost and good energy density promising 100+ hours of storage.

25 ENERGY STORAGE↗

Methods for Measuring and Computing the Reference Temperature in Newton’s Law of Cooling for External Flows

Newton’s law of cooling requires a reference temperature (𝑇 𝑟𝑒𝑓 ) to define the heat-transfer coefficient (ℎ). For external flows with multiple temperatures in the freestream, obtaining 𝑇 𝑟𝑒𝑓 is a challenge. One widely used method, referred to as the adiabatic-wall (AW) method, obtains 𝑇 𝑟𝑒𝑓 by requiring the surface of the solid exposed to convective heat transfer to be adiabatic. Another widely used method, referred to as the linear-extrapolation (LE) method, obtains 𝑇 𝑟𝑒𝑓 by measuring/computing the heat flux ($𝑞^{′′}_𝑠$) on the solid surface at two different surface temperatures (𝑇 𝑠 ) and then linearly extrapolating to $𝑞^{′′}_𝑠$ = 0. A third recently developed method, referred to as the state-space (SS) method, obtains 𝑇 𝑟𝑒𝑓 by probing the temperature space between the highest and lowest in the flow to account for the effects of 𝑇 𝑠 or $𝑞^{′′}_𝑠$ on 𝑇 𝑟𝑒𝑓 . This study examines the foundation and accuracy of these methods via a test problem involving film cooling of a flat plate where $𝑞^{′′}_𝑠$ switches signs on the plate’s surface. Results obtained show that only the SS method could guarantee a unique and physically meaningful 𝑇 𝑟𝑒𝑓 where 𝑇 𝑠 =𝑇 𝑟𝑒𝑓 on a nonadiabatic surface $𝑞^{′′}_𝑠$ = 0. The AW and LE methods both assume 𝑇 𝑟𝑒𝑓 to be independent of 𝑇 𝑠 , which the SS method shows to be incorrect. Though this study also showed the adiabatic-wall temperature, 𝑇 𝐴𝑊 , to be a good approximation of 𝑇 𝑟𝑒𝑓 (<10% relative error), huge errors can occur in ℎ about the solid surface where |𝑇 𝑠 −𝑇 𝐴𝑊 | is near zero because where 𝑇 𝑠 =𝑇 𝐴𝑊 , $𝑞^{′′}_𝑠$ ≠ 0.

Newton's law of cooling↗

Single‐Step Deformation Processing of Ultrathin Lithium Foil and Strip

Abstract Next‐generation, high‐efficiency energy storage and conversion systems require development of lithium metal batteries. But the high cost of production and constraints on thickness of lithium (anode) foils continue to limit adoption for integration into battery cell architectures. Here, a novel lithium anode manufacturing solution is demonstrated – single‐step production of ultrathin gauge foil formats directly from solid ingot. Hybrid cutting‐based deformation processes, involving large plastic strains and strain rates, produce foil to sub‐10 µm thickness, with surface quality even superior to present Li anode processing routes. Energy analysis shows the single‐stage processing is ≈50% more efficient than conventional processing by extrusion‐rolling. Through in situ force measurements and high‐speed imaging of the cutting it also characterize – for the first time – the flow stress of Li to strain rates of 800 sec −1 , revealing a power‐law relationship. The results present a paradigm shift in manufacturing and integration of solid lithium anodes for energy applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Flow Regimes for Gas Migration in Liquid Stemming Columns

Some underground engineering projects require the design of an effective stemming column to confine pressurized gases within boreholes. While solid materials such as grout plugs are commonly used, their exothermic curing can damage diagnostic tools, particularly under harsh in situ conditions involving elevated pressure and temperature.

42 ENGINEERING↗

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.

97 MATHEMATICS AND COMPUTING↗

Periodic boundary conditions for arbitrary deformations in molecular dynamics simulations

A generalization of the Lees-Edwards periodic boundary conditions (gLE-PBC) for molecular dynamics (MD) simulations is developed to allow for arbitrary deformations to be applied to the domain. The gLE-PBC domain remains a rectangular cuboid regardless of the applied deformation in contrast with the Lagrangian-rhomboid periodic boundary conditions (LR-PBC) where the domain deforms according to the applied deformation. Furthermore, the kinematics of gLE-PBC are validated against pure shear. The gLE-PBC method for interacting systems is then validated against the LR-PBC method and analytical solutions for a solid under isotropic compression, one-dimensional shearing, three-dimensional extension and shearing and for a liquid under Couette flow. Bulk physical properties extracted from the gLE-PBC simulations agree well with values calculated from equilibrium MD simulations. Three dimensional shearing and a deformation with a full velocity gradient matrix are also simulated, showing the range of problems gLE-PBC can explore.

74 ATOMIC AND MOLECULAR PHYSICS↗