Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “moving boundary”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A weighted shifted boundary method for immersed moving boundary simulations of Stokes' flow

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. The surrogate domain is constructed so as to avoid cut cells and the associated problematic implementation and numerical integration issues. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions: hence the name of the method, that shifts the location and values of the boundary conditions. Here, in this article, we extend the SBM to the simulation of incompressible Stokes flow, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach allows to drastically reduce spurious pressure oscillations in time, which are produced if the total volume of active fluid were to change abruptly over a time step. The proposed Weighted SBM (W-SBM) exactly preserves states of hydrostatic equilibrium, and induces small mass and momentum conservation errors, which converge as the grid is refined. This is in analogy to cutFEMs and related unfitted approaches, which rely on an affine representation of cut boundaries. We demonstrate the robustness and accuracy of the proposed method with an extensive suite of two-dimensional tests.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Instability and treatments of the coupled discrete element and lattice Boltzmann method by the immersed moving boundary scheme

The immersed moving boundary (IMB) scheme has been extensively used to couple the discrete element method (DEM) with the lattice Boltzmann method (LBM). In the literature, only the formulation of IMB for lattice nodal cells covered by a single-solid particle was given. The treatment of situations where a nodal cell is covered by two or more solid particles is seldom discussed. It is found that some numerical instability can occur for such situations due to an inappropriate computation of the weighting function in the IMB formulation. This work presents an enhanced treatment that can resolve the issue and validates it using some benchmark tests. Furthermore, to avoid the extra costs associated with the treatment and simplify the complicated procedure introduced, a simplified IMB scheme is proposed. The accuracy of both enhanced and simplified IMB schemes are validated by test cases including single-particle sedimentation, two-particle drafting-kissing-tumbling phenomenon, and multiple-particle sedimentation. Then, the robustness of both schemes is examined and discussed using a specially designed flow past cylinders test. Overall, the simplified IMB scheme is proved to be robust and sufficiently accurate and simpler and more effective than the enhanced scheme.

42 ENGINEERING↗
DOE OSTI2017

On mechanical waves and Doppler shifts from moving boundaries

We investigate the propagation of infinitesimal harmonic mechanical waves emitted from a boundary with variable velocity and arriving at a stationary observer. In the classical Doppler effect, X s (t)=vt is the location of the source with constant velocity v. In the present work, however, we consider a source co-located with a moving boundary x=X s (t), where X s (t) can have an arbitrary functional form. For ‘slowly moving’ boundaries (i.e., ones for which the timescale set by the mechanical motion is large in comparison to the inverse of the frequency of the emitted wave), we present a multiple-scale asymptotic analysis of the moving boundary problem for the linear wave equation. Here, we obtain a closed-form leading-order (with respect to the latter small parameter) solution and show that the variable velocity of the boundary results not only in frequency modulation but also in amplitude modulation of the received signal. Consequently, our results extend the applicability of two basic tenets of the theory of a moving source on a stationary domain, specifically that (i) $.\atop{x}_s$ for non-uniform boundary motion can be inserted in place of the constant velocity v in the classical Doppler formula and (ii) that the non-uniform boundary motion introduces variability in the amplitude of the wave. The specific examples of decelerating and oscillatory boundary motion are worked out and illustrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

A fourth order sharp immersed method for the incompressible Navier-Stokes equations with stationary and moving boundaries and interfaces

We propose a fourth order Navier-Stokes solver based on the immersed interface method (IIM), for flow problems with stationary and one-way coupled moving boundaries and interfaces. Our algorithm employs a Runge-Kutta-based projection method that maintains high-order temporal accuracy in both velocity and pressure for steady and unsteady velocity boundary conditions. Fourth order spatial accuracy is achieved through a novel fifth order IIM discretization scheme for the advection term, as well as existing high-order interface-corrected finite difference schemes for the other differential operators. Using a set of manufactured flow problems with stationary and moving boundaries, we demonstrate fourth order convergence of velocity and pressure in the infinity norm, both inside the domain and on the immersed boundaries. The solver’s performance is further validated through a range of practical flow simulations, highlighting its efficiency over a second order scheme. Finally, we showcase the ability of our immersed discretization scheme to handle interface-coupled multiphysics problems by solving a conjugate heat transfer problem with multiple immersed solids. Overall, the proposed approach robustly combines the efficiency of high order discretization schemes with the flexibility of immersed discretizations for flow problems with complex, moving boundaries and interfaces.

