Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fluid Structure Interaction”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

kynema-fmb [SWR-23-07]

Kynema-FMB (FKA: Kynema) is an open-source performance portable flexible multibody (FMB) dynamics solver designed for time-domain simulations. While originally tailored for wind turbine structural dynamics, the formulation and implementation are those of a general flexible-multidbody dynamics solver that can readily be applied to a wide range of systems. Kynema was designed with a narrow focus, namely to provide a lightweight, fast, accurate FMD solver for coupling to computational-fluid-dynamics (CFD) codes, especially the CFD codes in the Kynema suite, for fluid-structure-interaction (FSI) simulations. Kynema-FMB is equipped to model systems that can be represented as a collection of beams and rigid bodies that are connected through constraints. Degrees of freedom are defined in the inertial/global frame of reference and include displacements and rotations (formally as rotation matrices, but stored as quaternions). The underlying formulation is built on a Lie-group time integrator designed for index-3 differential-algebraic equations, which is second-order accurate in time (Bruls et al., 2012). Beam models are based on geometrically exact beam theory and are discretized as high-order spectral finite elements similar to those in BeamDyn (Wang et al., 2017). The governing equations for a FMD system like a wind turbine constitute a highly nonlinear system of constrained partial-differential equations. Kynema-FMB uses analytical Jacobians in the nonlinear-system solves in each time step. Linear systems use sparse storage and several third-party sparse-linear-system solvers are enabled. Ill conditioning of linear systems is mitigated with preconditioning described in Bottasso et al, 2008. Kynema-FMB is integrated with a simple open-source controller (ROSCO). There is an application programming interface (API) for coupling to geometry-resolved CFD (like that in Sharma et al., 2023) and actuator-force CFD (like that in Kuhn et al., 2025). In the latter, for actuator-line models, Kynema-FMB includes an internal blade-element solver that depends on user-provided lookup tables for coefficients of lift and drag, i.e., aerodynamic polars. Kynema-FMB is written in C++ and leverages Kokkos and Kokkos-Kernels (KokkosEcosystem) as its performance portability layer enabling simulations on both CPU and GPU systems. The repository is equipped with extensive automated testing at the unit and regression/system levels. The following describes the high-level development objectives conceived for Kynema: *Kynema will follow modern software development best practices, including test-driven development (TDD), version control, hierarchical automated testing, and continuous integration (CI) for a robust development environment. *The core data structures are memory efficient and enable vectorization and parallelization at multiple levels. *Data structures are data-oriented to exploit methods for accelerated computing including high utilization of chip resources (e.g., single instruction multiple data (SIMD) instruction sets) and parallelization using GP-GPUs. *The computational algorithms incorporate robust open-source libraries for mathematical operations, resource allocation, and data management. *The API design considers multiple stakeholder needs and ensure integration with existing and future ecosystems for data science, machine learning, and AI. *Kynema-FMB is written in modern C++ and leverages Kokkos as its performance-portability library with inspiration from the kynema stack.

Sprague, MichaelA.↗
DOE OSTI2020

An energy-based coupling approach to nonlocal interface problems

Nonlocal models provide accurate representations of physical phenomena ranging from fracture mechanics to complex subsurface flows, where traditional partial differential equations fail to capture effects caused by long-range forces at the microscale and mesoscale. However, the application of nonlocal models to problems involving interfaces such as multimaterial simulations and fluid-structure interaction, is hampered by the lack of a rigorous nonlocal interface theory needed to support numerical developments. In this paper, we use an energy-based approach to develop a mathematically rigorous nonlocal interface theory which provides a physically consistent extension of the classical perfect interface PDE formulation. Numerical examples validate the proposed framework and demonstrate the scope of our theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

An energy-based coupling approach to nonlocal interface problems.

Nonlocal models provide accurate representations of physical phenomena ranging from fracture mechanics to complex subsurface flows, settings in which traditional partial differential equation models fail to capture effects caused by long-range forces at the microscale and mesoscale. However, the application of nonlocal models to problems involving interfaces, such as multimaterial simulations and fluid-structure interaction, is hampered by the lack of a physically consistent interface theory which is needed to support numerical developments and, among other features, reduces to classical models in the limit as the extent of nonlocal interactions vanish. In this paper, we use an energy-based approach to develop a formulation of a nonlocal interface problem which provides a physically consistent extension of the classical perfect interface formulation for partial differential equations. Numerical examples in one and two dimensions validate the proposed framework and demonstrate the scope of our theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Involute Working Group – FSI Analysis of Fuel Plates Using Finite Volume and Finite Element Methods

