Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Refining”

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 829 records · Page 46

Statistical Analysis of CFD Solutions from the Third AIAA Drag Prediction Workshop

The first AIAA Drag Prediction Workshop, held in June 2001, evaluated the results from an extensive N-version test of a collection of Reynolds-Averaged Navier-Stokes CFD codes. The code-to-code scatter was more than an order of magnitude larger than desired for design and experimental validation of cruise conditions for a subsonic transport configuration. The second AIAA Drag Prediction Workshop, held in June 2003, emphasized the determination of installed pylon-nacelle drag increments and grid refinement studies. The code-to-code scatter was significantly reduced compared to the first DPW, but still larger than desired. However, grid refinement studies showed no significant improvement in code-to-code scatter with increasing grid refinement. The third Drag Prediction Workshop focused on the determination of installed side-of-body fairing drag increments and grid refinement studies for clean attached flow on wing alone configurations and for separated flow on the DLR-F6 subsonic transport model. This work evaluated the effect of grid refinement on the code-to-code scatter for the clean attached flow test cases and the separated flow test cases.

Statistical analysis↗

Simulation of an Isolated Tiltrotor in Hover with an Unstructured Overset-Grid RANS Solver

An unstructured overset-grid Reynolds Averaged Navier-Stokes (RANS) solver, FUN3D, is used to simulate an isolated tiltrotor in hover. An overview of the computational method is presented as well as the details of the overset-grid systems. Steady-state computations within a noninertial reference frame define the performance trends of the rotor across a range of the experimental collective settings. Results are presented to show the effects of off-body grid refinement and blade grid refinement. The computed performance and blade loading trends show good agreement with experimental results and previously published structured overset-grid computations. Off-body flow features indicate a significant improvement in the resolution of the first perpendicular blade vortex interaction with background grid refinement across the collective range. Considering experimental data uncertainty and effects of transition, the prediction of figure of merit on the baseline and refined grid is reasonable at the higher collective range- within 3 percent of the measured values. At the lower collective settings, the computed figure of merit is approximately 6 percent lower than the experimental data. A comparison of steady and unsteady results show that with temporal refinement, the dynamic results closely match the steady-state noninertial results which gives confidence in the accuracy of the dynamic overset-grid approach.

Lee-Rausch, Elizabeth M.↗

Best Merge Region Growing Segmentation with Integrated Non-Adjacent Region Object Aggregation

Best merge region growing normally produces segmentations with closed connected region objects. Recognizing that spectrally similar objects often appear in spatially separate locations, we present an approach for tightly integrating best merge region growing with non-adjacent region object aggregation, which we call Hierarchical Segmentation or HSeg. However, the original implementation of non-adjacent region object aggregation in HSeg required excessive computing time even for moderately sized images because of the required intercomparison of each region with all other regions. This problem was previously addressed by a recursive approximation of HSeg, called RHSeg. In this paper we introduce a refined implementation of non-adjacent region object aggregation in HSeg that reduces the computational requirements of HSeg without resorting to the recursive approximation. In this refinement, HSeg s region inter-comparisons among non-adjacent regions are limited to regions of a dynamically determined minimum size. We show that this refined version of HSeg can process moderately sized images in about the same amount of time as RHSeg incorporating the original HSeg. Nonetheless, RHSeg is still required for processing very large images due to its lower computer memory requirements and amenability to parallel processing. We then note a limitation of RHSeg with the original HSeg for high spatial resolution images, and show how incorporating the refined HSeg into RHSeg overcomes this limitation. The quality of the image segmentations produced by the refined HSeg is then compared with other available best merge segmentation approaches. Finally, we comment on the unique nature of the hierarchical segmentations produced by HSeg.

Tilton, James C.↗

Numerical Simulations of a Quiet SuperSonic Technology (QueSST) Aircraft Preliminary Design