42 ENGINEERING↗
DOE OSTI2019

Distributary Channel Networks as Moving Boundaries: Causes and Morphodynamic Effects

In this study, we propose an exploratory model to describe the morphodynamics of distributary channel network growth on river deltas. The interface between deep channels and the shallow, unchannelized delta front deposits is modeled as a moving boundary. Steady flow over the unchannelized delta front is friction dominated and modeled by Laplace's equation. Shear stress along the network boundary produces nonlinear erosion rates at the interface, causing the boundary to move and network elements (channels and branches) to form. The model was run for boundary conditions resembling the Wax Lake Delta in coastal Louisiana, 20 parameterizations of sediment transport, and 3 parameterizations of discharge. In each case, the model produced a complex channel network with channel number, width, bifurcation angle, and channel shape depending on the sediment transport formula. For reasonable sediment transport parameters and gradually increasing water discharge, the model produced network characteristics and progradation rates similar to theWax Lake Delta. Lastly, this suggests that the model contains the processes responsible for network growth, despite its abstract formulation.

58 GEOSCIENCES↗
DOE OSTI2018

A novel algorithm of immersed moving boundary scheme for fluid-particle interactions in DEMLBM

This paper presents an efficient and accurate Immersed Moving Boundary (IMB) algorithm for solving fluid-particle interactions in the framework of the Lattice Boltzmann Method (LBM). Although the IMB scheme has been widely employed in many fluid-particle coupling problems in a wide range of applications, the algorithm of its implementation, especially in identifying both fluid and solid boundary nodes for particles, is seldom reported. Besides, the computational cost of handling fluid-particle coupling is very expensive. To provide a bridge between theory and application and improve the computing efficiency of IMB, a novel boundary tracing procedure and an efficient method for computing the solid nodal ratio using Gaussian quadrature are proposed in this paper. Both accuracy and efficiency of the proposed algorithm are examined by two benchmark tests. It is also found that the IMB scheme are more efficient and stable compared to another widely used the Immersed Boundary Method (IBM) in LBM.

42 ENGINEERING↗
DOE OSTI2022

The Importance of a Moving Boundary Approach for Modeling the SEI Layer Growth to Predict Capacity Fade

One of the contributing factors to the aging of lithium-ion batteries is the growth of the solid-electrolyte interphase (SEI) layer. The growth of the SEI layer leads to the irreversible loss of lithium available for cycling and increases the resistance of the battery. Physics-based models in literature model the kinetically limited or solvent diffusion-limited growth. In such models, the interface resistance is a constant, and the contribution to the overpotential of the intercalation reaction from the SEI layer is considered to be ohmic. In this study, we propose a model that describes the growth of the SEI layer on the electrode surface as a moving interface. The transport of lithium ions and the solvent in the electrolyte are affected by this moving interface. The equations that govern the species transport and the potential drop across the SEI layer are derived from dilute solution theory and solved by transforming the coordinates of the moving boundary. The ion transport induces changes in the conductivity across the SEI layer, which affects the potential drop that arises due to its growth. The effects of this potential on capacity fade are studied over cycling the battery.

25 ENERGY STORAGE↗
DOE OSTI2016

Gradient elution moving boundary electrophoresis enables rapid analysis of acids in complex biomass-derived streams

Biomass conversion processes such as pretreatment, liquefaction, and pyrolysis often produce complex mixtures of intermediates that are a substantial challenge to analyze rapidly and reliably. To characterize these streams more comprehensively and efficiently, new techniques are needed to track species through biomass deconstruction and conversion processes. Here, we present the application of an emerging analytical method, gradient elution moving boundary electrophoresis (GEMBE), to quantify a suite of acids in a complex, biomass-derived streams from alkaline pretreatment of corn stover. GEMBE offers distinct advantages over common chromatography-spectrometry analytical approaches in terms of analysis time, sample preparation requirements, and cost of equipment. As demonstrated here, GEMBE is able to track 17 distinct compounds (oxalate, formate, succinate, malate, acetate, glycolate, protocatechuate, 3-hydroxypropanoate, lactate, glycerate, 2-hydroxybutanoate, 4-hydroxybenzoate, vanillate, p-coumarate, ferulate, sinapate, and acetovanillone). The lower limit of detection was compound dependent and ranged between 0.9 and 3.5 umol/L. Results from GEMBE were similar to recent results from an orthogonal method based on GCxGC-TOF/MS. Altogether, GEMBE offers a rapid, robust approach to analyze complex biomass-derived samples, and given the ease and convenience of deployment, may offer an analytical solution for online tracking of multiple types of biomass streams.