The three involute plate research reactors RHF, HFIR, and FRM II have expressed an interest in using computational software to carry their steady-state safety analysis. Since these tools represent a significant departure from the methods used currently (one-dimensional), the acceptability of the new approach by regulators requires thorough verification and validation of these tools. Therefore, Argonne National Laboratory and the three involute-plate reactors formed an informal group called the Involute Working Group aiming at qualifying computational tools to perform steady-state safety analysis. The present report focuses on a comparison of finite volume and finite element methods to model solids in fluid-structure interaction problems with the goal to estimate the coolant flow-induced fuel plate deflections obtained with the two methods. The finite volume method will be obsoleted in STARCCM+ by the end of 2021, nevertheless, this evaluation is important because the method was used by ANL researchers to model the response of the fuel plates, despite its drawbacks, which are discussed in the report. It was essential to check how those estimates compare to the results obtained with the finite element method that is considered superior for structural analysis. Various geometries, i.e., flat, cylindrical and circle-involute fuel plates, as well as coolant flow speed, were considered. The comparison shows that, independently of the plate geometry, the finite volume method significantly underestimates the deflection as compared to finite element method for coarser meshes. When the discretization is developed as a result of a mesh sensitivity study using finite element method, the result obtained using finite volume method can be a few times smaller than the corresponding finite element method solution. A code-to-code comparison , between STAR-CCM+ and LS-DYNA was included in the analysis. Within the LS-DYNA models, two types of finite element formulations were used: solid and shell finite elements. Mesh sensitivity study showed that both approaches converge to a similar value that was obtained with STAR-CCM+ finite element solver. The evaluation of the computational solvers was extended by adding two benchmark cases from the STAR-CCM+ Verification Suite and presented in the Appendix A. The selected cases are: (1) bending of a cantilever beam under external load, and (2) cylindrical shell deformation analysis, known in the literature as ‘Scordelis-Lo roof’. The problems were solved with finite volume, and finite element methods, and the results confirmed the previously discussed findings. The analysis shows that the finite element solver is superior to the finite volume solver in terms of representation of model geometry and estimating the structural behavior of fuel plates. Depending on the ratio of the load to the flexibility of the plate, the finite volume solver can greatly under- or overestimate the structural response if a very carefully selected mesh is not used.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Jan-Dec 2021 IC Project: w20_earthquake (Yearly Report)

The simulations conducted under this IC project used an in-house code developed at Los Alamos National Laboratory called HOSS (Hybrid Optimization Software Suite) and its extended version called HOSS-FSIS (fluid-structure interaction solver). HOSS has been under development in LANL’s Geophysics Group since 2004 and was used to simulate the signal propagation in pieces of material with pre-existing defects and also to simulate the near source effects of above ground nuclear explosions. One of the key feature of HOSS is that it is superior to pure finite element method and discrete element method implementations because it allows for a natural transition from continuum to discontinuum when the material breaks or fractures.

58 GEOSCIENCES↗
DOE OSTI2022

Design of High-deflection Foils for Marine Hydrokinetic Applications (Final Report)

This project used model-scale tank testing and fluid-structure-interaction (FSI) simulations to investigate the behavior of foils with large deflections and the effect of these deflections on crossflow turbine performance with the goal of determining the maximum allowable deflections consonant with efficiency and a robust, durable structure. A validated modeling and simulation approach was used in the design of an Ocean Renewable Power Company (ORPC) full-scale turbine. The methodology is applicable to the design of other marine hydrokinetic (MHK) devices, and to hydrokinetic turbines which experience significant deflections during operation.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

50 MW Segmented Ultralight Morphing Rotors for Wind Energy

