Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variable Geometry”

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,621 records · Page 90

Experimental Investigation of Distributed Sand-Grain Roughness Effects on Transition Onset and Turbulent Heating at Mach 6

An experimental investigation of distributed sand-grain surface roughness effects on hypersonic boundary-layer transition and turbulent heating has been performed. Two representative entry vehicle geometries, a sphere-cone aeroshell and a spherical-cap aeroshell, were considered. Cast ceramic models of each geometry were fabricated with distributed sand-grain roughness patterns of different heights that simulated an ablated thermal protection system. Wind tunnel testing was performed at Mach 6 over a range of Reynolds numbers sufficient to produce laminar, transitional, and turbulent flow. Aeroheating and boundary-layer transition onset data were obtained using global phosphor thermography. The experimental heating data are presented herein, as are comparisons to laminar and turbulent smooth-wall heat transfer distributions from computational flow field simulations. The boundary-layer transition data were found to correlate with a functional representation developed in prior roughness studies, although the data scatter was greater owing to the height variability of distributed sand-grain roughness.

Brian R Hollis↗

An Overset Grid Motion Driver for Resolving Entry Vehicle CFD Simulations with Large Changes in Wake Orientation

Time-accurate computational fluid dynamics (CFD) simulations using overset meshes can enable reduced computational cost compared to a rigid, monolithic mesh for blunt-body atmospheric entry vehicles. If the vehicle has a significant change in attitude, the high-density region of the overset mesh designed for resolving the subsonic wake can be kept in position as the wake orientation changes relative to the body. Current simulations generally keep the wake mesh fixed in attitude while the body rotates, assuming the wake orientation remains relatively constant. However, changes in the relative velocity over a long-duration simulation or effects of the body geometry can invalidate this assumption. This work presents a method for calculating the desired mesh orientation to resolve a blunt body wake using the atmosphere-relative velocity of the vehicle. A motion driver using these calculations is implemented in a CFD-in-the-loop flight model bridging the FUN3D flow solver and POST2 trajectory propagator. Overset meshing is performed using the Yoga extension for FUN3D. Results are presented for simulations with both overset and monolithic meshes under large changes in relative velocity (and therefore wake orientation). The results demonstrate that the motion driver successfully orients an overset wake mesh to capture the subsonic wake with sufficient resolution despite variable body rotation and body-relative freestream velocity.

CFD↗

An Overset Grid Motion Driver for Resolving Entry Vehicle CFD Simulations with Large Changes in Wake Orientation

Time-accurate computational fluid dynamics (CFD) simulations using overset meshes can enable reduced computational cost compared to a rigid, monolithic mesh for blunt-body atmospheric entry vehicles. If the vehicle has a significant change in attitude, the high-density region of the overset mesh designed for resolving the subsonic wake can be kept in position as the wake orientation changes relative to the body. Current simulations generally keep the wake mesh fixed in attitude while the body rotates, assuming the wake orientation remains relatively constant. However, changes in the relative velocity over a long-duration simulation or effects of the body geometry can invalidate this assumption. This work presents a method for calculating the desired mesh orientation to resolve a blunt body wake using the atmosphere-relative velocity of the vehicle. A motion driver using these calculations is implemented in a CFD-in-the-loop flight model bridging the FUN3D flow solver and POST2 trajectory propagator. Overset meshing is performed using the Yoga extension for FUN3D. Results are presented for simulations with both overset and monolithic meshes under large changes in relative velocity (and therefore wake orientation). The results demonstrate that the motion driver successfully orients an overset wake mesh to capture the subsonic wake with sufficient resolution despite variable body rotation and body-relative freestream velocity.

Computational Fluid Dynamics↗

Superstratum symbiosis

Superstrata are smooth horizonless microstate geometries for the supersymmetric D1-D5-P black hole in type IIB supergravity. In the CFT, “superstratum states” are defined to be the component of the supergraviton gas that is obtained by breaking the CFT into “|00⟩-strands” and acting on each strand with the “small,” anomaly-free superconformal generators. We show that the recently-constructed supercharged superstrata represent a final and crucial component for the construction of the supergravity dual of a generic superstratum state and how the supergravity solution faithfully represents all the coherent superstratum states of the CFT. For the supergravity alone, this shows that generic superstrata do indeed fluctuate as functions of three independent variables. Smoothness of the complete supergravity solution also involves “coiffuring constraints” at second-order in the fluctuations and we describe how these lead to new predictions for three-point functions in the dual CFT. We use a hybrid of the original and supercharged superstrata to construct families of single-mode superstrata that still have free moduli after one has fixed the asymptotic charges of the system. We also study scalar wave perturbations in a particular family of such solutions and show that the mass gap depends on the free moduli. This can have interesting implications for superstrata at non-zero temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Metal additively manufactured wavy fin cold-plate architecture for improved thermal-hydraulic performance

