Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “INCLINATION”

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,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Finite element model for piles in liquefiable ground

This paper develops a three dimensional finite element modelling method for piles in liquefiable ground and applies it to the analysis of seismic pile responses. A unified plasticity model for large post-liquefaction shear deformation of sand provides the basis for the effective and efficient modelling of liquefiable ground. Special attention is dedicated towards the modelling of piles and soil-pile interface to accurately reflect the behaviour of piles. A staged modelling procedure is adopted to appropriately generate the initial conditions for the soil and piles and achieve hydrostatic pore pressure prior to seismic loading. Three centrifuge shaking table tests on single piles, both with and without pile cap and superstructure, in level and inclined liquefiable ground are conducted and simulated in validation and application of the proposed method. Further studies to investigate the effects of pile cap, lateral spreading, and non-liquefiable surface layer are undertaken numerically using the validated method. The results show these aforementioned factors to be influential in the dynamic and residual response of piles in liquefiable ground.

Centrifuge test↗
DOE OSTI2014

Quasi-static energy absorption of hollow microlattice structures

We present a comprehensive modeling and numerical study focusing on the energy quasi-static crushing behavior and energy absorption characteristics of hollow tube microlattice structures. The peak stress and effective plateau stress of the hollow microlattice structures are deduced for different geometrical parameters which gives volume and mass densities of energy absorption, D-v and D-m, scale with the relative density, (rho) over bar, as D-v similar to (rho) over bar (1) (5) and D-m similar to (rho) over bar (0 5), respectively, fitting very well to the experimental results of both 60 degrees inclined and 90 degrees predominately microlattices. Then the strategies for energy absorption enhancement are proposed for the engineering design of microlattice structures. By introducing a gradient in the thickness or radius of the lattice members, the buckle propagation can be modulated resulting in an increase in energy absorption density that can exceed 40%. Liquid filler is another approach to improve energy absorption by strengthening the microtruss via circumference expansion, and the gain may be over 100% in terms of volume density. Insight into the correlations between microlattice architecture and energy absorption performance combined with the high degree of architecture control paves the way for designing high performance microlattice structures for a range of impact and impulse mitigation applications for vehicles and structures. (C) 2014 Elsevier Ltd. All rights reserved.

Liu, YL↗
DOE OSTI2021

Fiber/matrix debonding evaluation of SiCf/SiC composites using micropillar compression technique

Third-generation silicon carbide (SiC) composites reinforced by SiC fibers (Hi-Nicalon S [HNS] and Tyranno SA3 [SA3]) are attractive for use in next generation reactors owing to their high strength and chemical inertness at high temperatures, as well as enhanced radiation tolerance under neutron irradiation environments. To optimize composite performance, the interfacial mechanical properties of chemical vapor–infiltrated (CVI) SiC f /SiC composites are investigated in this effort by using a slant interface micropillar compression testing procedure. The micropillar test specimens, containing an inclined pyrolytic carbon (PyC) interphase, are prepared using a focused ion beam. The novel microcompression testing successfully quantifies the debond shear strength and internal friction coefficient of micropillar test samples by using the Mohr-Coulomb formulation. Furthermore, according to four types of SiCf/SiC composite microcompression test results, interfacial properties and debond mechanisms are significantly affected by the PyC layer thickness, the local bonding mechanism of the PyC interphase on the SiC fiber surface, and the surface roughness of fibers. Regardless of PyC thicknesses, SA3-reinforced CVI SiC f /SiC composites are found to have much higher debond shear strengths than HNS-reinforced SiC f /SiC CVI composites. By using this micropillar compression technique alongside analytical methods, we uncover new understandings of PyC interface properties. Additionally, the micropillar test results obtained are correlated with macroscopic mechanical properties of neutron-irradiated CVI SiCf/SiC composites.

42 ENGINEERING↗
DOE OSTI2015

Deformation fields near a steady fatigue crack with anisotropic plasticity