Reynolds Averaged Navier-Stokes (RANS) simulations were performed on a Lockheed Martin Quiet SuperSonic Technology (QueSST) aircraft preliminary design to assess inlet performance. The FUN3D flow solver and its adjoint-based grid refinement capability were used for the simulations in hopes of determining internal "best practices" for predicting inlet performance on top-aft-mounted inlets. Several parameters were explored including tetrahedral vs. pentahedral cells in/around the boundary-layer regions, an engine axis-aligned linear pressure sensor vs. a pressure box objective as the grid adaptation metric, and the number of grid adaptation cycles performed. Additional simulations were performed on manually refined grids for comparison with the adjoint-based adapted grids. Results showed poor agreement in predicted inlet performance on the refined grids compared to experimental data. This was true regardless of whether the refinement was adjoint-based or manual, the cell type in/near the boundary-layer regions, or the grid adaptation metric used. In addition, the 40-probe total pressure recovery was shown to decrease asymptotically as the number of adaptation cycles is increased. Solutions on the unadapted grids generally had better agreement with experimental data than their refined grid counterparts.

Inlet Performance↗

MPAS - Ocean Simulation Quality for Variable-Resolution North American Coastal Meshes

Climate model components utilizing unstructured meshes enable variableresolution, regionally enhanced simulations within global domains. Here we investigate the relationship between mesh quality and simulation statistics using the JIGSAW unstructured meshing library and the Model for Prediction Across ScalesOcean (MPASOcean) with a focus on Gulf Stream dynamics. In the base configuration, the refined region employs 8 km cells that extend 400 km from the coast of North America. This coastal refined region is embedded within a lowresolution global domain, with cell size varying latitudinally between 30 and 60 km. The resolution transition region between the refined region and background mesh is 600 km wide. Three sensitivity tests are conducted: 1) the quality of meshes is intentionally degraded so that horizontal cells are progressively more distorted; 2) the transition region from high to low resolution is steepened; and 3) resolution of the coastal refinement region is varied from 30 km to 8 km. Overall, the ocean simulations are shown to be robust to mesh resolution and quality alterations. Meshes that are substantially degraded still produce realistic currents, with Southern Ocean transports within 0.4% and Gulf Stream transports within 12% of highquality mesh results. The narrowest transition case of 100 km did not produce any spurious effects. Refined regions with high resolution produce eddy kinetic energy and sea surface height variability that are similar to the highresolution reference simulation. These results provide heuristics for the design criteria of variableresolution climate model domains.

Hoch, Kristin E.↗

TPSAS-NF1676L-13579-DND

During aerospace vehicle conceptual and preliminary design, empirical non-optimum factors are typically applied to predicted structural component weights to account for undefined manufacturing and design details. Non-optimum factors are developed here for 32 aluminum-lithium 2195 orthogrid panels comprising the liquid hydrogen tank barrel of the Space Shuttle External Tank using measured panel weights and manufacturing drawings. Minimum values for skin thickness, axial and circumferential blade stiffener thickness and spacing, and overall panel thickness are used to estimate individual panel weights. Panel non-optimum factors computed using a coarse weights model range from 1.21 to 1.77, and a refined weights model (including weld lands and skin and stiffener transition details) yields non-optimum factors of between 1.02 and 1.54. Acreage panels have an average 1.24 non-optimum factor using the coarse model, and 1.03 with the refined version. Non-optimum factors are also calculated for the aluminum-alloy gores that compose the ogive and ellipsoidal bulkheads of the Space Shuttle Super-Lightweight Tank propellant tanks. Minimum values for actual gore skin thicknesses and weld land dimensions are extracted from component production drawings, and are used to predict reference gore weights. These actual skin thicknesses are also compared to skin thicknesses predicted using classical structural mechanics and tank proof-test pressures. Both coarse and refined weights models are developed for the gores. The coarse model is based on the proof-test skin thicknesses, and the refined model uses the actual skin thicknesses and weld land details. Both models include manufacturing tolerance weight estimates. To determine the gore non-optimum factors, these reference weights are then compared to flight hardware weights reported in a mass properties database. The gore non-optimum factors computed using the coarse weights model range from 1.18 to 2.36, with an average non-optimum factor of 1.68. Application of the refined weights model yields non-optimum factors between 1.00 and 1.45, with an average non-optimum factor of 1.13. To demonstrate their utility, these calculated non-optimum factors are then used to predict more realistic (i.e., heavier) gore weights for a proposed heavy-lift launch vehicle's propellant tank bulkheads.