Rapid growth in artificial intelligence and data center workloads demands high-performance liquid cooling to manage increasing chip power. This study presents two metal-additive-manufactured cold plates with sinusoidal fins, constant-amplitude wavy fins and linearly variable-amplitude wavy fins and compares them against metal-additive-manufactured straight fins using experiments conducted at 1 kW heat dissipation as well as high-fidelity 3D conjugate computational fluid dynamic simulations. The cold plates were printed in AlSi10Mg material and underwent design using a Python-automated workflow prior to manufacture and testing. The experiments show that wavy fins reduce the normalized thermal resistance by 35 to 45 % at water flow rates from 1 to 4 LPM. At a fixed 20 kPa pressure drop, the variable-waviness design lowered peak surface temperature by 9 °C and thermal resistance by 51 %, while edge-channel maldistribution in the constant wavy fin design limited gains. A thermal resistance breakdown revealed that 55–63 % of the total thermal resistance in wavy designs comes from base heat conduction, 27–33 % from fin heat conduction, and 9–13 % from fin heat convection, indicating the need to address conduction bottlenecks. Parametric sweeps identify a 3 mm fin pitch as optimal, and that horizontal inlet/outlet manifolds further reduce pressure drop by 30–60 % and thermal resistance by 9–16 % relative to vertical inlet-outlet manifolds. The results yield comprehensive guidelines for fin geometry, manifold alignment, material selection and additive-manufacturing constraints to realize high-performance liquid-cooled cold plates for power-dense electronics.

3d printing↗

Propagation of plane waves in a variable area duct carrying a compressible subsonic flow

A complete numerical solution is obtained for plane wave propagation in a variable-area duct which carries a mean compressible flow. The equation solved is a matrix formulation of the equations first derived by Tsien in 1952. Sample calculations are presented for four particular flow configurations: (1) a duct-nozzle-duct combination with the acoustic wave traveling with the flow, (2) a duct-nozzle-duct combination with the acoustic wave traveling against the flow, (3) an inlet-nozzle-duct configuration with the acoustic wave traveling against the flow, and (4) a duct-nozzle-exit configuration with the acoustic wave traveling with the flow. For each of these cases, the variation of energy flux with frequency is shown for particular geometry and flow condition.

Davis, S. S.↗

Numerical Investigation of Butterfly Valve Performance in Variable Valve Sizes, Positions and Flow Regimes

Reliability and efficiency of valves are necessary for precise control and sufficient heat-flow to heat application plants for the integrated energy systems of nuclear power plants (NPPs). Strategic Management Analysis Requirement and Technology (SMART) valves’ ability to control flow and assess environmental parameters stands out for these requirements. Their ability to sustain the downstream flow rate, prevent reverse flow, and maintain pressure in the heat transport loop is much more efficient with the integration of sensors and intelligent algorithms. For assessing valve performance and monitoring, mechanical design and operating conditions are two important parameters. In this study, the butterfly valves of three different sizes are simulated with water and steam using STAR-CCM+ in various flow regimes and positions to analyze performance parameters to strategize an automated control system for efficiently balancing the heat–transport network. Also, flow behavior is studied using velocity and pressure fields for valve–body geometry optimization. It can be observed, through performance parameters, that the valves are suitable for operation between 30° and 90° positions with significantly low loss coefficients and high flow coefficients, and the performance parameters follow a certain pattern in both water and steam flow in each scenario.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Design principles from multiscale simulations to predict nanostructure in self-assembling ionic liquids

