Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Interface stability”

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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Nanoionics Drastically Accelerating Mass Transfer at Elevated Temperatures over 750 °C

Nanoionics were previously considered thermally unstable and infeasible for devices operating above 500 °C. Here, we elucidate the design principle for establishing stable nanoionics from various oxides. We utilized reversible solid oxide cells (SOCs) as the test bed and implemented nanoionics using atomic layer deposition (ALD). We demonstrate a straightforward, interface-controlled, practical approach to render a conformal, ∼15 nm thick ALD film, which initially thermodynamically favors the formation of a solid solution with the substrate into surface nanoionics with single or double layers of nanograins with random crystal orientations. The nanoionics exhibited conductivity estimated to be 7 orders of magnitude higher than that of their bulk-scale counterpart. They demonstrated conformability with uniform grain sizes of ∼15 nm, even after electrochemical operation for ∼500 h at 750 °C and 1000 h at 850 °C. The thermal stability and conductivity of such nanoionics represent a conceptual and technological framework in nanoionics.

atomic layer deposition↗
DOE OSTI2023

Stabilizing Liquids Using Interfacial Supramolecular Assemblies

Abstract Stabilizing liquids based on supramolecular assembly (non‐covalent intermolecular interactions) has attracted significant interest, due to the increasing demand for soft, liquid‐based devices where the shape of the liquid is far from the equilibrium spherical shape. The components comprising these interfacial assemblies must have sufficient binding energies to the interface to prevent their ejection from the interface when the assemblies are compressed. Here, we highlight recent advances in structuring liquids based on non‐covalent intermolecular interactions. We describe some of the progress made that reveals structure–property relationships. In addition to treating advances, we discuss some of the limitations and provide a perspective on future directions to inspire further studies on structured liquids based on supramolecular assembly.

Gu, Peiyang↗
DOE OSTI2023

Highly soluble copper(ι) iodide-based hybrid luminescent semiconductors containing molecular and one-dimensional coordinated anionic inorganic motifs

Copper(I) iodide-based inorganic–organic hybrid semiconductors are considered promising materials for various optoelectronic applications. Here, using imidazolium bridged benzotriazoles as multidentate ligands, we designed and synthesized a series of all-in-one (AIO)-type CuI hybrid materials that hold both ionic and coordinate bonds at the interface of organic and inorganic modules. Their structures range from 0D clusters to 2D extended networks built on various molecular (0D) and chain-like (1D) anionic inorganic motifs that are interconnected through cationic ligands via multiple Cu–N bonds. Benefiting from their unique bonding nature, these compounds exhibit high stability towards heat and moisture and can be well dissolved in polar aprotic solvents. They emit low-energy light spanning from yellow to red color (550–625 nm). Finally, the electronic structure, emission mechanism, and effect of different coordination modes on their photophysical properties were studied using both experimental and theoretical methods, which has provided insight into the structure–property relationship of these inorganic–organic semiconductors.

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

Optical N-plasmon: topological hydrodynamic excitations in graphene from repulsive Hall viscosity

Abstract Edge states occurring in Chern and quantum spin-Hall phases are signatures of the topological electronic band structure in two-dimensional (2D) materials. Recently, a new topological electromagnetic phase of graphene characterized by the optical N-invariant was proposed. Optical N-invariant arises from repulsive Hall viscosity in hydrodynamic many-body electron systems, distinct from the Chern and Z 2 invariants. In this paper, we introduce the topologically protected edge excitation—optical N-plasmon of interacting many-body electron systems in the topological optical N-phase. These optical N-plasmons are signatures of the topological plasmonic band structure in 2D materials. We demonstrate that optical N-plasmons exhibit unique dispersion relations, stability against various boundary conditions, and edge profiles when compared with the topologically trivial edge magneto plasmons. Based on the optical N-plasmon, we design an ultra sub-wavelength broadband topological hydrodynamic circulator, which is a chiral quantum radio-frequency circuit component crucial for information routing and interfacing quantum–classical computing systems. Furthermore, we reveal that optical N-plasmons can be effectively tuned by the neighboring dielectric environment without breaking the topological properties. Our work provides a smoking gun signature of topological electromagnetic phases occurring in 2D materials arising from repulsive Hall viscosity.

Physics↗
DOE OSTI2020

Engineered Barrier System R&D and International Collaborations – LANL (FY20): Spent Fuel and Waste Disposition