A multi-institutional team designed a 50 Megawatt (MW) rated wind turbine featuring downwind aeroelastic morphing to reduce blade loads and allow an ultralight segmented rotor (mass reduction of about 25% compared to a conventional upwind rotor). The team used a control co-design approach with state-of-the art simulations including for the rotor and tower design using non-linear fluid-structure interactions and control algorithms, and the team also designed, built, and field-tested an aeroelastically-scaled downwind rotor to demonstrate this novel technology and validate the design tool fidelity. In a follow-on phase, these results were used along with an updated Levelized Cost Of Energy (LCOE) methodology to co-design a more detailed set of 25 MW rated turbine designs (including individual pitch control) based on minimum LCOE with highly flexible blades for an Atlantic Ocean offshore fixed-bottom design targeted towards market technology evaluation. A set of upwind and downwind designs at 25 MW rated scales were found to provide the best LCOE, with strong improvements over all previous offshore reference turbines (NREL 5MW, DTU 10 MW and IEA 15 MW). The optimized design would represent the world’s largest offshore turbine design, which combines several state-of-the-art structural, aerodynamic, and control technologies into a new and optimized system concept.

17 WIND ENERGY↗
DOE OSTI2024

Impacts of Irradiation Structural Behavior on Thermal Hydraulics Safety Analysis to Support MURR LEU Conversion

The University of Missouri Research Reactor (MURR) located in Columbia, Missouri is one of six U.S. High Performance Research Reactors (USHPRR), including one critical facility, that is actively collaborating with U. S. Department of Energy (DOE) National Nuclear Security Administration (NNSA) Material Management and Minimization (M3) Office of Reactor Conversion and Uranium Supply to convert from the use of highly enriched uranium (HEU; ≥20 wt% U-235) to low-enriched uranium (LEU; <20 wt% U-235) fuel. A new type of very high-density LEU fuel based on an alloy of uranium and 10 wt% molybdenum (U-10Mo) is expected to allow the conversion to LEU of MURR, as well as four other USHPRR. MURR has been working with the Reactor Conversion Pillar at Argonne to perform fuel element design and fuel cycle performance analyses, steady-state thermal hydraulics safety analysis, and accident safety analyses in preparation for the conversion of MURR and to support a preliminary safety analysis report for conversion to LEU fuels. Subsequent analyses have also been performed, including transition cycles where all-fresh LEU fuel elements are introduced upon conversion and progressing through reactor operations the core is brought to equilibrium. Thermal hydraulics safety analyses performed as part of the above have employed an assumption on channel gap reduction due to burnup-related phenomena including fuel swelling, irradiation creep, and oxide layer buildup. Recently, a series of structural analyses have been performed on the MURR LEU fuel plates and an element due to significant differences between the plate and element designs of the MURR HEU and LEU fuels. In addition, NUREG-1537 indicates that structural phenomena are to be evaluated. Two separate types of structural analyses were performed for the MURR LEU fuel element: fluid-structure interaction (FSI) and irradiation thermo-mechanical. The FSI analysis evaluated the effects of hydraulic forces on the MURR LEU fuel element to quantify the flow-induced plate deflection, and a minimal impact to the channel gap thickness was predicted under prototypic and bounding conditions. The irradiation thermo-mechanical analysis evaluated the effects of fuel swelling, irradiation creep, and thermal expansion for the MURR LEU plates and the element for prototypic thermal and irradiation conditions based on a high-fidelity approach multiphysics approach. Overall, this thermo-mechanical analysis predicts smaller gap thickness changes in previously limiting regions. Larger changes are predicted in the middle of channels, and for end channels where power density is not typically a maximum. An additional thermo-mechanical analysis was performed for the outermost HEU fuel plate, which showed a similar magnitude of deflection as the outermost LEU plate. Due to substantial differences between the channel gap reductions assumed for the previous safety analyses and those predicted by the irradiation thermo-mechanical analysis, a need to evaluate their impact on the thermal hydraulics safety analyses arose. This report presents the results from the steady-state safety analyses for normal operation as well as the accident analyses for the two most limiting accident scenarios.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Open-Source Tidal Energy Converter (OSTEC) Testbed: Design Basis Report