In this work, from finite element simulations based on an irreversible, hysteretic cohesive interface model, a steady fatigue crack can be realized if the crack extension exceeds about twice the plastic zone size, and both the crack increment per loading cycle and the crack bridging zone size are smaller than the plastic zone size. The corresponding deformation fields develop a plastic wake behind the crack tip and a compressive residual stress field ahead of the crack tip. In addition, the Hill’s plasticity model is used to study the role of plastic anisotropy on the retardation of fatigue crack growth and the elastic strain fields. It is found that for Mode-I cyclic loading, an enhanced yield stress in directions that are inclined from the crack plane will lead to slower crack growth rate, but this retardation is insignificant for typical degrees of plastic anisotropy. Furthermore, these results provide key inputs for future comparisons to neutron and synchrotron diffraction measurements that provide full-field lattice strain mapping near fracture and fatigue crack tips, especially in textured materials such as wrought or rolled Mg alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Unusual size effects from tilted twin boundaries in nano-twinned metals

Tilted twin boundaries (TBs), whose plane normals are rotated at an angle from the parent grain axis, naturally occur in columnar-grained, nano-twinned (NT) metals. Here, using a combination of atomistic simulations and analytical modeling, we reveal that NT metals with the ideal, non-tilted TBs exhibit continuously increasing strength with decreasing twin thickness, and hence, no inverse twin thickness size effect on strength. In contrast, NT metals with tilted TBs exhibit an inverse size effect, and the critical twin thickness, below which strength decreases, increases as the TB-tilt angle increases. The analysis also identifies a critical value of TB tilt, for which strength becomes independent of twin thickness and is the weakest. The transition arises from a change in dislocation activity prevailing mostly on planes inclined to TBs to planes parallel to the TBs. Overall, these findings reveal a profound influence of TB tilt angle that could redirect the analysis and engineering of nano-twin structures.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Rare earth element components in atmospheric particulates in the Bayan Obo mine region

The Bayan Obo mine, located in Inner Mongolia, China, is the largest light rare earth body ever found in the world. The research for rare earth elements (REEs) enrichment in atmospheric particulates caused by mining and ore processing is fairly limited so far. In this paper, atmospheric particulates including total suspended particulate (TSP) matter and particles with an equivalent aerodynamic diameter less than 10 μm (PM{sub 10}) were collected around the Bayan Obo mine region, in August 2012 and March 2013, to analyze the levels and distributions of REEs in particles. The total concentrations of REEs for TSP were 149.8 and 239.6 ng/m{sup 3}, and those for PM{sub 10} were 42.8 and 68.9 ng/m{sup 3}, in August 2012 and March 2013, respectively. Enrichment factor was calculated for all 14 REEs in the TSP and PM{sub 10} and the results indicated that REEs enrichment in atmosphere particulates was caused by anthropogenic sources and influenced by the strong wind in springtime. The spatial distribution of REEs in TSP showed a strong gradient concentration in the prevailing wind direction. REE chondrite normalized patterns of TSP and PM{sub 10} were similar and the normalized curves inclined to the right side, showing the conspicuous fractionation between the light REEs and heavy REE, which supported by the chondrite normalized concentration ratios calculated for selected elements (La{sub N}/Yb{sub N}, La{sub N}/Sm{sub N}, Gd{sub N}/Yb{sub N}). - Highlights: • TSP and PM{sub 10} samples were collected to analyze the levels and distributions of REE. • Enrichment factors indicated that REE enrichment was caused by anthropogenic sources. • The distribution of REEs showed a strong gradient in the prevailing wind direction. • Obvious fractionation between LREEs and HREEs is observed in atmospheric particulates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Impact of molten salt inflow on the temperature distribution in thermal energy storage tanks at startup for central receiver concentrating solar power plants

Concentrating Solar Power (CSP) systems with molten salt thermal energy storage (TES) tanks are one of the most promising, renewable-based energy conversion technologies for larger-scale power generation. The TES tank is one of the most critical components in CSP plants due to its high-temperature operation (up to 565 °C), daily thermal cycling, and intermittent solar radiation conditions. The plant startup is one of the most challenging operation conditions that could lead to damaging thermal gradients due to low salt inventory levels. In this study an analytical model for the sparger ring was developed and integrated with a detailed computational fluid dynamics model of a commercial-scaled molten salt tank. The integrated model allows an accurate representation of the tank operation to evaluate the effect of molten salt inflow on the mixing process. The tank filling process during plant startup was analyzed considering sparger rings with variations in design features, including inlet orifice configurations, number of orifices, direction of the inlets, and orifice diameter. The results demonstrated that higher temperature gradients are obtained in the tank floor during the plant startup. A sparger ring configuration with a predetermined orifice inlet inclination (30°, 45° and 60°) leads to significant temperature differences in the floor, between 57 °C and 62 °C, but better homogeneity in the temperature of the salt inventory. Lower salt inflow velocities result in a more homogeneous floor temperature, with maximum temperature differences under 38 °C. The sparger ring configuration with 52 orifices of 1-in. diameter and vertical flow showed better homogeneity in the temperature differences as a function of the salt level and lower temperature gradients in the tank floor. Because large temperature gradients in the tank's floor have been identified as one of the main factors contributing to tank failures, assessing various sparger ring design features is fundamental to determining proper inflow conditions that lead to low-temperature gradients and reducing failure susceptibility.