The United States Department of Energy’s Spent Fuel and Waste Disposition (SFWD) program is investigating the design and safety function of generic nuclear geologic repositories in a variety of geologic settings (salt, argillite, and crystalline rock). Different configurations and loadings of spent nuclear fuel and waste within disposal canisters are also being investigated, some of which have the potential to generate repository temperatures higher than previously considered (i.e., temperatures >100ºC) by foreign and domestic concepts. This report expands on engineered barrier material stability in a high temperature crystalline rock repository through high temperature hydrothermal experiments. Experiments were designed to develop engineered barrier system (EBS) concepts in a hightemperature crystalline environment in 1) bentonite-Grimsel Granodiorite interactions, 2) bentonite-cement reactions, and 3) interaction between waste canister materials and bentonite. Experiment results are applied to understanding long-term repository performance in terms of radionuclide isolation. One hydrothermal experiment was completed in the rocking autoclaves at LANL in FY-20: IEBS-6 (Grimsel Granodiorite + Wyoming bentonite + cured ordinary Portland cement + Grimsel Granodiorite synthetic groundwater, 250ºC/150 bar, 8 weeks). Several other experiments were planned but were delayed due to the COVID-19 pause in laboratory work. In addition, some characterization of the reaction products of experiments conducted in FY-20 was prevented by COVID-19. Mineral phase chemistry was not measured via electron microprobe analyses for IEBS-6. Quantitative X-ray diffraction results from IEBS-6 and HBT-1 were not completed. The missing analyses and discussion of the results will be included in next year’s report. New characterization that was completed includes scanning electron microscopy of reaction products from IEBS-6, quantitative X-ray diffraction results are presented from IEBS-1 through IBES-5, preliminary scanning electron microscope images and chemical analyses for IEBS-6, X-ray diffraction of the clay fraction from all experiments, measurement/imaging of mineral growth on the surface of steel coupons. Major observations pertaining to bentonite stability in a Grimsel Granodiorite environment include the stability of Na-montmorillonite at 250°C, the formation of trace CSH phases, and the formation of bentonite colloids on experiment cooling. The addition of a cured chip of Portland cement to the bentonite-Grimsel system results in slightly higher pH values and the formation of diverse secondary mineral phases that were not observed in the previous experiments (e.g., analcime, garronite, CSH phases). The new characterization efforts related to the interaction of stainless-steel coupons and bentonite clay focused on thickness and mineralogy of phases that formed at the steel surface. In Wyoming bentonite + Grimsel Granodiorite systems, newly formed minerals at the bentonitesteel coupon interface included alteration of the outermost steel edge to Fe,Cr-oxide phases, followed by Fe-rich phyllosilicates (Fe-saponite, chlorite) and interbedded Fe,Cr,Ni-sulfide phases (pentlandite). Hydrothermal experiments were completed to assess uranium-carbonate complexation at conditions relevant to high-temperature disposal. Autoclave solubility experiments were conducted at 150 to 250ºC with a range of carbonate and uranium concentrations. The experiment results were characterized via situ UV-Visible spectroscopy and synchrotron-based in situ XAS techniques. Results show a significant decrease in the stability of uranyl-carbonate complexes at temperatures above 100ºC along with the precipitation of uranium oxides. Further, at T > 200ºC, results show that uranyl-hydroxyl complexes control solubility of uranium instead of the previously predicted uranyl-carbonate species. These results are significant for understanding the mobility of uranium in the EBS, which will likely contain carbonate-rich fluids. International research efforts focused on three main areas: 1) participation in international conferences, 2) building collaborations with foreign repository programs, and 3) the initiation of an experimental program to complement the full-scale HotBENT test at the Grimsel test site. This experiment included Wyoming bentonite + low carbon steel + Grimsel Granodiorite synthetic groundwater and was run at the planned maximum temperature of the HotBENT test (200°C). Complete characterization of reaction products was hindered by disruptions to laboratory work but will be reported in the next FY. The experimental results obtained in FY-20 continue to document the wide-ranging effects of bulk composition and pressure-temperature conditions in the mineralogical and geochemical evolution of a high-temperature repository environment. Concepts developed will be used to inform models of long-term material stability in a generic crystalline rock-hosted repository.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Electrolyte Reactivity on the MgV 2 O 4 Cathode Surface

