Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Domain Overlapping Coupling”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Coupling of CTF and TRACE for Modeling of Transients

This report documents the improvements that have been made to the capabilities for coupling CTF to systems codes-specifically, the US Nuclear Regulatory Commission (NRC) TRACE code. An initial systems coupling capability had been set up previously using a nonoverlapping domain approach with the codes exchanging data at the core boundaries. The present work adds a new approach using overlapping domains, in which the system code models the core as well. A new input format has been added to allow the user to specify the physical quantities to be exchanged and their location in the system model, which gives the flexibility of applying one-way or two-way coupling between the codes using the desired data exchanges. In addition to applying thermal-hydraulics (T/H) boundary condition (BC) values obtained from TRACE, a capability was added to allow CTF to apply flow resistance feedback to TRACE to match the CTF core pressure drop. Support was added for executing parallel CTF models within the CTF-systems coupling. The system coupling capability was successfully applied to a parallel MSLB transient, demonstrating that both the one-way and two-way coupling behaved as expected and provided substantial improvements to numerical stability and routine compared to the previous nonoverlapping domain coupling. An initial capability was also developed for performing restart calculations in CTF which will be used in the future for restarting CTF-TRACE simulations at specific points in the transient simulation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

An optimization-based strategy for peridynamic-FEM coupling and for the prescription of nonlocal boundary conditions

We develop and analyze an optimization-based method for the coupling of a static peridynamic (PD) model and a static classical elasticity model. The approach formulates the coupling as a control problem in which the states are the solutions of the PD and classical equations, the objective is to minimize their mismatch on an overlap of the PD and classical domains, and the controls are virtual volume constraints and boundary conditions applied at the local-nonlocal interface. Our numerical tests performed on three-dimensional geometries illustrate the consistency and accuracy of our method, its numerical convergence, and its applicability to realistic engineering geometries. We demonstrate the coupling strategy as a means to reduce computational expense by confining the nonlocal model to a subdomain of interest, and as a means to transmit local (e.g., traction) boundary conditions applied at a surface to a nonlocal model in the bulk of the domain.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Time resolved mechanism of the isotope selectivity in the ultrafast light induced dissociation in N 2

The time evolution of a vacuum ultraviolet excited N2 molecule is followed all the way from an ultrafast excitation to dissociation by a quantum mechanical simulation. The primary aim is to discern the role of the excitation by a pulse short compared to the vibrational period, to discern the different coupling mechanisms between different electronic states, nonadiabatic, spin orbit, and to analyze the origin of any isotopic effect. We compare the picture in the time and energy domains. The initial ultrafast excitation pumps the molecule to a coherent electronic wave packet to which several singlet bound electronic states contribute. The total nonstationary wave function is given as a coherent sum of nuclear wave packets on each electronic state times the stationary electronic wave function. When the wave packets on different electronic states overlap, they are coupled in a mass-dependent manner whether one uses an adiabatic or a diabatic electronic basis. A weak spin-orbit coupling acts as a bottleneck between the bound singlet part of phase space and the triplet manifold of states in which dissociation takes place. To describe the spin-orbit perturbation that is ongoing in time, an energy-resolved eigenstate representation appears to be more intuitive. In the eigenstate basis, the singlet-to-triplet population transfer is large only between those vibronic eigenstates that are quasiresonant in energy. The states in resonance are different for different excitation energy ranges. Finally, the resonances are mass dependent, which explains the control of the isotope effect through the profile of the pulse.

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

A Code-Agnostic Driver Application for Coupled Neutronics and Thermal-Hydraulic Simulations