09 BIOMASS FUELS↗
DOE OSTI2022

A Moving Embedded Boundary Approach for the Compressible Navier-Stokes Equations in a Block-Structured Adaptive Refinement Framework

A computational technique has been developed to perform compressible flow simulations involving moving boundaries using an embedded boundary approach within the block-structured adaptive mesh refinement (SAMR) framework of AMReX [1], [91], [92]. We leverage the SAMR capability to obtain quantitatively accurate results whilst using robust, second-order finite volume schemes. A conservative, unsplit, cut-cell approach is utilized and a ghost-cell approach is developed for computing the flux on the moving, embedded boundary faces. A third-order least-squares formulation has been developed to compute the wall velocity gradients, and was found to significantly improve the performance of the solver in terms of the quantitative comparison of surface quantities such as the skin friction coefficient. Various test cases are performed to validate the method, and compared with analytical, experimental, and other numerical results in literature. Inviscid and viscous test cases are performed that span a wide regime of flow speeds - acoustic (harmonically pulsating sphere), smooth flows (expansion fan created by a receding piston) and flows with shocks (shock-cylinder interaction, shock-wedge interaction, pitching NACA 0012 airfoil and shock-cone interaction). A closed system with moving boundaries - an oscillating piston in a cylinder, showed that the percentage error in mass within the system decreases with refinement, demonstrating that the numerical scheme is conservative with grid refinement, but is not discretely conservative. Viscous test cases involve that of a horizontally moving cylinder at Re = 40, an inline oscillating cylinder at Re = 100, and a transversely oscillating cylinder at Re = 185. The judicious use of adaptive mesh refinement with appropriate refinement criteria to capture the regions of interest leads to well-resolved flow features, and good quantitative comparison is observed with the results available in literature.

adaptive refinement↗
DOE OSTI2017

Numerical Solution of Moving Phase Boundary and Diffusion-Induced Stress of Sn Anode in the Lithium-Ion Battery

Here, we have previously observed a large transient stress in Sn film anodes at the beginning of the Sn-Li 2 Sn 5 phase transformation. To understand this behavior, we use numerical modeling to simulate the kinetics of the 1-D moving boundary and Li diffusion in the Sn anodes. A mixture of diffusion-controlled and interface-controlled kinetics is found. The Li concentration in the Li 2 Sn 5 phase remains near a steady-state profile as the phase boundary propagates, whereas the Li diffusion in Sn is more complicated. Li continuously diffuses into the Sn layer and produces a supersaturation; the Li can then diffuse toward the Sn/Li 2 Sn 5 interface and contribute to further phase transformation. Finally, the evolution of Li concentration in the Sn induces strain which involves rate-dependent plasticity and elastic unloading, resulting in the complex stress evolution that is observed. In the long term, the measured stress is dominated by the stress in the growing Li 2 Sn 5 phase.

25 ENERGY STORAGE↗
DOE OSTI2020

Extended Finite Element Based Approach in Additive Manufacturing Modeling for Optimizing Highly Complex Manifold in Protonic Ceramic Electrochemical Cells

The objective of this project is to develop a novel numerical approach based on the extended finite element method (XFEM) for modeling moving boundaries of material deposited in additive manufacturing (AM) process with improved accuracy and reduced computational cost. One of the major challenges in simulating AM processes is that the boundaries of the computational domain need to change as material is deposited. Previous approaches have achieved this by activating new finite elements, but this requires a high level of mesh refinement to capture small movement of the boundary. XFEM allows the solution boundary to move independently of the mesh, permitting smooth representation of the evolution of the boundary as material is deposited. As a result, this new method will significantly improve the simulation accuracy while reducing the computational cost compared with existing approaches. Furthermore, this project will bring significant improvements for the design and AM process optimization in an iterative fashion among numerical simuation, parameter optimization, and experimental validation. Specifically, this project will support the protonic ceramic electrochemical cells (PCEC) stack development at INL by providing insights into the density, residual stress, thermal and mechanical properties of the PCEC manifold and interconnect, which is critical in enhancing the system lifetime and reducing the overall cost. Upon the success of proposed development and validation of the proposed approach, simulation will be applied to determine an optimal set of AM process parameters for the PCEC manifold production. This project will consolidate Idaho National Laboratory (INL)’s simulation capabilities provided by the MOOSE framework and the Valhalla AM simulation application, with encouraging expansion to emerging PCEC applications and AM technologies at INL and beyond.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Probing the decelerating trajectory of a Raman soliton using temporal reflection