14 SOLAR ENERGY↗
DOE OSTI2024

Elucidate the molecular basis of ampholytic chitosan as a high-performance cryoprotectant to myosin denaturation: The importance of saccharide charges

The uses of charged poly/oligo-saccharides enable the retarding of protein denaturation against various environmental stresses during food storage and manufacturing. However, at subzero temperatures, the molecular basis of such stabilization behaviors, i.e., cryoprotections, remain less explored. Here, in this study, we introduced an ampholytic saccharide, carboxymethyl chitooligosaccharide (CMCO) that effectively inhibited the freezing-induced myosin denaturation. The in-depth cryoprotective mechanism was systematically investigated by using molecular dynamic simulation and multispectral characterizations. Results showed that CMCO may interact with myosin through hydrogen bonding and electrostatic interactions, which caused the expelling of water at protein surfaces and the reduced conformational flexibility of myosin molecules. Due to this water replacement event, both secondary and ternary structures of myosin became freezing-resistant, leading to the inhibited protein aggregations and retained functionalities, such as solubility, Ca 2+ -ATPase activity, and gelling properties. Moreover, cryoprotective behaviors of CMCO were charge-dependent. CMCO with a higher degree of carboxymethyl substitution (DS: 1.2) was inclined to bind and stabilize myosin molecules better than the low-DS (DS: 0.8) one, even though both outperformed other cryoprotective saccharides. Therefore, this investigation not only introducing a high-performance myosin cryoprotectant, but also elaborated the cryoprotective mechanism of ampholytic saccharides.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Perspective on the oil-dominated gas hydrate plugging conceptual picture as applied to transient Shut-In/Restart

Gas hydrate formation is rarely an issue during continuous production conditions. Unplanned and/or extended shut-downs comprise the majority of scenarios where overlapping concern and uncertainty for gas hydrate blockages exist. To date, most experimental studies and theoretical analysis of gas hydrate formation derive from the steady-state conceptual picture. Yet, significant distinctions induced by phase separation and prolonged cooling suggest an altered conceptual picture is required for transient shut-in/restart operations. To address this outstanding issue, this work presents a conceptual picture for the base case of a non-surface active hydrocarbon undergoing cold restart without enacted mitigation efforts. First, phase separation upon shut-in leads to the localization of gas hydrate ingredients due to fluids density differences and pipeline inclination. Cooling to the seafloor temperature, then, induces gas hydrate film formation at the hydrocarbon interface and maximizes subcooling, resulting in rapid and confined gas hydrate growth occurring upon restart. Furthermore, this ultimately leads to plugging dominated by the deposition of porous, low conversion gas hydrate masses. Direct visual observations and results from Raman spectroscopy, hydrate contact angle and film growth, micromechanical force apparatus, rocking cell, and flow loop studies support the mechanisms outlined within each step. Similar to the influence provided by the original Colorado School of Mines steady-state conceptual picture that was introduced about 15 years ago, this transient counterpart builds a foundation that could inspire new directions of hydrate research and analysis in areas where uncertainty abounds.

02 PETROLEUM↗
DOE OSTI2018

Scrape-off layer plasma and neutral characteristics and their interactions with walls for FNSF

Simulations of the heat flux on plasma facing components from exhausting core plasma are reported for two possible Fusion Nuclear Science Facility (FNSF) divertor configurations. One configuration utilizes divertor plates strongly inclined with respect to the poloidal magnetic flux surfaces like that planned for ITER and results in a partially detached divertor-plasma. The second configuration has divertor plates orthogonal to the flux surfaces, which leads to a fully detached divertor-plasma if the width of the divertor region is sufficient. Both configurations use scrape-off layer impurity seeding to yield an acceptable peak heat flux of ~10 MW/m 2 or smaller on the divertor plates and chamber walls. The roles of recycled hydrogenic atoms and molecules are investigated and distribution of sputtering tungsten throughout the edge region modeled. Furthermore, the simulations are performed with the UEDGE 2D transport code to model both plasma and neutral components with supplementary neutral modeling performed with the DEGAS 2 Monte Carlo code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