This report describes the design basis and design details for an instrumented marine turbine system intended to serve as the DOE’s marine tidal turbine test bed for foundational open-source R&D and data generation to advance our understanding and to identify knowledge gaps on the techno-economic performance of tidal energy converters (TEC) under real tidal flow conditions and at sufficiently large scale to enable upscaling of fluid-structure-interactions and component and material load responses

16 TIDAL AND WAVE POWER↗
DOE OSTI2026

Dynamic Wind Loading on Compliant, Floating PV Arrays: Cooperative Research and Development (Final Report)

PV panels are subjected to wind loads during normal outdoor operation, where th wind speed, wind direction, panel angle, and array layout play a large role in the overall loading magnitude. For floating PV systems, where panels are installed on floating rafts, these forces can lead to a dynamic displacement of the raft and mounted PV hardware. The Contractor and Participant will perform fluid-structure interaction simulations of this phenomenon for a variety of wind speeds, directions, and panel angles to characterize these forces and help design mooring/tethering lines to resist and anchor raft movement.

14 SOLAR ENERGY↗
DOE OSTI2024

A hybrid finite volume method and smoothed particle hydrodynamics approach for efficient and accurate blast simulations

Modeling strong shock waves in fluids remains a persistent challenge in computational physics. Essential to research efforts in industry and defense, numerous methods have been devised to improve the accuracy and efficiency of shock simulations. A novel, hybrid Finite Volume Method (FVM)-Smoothed Particle Hydrodynamics (SPH) approach is capable of further improving efficiency and retaining accuracy by exploiting the favorable characteristics of each respective method. This hybrid approach is presented for shock capturing in compressible fluids. The Python framework Pyro2 is employed to simulate a coarse FVM mesh, while the Python framework PySPH is utilized to model the fluid in regions with high gradients through SPH particles. The performance of the hybrid FVM-SPH scheme, compared to the individual FVM and SPH methods, is assessed in 1 kt and 10 kt blast simulations. Our results indicate that the hybrid approach offers higher computational efficiency than SPH while preserving its accuracy and characteristics. The hybrid approach had a relative speedup of 11.3x and 22.3x over the FVM and SPH approaches for the 1 kt simulation and a relative speedup of 14.7x and 20.9x over the FVM and SPH approaches for the 10 kt simulation. The hybrid SPH algorithm enables future compressible fluid simulations with more extensive capabilities than grid-based methods alone, presenting potential applications in modeling fluid-structure interactions and solid deformation and fracturing in blast simulations.

Myers, Conner↗
DOE OSTI2022

Removal of Non-Specifically Bound Proteins Using Rayleigh Waves Generated on ST-Quartz Substrates

Label-free biosensors are plagued by the issue of non-specific protein binding which negatively affects sensing parameters such as sensitivity, selectivity, and limit-of-detection. In the current work, we explore the possibility of using the Rayleigh waves in ST-Quartz devices to efficiently remove non-specifically bound proteins via acoustic streaming. A coupled-field finite element (FE) fluid structure interaction (FSI) model of a surface acoustic wave (SAW) device based on ST-Quartz substrate in contact with a liquid loading was first used to predict trends in forces related to SAW-induced acoustic streaming. Based on model predictions, it is found that the computed SAW body force is sufficient to overcome adhesive forces between particles and a surface while lift and drag forces prevent reattachment for a range of SAW frequencies. We further performed experiments to validate the model predictions and observe that the excitation of Rayleigh SAWs removed non-specifically bound (NSB) antigens and antibodies from sensing and non-sensing regions, while rinsing and blocking agents were ineffective. An amplified RF signal applied to the device input disrupted the specific interactions between antigens and their capture antibody as well. ST-quartz allows propagation of Rayleigh and leaky SH-SAW waves in orthogonal directions. Thus, the results reported here could allow integration of three important biosensor functions on a single chip, i.e., removal of non-specific binding, mixing, and sensing in the liquid phase.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Aeroelastic Modeling and Full-Scale Loads Measurements for Investigation of Single-Axis PV Tracker Wind-Driven Dynamic Instabilities