While the literature has numerous examples of Monte Carlo and computational fluid dynamics (CFD) coupling, most are hard-wired codes intended primarily for research rather than as standalone, general-purpose applications. In this work, we describe an open source application, ENRICO, that enables coupled neutronic and thermal-hydraulic simulations between multiple codes that can be chosen at runtime (as opposed to a coupling between two specific codes). The application has been designed such that the control flow logic, domain mapping, nonlinear fixed-point iteration, solution transfers, and convergence checks are all agnostic to the underlying physics solvers used. Special emphasis has also been placed on enabling efficient execution on distributed-memory computing environments. The transfer of solution fields between solvers is performed in memory rather than through filesystem I/O. Additionally, solvers can be configured to run on overlapping or disjoint sets of processes. To date, coupling with the OpenMC and Shift Monte Carlo codes, the Nek5000 CFD code, and a simplified heat diffusion and subchannel solver has been implemented in ENRICO. We present results for coupled simulations of a single light-water reactor fuel assembly based on the NuScale reactor using various combinations of the physics solvers. For this problem, the coupled simulations are shown to converge in about four Picard iterations. A comparison of the heat source and temperature distributions computed by ENRICO using OpenMC coupled with Nek5000 and Shift coupled with Nek5000 illustrates remarkable agreement between the codes.

42 ENGINEERING↗
DOE OSTI2023

Coupling of CTF and RELAP5-3D Within an Enhanced Fidelity Nuclear Power Plant Simulator

A robust and accurate multiphysics engineering simulator is being developed to model the core behavior and system response of pressurized water reactors. This simulator relies on the NESTLE and CTF computer codes to model the neutronics and thermal hydraulics (TH), respectively, inside the core on a nodal scale and on the Reactor Excursion and Leak Analysis Program—Three Dimensional (RELAP5-3D) to model the entire nuclear steam supply system. The RELAP5-3D model includes highly detailed nodalization and multidimensional flow modeling throughout the vessel. Previously, pin-resolved data generated via the Virtual Environment for Reactor Analysis core simulator were used to improve the accuracy of the NESTLE core predictions. The engineering simulator being developed as part of this work uses the 3KEYMASTER platform to couple the enhanced NESTLE model to a nodal-fidelity CTF model to balance run time with accuracy; NESTLE provides node-dependent powers to CTF, and CTF provides node-dependent coolant densities and fuel temperatures to NESTLE.An overlapping domain approach is used for the core TH in which RELAP5-3D provides core boundary conditions based on the system response and CTF provides a node-dependent coolant heating rate to the RELAP5-3D core solution. In the preliminary TH demonstration discussed in this paper, CTF and RELAP5-3D provided similar steady-state core predictions, indicating the hydraulic compatibility between the codes, as well as reasonable and expected behavior under hypothetical transient conditions. This provides an initial step in ongoing efforts toward a robust, multiscale TH/neutronics engineering simulator capability.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Concurrent two-way coupling of global and local models across internal boundaries with non-matching discretizations

Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.