CO 2 storage and potential fault instability in the St. Lawrence Lowlands sedimentary basin (Quebec, Canada): Insights from coupled reservoir-geomechanical modeling

In this paper, coupled reservoir-geomechanical (TOUGH-FLAC) modeling is applied for the first time to the St. Lawrence Lowlands region to evaluate the potential for shear failure along pre-existing high-angle normal faults, as well as the potential for tensile failure in the caprock units (Utica Shale and Lorraine Group). This activity is part of a general assessment of the potential for safe CO 2 injection into a sandstone reservoir (the Covey Hill Formation) within an Early Paleozoic sedimentary basin. Field and subsurface data are used to estimate the sealing properties of two reservoir-bounding faults (Yamaska and Champlain faults). The spatial variations in fluid pressure, effective minimum horizontal stress, and shear strain are calculated for different injection rates, using a simplified 2D geological model of the Becancour area, located ~110 km southwest of Quebec City. The simulation results show that initial fault permeability affects the timing, localization, rate, and length of fault shear slip. Contrary to the conventional view, our results suggest that shear failure may start earlier for a permeable fault than for a sealing fault, depending on the site-specific geologic setting. In simulations of a permeable fault, shear slip is nucleated along a 60 m long fault segment in a thin and brittle caprock unit (Utica Shale) trapped below a thicker and more ductile caprock unit (Lorraine Group) – and then subsequently progresses up to the surface. In the case of a sealing fault, shear failure occurs later in time and is localized along a fault segment (300 m) below the caprock units. The presence of the inclined low-permeable Yamaska Fault close to the injection well causes asymmetric fluid-pressure buildup and lateral migration of the CO 2 plume away from the fault, reducing the overall risk of CO 2 leakage along faults. Finally, fluid-pressure-induced tensile fracturing occurs only under extremely high injection rates and is localized below the caprock units, which remain intact, preventing upward CO 2 migration.

58 GEOSCIENCES↗
DOE OSTI2020

Aluminum-doped calcium manganite particles for solar thermochemical energy storage: Reactor design, particle characterization, and heat and mass transfer modeling

A two-step cycle is discussed for solar thermochemical energy storage based on particulate aluminum-doped calcium manganite reduction/oxidation reactions for direct integration into Air Brayton cycles. The two steps encompass 1) the storage of concentrated solar irradiation within endothermic reduction of aluminum-doped calcium manganite and 2) the delivery of heat to an Air-Brayton cycle via exothermic re-oxidation of oxygen-deficient aluminum-doped calcium magnanite. A 5 kW th scale solar thermochemical reactor operating under vacuum was designed, modeled, and optimized to thermally reduce a continuous, gravity-driven flow of aluminum-doped calcium manganite particles. The granular flows were characterized in a tilt-flow rig, and particle image velocimetry was used to determine flow properties via frictional and velocity scaling relationships. Flow properties were integrated into a detailed heat and mass transfer model of the solar thermochemical reactor. A reactor design with 31° inclination angle, 230 g/min of particles, and 5.2 kW th radiative input from the high-flux solar simulator was discovered to produce an outlet flow temperature of 1158 K, with stoichiometric deviations of 0.076 and a storage efficiency of 0.628 while avoiding particle overheating and promoting longer particle residence times.

14 SOLAR ENERGY↗
DOE OSTI2025

Fuel cells for single-aisle regional aircraft: System configuration, performance and cost