Predictive understanding of the solvation-dependent reactivity and molecular interaction of electrolyte ions and solvent molecules on reactive electrodes has been a major challenge but is essential for addressing instabilities and surface passivation that occur at electrode-electrolyte interface (EEI) of multivalent Mg batteries. In this work, the isolated intrinsic reactivities of prominent chemical species present in magnesium bis(trifluoromethanesulfonimide) (Mg(TFSI) 2 ) in diglyme (G2) electrolytes, including ionic (TFSI - , [Mg(TFSI)] + , [Mg(TFSI):G2] + , [Mg(TFSI):2G2] + ) as well as neutral molecules (G2) on magnesium vanadate cathode (MgV 2 O 4 ) surface has been studied using a combination of first-principles calculations and multimodal analysis of well-defined cathode electrolyte interphase (CEI) layers. Here, our calculations show that non-solvated [Mg(TFSI)] + is the strongest adsorbing species on the MgV 2 O 4 surface compared to all other ions while fully solvated [Mg(TFSI):2G2] + are least favorable to decomposition. The cleavage of C-S bonds in TFSI - to form CF 3 - is predicted to be most desired pathway for all ionic species, which is followed by the cleavage of C-O bonds of G2 to yield CH 3 + or OCH 3 - species. The strong stabilization and electron transfer between ionic electrolyte species and MgV 2 O 4 is found to significantly favor these decomposition reactions on the surface compared to intrinsic gas phase dissociation. Experimentally, we used state-of-the-art ion soft landing to selectively deposit mass-selected TFSI - , [Mg(TFSI):G2] + and [Mg(TFSI):2G2] + on MgV 2 O 4 thin film to form well-defined electrolyte-MgV 2 O 4 interface. Analysis of soft-landed interphase using X-ray photoelectron, X-ray absorption near edge structure, electron energy-loss spectroscopies as well as transmission electron microscopy confirmed the presence of decomposition species (e.g., MgFx, carbonates) formed in the interfacial region and the higher amount of MgFx with [Mg(TFSI):G2] + , which corroborates the theoretical observation. Overall, we established the mechanistic pathway for the electrolyte-induced formation of passivating fluorides on MgV 2 O 4 cathode facilitated by the surface adsorption and charge transfer, which provided essential knowledge for rational design of stable electrolytes for multivalent cathodes.

Cathode-electrolyte interphase formation↗
DOE OSTI2026

Pressure Heterogeneity and Material Utilization in Thin-Film Solid-State Cathodes

Intimate interfacial contact between the solid electrolyte and the cathode active material is critical for maximizing cathode utilization in solid-state batteries. However, volume changes during electrochemical cycling induce internal stresses that drive interfacial degradation, particularly under nonuniform stack pressure. In this study, we employ in situ energy-dispersive X-ray diffraction tomography to visualize and quantify reaction heterogeneities across a 3 mm-diameter solid-state cathode with a well-defined interface. Our results reveal that regions under a lower stack pressure exhibit reduced material utilization and reversibility, which negatively affect the high-pressure regions. Interfacial degradation further impedes lithium-ion transport and amplifies microscale reaction heterogeneity. These findings highlight the critical role of stack pressure distribution in governing interfacial stability and electrochemical performance, offering important design insights into practical solid-state battery systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Immunogenicity and therapeutic targeting of a public neoantigen derived from mutated PIK3CA

Public neoantigens (NeoAgs) represent an elite class of shared cancer-specific epitopes derived from recurrently mutated driver genes. Here we describe a high-throughput platform combining single-cell transcriptomic and T cell receptor (TCR) sequencing to establish whether mutant PIK3CA, among the most frequently genomically altered driver oncogenes, generates an immunogenic public NeoAg. Using this strategy, we developed a panel of TCRs that recognize an endogenously processed neopeptide encompassing a common PIK3CA hotspot mutation restricted by the prevalent human leukocyte antigen (HLA)-A*03:01 allele. Mechanistically, immunogenicity to this public NeoAg arises from enhanced neopeptide/HLA complex stability caused by a preferred HLA anchor substitution. Structural studies indicated that the HLA-bound neopeptide presents a comparatively ‘featureless’ surface dominated by the peptide’s backbone. To bind this epitope with high specificity and affinity, we discovered that a lead TCR clinical candidate engages the neopeptide through an extended interface facilitated by an unusually long CDR3β loop. In patients with diverse malignancies, we observed NeoAg clonal conservation and spontaneous immunogenicity to the neoepitope. Finally, adoptive transfer of TCR-engineered T cells led to tumor regression in vivo in mice bearing PIK3CA-mutant tumors but not wild-type PIK3CA tumors. Together, these findings establish the immunogenicity and therapeutic potential of a mutant PIK3CA-derived public NeoAg.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Structure of metallochaperone in complex with the cobalamin-binding domain of its target mutase provides insight into cofactor delivery

