Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Immersed 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 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 OSTI2025

A high-order computational framework for particle-resolved simulations of disperse multiphase flows

This work presents a high-order numerical approach for particle-resolved simulations of disperse multiphase flows, where the Navier-Stokes equations for fluid flow are solved using a high-order spectral element method in the Eulerian framework, and the particle phase is directly simulated with a discrete element method. The coupling between particles and fluids is explicitly handled using an adapted direct-forcing immersed boundary method. Unlike the conventional schemes, a high-order barycentric Lagrange interpolation method and a Gaussian projection kernel are used to ensure accurate momentum exchange between local boundary points and surrounding fluid nodes in the framework of high-order fluid solver. Benchmark tests of increasing complexity are conducted to demonstrate the accuracy and efficiency of our method. Here, it is found that our approach exhibits an excellent convergence performance, as the fluid element/grid is refined and the number of boundary points increases. Compared to conventional low-order methods, the proposed high-order framework enables the use of substantially larger fluid elements while maintaining high accuracy in modeling fluid-particle interactions, owing to the enhanced resolution of high-order basis functions. Moreover, since the primary unknowns are stored at element or grid nodes, the high-order approach offers improved efficiency in both CPU memory usage and total computational cost.

42 ENGINEERING↗
DOE OSTI2025

Explicit block encodings of boundary value problems for many-body elliptic operators

Simulation of physical systems is one of the most promising use cases of future digital quantum computers. In this work we systematically analyze the quantum circuit complexities of block encoding the discretized elliptic operators that arise extensively in numerical simulations for partial differential equations, including high-dimensional instances for many-body simulations. When restricted to rectangular domains with separable boundary conditions, we provide explicit circuits to block encode the many-body Laplacian with separable periodic, Dirichlet, Neumann, and Robin boundary conditions, using standard discretization techniques from low-order finite difference methods. To obtain high-precision, we introduce a scheme based on periodic extensions to solve Dirichlet and Neumann boundary value problems using a high-order finite difference method, with only a constant increase in total circuit depth and subnormalization factor. We then present a scheme to implement block encodings of differential operators acting on more arbitrary domains, inspired by Cartesian immersed boundary methods. We then block encode the many-body convective operator, which describes interacting particles experiencing a force generated by a pair-wise potential given as an inverse power law of the interparticle distance. This work provides concrete recipes that are readily translated into quantum circuits, with depth logarithmic in the total Hilbert space dimension, that block encode operators arising broadly in applications involving the quantum simulation of quantum and classical many-body mechanics.