K Chauncey Wu↗

Grid-Stretching Capability for the GEOS-Chem 13.0.0 Atmospheric Chemistry Model

Modeling atmospheric chemistry at fine resolution globally is computationally expensive; the capability to focus on specific geographic regions using a multi scale grid is desirable. Here, we develop, validate, and demonstrate stretched grids in the GEOS-Chem atmospheric chemistry model in its high-performance implementation (GCHP). These multiscale grids are specified at runtime by four parameters that offer users nimble control of the region that is refined and the resolution of the refinement. We validate the stretched-grid simulation versus global cubed-sphere simulations. We demonstrate the operation and flexibility of stretched-grid simulations with two case studies that compare simulated tropospheric NO2column densities from stretched-grid and cubed-sphere simulations to retrieved column densities from the TROPOspheric Monitoring Instrument (TROPOMI). The first case study uses a stretched grid with a broad refinement covering the contiguous US to produces imulated columns that perform similarly to a C180 (∼50 km) cubed-sphere simulation at less than one-ninth the computational expense. The second case study experiments with a large stretch-factor for a global stretched-grid simulation with a highly localized refinement with∼10 km resolution for California. We find that the refinement improves spatial agreement with TROPOMI columns compared to a C90 cubed-sphere simulation of comparable computational demands. Overall we find that stretched grids in GEOS-Chem are a practical tool for fine resolution regional- or continental-scale simulations of atmospheric chemistry. Stretched grids are available in GEOS-Chem version 13.0.0

GEOS-Chem↗

VIPER Lunar Rover Agile Mission Systems

Agile development methods, which have gone from outlier to mainstream in software development, are poised to expand into all aspects of space mission development. Modern software development operates on a principle of continuous deployment, where progress is verified not with conventional metrics, but with a continuous build, available to key stakeholders, enabling direct examination of the state of the code base, and assessment of progress through demonstration of capability. Delivery times are measured in weeks, not months. Stakeholders are part of the process on an ongoing basis. The cost of change is comparatively low and requirements, which often are not precisely defined at the start of a project, may be iteratively refined in a series of agile development cycles. Agile methods are compatible with traditional system engineering methods and may be tailored to the space operations environment. The low cost of change and iterative development cycles of agile enable requirements to be defined as outcomes and constraints, with design details to be refined during the development cycle. We are now at a point where agile methods may be extended beyond software, to Mission Systems, including the Mission Operations System and the Ground Data System. For NASA’s VIPER Lunar Rover Mission, scheduled to land at a lunar pole in late 2023, we are developing the Mission System using agile methods. As in agile software, where the measure of progress is working code, in agile mission system development, the measure of capability is what we can demonstrate. Demonstrations over presentations. We demonstrate mission system capability using simulations. The concept of operations, from commanding, to driving the rover, to how we downlink images for evaluation for a near-real time command cycle, will be tested and proven in simulation, years before we begin the traditional simulation cycle for training. “Say it then simulate it.” We develop and refine our designs using simulations, with an emphasis on new components of the system that are not well known early. For example, the required duration of a mission planning cycle for a lunar surface asset such as VIPER, that operates twenty-four hours a day, seven days a week, with continuous communications and a unique set of constraints based on the physics of the lunar poles and the line of site to Earth, is a unique problem in mission planning that is unlikely to be solved in a series of meetings. A small number of requirements specifying the outcomes may serve as the jumping off point to an agile development cycle, with demonstration in simulations. We have already demonstrated this process with simulations of rover driver decision time. VIPER is driven using near-real time command and control to waypoints. The driver decision time between waypoints is a fundamental enabling unit of productivity to accomplish the mission timeline. We have validated driver decision time in simulations of rover driving at the lunar South Pole, using the prototype mission tools for driving, command and control. The capability to develop and refine designs using simulations as part of agile Mission System development cycle changes the nature of team interactions, creating a focus on doing, rather than analyzing and documenting. Waterfall development cycles were, in part, a product of the significant cost of change in the early days of spaceflight. When the cost of change is high, it is vital to get your requirements right at the outset, because the system will be built to those specifications, and, when change is expensive, you better get it right early. However, modern technology has greatly lowered the cost of change, enabling iterative, rapid development cycles, in which key operations concepts may be tested and refined during development. Extending agile development to the Mission System for VIPER is a significant step in moving agile development methods for space operations beyond software, to the Mission System.