Molecular dynamics simulations (up to the nanoscale) were performed on the 3-methyl-1-pentylimidazolium ionic liquid cation paired with three anions; chloride, nitrate, and thiocyanate as aqueous mixtures, using the effective fragment potential (EFP) method, a computationally inexpensive way of modeling intermolecular interactions. The simulations provided insight (preferred geometries, radial distribution functions and theoretical proton NMR resonances) into the interactions within the ionic domain and are validated against 1 H NMR spectroscopy and small- and wide-angle X-ray scattering experiments on 1-decyl-3-methylimidazolium. Ionic liquids containing thiocyanate typically resist gelation and form poorly ordered lamellar structures upon mixing with water. Conversely, chloride, a strongly coordinating anion, normally forms strong physical gels and produces well-ordered nanostructures adopting a variety of structural motifs over a very wide range of water compositions. Nitrate is intermediate in character, whereby upon dispersal in water it displays a range of viscosities and self-assembles into nanostructures with considerable variability in the fidelity of ordering and symmetry, as a function of water content in the binary mixtures. The observed changes in the macro and nanoscale characteristics were directly correlated to ionic domain structures and intermolecular interactions as theoretically predicted by the analysis of MD trajectories and calculated RDFs. Specifically, both chloride and nitrate are positioned in the plane of the cation. Anion to cation proximity is dependent on water content. Thiocyanate is more susceptible to water insertion into the second solvent shell. Experimental 1 H NMR chemical shifts monitor the site-specific competition dependence with water content in the binary mixtures. As a result, thiocyanate preferentially sits above and below the aromatic ring plane, a state disallowing interaction with the protons on the imidazolium ring.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Combustion instability analysis for liquid propellant rocket engines

The multi-dimensional numerical model has been developed to analyze the nonlinear combustion instabilities in liquid-fueled engines. The present pressure-based approach can handle the implicit pressure-velocity coupling in a non-iterative way. The additional scalar conservation equations for the chemical species, the energy, and the turbulent transport quantities can be handled by the same predictor-corrector sequences. This method is time-accurate and it can be applicable to the all-speed, transient, multi-phase, and reacting flows. Special emphasis is given to the acoustic/vaporization interaction which may act as the crucial rate-controlling mechanism in the liquid-fueled rocket engines. The subcritical vaporization is modeled to account for the effects of variable thermophysical properties, non-unitary Lewis number in the gas-film, the Stefan flow effect, and the effect of transient liquid heating. The test cases include the one-dimenisonal fast transient non-reacting and reacting flows, and the multi-dimensional combustion instabilities encountered in the liquid-fueled rocket thrust chamber. The present numerical model successfully demonstrated the capability to simulate the fast transient spray-combusting flows in terms of the limiting-cycle amplitude phenomena, correspondence between combustion and acoustics, and the steep-fronted wave and flame propagation. The investigated parameters include the spray initial conditions, air-fuel mixture ratios, and the engine geometry. Stable and unstable operating conditions are found for the liquid-fueled combustors. Under certain conditions, the limiting cycle behavior of the combusting flowfields is obtained. The numerical results indicate that the spray vaporization processes play an important role in releasing thermal energy and driving the combustion instability.

Kim, Y. M.↗

Magellan stereo images and Venusian geology

Areas of Venus imaged by Magellan radar with multiple viewing conditions provide unique data that will contribute to the solution of venusian geologic problems and provide a basis for quantitative comparison of venusian landforms with those on other planetary bodies. Three sets of images with different viewing conditions have been acquired: (1) left-looking with variable incidence angles (cycle 1 profile), (2) right-looking with nearly constant incidence angles (cycle 2 profile), and (3) left-looking with variable incidence angles that are almost always smaller than those in (1) (cycle 3 profiles). The unique data provided by paired images of the same scene with different incidence angles arises from image displacements caused by the relief of individual landforms at scales comparable to the ground-range and azimuth resolutions of the images. There are two aspects of the data: (1) Stereopsis achieved by simultaneous viewing of paired left-looking images of the same scene permits three-dimensional perception and interpretation of the morphologies of landforms at resolutions much finer than the altimetry footprints. (2) Measurements of differences of image displacements (parallax) on paired images with known imaging geometries provide quantitative estimates of the relief and shapes of landforms. The potential scientific contributions of the data can be grouped into two interrelated classes: (A) geologic mapping, analysis, and interpretation and (B) topical studies that involve topographic measurements. Stereopsis, without quantitative measurements, enhances geologic mapping, analysis, and interpretation of the rock units of Venus to a degree that cannot be overestimated. In geologic mapping, assemblages of landforms, assessments of backscatter and variations in backscatter, and fine-scale topography are used to define and characterize geologic map units that represent laterally continuous deposits or rock units. Stereopsis adds the important dimension of local relief for characterization of geologic units at a scale that is not possible with Magellan altimetry or products derived from the altimetry. Relative ages of the geologic units are determined using the well-known principles of superposition and intersection. Here, the perception of relief is invaluable because superposition relations among the geological units are more readily and clearly established. The recognition of folds, faults, and fault systems, regardless of their orientations, is facilitated with stereopsis so that sequences of deformation of the geologic units can be determined and structural analyses vastly improved. Shapes of landforms are readily perceived so that they can be properly interpreted. The end result of the mapping, analyses, and interpretations is a geologic history of Venus that includes the sequences of formation and deformation of various geologic units. Measurements of relief at the finest scale possible are necessary for numerous topical studies. Standard altimetry will provide the necessary information on the relief of most large landforms, but it tends to underestimate the relief of small landforms and distorts their shapes. Although special processing of the altimeter echoes improves the estimates of the relief and shapes of some landforms, there are uncertainties in the interpretations of the echoes. Examples of topical studies requiring measurements of relief are given.