Temporal reflection is a process where an optical pulse reflects off a moving boundary with different refractive indices across it. In a dispersive medium, this process creates a reflected pulse with a frequency shift that changes its speed. Such frequency shifts depend on the speed of the moving boundary. In this work, we propose and experimentally show that it is possible to probe the trajectory of the boundary by measuring the frequency shifts while changing the initial delay between the incident pulse and the boundary. We demonstrate this effect by reflecting a probe pulse off a short soliton, acting as a moving boundary that decelerates inside a photonic crystal fiber because of intrapulse Raman scattering. We deduce trajectory of the soliton from the measured spectral data for the reflected pulse.

47 OTHER INSTRUMENTATION↗
DOE OSTI2019

On the implicit immersed boundary method in coupled discrete element and lattice Boltzmann method

The coupled discrete element method and lattice Boltzmann method (DEMLBM) has increasingly drawn attention of researchers in geomechanics due to its mesoscopic nature since 2000. Immersed boundary method (IBM) and immersed moving boundary (IMB) are two popular schemes for coupling fluid particle in DEMLBM. This work aims at coupling DEM and LBM using the latest IBM algorithm and investigating its accuracy, computational efficiency, and applicability. Two benchmark tests, interstitial fluid flow in an ideal packing and single particle sedimentation in viscous fluid, are carried out to demonstrate the accuracy of IBM through semi-empirical Ergun equation, finite element method (FEM), and IMB. Then, simulations of particle migration with relatively large velocity in Poiseuille flow are utilized to address limitations of IBM in DEMLBM modeling. In addition, advantages and deficiencies of IBM are discussed and compared with IMB. It is found that the accuracy of IBM can be only guaranteed when sufficient boundary points are used and it is not suitable for geomechanical problems involving large fluid or particle velocity.

58 GEOSCIENCES↗
DOE OSTI2022

A coupled DEM-IMB-LBM model for simulating methane hydrate exploitation involving particle dissolution

The coupled discrete element and lattice Boltzmann method using an immersed moving boundary scheme was extended to simulate methane hydrate exploitation involving mass transport and particle dissolution. In this coupled DEM-IMB-LBM model, a new Dirichlet-type thermal boundary condition is extended to simulate moving curved boundaries with constant concentration. A novel periodic boundary including an efficient searching algorithm for particle contact is proposed to reduce the computational cost and boundary effect. So, this model is validated by two numerical examples: a circular particle with concentration convection-diffusion moving in a horizontal channel and mass transport from a cylinder particle in a simple shear flow. The numerical results obtained from the proposed model agree well with previous studies. To further demonstrate the capacity of the proposed model, simulations of methane hydrate exploitation including two formations in marine sediments are carried out. The numerical results indicate that the coupled DEM-IMB-LBM is not only capable of simulating the dissolution of hydrate particles at the grain level, but also recover the sand erosion and migration process in a fundamental perspective during the methane hydrate exploitation process.

42 ENGINEERING↗
DOE OSTI2019

A posteriori error analysis of fluid–structure interactions: Time dependent error

A posteriori error analysis is a technique to quantify the error in particular simulations of a numerical approximation method. In this article, we use such an approach to analyze how various error components propagate in certain moving boundary problems. We study quasi-steady state simulations where slowly moving boundaries remain in mechanical equilibrium with a surrounding fluid. Such problems can be numerically approximated with the Method of Regularized Stokelets (MRS), a popular method used for studying viscous fluid–structure interactions, especially in biological applications. Our approach to monitoring the regularization error of the MRS is novel, along with the derivation of linearized adjoint equations to the governing equations of the MRS with a elastic elements. Our main numerical results provide a clear illustration of how the error evolves over time in several MRS simulations.

97 MATHEMATICS AND COMPUTING↗