Agile↗

Measurement of Surface Deformation Related to the December 2018 Mt. Etna Eruption Using Time-Series Interferometry and Magma Modeling for Hazard Zone Mapping

Mount Etna has erupted several times since it was first formed. Recently, Mount Etna began erupting again over 24–27 December 2018. Because it erupts frequently, Mount Etna should be observed on a frequent basis. From June 2018 to October 2019, 34 and 56 interferometric synthetic aperture radar (InSAR) images were acquired from the ascending and descending tracks of the Sentinel-1 satellite, respectively. We employed the Stanford Method for Persistent Scatterers (StaMPS) and a refined small baseline subset (SBAS) InSAR method to produce a surface deformation time-series map. In the time-series analysis, the phase signal remained unaltered with time. The Okada model was then applied to the result to generate a modeled interferogram, and the Q-LavHA program was run to generate a lava flow prediction model. A direct comparison of the results showed that Persistent Scatterers Interferometry (PSI)-StaMPS and the refined SBAS technique were comparable in terms of the displacement pattern, with slightly different velocity values obtained for individual points. In particular, a velocity range of −25 to 21 cm/yr was obtained from PSI-StaMPS, whereas a range of −30 to 25 cm/yr was obtained from the refined SBAS method. Upon computation of the vertical and east-west displacement components based on ascending and descending track data using both methods, deformation velocities of 51.5 and 52.5 cm/yr in the westerly direction on the western flank of Mount Etna were obtained from PSI-StaMPS and the refined SBAS method, respectively, whereas on the eastern flank, deformation toward the east was estimated to occur at a velocity of 50.1 or 54.2 cm/yr, respectively. PSIStaMPS estimated a vertical deformation velocity of −5.3 to 18.3 cm/yr, whereas the refined SBAS method produced a velocity range of approximately −7 to 19 cm/yr. The interferogram obtained via Okada modeling showed two fault sources in the 2018 Mount Etna eruption and a total volume change of approximately 12.39 × 10 6 m 3 . From the modeling results, a lava flow prediction model was generated using the Q-LavHA program. The approaches described in this study can be used by government officials, authorities, and other decision-makers to monitor and assess the risk of volcanic activity in the region.

Suci Ramayanti↗

Investigation of CAD-based Geometry Workflows for Multiphysics Fusion Problems Using OpenMC and MOOSE

Fusion system designs are complex and require intricate and accurate meshes to be properly modeled. In this study, we investigate the use of CAD-based geometry workflows in fusion systems multiphysics problems. A simplified tokamak was introduced and modeled in CAD using a multiphysics coupling of OpenMC Monte Carlo transport and MOOSE heat conduction. The meshed geometry was prepared using direct accelerated geometry Monte Carlo (DAGMC) for particle transport, and a volumetric mesh was also prepared to be used in MOOSE and to tally OpenMC results. Cardinal was used to run OpenMC Monte Carlo particle transport within MOOSE framework. The heat source distribution and tritium production were calculated in OpenMC. The data transfer system was used to transfer heat source and temperature distribution between OpenMC and MOOSE. Two computational studies related to mesh refinement were performed: (1) refining the DAGMC and volumetric meshes used for tallying results and solving heat conduction and (2) only refining the DAGMC particle transport mesh. The refinement of the tally mesh has a much larger effect on the runtime compared to the refinement of the DAGMC particle transport surface mesh.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dyn$\mathrm{AMO}$: Multi-agent reinforcement learning for dynamic anticipatory mesh optimization with applications to hyperbolic conservation laws