Moore, H. J.↗

Bayesian Vision for Shape Recovery

We present a new Bayesian vision technique that aims at recovering a shape from two or more noisy observations taken under similar lighting conditions. The shape is parametrized by a piecewise linear height field, textured by a piecewise linear irradiance field, and we assume Gaussian Markovian priors for both shape vertices and irradiance variables. The observation process. also known as rendering, is modeled by a non-affine projection (e.g. perspective projection) followed by a convolution with a piecewise linear point spread function. and contamination by additive Gaussian noise. We assume that the observation parameters are calibrated beforehand. The major novelty of the proposed method consists of marginalizing out the irradiances considered as nuisance parameters, which is achieved by Laplace approximations. This reduces the inference to minimizing an energy that only depends on the shape vertices, and therefore allows an efficient Iterated Conditional Mode (ICM) optimization scheme to be implemented. A Gaussian approximation of the posterior shape density is computed, thus providing estimates both the geometry and its uncertainty. We illustrate the effectiveness of the new method by shape reconstruction results in a 2D case. A 3D version is currently under development and aims at recovering a surface from multiple images, reconstructing the topography by marginalizing out both albedo and shading.

Jalobeanu, Andre↗

Characterization of the Modal Characteristics of Structures Operating in Dense Liquid Oxygen Turbopumps

A number of valuable conclusions can be drawn from this study. First, knockdown factors for a specific fluid are not constant but instead are dependent on the mode shape, although the largest this variability gets is about 10% for LOX, the densest fluid. The factors decrease the most for lower frequency shapes and less for higher ones. It follows, therefore, that mode number mismatch between air and fluid operation becomes not only possible, but common, as a knockdown factor for a particular mode shape may be higher than for another mode shape. Since this is a function of added mass, the mismatch is more prevalent for higher density fluids, but it initiates even for very low density ones. Another important conclusion reached is that it appears that the basic mode shapes of a structure do not change if it is fully symmetric, which includes its geometry and boundary conditions. There is some indication of small changes in the relative magnitudes within the mode shape. This conclusion is evident in the results from the cantilever rectangular plate and the inducer, which are not symmetric, and the fixed-fixed plate and the annular disk, which are. For non-symmetric structures, though, the mode shapes almost universally change for dense fluids, as shown by the very low MAC calculations. For the inducer in particular, the changes follow a trend of reduced parabolic and sine wavelengths with increasing density. It is critical to recognize the change in mode shape for several reasons. First, model updating with modal test becomes problematic if the shapes change. Second, design to avoid resonance is highly critical on the mode shape for modes other than the primary ones, as resonance is only a factor when the excitation shape matches the mode shape. Finally, application of the modal superposition method of forced response analysis is dependent on the use of accurate mode shapes. A more-refined assessment of the "knockdown" factor values and ranges than any previously reported in the literature for a realistic engineering structure is also presented in this paper. This data is of tremendous benefit for preliminary analysis and design, where a quick estimate is necessary. These results are important not just for rocket engine turbomachinery, but for water pumps and turbines, propellers, and any other structure operating in a heavy fluid with dynamic excitation. The clear avenue for future work for this endeavor is to expand the analytical techniques discussed in the literature to develop analytical expressions and justification for the mode shape changes and associated frequency knockdowns. These expressions must be able to accurately predict the functional relationship to the shapes, which will enable accurate tracing of the mode number from vacuum analysis (or testing in air) to analysis and operation in the intended fluid environment.

Chiu, Joseph↗