Variational Multiscale Discontinuous Galerkin (VMD↗
DOE OSTI2023

Quantum and anharmonic effects in non-adiabatic transition state theory

Quantitative descriptions of non-adiabatic transition rates at intermediate temperatures are challenging due to the simultaneous importance of quantum and anharmonic effects. Here, in this paper, the interplay between quantum effects—for motion across or along the seam of crossing—and anharmonicity in the seam potential is considered within the weak coupling limit. The well-known expression for quantized 1-D motion across the seam (i.e., tunneling) in the linear terms approximation is derived in the thermal domain using the Lagrangian formalism, which is then applied to the case when tunneling is distributed along the seam of crossing (treating motion along the seam classically). For high-frequency quantum modes, a vibrationally adiabatic (VA) approach is developed that introduces to the non-adiabatic rate constant a factor associated with high-frequency wavefunction overlap; this approach treats the high-frequency motion along the seam quantum mechanically. To test these methodologies, the reaction N 2 O ↔ N 2 + O( 3 P) was chosen. CCSD(T)-F12b/cc-pVTZ-F12 explorations of the 3 A'- 1 A' seam of N 2 O revealed that seam anharmonicity has a strong effect on the rate constant (a factor of ~20 at 2000 K). Several quantum effects were found to be significant at intermediate/lower temperatures, including the quantum N–N vibration that was coupled with seam anharmonicity using the VA approach. Finally, a 1-D approximation to non-adiabatic instanton theory is presented to estimate the validity limit of the linear terms model at low temperatures (⁠~ 250 K for N 2 O). We recommend that the assumptions built into many statistical theories for non-adiabatic reactions—harmonic behavior, classical motion, linear terms, and weak coupling—should be verified on a case-by-case basis.

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

A Fundamentally New Coupled Approach to Contact Mechanics via the Dirichlet‐Neumann Schwarz Alternating Method

Contact phenomena are crucial for understanding the behavior of mechanical systems. However, existing computational approaches for simulating mechanical contact often face numerical challenges, such as inaccurate physical predictions, energy conservation errors, and unwanted oscillations. Here, we introduce an alternative technique for simulating dynamic contact based on the non‐overlapping Schwarz alternating method, originally developed for domain decomposition. In multibody contact scenarios, this method treats each body as a separate, non‐overlapping domain and prevents interpenetration using an alternating Dirichlet–Neumann iterative process. This approach has a strong theoretical foundation, eliminates the need for contact constraints, and offers flexibility, making it ideal for multiscale and multiphysics applications. We conducted a numerical comparison between the Schwarz method and traditional methods, such as the Lagrange multiplier and penalty methods, focusing on a benchmark impact problem. Our results indicate that the Schwarz alternating method outperforms traditional methods in several key areas: it provides more accurate predictions for various measurable quantities and demonstrates exceptional energy conservation capabilities. To address unwanted oscillations in contact velocities and forces, we explored various algorithms and stabilization techniques, ultimately opting for the naïve‐stabilized Newmark scheme for its simplicity and effectiveness. Additionally, we validated the efficiency of the Schwarz method in a three‐dimensional impact problem, highlighting its inherent capacity to accommodate different mesh topologies, time‐integration schemes, and time steps for each interacting body.

Schwarz alternating method↗
DOE OSTI2025

SICAS, a new integrated divertor-SOL-pedestal-core framework for self-consistent modeling of ion and impurity transport in fusion devices

We present a new framework for core-edge integration studies named SICAS (SOLPS-ITER coupled to ASTRA-STRAHL) which enables high fidelity simulations of the core, edge, and divertor regions encompassing the transport of ions as well as the impurities through the entire plasma domain. SICAS handles the exchanging of the particle and power fluxes as well as transport coefficients to ensure consistency through the codes. An overlapping region is defined from the inside of the separatrix of the plasma to the SOLPS-ITER core boundary allowing for matching profiles and fluxes between the two codes for a self-consistent approach. The results presented here demonstrated the flexibility of SICAS to simulate different configurations, scenarios, divertor geometries, and plasma species with good agreement with DIII-D experimental data. This tool opens new possibilities in integrated modeling of fusion devices integrating all relevant phenomena in the core and the divertor plasmas. These capabilities are required for the interpretation of current experiments as well as the design of new devices.

core-edge integration↗
DOE OSTI2022

A micromorphic filter for determining stress and deformation measures from direct numerical simulations of lower length scale behavior

A micromorphic filter is presented for the extraction of stress and deformation measures from underlying direct numerical simulation (DNS) of microstructural mechanical response. The filter is consistent with the micromorphic continuum theory of Eringen and Suhubi (1964), and allows the interrogation of a DNS to generate higher order constitutive models. The filter adopts aspects of the overlap coupling techniques for atomistic-continuum scale-bridging in order to update the nodal degrees of freedom of the finite-element-based filter. Here, two methods for the extraction of quantities of interest from the underlying DNS are presented. The first is a least squares fitting technique which is appropriate for finite-element-style macro-scale domains, and the second is a variationally based approach which generalizes the method to any macro-scale domain which has a variational form. No restrictions on the form of the micro-scale simulation are inherent to the approach, provided that the required quantities (stress, position, etc.) can be provided. The filter is applied to a simple homogeneous medium undergoing homogeneous deformation to demonstrate the recovery of a classical continuum response, what we call the null case. We then apply the framework to the same homogeneous medium with a central and edge penny-shaped crack to show micromorphic effects. We demonstrate that varying the size of the filter causes the micromorphic behaviors to arise and then dissipate as the filter becomes smaller than the crack's length scale. Additionally we show the effect of filter size on a chiral element network along with a plate perforated with evenly spaced and perturbed holes through the thickness. The micromorphic effects are shown to be present for the chiral elements under any loading but become most pronounced when a single element is located within each averaging domain. Micromorphic effects become pronounced for the plate upon the onset of localization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Charmed and ϕ meson decay constants from 2+1-flavor lattice QCD *

Abstract On a lattice with 2+1-flavor dynamical domain-wall fermions at the physical pion mass, we calculate the decay constants of , , and . The lattice size is , which corresponds to a spatial extension of fm, with a lattice spacing of fm. For the valence light, strange, and charm quarks, we use overlap fermions at several mass points close to their physical values. Our results at the physical point are MeV, MeV, MeV, MeV, and MeV. The couplings of and to the tensor current ( ) can be derived from ratios and , respectively, which are the first lattice quantum chromodynamics (QCD) results. We also obtain ratios and , which reflect the size of heavy quark symmetry breaking in charmed mesons. Ratios and can be taken as a measure of SU (3) flavor symmetry breaking.

Chen, Ying↗
DOE OSTI2022

Selective posttranslational inhibition of Ca V β 1 -associated voltage-dependent calcium channels with a functionalized nanobody

Ca 2+ influx through high-voltage-activated calcium channels (HVACCs) controls diverse cellular functions. A critical feature enabling a singular signal, Ca 2+ influx, to mediate disparate functions is diversity of HVACC pore-forming α1 and auxiliary Ca V β 1 –Ca V β 4 subunits. Selective Ca V α 1 blockers have enabled deciphering their unique physiological roles. By contrast, the capacity to post-translationally inhibit HVACCs based on Ca V β isoform is non-existent. Conventional gene knockout/shRNA approaches do not adequately address this deficit owing to subunit reshuffling and partially overlapping functions of Ca V β isoforms. Here, we identify a nanobody (nb.E8) that selectively binds Ca V β 1 SH3 domain and inhibits Ca V β 1 -associated HVACCs by reducing channel surface density, decreasing open probability, and speeding inactivation. Functionalizing nb.E8 with Nedd4L HECT domain yielded Chisel-1 which eliminated current through Ca V β 1 -reconstituted Ca V 1/Ca V 2 and native Ca V 1.1 channels in skeletal muscle, strongly suppressed depolarization-evoked Ca 2+ influx and excitation-transcription coupling in hippocampal neurons, but was inert against Ca V β 2 -associated Ca V 1.2 in cardiomyocytes. The results introduce an original method for probing distinctive functions of ion channel auxiliary subunit isoforms, reveal additional dimensions of Ca V β 1 signaling in neurons, and describe a genetically-encoded HVACC inhibitor with unique properties.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Electric Utility Industry Standards Landscape

The electric utility industry relies on robust communication protocols to manage complex electrical grid data. The inherent networked nature of electrical grids, coupled with the radial structure of the “last mile” portion delivering power to end-use customers, presents difficulties in describing electrical models using simple data constructs. The paper provides an overview of communication protocols that address electric utility data, including grid data, and classifies these protocols through identification key characteristics that make them suitable for different electric grid data domains. Because of their variety, grid edge devices and their associated dedicated protocols are assessed by groups: those primarily designed for energy production and storage, those related to flexible loads, and those related to electric vehicles. The intent of this report is to provide guidance for stakeholders to navigate the challenges posed by the numerous overlapping protocols available to address electric grid data. Although further industry review, refinement, and validation of the categorization of these protocols is recommended, the authors propose this categorization as a start to improve electric grid awareness and understanding. In addition, this report includes three recommended industry actions regarding protocols to improve communications on the electric grid.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Watching Polarons Move in the Energy and Frequency Domains Using Color Impedance Spectroscopy

The hybrid electronic–ionic transport property of π-conjugated polymers enables new (opto)electrochemical device constructs for energy conversion and storage and biosensing applications. One major challenge is separating the energy and frequency dependence of Faradaic events–those involving charge transfer and the redox processes of the conjugated backbone–from the non-Faradaic components, such as ionic motion. Herein, we combine optical spectroscopy with electrochemical impedance spectroscopy (EIS) to resolve the frequency response of ionic–electronic coupling as a function of electrochemical doping potential. First, using EIS, we identify two different frequency regimes resulting in potential-dependent capacitive elements on the order of ~10 μF/cm 2 in a high-frequency regime and ~50–150 μF/cm 2 in a low-frequency regime. Given the larger magnitude and greater potential dependence, we posit that polaronic motion is more likely to occur at low frequencies (<1 kHz) and overlaps with ionic motion. The use of color impedance spectroscopy (CIS) enables observation of polaronic motion with frequency modulation. We observe that higher doping potentials show a greater motion of polarons above the DC-bias baseline concentration for onset in electrochemical doping, but all potentials considered demonstrate a critical frequency at which the polaronic motion is “frozen” (~40 Hz). Furthermore, this critical information obtained from CIS in highly dielectric environments offers a unique figure of merit for future studies on electronic–ionic coupling by which to compare across polymer/electrolyte interfaces, including the role of a charge-supporting electrolyte, a solvent, and alternative Faradaic processes (e.g., electrocatalysis).

36 MATERIALS SCIENCE↗
DOE OSTI2024

Two-stage crash process in resistive drift ballooning mode driven ELM crash

We report a two-stage crash process in edge localized mode (ELM) driven by resistive drift-ballooning modes (RDBMs) numerically simulated in a full annular torus domain with a scale-separated four-field reduced MHD (RMHD) model using the BOUT++ framework. In the early nonlinear phase, the small first crash is triggered by linearly unstable RDBMs, and m/n=2/1 magnetic islands are nonlinearly excited by nonlinear coupling of RDBMs as well as their higher harmonics. Here, m is the poloidal mode number, n is the toroidal mode number, the q = 2 rational surface exists near the pressure gradient peak, and q is the safety factor. Simultaneously, middle-n RDBM turbulence develops but is poloidally localized around X-points of the magnetic islands, leading to the small energy loss. The second large crash occurs in the late nonlinear phase. Higher harmonics of m/n=2/1 magnetic islands well develop around the q = 2 surface via nonlinear coupling and make the magnetic field stochastic by magnetic island overlapping. Turbulence heat transport develops at X-points of higher harmonics of m/n=2/1 magnetic islands, resulting in the turbulence spreading in the poloidal direction. The large second crash is triggered when the turbulence covers the whole poloidal region so that the magnetic island generation and magnetic field stochastization before the large crash can be interpreted as ELM precursors. It is concluded that the ELM trigger is attributed to the turbulent spreading in the poloidal direction in synchronization with the magnetic field stochastization and the crash is driven by E × B convection rather than the conventional Rechester–Rosenbluth anomalous electron heat transport.

Seto, H. (ORCID:0000000259278643)↗
DOE OSTI2020

Layer-resolved magnetic proximity effect in van der Waals heterostructures

Magnetic proximity effects are integral to manipulating spintronic, superconducting, excitonic and topological phenomena in heterostructures. These effects are highly sensitive to the interfacial electronic properties, such as electron wavefunction overlap and band alignment. The recent emergence of magnetic two-dimensional materials opens new possibilities for exploring proximity effects in van der Waals heterostructures. In particular, atomically thin CrI 3 exhibits layered antiferromagnetism, in which adjacent ferromagnetic monolayers are antiferromagnetically coupled. In this paper, we report a layer-resolved magnetic proximity effect in heterostructures formed by monolayer WSe 2 and bi/trilayer CrI 3 . By controlling the individual layer magnetization in CrI 3 with a magnetic field, we show that the spin-dependent charge transfer between WSe 2 and CrI 3 is dominated by the interfacial CrI 3 layer, while the proximity exchange field is highly sensitive to the layered magnetic structure as a whole. In combination with reflective magnetic circular dichroism measurements, these properties allow the use of monolayer WSe 2 as a spatially sensitive magnetic sensor to map out layered antiferromagnetic domain structures at zero magnetic field as well as antiferromagnetic/ferromagnetic domains at finite magnetic fields. Our work reveals a way to control proximity effects and probe interfacial magnetic order via van der Waals engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Enabling quantitative analysis of complex polymer blends by infrared nanospectroscopy and isotopic deuteration

Atomic-force microscopy coupled with infrared spectroscopy (AFM-IR) deciphers surface morphology of thin-film polymer blends and composites by simultaneously mapping physical topography and chemical composition. However, acquiring quantitative phase and composition information from multi-component blends can be challenging using AFM-IR due to the possible overlapping infrared absorption bands between different species. Isotope labeling one of the blend components introduces a new type of bond (carbon-deuterium vibration) that can be targeted using AFM-IR and responds at wavelengths sufficiently shifted toward unoccupied regions (around 2200 cm –1 ). In this project, AFM-IR was used to probe the surface morphology and chemical composition of three polymer blends containing deuterated polystyrene; each blend is expected to exhibit various degrees of miscibility. AFM-IR results successfully demonstrated that deuterium labeling prevents infrared spectral overlap and enables the visualization of blend phases that could not normally be distinguished by other scanning probe techniques. The nanoscale domain composition was resolved by fast infrared spectrum analysis. Overall, we presented isotope labeling as a robust approach for circumventing obstacles preventing the quantitative analysis of multiphase systems by AFM-IR.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A massively-parallel, unstructured overset method to simulate moving bodies in turbulent flows

An unstructured overset method capable of performing direct numerical simulation (DNS) and large eddy simulation (LES) of many (O(10 5 )) moving bodies, utilizing many computational cores (O(10 5 )), in turbulent, incompressible fluid flow is presented. Unstructured meshes are attached to bodies and placed within a fixed background domain. Body meshes are allowed to arbitrarily overlap and move throughout the domain. Within each mesh a high resolution, unstructured, non-dissipative finite volume method is used to solve for the flow field. Boundary conditions for each mesh are provided by interpolation from flow solutions on overlapping meshes. When many unstructured meshes of different resolution overlap, care is required in the connection between the different flow solutions. An interpolant is created which seeks to preserve volume conservation of flow quantities between meshes regardless of overlapping mesh differences. An implicit fractional step method is used for time advancement, requiring the calculation of a predicted fluid velocity and corrector pressure field. For the predictor step, the resulting interpolation is directly introduced into the implicit equations for the predicted flow field. For the corrected pressure field, the continuity between meshes is weakly enforced using a penalty formulation. The pressure formulation is symmetric, positive-definite and non-singular resulting in a formulation which is readily solvable using traditional iterative matrix inversion techniques. An Arbitrary Euler-Lagrangian (ALE) method coupled to a 6 degrees of freedom rigid body equation system (6-DOF) is used for body motion. For rotation, a quaternion representation is used to solve Euler's equations of rigid body motion. A linear spring damper model, which uses geometry information readily available from the overset assembly process, is used for collisions. Validation of the method for canonical flow fields is presented including assessment of order of accuracy and kinetic energy conservation properties. Particle-resolved direct numerical simulation (PR-DNS) of single particles in various flow fields are presented for validation. PR-DNS results of 50,000 spherical particles freely moving within turbulent channel flow are shown as a demonstration of the method at full scale. LES results of a marine propeller under crashback conditions are shown to demonstrate the ability to simulate highly unsteady turbulent flows over complex, moving geometries.

97 MATHEMATICS AND COMPUTING↗