A hydrogen fuel cell propelled electric aircraft can compete with incumbent turbofan technologies for single-aisle regional aircraft by coupling design of stack, air handling, thermal management, propulsion, and airframe to optimize performance. The stack operates at 95°C to facilitate heat rejection during take-off and below 75°C during cruise to extend lifetime and is oversized to satisfy power requirements at end of life. A multi-stage turbocompressor with a compression ratio >10 is selected to reach high stack power density at 11,300-m cruise altitude. The propulsion system is configured to accommodate air handling within the core duct, an inclined heat exchanger in the outer duct to limit the nacelle size, and variable area nozzles to independently control mass flows through the core and bypass ducts. The airframe is modified for maximum lift coefficient and longer balanced field length for dramatically reduced thrust during take-off, and the fuselage is stretched by 20% to store liquid hydrogen (LH 2 ). Modularization of power systems promotes safety in one engine inoperative scenarios and allows reaching specific power metrics for stack, balance-of-plant and fuel cell system (FCS), necessary for acceptable take-off weight. In conclusion, cost parity requires increase in FCS lifetime, LH 2 cost reduction, and improved FCS specific power.

Catalyst durability↗
DOE OSTI2017

A broadly-applicable unified closure relation for Taylor bubble rise velocity in pipes with stagnant liquid