FUN3D Analyses in Support of the 1st AIAA Stability and Control Prediction Workshop

The 1st AIAA Stability and Control Prediction Workshop was created to establish best practices for the prediction of stability & control derivatives using computational fluid dynamics and assess the limitations of these computational methods when those best practices are applied. The inaugural workshop considers the ONERA version of the NASA/Boeing Common Research Model(CRM), which includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. Wind tunnel tests have been conducted for this configuration with longitudinal tests having been previously published, in addition to unpublished data at small sideslip angles that will serve as ‘blind’ data for workshop data comparisons. Participants were provided a ‘family’ of unstructured grids for the full-span ONERA CRM model with the wind-tunnel sting included. This family of mixed-element grids consists of 5 levels of refinement (tiny, coarse, medium, fine, and extra fine) with surface and volume mesh scaling, resulting in a size range of 14.6 to 53.4 million nodes. In addition, a medium refinement mesh has been provided for the ONERA CRM configuration without a sting to evaluate the sting’s impact on static longitudinal stability characteristics. In addition to these workshop-provided grids, the present work also considers an equivalent ‘family’ of computational grids generated using Heldenmesh™, a rapid grid generation software by Helden Aerospace Corporation for creating high-quality, three-dimensional, mixed-element unstructured meshes. Because of the authors’ familiarity with this software, these additional grids were generated as a comparison to the workshop-provided grids and to better understand the implications of using volume-mirrored grids for stability and control predictions. The present work will contribute to the workshop with test case data generated using the NASA FUN3D code, a parallelized, unstructured, node-based, finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. Numerical simulations will be conducted using the Quadratic Constitutive Relationship (QCR) version of the Spalart-Allmaras (SA) turbulence model with negative turbulence variable provisions. Both steady and2nd-order, time-accurate simulation results are to be generated and compared for select test cases, as time permits, to investigate their impact on FUN3D predictions. The present work will consider the three primary workshop test cases: (1) grid convergence study, (2) Mach number effect on static stability, and (3) wind tunnel sting increments. Additionally, data will be provided for the two optional test cases, which include:(1) static stability derivative calculations and (2)sideslip angle sweeps. In each of the test cases, the vehicle is stationary, and the body is assumed to be rigid, where vehicle deformation has been accounted for in the model configuration geometry. For all test cases, longitudinal and lateral force and moment aerodynamic coefficients will be provided for the total configuration, in addition to a component-level breakdown that includes the port wing, starboard wing, fuselage, and tail.

CFD↗

Feature Based Qualification of 17-4PH Stainless Steel to Evaluate Location-Specific Variability in Wire Arc Additive Manufacturing

Qualifying large-scale metal additive manufacturing (M-AM) technologies such as wire arc additive manufacturing (WAAM) can be challenging. This is especially significant in precipitation hardened martensitic stainless steels like SS 17-4PH, where thermal histories induce location-specific microstructural variability and property anisotropy. The Department of Defense (DOD) and the United States Army Combat Capabilities Development Command Ground Vehicle Systems Center (GVSC) Ground Vehicle Materials Engineering (GVME) aim to build robust and qualified large-scale M-AM workflows that could reduce the time and cost through quick and informed evaluation, testing, and development of feedstock, processes, and parts. The report presents the findings from the collaborative efforts between Oak Ridge National Laboratory (ORNL) and the U.S. Army GVSC GVME. The aim of this project was to develop a geometric feature-based qualification framework for WAAM of SS 17-4PH components. This report outlines selection methodology of representative build geometries, optimization of WAAM process parameters, in-situ monitoring, microstructure-property evaluation, thermal simulations, as well as data visualization techniques incorporated in this project. The results from this project demonstrate a clear understanding of thermal history dependent phase evolution and consequent location-specific property variations in WAAM of SS 17-4PH. These results in conjunction with the data-driven methodologies used in this project are expected to reduce qualification timelines, improve predictability, and accelerate the development of reliable feature-based qualification strategies for part production via large-scale M-AM technologies.

36 MATERIALS SCIENCE↗

Assessment of a parabolic analysis for axisymmetric internal flows in rocket and turbomachinery ducts