Kharazi, Tyler [University of California, Berkeley↗
DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗
DOE OSTI2026

Two-phase flow numerical analysis of electrode geometry for alkaline water electrolyzers

Hydrogen is a promising component of a future energy-secure and efficient economy, but its competitiveness depends on reducing production costs. One strategy is to operate alkaline water electrolyzers at higher current densities to increase output. However, this intensifies performance losses due to gas bubble accumulation, which blocks transport pathways and deactivates electrochemically active surfaces. Enhancing bubble evacuation through electrode design is therefore essential. Previous studies have explored various approaches — such as modifying surface morphology, applying sonication or pressure modulation, and introducing surfactants — but these efforts have addressed a limited range of conditions due to the complexity of two-phase flow and electrode geometries. Experiments have also largely been focused on either cell level improvements, which lack the information necessary to isolate each contributing factor, or on modified geometries that are not relevant to practical cell operation. From a modeling perspective, conventional Eulerian multiphase models do not track the complex gas–liquid interfacial dynamics and often neglect surface tension and contact angle effects, reducing their predictive accuracy. To provide insights on the effects of different electrode geometries on the performance of alklaine water electrolyzers this work employs an immersed boundary volume-of-fluid method to simulate bubble behavior in 3D porous electrodes. Multiple base electrode geometries, typically used in practice, with varying porosity are evaluated under a constant surface gas generation rate. Simulation data is analyzed to quantify electrode gas coverage, bubble size dynamics and other relevant metrics. Results show that porosity strongly influences bubble accumulation on electrode surfaces, with higher porosity reducing gas coverage, and its not strictly dependent on the electrode geometry. However, the electrode’s base geometry significantly affects gas accumulation at the separator gap, independent of porosity. A foam electrode geometry resulted in the lowest gas coverage of all electrodes with a median volumetric gas coverage of 11%, but at the cost of a 70% reduction in active area compared with the largest surface area electrode, while gyroid electrodes showed the best trade-off between gas coverage, particularly at the separator surface, and electrochemically active area. In conclusion, the results highlight the need for holistic electrode design strategies.

Alkaline↗
DOE OSTI2025

Numerical simulation of frost formation and heat transfer on fin-and-tube heat exchangers in turbulent cross-flow

Frost formation in fin-and-tube heat exchangers in turbulent cross-flow presents significant challenges in industrial refrigeration applications, affecting heat transfer efficiency and operational reliability. The purpose of this work is to investigate frost deposition and growth on a staggered bank of a fin-and-tube freezer coil under turbulent forced convection conditions. The focus here is on investigating conditions that closely replicate real-world scenarios in large walk-in industrial freezers. Using a direct numerical simulation approach, we examine the flow dynamics and thermal behaviour in the presence of frost, considering turbulent regimes characterized by a Reynolds number in the range 1050 ≤ R e D , avg ≤ 4800 , with the characteristic length being the outer diameter of the tube and the velocity being the bulk fluid velocity between the plates (fins). Computational fluid dynamics simulations are employed to resolve the interactions between turbulent airflow and the frost layer. Our approach incorporates a modified immersed boundary method and a slow-time acceleration technique to address the complex dynamic interface between the continuously evolving frost layer and the flowing air stream. Our findings indicate that frost forms more on the sides of the finned surfaces (plates) and less on the tubes themselves. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics↗
DOE OSTI2025

Q2 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

This report describes the work and activities carried out towards the completion of each of the following milestones in FY25 Q2: 1. Demonstrate workflow for generating self-consistent CESOL plasma profiles + first wall and divertor loading prediction and generate the CAT plasma and neutron loading needed for further engineering analysis: $\circ$ Run CESOL with BOUT++/Hermes-3 and immersed boundary condition to directly map to wall: • Run BOUT++/Hermes-3 through the IPS workflow to find radial particle and energy diffusivities to match either the Eich or the physics-based scaling of the SOL heat flux width, and • Expand source of first wall heat flux to include charged particles, neutrals, and radiation from the core+edge. 2. Generate medium fidelity parametrized CAD: $\circ$ Develop the TRACER tool to read an existing CAD, regenerate the geometry based on vertex location and connectivity information, define vertex translation and parameters needed for scaling the CAD, and $\circ$ Utilize the FreeGS code to determine CAT PF coil placement, including minimizing the number of coils, coil current, and electromechanical stresses. 3. Utilize plasma loading for engineering analysis: $\circ$ Couple the plasma loading to input for OpenFOAM and demonstrate initial test of thermal analysis of CAT first wall loading with typical DCLL blanket component cooling boundary conditions. 4. Demonstrate nuclear analysis: $\circ$ Apply initial analysis of tritium transport in DCLL blanket by evaluating spatially resolved tritium generation rates, tritium diffusion and convection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2025

Analysis of the weighted shifted boundary method for the Poisson and Stokes problems

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. Accuracy is maintained by properly shifting the location and values of the boundary conditions. This avoids integration over cut cells and the associated implementation issues. Recently, the Weighted SBM (WSBM) was proposed for the Navier-Stokes equations with free surfaces and the Stokes flow with moving boundaries. The attribute “weighted” in the name WSBM stems from the fact that its variational form is weighted with the elemental volume fraction of active fluid. The motivation for the development of the WSBM was the preservation of the volume of active fluid to a higher degree of accuracy, which in turn resulted in improved stability and robustness characteristics in moving-boundary, time-dependent simulations. In this article, we present the numerical analysis of the WSBM formulations for the Poisson and Stokes problems. We give mathematical conditions under which the bilinear forms defining the discrete variational formulations are uniformly coercive (Poisson problem) or inf-sup stable (Stokes problem). By these results, stability and optimal convergence is proven in the natural norm; L2-error estimates can also be derived.

Approximate domain boundaries↗
DOE OSTI2024

Adaptive immersed isogeometric level-set topology optimization

Here, this paper presents for the first time an adaptive immersed approach for level-set topology optimization using higher-order truncated hierarchical B-spline discretizations for design and state variable fields. Boundaries and interfaces are represented implicitly by the iso-contour of one or multiple level-set functions. An immersed finite element method, the eXtended IsoGeometric Analysis, is used to predict the physical response. The proposed optimization framework affords different adaptively refined higher-order B-spline discretizations for individual design and state variable fields. The increased continuity of higher-order B-spline discretizations together with local refinement enables direct control over the accuracy of the representation of each field while simultaneously reducing computational cost compared to uniformly refined discretizations. A flexible mesh adaptation strategy enables local refinement based on geometric measures or physics-based error indicators. These adaptive discretization and analysis approaches are integrated into gradient-based optimization schemes, evaluating the design sensitivities using the adjoint method. Numerical studies illustrate the features of the proposed framework with static, linear elastic, multi-material, two- and three-dimensional problems. The examples provide insight into the effect of refining the design variable field on the optimization result and the convergence rate of the optimization process. Using coarse higher-order B-spline discretizations for level-set fields promotes the development of smooth designs and suppresses the emergence of small features. Moreover, adaptive mesh refinement for state variable fields results in a reduction of overall computational cost. Higher-order B-spline discretizations are especially interesting when evaluating gradients of state variable fields due to their higher inter-element continuity.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Corrosion evaluation of Al-Cu-Mn-Zr cast alloys in 3.5% NaCl solution

Corrosion behavior of cast Al-Cu-Mn-Zr (ACMZ) and RR350 alloys was compared to a cast 319 alloy in 3.5 wt.% NaCl. After 168 h immersion, ACMZ and RR350 alloys suffered from preferential attack adjacent to intermetallic particles decorated at grain boundaries while the attack in 319 occurred in eutectic Al-Si dendritic boundaries. Electrochemical data allowed semiquantitative comparison of alloy resistance to corrosion initiation, and ACMZ type alloys, including RR350 and three alloys with higher Cu, were considered more resistant than 319 due to the absence of deleterious Si particles. In case of 319, such Si particles presumably drove higher micro-galvanic influence to initiate and sustain Al corrosion. With lower susceptibility to corrosion initiation, ACMZ alloys should exhibit higher or at minimum similar resistance compared to cast 319.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Variationally consistent Maxwell stress in flexoelectric structures under finite deformation and immersed in free space

Maxwell stress refers to the mechanical stress exerted on a dielectric material due to the presence of electric fields. It plays a significant role in the interaction between a dielectric material and the surrounding free space under finite deformation. Previous research on finite deformation of flexoelectricity mainly adopted a modified form of Maxwell stress, potentially not able to correctly capture some physical phenomena, such as the compression of a dielectric droplet in an electric field. In this work, we propose a consistent and complete variational principle for flexoelectricity, in which the Maxwell stress emerges naturally from the derivation, without introducing additional assumptions. An Isogeometric analysis-based numerical framework is developed accordingly and verified by both linear and nonlinear benchmark cases compared with experimental results. The present framework successfully captures and quantifies the behaviors of conductive liquids and soft dielectric solids subjected to an external electric field. Finally, a novel scenario is investigated in which a flexoelectric beam immersed in free space is analyzed, showing the interesting distribution of Maxwell stress-induced tractions at opposing boundaries. The test demonstrates that a higher dielectric constant can effectively enhance the material's stiffness in response to the external electric loading.

36 MATERIALS SCIENCE↗
DOE OSTI2025

On the motion of compact objects in relativistic viscous fluids

We present a world-line effective field theory of compact objects moving relativistically through a viscous fluid. The theory is valid when velocity gradients are small compared to the inverse size of the object. Working within the EFT eliminates the need to solve a boundary value problem by turning all interactions between the fluid and the object into a source term in the action. We use the EFT to derive the relativistic equations of motion for a compact object immersed in a viscous fluid in a curved background, when the relative velocity of the object and the fluid is small compared to the speed of light.

astrophysical black holes↗
DOE OSTI2025

Spherical and sessile droplet dynamics by fluctuating hydrodynamics

We simulate the mesoscopic dynamics of droplets formed by phase-separated fluids at nanometer scales where thermal fluctuations are significant. Both spherical droplets fully immersed in a second fluid and sessile droplets which are also in contact with a solid surface are studied. Our model combines a Cahn–Hilliard formulation with incompressible fluctuating hydrodynamics; for sessile droplets, the fluid–solid contact angle is specified as a boundary condition. Deterministic simulations with an applied body force are used to measure the droplets' mobility from which a diffusion coefficient is obtained using the Einstein relation. Stochastic simulations are independently used to obtain a diffusion coefficient from a linear fit of the variance of a droplet's position with time. In some scenarios, these two measurements give the same value but not in the case of a spherical droplet initialized near a slip wall or in the case of sessile droplets with large contact angles (≥90°) on both slip and no-slip surfaces.

Bell, John B↗
DOE OSTI2025

Resistive hose modes in tokamak runaway electron beams II

Resistive hose instabilities of runaway electron (RE) beams immersed in resistive background plasmas are examined with analytic and numerical calculations. The RE beam-plasma equilibria considered are characteristic of the situation observed post-thermal quench in a tokamak disruption. An analytic linear dispersion relation is presented for the case of a uniform RE current density profile with a sharp boundary in cylindrical geometry. Initial value linear calculations for a more general profile in toroidal geometry find that reducing aspect ratio increases the resistive hose mode growth rates with fixed safety factor profile. Nonlinear calculations in cylindrical and toroidal geometry find that the resistive hose instability-driven fluctuations relax the gradient of the current density profile. In toroidal geometry, changes to the magnetic topology are observed as a result of the resistive hose activity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Irreversible thermodynamics of curved lipid membranes. II. Permeability and osmosis

We present a theory that combines the framework of irreversible thermodynamics with modified integral theorems to model arbitrarily curved and deforming membranes immersed in bulk fluid solutions. We study the coupling between the mechanics and permeability of a viscous and elastically bendable membrane, and a multicomponent bulk fluid solution. An equation for the internal entropy production for irreversibilities at the membrane is derived, determining the generalized thermodynamic forces and fluxes, from which we identify the deviatoric stress as a novel driving force for permeability. A complete set of equations of motion, constitutive laws, and boundary conditions to model the lipid membrane and bulk fluid system are provided.

Alkadri, Ahmad M↗
DOE OSTI2025

Underwater thermomagnetic generator for remote marine thermal energy harvesting and sensing

Thermomagnetic generators offer a promising approach for sustainable power generation in remote marine environments. Here, this study presents the design, thermal modeling, and experimental validation of a passively driven underwater thermomagnetic generator developed for powering ocean observation and monitoring sensors. The generator was evaluated under varying working fluids, thermal boundary conditions, and extended operation to assess real-world applicability. Two fluids, deionized water and silicone-based Thermal C5, were tested under simulated shallow- and deep-ocean conditions. Deionized water outperformed Thermal C5, especially in colder environments (~5°C), achieving a peak output of 2.7 mW due to larger temperature gradients and enhanced convective-evaporative heat transfer. Long-duration tests revealed a transient evaporation-condensation cycle that temporarily reduced rotor immersion and performance before stabilizing. The generator powered commercial marine sensors for over 6 h without external batteries, demonstrating the viability of compact, passively cooled thermomagnetic systems for autonomous, off-grid marine sensing.

13 HYDRO ENERGY↗
DOE OSTI2025

Unavoidable multilevel biaxial symmetry breaking in chiral hybrid liquid crystals

Chiral nematic or cholesteric liquid crystals (LCs) are fluid mesophases with long-ranged orientational order featuring a quasilayered periodicity imparted by a helical director configuration but lacking long-range positional order. Doping molecular cholesteric LCs with strongly anisotropic uniaxial colloidal particles adds another level of complexity because of the interplay between weak surface-anchoring boundary conditions and bulk-based elastic distortions near the particle-LC interface. Using cylindrical colloidal disks and rods with different geometric shapes and surface conditions, we demonstrate that these colloidal inclusions generically exhibit biaxial orientational probability distributions which may impart anomalously strong local biaxiality onto the hybrid cholesteric LC structure. Unlike nonchiral hybrid molecular-colloidal LCs, where biaxial order emerges only at critical colloid volume fractions exceeding some uniaxial-biaxial transition value, the orientational probability of the colloidal inclusions immersed in chiral nematic hosts is unambiguously biaxial even at infinite dilution. We demonstrate that the colloids induce local biaxial perturbations within the molecular orientational order of the LC host medium which strongly enhances the weak but native biaxial order of chiral nematic LC induced by the chiral symmetry breaking of the director field. With the help of analytical modeling and computer simulations based on the Landau–de Gennes free energy of the host LC around the colloids, we rationalize the observed multilevel biaxial order and conclude that it is not only unavoidable but also strongly enhanced compared to both achiral hybrid LCs and purely molecular cholesteric LCs.

biaxial liquid crystals↗