G-protein metallochaperone MeaB in bacteria [methylmalonic aciduria type A (MMAA) in humans] is responsible for facilitating the delivery of adenosylcobalamin (AdoCbl) to methylmalonyl-CoA mutase (MCM), the only AdoCbl-dependent enzyme in humans. Genetic defects in the switch III region of MMAA lead to the genetic disorder methylmalonic aciduria in which the body is unable to process certain lipids. Here, we present a crystal structure of Methylobacterium extorquens MeaB bound to a nonhydrolyzable guanosine triphosphate (GTP) analog guanosine-5'-[(β,γ)-methyleno]triphosphate (GMPPCP) with the Cbl-binding domain of its target mutase enzyme (MeMCM cbl ). This structure provides an explanation for the stimulation of the GTP hydrolyase activity of MeaB afforded by target protein binding. We find that upon MCM cbl association, one protomer of the MeaB dimer rotates ~180°, such that the inactive state of MeaB is converted to an active state in which the nucleotide substrate is now surrounded by catalytic residues. Importantly, it is the switch III region that undergoes the largest change, rearranging to make direct contacts with the terminal phosphate of GMPPCP. These structural data additionally provide insights into the molecular basis by which this metallochaperone contributes to AdoCbl delivery without directly binding the cofactor. Our data suggest a model in which GTP-bound MeaB stabilizes a conformation of MCM that is open for AdoCbl insertion, and GTP hydrolysis, as signaled by switch III residues, allows MCM to close and trap its cofactor. Substitutions of switch III residues destabilize the active state of MeaB through loss of protein:nucleotide and protein:protein interactions at the dimer interface, thus uncoupling GTP hydrolysis from AdoCbl delivery.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

In-plane Néel wall chirality and orientation of interfacial Dzyaloshinskii-Moriya vector in magnetic films

The interfacial Dzyaloshinskii-Moriya interaction (DMI) is of great interest as it can stabilize chiral spin structures in thin films. Experiments verifying the orientation of the interfacial DMI vector remain rare, in part due to the difficulty of separating vector components of DMI. In this study, Fe/Ni bilayers and Co/Ni multilayers were deposited epitaxially onto Cu(001) and Pt(111) substrates, respectively. By tailoring the effective anisotropy, spin reorientation transitions (SRTs) are employed to probe the orientation of the DMI vector by measuring the spin structure of domain walls on both sides of the SRTs. In our results, the interfacial DMI is found to be sufficiently strong to stabilize chiral Néel walls in the out-of-plane magnetized regimes, while achiral Néel walls are observed in the in-plane magnetized regimes. Furthermore, these findings experimentally confirm that the out-of-plane component of the DMI vector is insignificant in these fcc(001) and fcc(111) oriented interfaces, even in the presence of atomic steps.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

PID-Regulated Heating System for PIP-II Reference Line

The Proton Improvement Project-2 centers on building a new superconducting linear particle accelerator (Linac) at Fermilab. At the heart of the accelerator is the reference line, a critical system that defines the ideal path for the particle beam as it passes through magnets, RF cavities, and other beamline elements. Temperature stability is crucial for the reliable operation of RF components, such as mixers and filters. Fluctuations affect key performance parameters like conversion loss, isolation, and linearity. To mitigate any drift caused by ambient temperature changes, a heating plate assembly is utilized to maintain key components at a controlled temperature of 40°C. The system utilizes an aluminum 36”x36”x0.5” heat plate powered by a MOSFET-based control circuit, delivering approximately 460 W of thermal energy through a resistor array. Real-time temperature feedback is provided by a PT100 Resistance Temperature Detector (RTD), which interfaces with a Proportional–Integral–Derivative (PID) control algorithm to maintain closed-loop temperature regulation. The control signal actively modulates the gate voltage of an N channel MOSFET, dynamically adjusting power delivery in response to deviations from the temperature setpoint. Simulations and LTspice models validate the functionality and responsiveness of the circuit under varying conditions. The prototype has successfully demonstrated stable thermal control, paving the way for integration into the PIP-II infrastructure. The final design will feature an expanded resistor array, as well as communication with a PLC for continuous data acquisition and diagnostics. This work directly supports Fermilab’s broader mission by contributing to the stability and reliability of core accelerator systems, enhancing the precision of particle beam delivery for future physics experiments.