The flow paths in gas turbine passages encompass a wide range of flow properties such as Reynolds number and Mach number as well as many other variable flow conditions such as swirl, free-stream turbulence, and laminar/turbulent transition. An existing computer program, the Axisymmetric Diffuser Duct (ADD) code, which calculates compressible, turbulent, swirling flow through axisymmetric ducts has been modified to include the effects of free-stream turbulence and laminar/turbulent transition. The program has been evaluated on a matrix of test cases to determine its accuracy, robustness, and limits of applicability. This improved version of the ADD code calculates solutions which compare well with available data and can now be applied to a wider range of problems than previously possible. In addition a concept called local enhancement was developed and tested on a simple two-dimensional geometry in order to demonstrate a method to reduce computer time. In this concept the pressure distribution is calculated on a coarse grid using the ADD code and the viscous layer is locally enhanced using a boundary layer analysis. By applying this concept, an order of magnitude reduction in computer time was possible without any loss in accuracy.

Power, G. D.↗

On the Feasibility of Monitoring Carbon Monoxide in the Lower Troposphere from a Constellation of Northern Hemisphere Geostationary Satellites (PART 1)

By the end of the current decade, there are plans to deploy several geostationary Earth orbit (GEO) satellite missions for atmospheric composition over North America, East Asia and Europe with additional missions proposed. Together, these present the possibility of a constellation of geostationary platforms to achieve continuous time-resolved high-density observations over continental domains for mapping pollutant sources and variability at diurnal and local scales. In this paper, we use a novel approach to sample a very high global resolution model (GEOS-5 at 7 km horizontal resolution) to produce a dataset of synthetic carbon monoxide pollution observations representative of those potentially obtainable from a GEO satellite constellation with predicted measurement sensitivities based on current remote sensing capabilities. Part 1 of this study focuses on the production of simulated synthetic measurements for air quality OSSEs (Observing System Simulation Experiments). We simulate carbon monoxide nadir retrievals using a technique that provides realistic measurements with very low computational cost. We discuss the sampling methodology: the projection of footprints and areas of regard for geostationary geometries over each of the North America, East Asia and Europe regions; the regression method to simulate measurement sensitivity; and the measurement error simulation. A detailed analysis of the simulated observation sensitivity is performed, and limitations of the method are discussed. We also describe impacts from clouds, showing that the efficiency of an instrument making atmospheric composition measurements on a geostationary platform is dependent on the dominant weather regime over a given region and the pixel size resolution. These results demonstrate the viability of the "instrument simulator" step for an OSSE to assess the performance of a constellation of geostationary satellites for air quality measurements.

GEOS-5↗

A Precessing Be Disc as a Possible Model for Occultation Events in GX 304-1

We report on the RXTE detection of a sudden increase in the absorption column density, NH, during the 2011 May outburst of GX 304−1. The NH increased up to ∼16 × 10(exp 22) atoms cm(exp −2), which is a factor of 3–4 larger than what is usually measured during the outbursts of GX 304−1 as covered by RXTE. Additionally, an increase in the variability of the hardness ratio as calculated from the energy resolved RXTE-Proportional Counter Array light curves is measured during this time range. We interpret these facts as an occultation event of the neutron star by material in the line of sight. Using a simple 3D model of an inclined and precessing Be disc around the Be-type companion, we are able to qualitatively explain the NH evolution over time. We are able to constrain the Be disc density to be of the order of 10(exp −11)g cm(exp −3). Our model strengthens the idea of inclined Be discs as origin of double-peaked outbursts as the derived geometry allows accretion twice per orbit under certain conditions.

Kühnel, M.↗

Computational design of microarchitected porous electrodes for redox flow batteries

Porous electrodes are used as the core reactive component across electrochemical technologies. In flowing systems, controlling the fluid distribution, species transport, and reactive environment is critical to attaining high performance. However, conventional electrode materials like felts and papers provide few opportunities for precise engineering of the electrode and its microstructure. To address these limitations, architected electrodes composed of unit cells with spatially varying geometry determined via computational optimization are proposed. Resolved simulation is employed to develop a homogenized description of the constituent unit cells. These effective properties serve as inputs to a continuum model for the electrode when used in the negative half-cell of a vanadium redox flow battery. Porosity distributions minimizing power loss are then determined via computational design optimization to generate architected porosity electrodes. The architected electrodes are compared to bulk, uniform porosity electrodes and found to lead to increased power efficiency across operating flow rates and currents. The design methodology is further used to generate a scaled-up electrode with comparable power efficiency to the bench-scale systems. Finally, the variable porosity architecture and computational design methodology presented here thus offers a novel pathway for automatically generating spatially engineered electrode structures with improved power performance.

25 ENERGY STORAGE↗