Taylor bubble velocity in slug flow is a closure relation which significantly affects the prediction of liquid holdup (or void fraction) and pressure gradient in mechanistic models of slug flow for oil and gas pipe applications. In this work, we use a validated Computational Fluid Dynamics (CFD) approach to simulate the motion of Taylor bubbles in pipes; the interface is tracked with a Level-Set method implemented in a commercial code. A large numerical database is generated covering the most ample range of fluid properties and pipe inclination angles explored to date (Eo ∊ [10, 700],Mo ∊ [1 x 10 -6 , 5 x 10 3 , and θ ∊ [0°, 90°]). A unified Taylor bubble rise velocity correlation is extracted from that database. The new correlation predicts the numerical database with 8.6% absolute average relative error and a coefficient of determination R 2 = 0.97 and other available experimental data with 13.0% absolute average relative error and R 2 = 0.84 outperforming existing correlations and models.

Mechanics↗
DOE OSTI2018

A unified framework for polycrystal plasticity with grain boundary evolution

Plastic deformation in polycrystals is governed by the interplay between intra-granular slip and grain boundary-mediated plasticity. However, while the role played by bulk dislocations is relatively well-understood, the contribution of grain boundaries (GBs) has only recently begun to be studied. GB plasticity is known to play a key role along with bulk plasticity under a wide range of conditions, such as dynamic recovery, superplasticity, severe plastic deformation, etc., and developing models capable of simultaneously capturing GB and bulk plasticity has become a topic of high relevance. In this paper we develop a thermodynamically-consistent polycrystal plasticity model capable of simulating a variety of grain boundary-mediated plastic processes in conjunction with bulk dislocation slip. The model starts from the description of a single crystal and creates lattice strain-free polycrystalline configurations by using a specially-designed multiplicative decomposition developed by the authors. This leads to the introduction of a particular class of geometrically necessary dislocations (GND) that define fundamental GB features such as misorientation and inclination. The evolution of the system is based on an energy functional that uses a non-standard function of the GND tensor to account for the grain boundary energy, as well as for the standard elastic energy. Our implementation builds on smooth descriptions of GBs inspired on diffuse-interface models of grain evolution for numerical convenience. Here we demonstrate the generality and potential of the methodology by simulating a wide variety of phenomena such as shear-induced GB sliding, coupled GB motion, curvature-induced grain rotation and shrinkage, and polygonization via dislocation sub-grain formation.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Low- and high-order accurate boundary conditions: From Stokes to Darcy porous flow modeled with standard and improved Brinkman lattice Boltzmann schemes

The present contribution focuses on the accuracy of reflection-type boundary conditions in the Stokes–Brinkman–Darcy modeling of porous flows solved with the lattice Boltzmann method (LBM), which we operate with the two-relaxation-time (TRT) collision and the Brinkman-force based scheme (BF), called BF-TRT scheme. In parallel, we compare it with the Stokes–Brinkman–Darcy linear finite element method (FEM) where the Dirichlet boundary conditions are enforced on grid vertices. In bulk, both BF-TRT and FEM share the same defect: in their discretization a correction to the modeled Brinkman equation appears, given by the discrete Laplacian of the velocity-proportional resistance force. This correction modifies the effective Brinkman viscosity, playing a crucial role in the triggering of spurious oscillations in the bulk solution. While the exact form of this defect is available in lattice-aligned, straight or diagonal, flows; in arbitrary flow/lattice orientations its approximation is constructed. At boundaries, we verify that such a Brinkman viscosity correction has an even more harmful impact. Already at the first order, it shifts the location of the no-slip wall condition supported by traditional LBM boundary schemes, such as the bounce-back rule. For that reason, this work develops a new class of boundary schemes to prescribe the Dirichlet velocity condition at an arbitrary wall/boundary-node distance and that supports a higher order accuracy in the accommodation of the TRT-Brinkman solutions. For their modeling, we consider the standard BF scheme and its improved version, called IBF; this latter is generalized in this work to suppress or to reduce the viscosity correction in arbitrarily oriented flows. Our framework extends the one- and two-point families of linear and parabolic link-wise boundary schemes, respectively called B-LI and B-MLI, which avoid the interference of the Brinkman viscosity correction in their closure relations. The performance of LBM and FEM is thoroughly evaluated in three benchmark tests, which are run throughout three distinctive permeability regimes. The first configuration is a horizontal porous channel, studied with a symbolic approach, where we construct the exact solutions of FEM and BF/IBF with different boundary schemes. The second problem refers to an inclined porous channel flow, which brings in as new challenge the formation of spurious boundary layers in LBM; that is, numerical artefacts that arise due to a deficient accommodation of the bulk solution by the low-accurate boundary scheme. The third problem considers a porous flow past a periodic square array of solid cylinders, which intensifies the previous two tests with the simulation of a more complex flow pattern. The ensemble of numerical tests provides guidelines on the effect of grid resolution and the TRT free collision parameter over the accuracy and the quality of the velocity field, spanning from Stokes to Darcy permeability regimes. It is shown that, with the use of the high-order accurate boundary schemes, the simple, uniform-mesh-based TRT-LBM formulation can even surpass the accuracy of FEM employing hardworking body-fitted meshes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

An immersed boundary method for fluid–structure interaction with compressible multiphase flows

This paper presents a two-dimensional immersed boundary method for fluid–structure interaction with compressible multiphase flows involving large structure deformations. This method involves three important parts: flow solver, structure solver and fluid–structure interaction coupling. In the flow solver, the compressible multiphase Navier–Stokes equations for ideal gases are solved by a finite difference method based on a staggered Cartesian mesh, where a fifth-order accuracy Weighted Essentially Non-Oscillation (WENO) scheme is used to handle spatial discretization of the convective term, a fourth-order central difference scheme is employed to discretize the viscous term, the third-order TVD Runge–Kutta scheme is used to discretize the temporal term, and the level-set method is adopted to capture the multi-material interface. In this work, the structure considered is a geometrically non-linear beam which is solved by using a finite element method based on the absolute nodal coordinate formulation (ANCF). The fluid dynamics and the structure motion are coupled in a partitioned iterative manner with a feedback penalty immersed boundary method where the flow dynamics is defined on a fixed Lagrangian grid and the structure dynamics is described on a global coordinate. We perform several validation cases (including fluid over a cylinder, structure dynamics, flow induced vibration of a flexible plate, deformation of a flexible panel induced by shock waves in a shock tube, an inclined flexible plate in a hypersonic flow, and shock-induced collapse of a cylindrical helium cavity in the air), and compare the results with experimental and other numerical data. The present results agree well with the published data and the current experiment. Finally, we further demonstrate the versatility of the present method by applying it to a flexible plate interacting with multiphase flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A new re-redistribution scheme for weighted state redistribution with adaptive mesh refinement

State redistribution (SRD) is a recently developed technique for stabilizing cut cells that result from finite-volume embedded boundary methods. SRD has been successfully applied to a variety of compressible and incompressible flow problems. When used in conjunction with adaptive mesh refinement (AMR), additional steps are needed to preserve the accuracy and conservation properties of the solution if the embedded boundary is not restricted to a single level of the mesh hierarchy. In this work, we extend the weighted state redistribution algorithm to cases where cut cells live at or near a coarse-fine interface within the domain. Here, we present numerical results that demonstrate that the algorithm is conservative when the coarse-fine interface intersects the embedded boundary. Additionally we compare the numerical solution of the Sod shock tube problem in an inclined cylinder with the analytic solution, and we compare the simulation of a shock hitting a cylindrical obstacle with experimental data. Finally we demonstrate the methodology for simulation of the multicomponent compressible Navier-Stokes equations in a piston-bowl geometry, and discuss the computational efficiency gained by not requiring the entire embedded boundary to be defined at the finest level.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