While wind tunnel testing and proprietary industry modeling tools have been used for years to design and develop PV tracker systems, wind induced dynamic failures are becoming more prevalent and high visibility. Designers and manufactures have reacted by developing new systems or add-on products to improve system dynamics adding to overall system costs. To better understand the physics and aeroelastic behavior that leads to dynamic instabilities NREL researches have taken up a first of a kind study to both develop open source aeroelastic modeling tools and measure tracker loads on a single axis full-scale tracker system at a high wind site. These modeling tools can simulate the fluid-structure interaction driving torsional instabilities under a wide range of turbulent inflow conditions and stow angles. The results of these simulations reveal the flow features associated with panel rotation and the induced loads in the structure. The experimental load measurements can be used to validate models and provide a quantifiable understanding of tracker dynamics, critical loads paths, identify instability markers, inform resilient design, and suggest the most favorable stow approach.

14 SOLAR ENERGY↗
DOE OSTI2021

Prototype and Codesign of Nascent Flexible Wave Energy Converter Concepts

Prevailing mainstream ocean wave energy converter (WEC) archetypes are often relatively large monolithic rigid bodies with singular means of power take-off (e.g., relative motion between rigid bodies) and solo prime mover mechanisms (e.g. a rotary generator, a hydraulic piston). These mainstream WEC archetypes have existed for decades with some having histories spanning more than a century of significant investments and interests. As an alternative to the domain of mainstream WEC archetypes is the domain of distributed embedded energy converter technologies (also known as DEEC-Tec). Use of DEEC-Tec to create WECs - often called flexWECs - is a nascent domain of research and development that centers itself upon the interconnection of many relatively small energy converters - very small energy transducers - that interconnect and combine to create a much larger, overall, energy converting structure. The DEEC-Tec domain enables promising possibilities such as: direct in situ ocean wave energy harvesting and conversion throughout the volume of an entire flexWEC; inherent broad-banded ocean wave energy capture and conversion arising from near semi-continuous structural degrees of freedom; an innate absence of focusing and concentrating of ocean wave energy into a single prime mover and/or transmission system; and built-in implicit redundancy. As such, the DEEC-Tec domain could revolutionize how we conceptualize ocean wave energy conversion - both in terms of ocean wave energy converter design and what we view as viable environments and circumstances for ocean wave energy conversion. The research and development opportunities for the DEEC-Tec domain, therefore, are vast compared to the prevailing mainstream ocean wave energy converter domain. Correspondingly, there are ample opportunities to invest time and money into the DEEC-Tec domain with an aim to enable a paradigm shift in how ocean waves can be utilized as a more viable form of renewable energy. To that end, two seedling efforts have been enacted providing a very cursory investigation into the DEEC-Tec domain: (1) a preliminary simple fabrication/prototyping effort and (2) a preliminary co-design effort. The prototype seedling centers upon the design of a fabrication and evaluation rig for individual DEEC-Tec transducers; e.g., dielectric elastomer generators. The co-design seedling centers upon the analysis of those methods suitable for DEEC-Tec, ocean wave energy conversion, and fluid-structure interaction co-design efforts.