Mosher, Alexander [Fermilab]↗
DOE OSTI2024

Design of Star‐Shaped Trimer Acceptors for High‐Performance (Efficiency > 19%), Photostable, and Mechanically Robust Organic Solar Cells

Abstract High power conversion efficiency (PCE), long‐term stability, and mechanical robustness are prerequisites for the commercial applications of organic solar cells (OSCs). In this study, a new star‐shaped trimer acceptor (TYT‐S) is developed and high‐performance OSCs with a PCE of 19.0%, high photo‐stability (t 80% lifetime = 2600 h under 1‐sun illumination), and mechanical robustness with a crack‐onset strain (COS) of 21.6% are achieved. The isotropic molecular structure of TYT‐S affords efficient multi‐directional charge transport and high electron mobility. Furthermore, its amorphous structure prevents the formation of brittle crystal‐to‐crystal interfaces, significantly enhancing the mechanical properties of the OSC. As a result, the TYT‐S‐based OSCs demonstrate a significantly higher PCE (19.0%) and stretchability (COS = 21.6%) than the linear‐shaped trimer acceptor (TYT‐L)‐based OSCs (PCE = 17.5% and COS = 6.4%) and the small‐molecule acceptor (MYT)‐based OSCs (PCE = 16.5% and COS = 1.3%). In addition, the increased molecular size of TYT‐S, relative to that of MYT and dimer (DYT), suppresses the diffusion kinetics of the acceptor molecules, substantially improving the photostability of the OSCs. Finally, to effectively demonstrate the potential of TYT‐S, intrinsically stretchable (IS)‐OSCs are constructed. The TYT‐S‐based IS‐OSCs exhibit high device stretchability (strain at PCE 80% = 31%) and PCE of 14.4%.

Chemistry↗
DOE OSTI2023

Stabilizing Liquids Using Interfacial Supramolecular Assemblies

Stabilizing liquids based on supramolecular assembly (non-covalent intermolecular interactions) has attracted significant interest, due to the increasing demand for soft, liquid-based devices where the shape of the liquid is far from the equilibrium spherical shape. The components comprising these interfacial assemblies must have sufficient binding energies to the interface to prevent their ejection from the interface when the assemblies are compressed. In this report we highlight recent advances in structuring liquids based on non-covalent intermolecular interactions. We describe some of the progress made that reveals structure–property relationships. In addition to treating advances, we discuss some of the limitations and provide a perspective on future directions to inspire further studies on structured liquids based on supramolecular assembly.

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

Bipolar Redox‐Active Molecules in Non‐Aqueous Organic Redox Flow Batteries: Status and Challenges

Abstract In this Review, recent advances in the development of organic bipolar redox‐active molecules (BRMs) for redox‐flow batteries (RFBs) are discussed, with special emphasis on their utilization in non‐aqueous systems. BRMs are single molecules that can be used as both anolyte and catholyte. Their unique redox chemistry enables identical components to be present in RFBs, which has distinct advantages in storage, mitigating crossover and compound degradation. Recent scientific innovations have uncovered three main strategies in preparing these bipolar species: developing stable radicals that are capable of both positive and negative redox reactions, covalently bridging electron donor‐acceptor motifs on one single molecule, and forming biredox eutectic mixtures through physically mixing active components. These approaches not only enabled symmetric flow battery design, but also resulted in electroactive materials with improved properties such as solubility, stability, and cell voltage. Nevertheless, much is still to be done. Importantly, developing BRMs that are competitive with vanadium RFBs in all performance indicators (e. g. solubility >1.0 M, voltage output >2.0 V, and cycling longevity >10 000 cycles) remains a challenge. Mechanistic studies including molecular dynamics, plausible degradation processes, and electrode/solution interface are necessary to peer into key parameters that dominate their battery performance. Accordingly, we review the current status and suggest future opportunities in these areas.

Li, Min↗
DOE OSTI2021

Mechanistic understanding of support effect on the activity and selectivity of indium oxide catalysts for CO 2 hydrogenation