Here we introduce DynAMO, a reinforcement learning paradigm for Dynamic Anticipatory Mesh Optimization. Adaptive mesh refinement is an effective tool for optimizing computational cost and solution accuracy in numerical methods for partial differential equations. However, traditional adaptive mesh refinement approaches for time-dependent problems typically rely only on instantaneous error indicators to guide adaptivity. As a result, standard strategies often require frequent remeshing to maintain accuracy. In the DynAMO approach, multi-agent reinforcement learning is used to discover new local refinement policies that can anticipate and respond to future solution states by producing meshes that deliver more accurate solutions for longer time intervals. By applying DynAMO to discontinuous Galerkin methods for the linear advection and compressible Euler equations in two dimensions, we demonstrate that this new mesh refinement paradigm can outperform conventional threshold-based strategies while also generalizing to different mesh sizes, remeshing and simulation times, and initial conditions.

97 MATHEMATICS AND COMPUTING↗

Localization of proteins involved in the biogenesis and repair of the photosynthetic apparatus to thylakoid subdomains in Arabidopsis

Abstract Thylakoid membranes in chloroplasts and cyanobacteria harbor the multisubunit protein complexes that catalyze the light reactions of photosynthesis. In plant chloroplasts, the thylakoid membrane system comprises a highly organized network with several subcompartments that differ in composition and morphology: grana stacks, unstacked stromal lamellae, and grana margins at the interface between stacked and unstacked regions. The localization of components of the photosynthetic apparatus among these subcompartments has been well characterized. However, less is known about the localization of proteins involved in the biogenesis and repair of the photosynthetic apparatus, the partitioning of proteins between two recently resolved components of the traditional margin fraction (refined margins and curvature), and the effects of light on these features. In this study, we analyzed the partitioning of numerous thylakoid biogenesis and repair factors among grana, curvature, refined margin, and stromal lamellae fractions of Arabidopsis thylakoid membranes, comparing the results from illuminated and dark‐adapted plants. Several proteins previously shown to localize to a margin fraction partitioned in varying ways among the resolved curvature and refined margin fractions. For example, the ALB3 insertase and FtsH protease involved in photosystem II (PSII) repair were concentrated in the refined margin fraction, whereas TAT translocon subunits and proteins involved in early steps in photosystem assembly were concentrated in the curvature fraction. By contrast, two photosystem assembly factors that facilitate late assembly steps were depleted from the curvature fraction. The enrichment of the PSII subunit OE23/PsbP in the curvature fraction set it apart from other PSII subunits, supporting the previous conjecture that OE23/PsbP assists in PSII biogenesis and/or repair. The PSII assembly factor PAM68 partitioned differently among thylakoid fractions from dark‐adapted plants and illuminated plants and was the only analyzed protein to convincingly do so. These results demonstrate an unanticipated spatial heterogeneity of photosystem biogenesis and repair functions in thylakoid membranes and reveal the curvature fraction to be a focal point of early photosystem biogenesis.

Chotewutmontri, Prakitchai↗

hr-Adaptivity for nonconforming high-order meshes with the target matrix optimization paradigm

In this study, we present an \(hr\) -adaptivity framework for optimization of high-order meshes. This work extends the r -adaptivity method by Dobrev et al. (Comput Fluids, 2020), where we utilized the Target-Matrix Optimization Paradigm (TMOP) to minimize a functional that depends on each element’s current and target geometric parameters: element aspect-ratio , size , skew , and rotation . Since fixed mesh topology limits the ability to achieve the target size and aspect-ratio at each position, in this paper, we augment the r -adaptivity framework with nonconforming adaptive mesh refinement to further reduce the error with respect to the target geometric parameters. The proposed formulation, referred to as \(hr\) -adaptivity, introduces TMOP-based quality estimators to satisfy the aspect-ratio target via anisotropic refinements and size target via isotropic refinements in each element of the mesh. The methodology presented is purely algebraic, extends to both simplices and hexahedra/quadrilaterals of any order, and supports nonconforming isotropic and anisotropic refinements in 2D and 3D. Using a problem with a known exact solution, we demonstrate the effectiveness of \(hr\) -adaptivity over both r - and \(h\) -adaptivity in obtaining similar accuracy in the solution with significantly fewer mesh nodes. We also present several examples that show that \(hr\) -adaptivity can help satisfy geometric targets even when \(r\) -adaptivity fails to do so, due to the topology of the initial mesh.

42 ENGINEERING↗

Extended isogeometric analysis of multi-material and multi-physics problems using hierarchical B-splines