50 EE - Wind and Water Power Program - Water (EE-4↗
DOE OSTI2021

Application of Machine Learning to Predict the Response of the Liquid Mercury Target at the Spallation Neutron Source

The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory is currently the most powerful accelerator-driven neutron source in the world. The intense proton pulses strike on SNS’s mercury target to provide bright neutron beams, which also leads to severe fluid-structure interactions inside the target. Prediction of resultant loading on the target is difficult particularly when helium gas is intentionally injected into mercury to reduce the loading and mitigate the pitting damage on the target’s internal walls. Leveraging the power of machine learning and the measured target strain, we have developed machine learning surrogates for modeling the discrepancy between simulations and experimental strain data. We then employ these surrogates to guide the refinement of the high-fidelity mercury/helium mixture model to predict a better match of target strain response.

Lin, Lianshan↗
DOE OSTI2024

Effects of Interactions Between Produced Formation Fluid and Rock Matrix on Pore Structure of Caney Shale, Southern Oklahoma

ABSTRACT: Rock-fluid interactions change properties of shales during exploitation. To investigate effects of rock-fluid interactions on pore structure of shales matrix after hydraulic fracturing, powder samples from two late Mississippian Caney Shale cores in the Ardmore Basin, southern Oklahoma, were used to react with formation produced fluid from the field in the batch reactor analysis. X-ray diffraction for mineralogy and Low-pressure nitrogen adsorption isotherms for pore structure were measured for original, after-7days, and after-30days samples. Results show that the samples consist mainly of quartz, followed by clay minerals, carbonates, and feldspar. The pore sizes of micropore (<2 nm) and mesopore (2-50 nm) increase 14%-233% due to dissolution of pyrite, feldspar, and carbonates after 7 days. Due to the transformation from smectite to illite and the increase of pore size, the specific surface area (SSA) decreases after 7-days interactions. After 30-days interactions, the micropore volume slightly increases and the mesopore and macropore volume decreases. Due to the decrease of pore size, the SSA of 30-days reacted samples increases correspondingly and is lower (for the clay-rich sample) or higher (for the calcareous sample) than that of the unreacted samples. Findings improve our understanding of dynamic alteration of shale properties during production. 1. INTRODUCTION Energy demand will continuously grow owing to the increasing global population as well as energy consumption (EIA, 2023). On the other hand, shale gas and oil reshaped the energy market in the United States, enabling the United States to become a net-export of natural gas country in 2017 (EIA, 2023). However, shale reservoirs are challenging tight formations that are still poorly understood in the extraction and production of hydrocarbons (Ross and Bustin, 2009; Curtis et al., 2012; Xiong et al., 2015, 2021a; Li Y. et al., 2016; Gong et al., 2019a; Benge et al., 2021; Awejori et al., 2022; Huang et al., 2022). One of the most challenging topics is the rock-fluid interactions post hydraulic fracturing and its subsequent impacts on the pore structures of fractured formation matrix.

Xiong, Fengyang↗
DOE OSTI2022

Modeling and Analysis of a Novel Offshore Binary Species Free-Floating Longline Macroalgal Farming System

The investigation of innovative macroalgal cultivation is important and needed to optimize farming operations, increase biomass production, reduce the impact on the ecosystem, and lower system and operational costs. However, most macroalgal farming systems (MFSs) are stationary, which need to occupy a substantial coastal area, require extensive investment in farm infrastructure, and cost high fertilizer and anchoring expenses. This study aims to model, analyze, and support a novel binary species free-floating longline macroalgal cultivation concept. The expected outcomes could provide a basis for the design and application of the novel MFS to improve biomass production, decrease costs, and reduce the impact on the local ecosystem. In this paper, Saccharina latissima and Nereocystis luetkeana were modeled and validated, and coupled with longline to simulate the binary species MFS free float in various growth periods and associated locations along the US west coast. Further, the numerical predictions indicated the possibility of failure on the longline and breakage at the kelp holdfasts is low. However, the large forces due to an instantaneous change in dynamic loads caused by loss of hydrostatic buoyancy when the longline stretches out of the water would damage the kelps. Buoy-longline contact interactions could damage the buoy, resulting in the loss of the system by sinking. Furthermore, the kelp-longline and kelp-kelp entanglements could potentially cause kelp damage.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A sharp immersed method for 2D flow-body interactions using the vorticity-velocity Navier-Stokes equations

Immersed methods discretize boundary conditions for complex geometries on background Cartesian grids. Here, this makes such methods especially suitable for two-way coupled flow-body problems, where the body mechanics are partially driven by hydrodynamic forces. However, for the vorticity-velocity form of the Navier-Stokes equations, existing immersed geometry discretizations for two-way coupled problems only achieve first order spatial accuracy near solid boundaries. Here we introduce a sharp-interface approach based on the immersed interface method to handle the one- and two-way coupling between an incompressible flow and one or more rigid bodies using the 2D vorticity-velocity Navier-Stokes equations. Our main contributions are three-fold. First, we develop and analyze a moving boundary treatment for sharp immersed methods that can be applied to PDEs with implicitly defined boundary conditions, such as those commonly imposed on the vorticity field. Second, we develop a two-way coupling methodology for the vorticity-velocity Navier-Stokes equations based on control-volume momentum balance that does not require the pressure field. Third, we show through extensive testing and validation that our resulting flow-body solver reaches second-order accuracy for most practical scenarios, and provides significant efficiency benefits compared to a representative first-order approach.

42 ENGINEERING↗