Herein we present a mechanistic study on the support effect (ZrO 2 and CeO 2 ) of In 2 O 3 catalysts in CO 2 hydrogenation by a combined experimental and computational approach. Kinetic experiments and surface characterization suggested that the activity of In 2 O 3 catalysts cannot be simply correlated with the abundance of surface oxygen vacancies (O v ) formed by either H 2 -reduction or thermal treatment, which has been frequently invoked in previous studies. The support effect should originate from the electronic interactions between In 2 O 3 and the support oxide, rather than geometric factors or the difference in the particle size of In 2 O 3 . Theoretical modelling revealed that surface O v facilitate the formation and stabilization of the formate (HCOO*) intermediate. While a carbonate-like structure is favored for CO 2 adsorption on CeO 2 -supported or unsupported In 2 O 3 catalysts, CO 2 tends to bind strongly in a bent configuration on the O v site at the In 2 O 3 -ZrO 2 interface. The distinct CO 2 adsorption structures on different supported In 2 O 3 catalysts may account for the different reaction energy profiles in the subsequent hydrogenation reactions, especially the rate-limiting step, i.e., hydrogenation of HCOO* to CH 2 O* and methoxy (CH 3 O*). The relatively higher methanol selectivity of In 2 O 3 catalyst supported on ZrO 2 with respect to that on CeO 2 are suggested to stem from the greater energy difference (Δ$E_a$) between the parallel hydrogenation and C-O bond cleavage of HCOO*, which leads to the formation of methanol and CO, respectively. We report this study underlines the important role of metal-oxide-interface in determining the catalytic behavior of oxide-supported In 2 O 3 catalysts in CO 2 conversion.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

Development of intelligent control module for the J-TEXT electron cyclotron emission imaging system

An Electron Cyclotron Emission Imaging (ECEI) system has been developed on the J-TEXT tokamak and employed for the multi-dimensional electron temperature fluctuation imaging measurement in 2019 campaign. The intelligent control module has been developed and been used to flexibly setup the 256-channel ECEI system configurations for different plasma scenarios. The adjustments of key parameters are merged in this module including size and position of the geometric observation window, and the dynamic range of the signal and image resolution. The user-friendly interface is also developed for the tokamak operator to set the ECEI system focusing on physics interested region. Furthermore, this module is installed close to the diagnostic system around the J-TEXT tokamak, which can be remotely controlled by authorized operators in the control room. It significantly improves ECEI flexibility and stability, and also keeps the operation technician away from the harsh radiation environment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

An analysis of the spatio-temporal resolution of the immersed boundary method with direct forcing

The immersed boundary method (IBM) with direct forcing is very popular in the simulation of rigid particulate flows. In the IBM, an interaction force is introduced at the interface between fluid and particle in order to approximate the no-slip boundary condition. The interaction force is calculated through dividing the velocity difference (or error) between fluid and particle at the interface by the time step. Here, a dynamic equation for the velocity difference is derived. Additionally, analyses on the dynamic equation provide a few new findings: (i) The interaction force is the solution of a least-squares error problem, with the direct implication that the Lagrangian marker distribution has no effect on the large scale flow structure once the distribution of Lagrangian markers become saturated along the interface (i.e., each marker remains properly correlated with all its neighbors); (ii) The Lagrangian volume-weight is a relaxation factor to control how fast the velocity error decays to the ideal value of zero; (iii) The optimal choice of the Lagrangian volume-weight is the largest value permissible by a stability condition. A comprehensive convergence analysis with regard to the spatial and temporal resolution is presented for the velocity error and also for the shear-stress and surface pressure. In three simple canonical problems, it is analytically and numerically shown that the IBM results converge to the theoretical solutions obtained with precise imposition of no-slip and no-penetration boundary conditions. It is observed that it is not necessary to match the Lagrangian marker volume-weight to that of the local Eulerian cell volume and in fact this matching leads to lower than optimal computational efficiency. However, it is found that extremely high Eulerian grid resolution and small time step have to be used to obtain high precision simulation results. Especially, the time step should be inversely proportional to the particle Reynolds number for low Reynolds number flows. For high frequency oscillation problems, the grid size needs to be reduced by a factor of the square root of the frequency, and the time step to be reduced by a factor of the frequency. The theoretical findings here can be used to alleviate the technical difficulties in simulating non-spherical particles by not requiring the Lagrangian marker distribution to match the Eulerian grids and also in the implementation of IBM on non-uniform Eulerian grids. The present work also provides simple practical guidance on the choice of temporal and spatial resolution so as to control the simulation error a priori.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