Here, this paper presents an immersed, isogeometric finite element framework to predict the response of multi-material, multi-physics problems with complex geometries using locally refined discretizations. To circumvent the need to generate conformal meshes, this work uses an extended finite element method (XFEM) to discretize the governing equations on non-conforming, embedding meshes. A flexible approach to create truncated hierarchical B-splines discretizations is presented. This approach enables the refinement of each state variable field individually to meet field-specific accuracy requirements. To obtain an immersed geometry representation that is consistent across all hierarchically refined B-spline discretizations, the geometry is immersed into a single mesh, the XFEM background mesh, which is constructed from the union of all hierarchical B-spline meshes. An extraction operator is introduced to represent the truncated hierarchical B-spline bases in terms of Lagrange shape functions on the XFEM background mesh without loss of accuracy. The truncated hierarchical B-spline bases are enriched using a generalized Heaviside enrichment strategy to accommodate small geometric features and multi-material problems. The governing equations are augmented by a formulation of the face-oriented ghost stabilization enhanced for locally refined B-spline bases. We present examples for two- and three-dimensional linear elastic and thermo-elastic problems. The numerical results validate the accuracy of our framework. The results also demonstrate the applicability of the proposed framework to large, geometrically complex problems.

42 ENGINEERING↗

HHL algorithm with mapping function and enhanced sampling for model predictive control in microgrids

Here, this paper presents a refined quantum Harrow Hassidim Lloyd (HHL) algorithm for microgrid control. The first novelty of the developed method is that a mapping shift function enables the original HHL algorithm to handle general linear equations with non-singular and indefinite matrix. Second, a method of Matrix Extension for Amplifying Sampling Probabilities of Intended Solution (ME-ASPI) is proposed to design the reformulated linear algebraic equations, allowing for improved sampling efficiency of the quantum tomography in the refined HHL algorithm. Then, we applied the method to solve the model predictive control (MPC) problem in nonlinear dynamical microgrids. Specifically, with the ME-ASPI method, the refined HHL algorithm can effectively obtain the intended partial optimal control inputs for MPC. The optimization of quadratic programming problem in each time step of MPC is transformed into a linear system problem, which is addressed by the proposed quantum solver through using only partial information, with the time complexity improved from $\mathscr{O}(\mathscr{N}^{2.37286})$ classically to $\mathscr{O}(\mathscr{N}^{2} log \mathscr{N}$ x $p$ log $p)$ in quantum. Numerical examples have validated the effectiveness of the refined HHL algorithm with the proposed mapping function and the ME-ASPI method. By leveraging quantum properties, the proposed method provides a hybrid quantum–classical framework for microgrid control. This generic method can also potentially tackle many other challenges in analyzing and controlling general complex engineered systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Mass transport in a highly immiscible alloy on extended shear deformation

Forced mixing to a single phase or supersaturated solid solution (SSS) and its prerequisite microstructure evolution in immiscible systems has been a focus of research for both fundamental science and a variety of applications. Controlling the formation of SSS by shear deformation assisted processing could enable a material design beyond conventional equilibrium microstructure in immiscible systems. Here, a highly immiscible (mixing enthalpy of ~ 20 kJ·mol -1 ) Cu-50 at. % Cr binary alloy was employed to investigate the microstructure evolution and localized tendencies of SSS during severe shear deformation. The present results demonstrate distinctive microstructure refinement process in each phase and how they lead to localized SSS as a function of shear strain. Preferential dynamic recrystallization occurs in the softer Cu phase due to strain localization, leading to substantial grain refinement. The refinement of Cr phase in the top-most layer, however, is enabled by the progressive evolution of grain lamination, splitting, and fragmentation as a function of shear strain. The eventual SSS is found to be strongly dependent on the local environments that affect the dislocation activity, including the level of microstructure refinement, the interfacial orientation relationship, hardness difference, and supposed stability of oxidation. Ab initio simulations comfirm that it is more favorable to oxidize Cr than Cu at incoherent Cu/Cr interfaces which then limits the mass transport on an incoherent boundary. Our results shed light on the underpinning mechanism for non-equilibrium mass transport in immiscible systems upon severe deformation that can be applicable to a variety of processing techniques aimed at producing immiscible alloys with superior mechanical properties.

36 MATERIALS SCIENCE↗

Synchrotron X-ray diffraction of pyrolusite (MnO 2 ) and rutile (TiO 2 ) during heating to ~1000 °C

To compare thermal expansion behaviors in isomorphic structures with different redox behaviors, we separately heated pyrolusite (Mn 4+ O 2 ) and rutile (Ti 4+ O 2 ) powders from 25 to ~1000 °C at a rate of 0.05 °C/s at ambient pressure and analyzed them using synchrotron X-ray diffraction and Rietveld refinement. The pyrolusite exhibited two reduction-induced phase decompositions. Between 531 and 583 °C, pyrolusite decomposed to bixbyite (Mn 3+ 2 O 3 ), and between 972 and 998 °C, bixbyite decomposed to hausmannite (Mn 2+ Mn 3+ 2 O 4 ). The rutile experienced no phase decompositions or transitions over this same temperature interval. The refined bond lengths and angles for pyrolusite showed that the Mn 4+ coordination octahedron became more distorted near the phase decomposition to bixbyite due to the approach of two coordinating oxygen atoms. The 3d electrons of each Mn 4+ ion in pyrolusite make π bonds with the 2p electrons of the surrounding O 2- ions, whereas the Ti 4+ of rutile has no 3d electrons. Thus, π bonding between octahedral Mn 4+ and the surrounding O 2- anions increases the strength of the Mn–O bonds of pyrolusite relative to the Ti–O bonds of rutile. However, we observed a small decrease in O–O distances in pyrolusite before the decomposition, suggesting that at high temperature, increased π bonding between adjacent O anions precedes the release of O 2 and the reduction of Mn 4+ to Mn 3+ . Analogous behavior was observed for bixbyite before its reductive phase decomposition to hausmannite. In contrast, no anomalous changes in the O–O distances occurred for rutile. Our X-ray diffraction analyses provided accurate thermal expansion coefficients for these materials over a broader temperature range than reported in previous studies. The Mn–O bond lengths in pyrolusite were shorter and stronger than the Ti–O bonds of rutile; for example, at 60 °C, bond distances in pyrolusite refined to 1.8854(6) Å whereas the bond lengths of rutile were 1.9575(4) Å. Accordingly, the axial and volumetric thermal expansion coefficients refined for rutile were greater than those for pyrolusite. Further, the axial coefficients of thermal expansion (CTE) for pyrolusite were 96% (a-axis) and 69% (c-axis) of the corresponding values for rutile, and the volumetric CTE for pyrolusite was 85% of that for rutile.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An airfoil-based synthetic actuator disk model for wind turbine aerodynamic and structural analysis

Here, this study introduces an airfoil-based refinement technique to enhance the Actuator Disk Model (ADM) for improved wind turbine aerodynamic load prediction and structural simulation in conjunction with Large Eddy Simulations of the wind flow. While ADM offers higher computational efficiency than the more detailed but resource-intensive Actuator Line Model (ALM), it traditionally lacks the resolution needed to capture the localized blade forces accurately. To address this limitation, we introduce a refinement technique that uses airfoil-specific data and employs interpolation-based grid point refinement, achieving ALM-comparable accuracy while preserving ADM's efficiency. Unlike conventional ADM that provides only rotor-disk averaged forces, our synthetic method tracks transient aerodynamic load variations over multiple blade revolutions, allowing us to calculate the distributions of maximum and minimum loads during typical cycles. Applied to the NREL 5 MW reference turbine, our enhanced ADM accurately predicts key aerodynamic parameters (angle of attack, axial velocity, lift, drag, axial and tangential forces along the blades) as well as structural responses (blade tip deflection, maximum stress, and stress concentration). Our results show that the tip deflection ranges from 2.33m (3.69 % of blade length) to 4.28m (6.79 %), with maximum stress concentration occurring near the blade root. This research demonstrates that a refined synthetic ADM approach can serve as a computationally efficient alternative for both aerodynamic analysis and structural simulation of wind turbine blades subjected to realistic wind fields.

17 WIND ENERGY↗