Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “embedded boundary”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

Second Annual Progress Report on Correlation Between Microstructure and Mechanical Properties of Neutron-Irradiated Ferritic-Martensitic and Austenitic Steels

Ferritic-martensitic steels G92-2b (an optimized Grade 92 heat), NF616 and T91, and austenitic stainless steel 800H and its Grain Boundary Engineering (GBE)-treated version 800H-TMP (ThermoMechanical Processing) were irradiated in the High Flux Isotope Reactor (HFIR) of Oak Ridge National Laboratory (ORNL) and the Advanced Test Reactor (ATR) of Idaho National Laboratory (INL). Selected G92-2b samples were irradiated up to 14.66 dpa in the HFIR at two temperature ranges: 400–496.7°C and 683.3– 720°C. NF616 and T91 were irradiated in the ATR up to 8.16 dpa with the irradiation temperatures ranged from 241°C to 447.5°C. Alloy 800H and 800H-TMP samples were irradiated in both the HFIR and the ATR. Selected 800H and 800H-TMP samples had HFIR irradiation to 1.28 dpa at 580°C and ATR irradiation up to 9.12 dpa at 359°C to 431°C. Vickers hardness measurements, fractography, and microstructural characterization were performed on the selected samples in the Low Activation Materials Design and Analysis (LAMDA) laboratory. Radiation-hardening of G92-2b was observed at the lower doses and lower irradiation temperatures (400- 496.7°C), with GB03 (0.52 dpa at 400°C) and GB04 (7.44 dpa at ~490°C) showing ~12% and ~8% hardening, respectively. Softening by ~14% was observed for GB05 (14.66 dpa at 496.7°C). Radiationsoftening of G92-2b was more prevalent at the higher irradiation temperatures (683.3-~720°C), with GB10 (0.46 dpa at 683.3°C), GB11 (7.44 dpa at ~720°C), and GB12 (14.63 dpa at ~720°C) showing ~8%, ~8%, and ~40% softening, respectively. Radiation-hardening of NF616 and T91 was observed with the hardness increased by ~37% to ~65% depending on the irradiation doses and irradiation temperatures. Within the studied irradiation conditions of NF616 and T91, samples with a higher dose had a larger hardness after irradiation. All the tested alloy 800H and 800H-TMP samples in this work showed radiation-hardening by ~96±7% to ~152±10%. Alloy 800H-TMP tended to have slightly smaller radiation-hardening than alloy 800H. The fractography results of G92-2b sample GB03, GB10, and GB11, together with the previously characterized fractography of GB04, GB05, and GB12, indicated that the ductility of G92-2b was maintained up to 14.66 dpa at the lower irradiation temperatures of 400-496.7°C, while some loss of ductility (less necking) was observed for higher doses at the higher irradiation temperatures of 683.3-720°C. This agrees with the previously reported tensile test results of G92-2b, where the elongation of G92-2b was reduced at higher doses at the higher irradiation temperatures. Dimple sizes increased at higher doses, which are more evident at the higher irradiation temperatures of 683.3-~720°C. The fractography of NF616 sample D2 (2.96 dpa at 291.5°C), D4 (5.91 dpa at 359°C), and D6 (8.16 dpa at 431°C) indicated loss of ductility with negligible necking for sample D2, while ductile failure for samples D4 and D6. Fractography of alloy 800H and 800H-TMP samples in various irradiation conditions showed ductile failure with obvious necking. Dimples were observed, with some of them containing large Ti-rich particles, in all the 800H/800H-TMP samples. Electron backscatter diffraction characterization of GB12 indicated the recovery of the lath structure, which was generally replaced by an equiaxed grain structure. Transmission electron microscopy (TEM) characterization showed the presence of frequent M 23 C 6 (M = primarily Cr), MX (M = primarily V), spherical Nb(C,N) precipitates, and occasional Laves phase precipitates in the G92-2b samples. MX precipitates with sizes of 20-30 nm were observed at boundaries of smaller grains, indicating the pinning effect of the V-rich precipitates. The lath structure recovery was more evident at the higher irradiation temperatures (683.3-~720°C), with decreased densities of line dislocations and M23C6 precipitates. The irradiated T91 (TA04) showed the growth of M23C6 precipitates to 101 ± 40 nm from the initial 68 ± 22 nm in the unirradiated condition. Dislocation loops of both {100} and {111} types were present in TA04. TEM characterization was also performed on the irradiated 800H (N4, N5, N6, and AR2) and 800H-TMP (P4, P5, P6, and HG1). Accumulation of M 23 C 6 precipitates at grain boundaries was observed in all the xii 800H/800H-TMP samples, and the presence of Ti(C,N) precipitates at grain boundaries and in the matrix was observed in the irradiated 800H-TMP. Some Ti(C,N) precipitates are embedded in the M 23 C 6 precipitates, maintaining specific orientation relationships between the precipitates and between the precipitate and matrix. In addition, nanoscale Si-rich clusters were observed in the matrix of all the 800H/800H-TMP samples, with EDS Si maps tending to have a lower contrast in 800H-TMP samples. Atom probe tomography was conducted on the same samples, supported by a Rapid Turnaround Examination project under Nuclear Science User Facilities. The results are being analyzed to be integrated with the TEM results for a confident description of the γ’ precipitates. Dislocation loops also formed in all the 800H/800H-TMP samples. The density and the average size of dislocation loops were quantified to be in the order of 10 22 – 10 23 m -3 and 11.7 – 15.9 nm, respectively, in the ATR-irradiated 800H/800H-TMP samples. The loop densities in irradiated 800H were higher than that in irradiated 800H-TMP under the same irradiation conditions. Further systematic data analyses, together with some complementary experiments, will be pursued for these samples to foster peer-reviewed journal article publications.

36 MATERIALS SCIENCE↗

Analysis of supersonic combustion flow fields with embedded subsonic regions

The viscous characteristic analysis for supersonic chemically reacting flows was extended to include provisions for analyzing embedded subsonic regions. The numerical method developed to analyze this mixed subsonic-supersonic flow fields is described. The boundary conditions are discussed related to the supersonic-subsonic and subsonic-supersonic transition, as well as a heuristic description of several other numerical schemes for analyzing this problem. An analysis of shock waves generated either by pressure mismatch between the injected fluid and surrounding flow or by chemical heat release is also described.

Dash, S.↗

Methods for analysis of cracks in three-dimensional solids

Analytical and numerical methods evaluating the stress-intensity factors for three-dimensional cracks in solids are presented, with reference to fatigue failure in aerospace structures. The exact solutions for embedded elliptical and circular cracks in infinite solids, and the approximate methods, including the finite-element, the boundary-integral equation, the line-spring models, and the mixed methods are discussed. Among the mixed methods, the superposition of analytical and finite element methods, the stress-difference, the discretization-error, the alternating, and the finite element-alternating methods are reviewed. Comparison of the stress-intensity factor solutions for some three-dimensional crack configurations showed good agreement. Thus, the choice of a particular method in evaluating the stress-intensity factor is limited only to the availability of resources and computer programs.

Raju, I. S.↗

Numerical investigation of multi-element airfoils

The flow over multi-element airfoils with flat-plate lift-enhancing tabs was numerically investigated. Tabs ranging in height from 0.25 percent to 1.25 percent of the reference airfoil chord were studied near the trailing edge of the main-element. This two-dimensional numerical simulation employed an incompressible Navier-Stokes solver on a structured, embedded grid topology. New grid refinements were used to improve the accuracy of the solution near the overlapping grid boundaries. The effects of various tabs were studied at a constant Reynolds number on a two-element airfoil with a slotted flap. Both computed and measured results indicated that a tab in the main-element cove improved the maximum lift and lift-to-drag ratio relative to the baseline airfoil without a tab. Computed streamlines revealed that the additional turning caused by the tab may reduce the amount of separated flow on the flap. A three-element airfoil was also studied over a range of Reynolds numbers. For the optimized flap rigging, the computed and measured Reynolds number effects were similar. When the flap was moved from the optimum position, numerical results indicated that a tab may help to reoptimize the airfoil to within 1 percent of the optimum flap case.

Cummings, Russell M.↗

Thermal Conductivity Estimation from Transient Test Data with Embedded Thermocouples using Genetic Algorithm Optimization

Inverse heat transfer methodology previously developed to estimate thermal properties of high temperature fibrous insulation from embedded thermocouples was applied to the thermoplastic polymer, polyether-ether-ketone (PEEK). A small experimental setup was utilized to test the PEEK sample between temperatures of 300 K and 525 K at atmospheric pressure. Cylindrical plugs of PEEK with three thermocouples embedded at various depths were incorporated in the test sample. The experimental PEEK thermocouple data were used as the boundary and initial conditions of a one-dimensional numerical thermal model to predict the internal temperatures of the material. The thermal conductivity of PEEK was estimated with the Continuous Genetic Algorithm optimization technique by searching for the coefficients of a functional form of thermal conductivity that minimized the difference between the experimentally measured and model predicted internal temperature values. The thermal testing, one-dimensional numerical thermal model, optimization algorithm, analysis, and results are presented.

Thermal Properties↗

Conceptual Design of a Counter-Rotating Fan System for Distributed Boundary Layer Ingesting Propulsion

The present paper details the design of the counter rotating fans for a Turboelectric Distributed Propulsion (TeDP) system. Sixteen propulsors installed in mail-slot-shape nacelles are embedded on an aerodynamically optimized hybrid wing-body configuration. The hybrid-wing/body (HWB) configuration which was previously designed to satisfy the conditions of trim, longitudinally static stability and specific cargo space is employed as the baseline configuration in pursuing an optimal distributed propulsion system. A set of distributed propulsors is conceptually designed and the collective performance is evaluated against the target thrust mandated by the mission requirements. The concept of the distributed propulsion allows the fan pressure ratio to be around 1.27~1.32 for the target thrust. In addition, further splitting of the fan pressure ratio by using the counter-rotating fans for each slot realizes the target pressure ratio with low tip speed. In the distributed propulsion system, the nature of the flow conditions and/or the thickness of the ingested boundary layer may differ and result in different propulsive reaction of each individual propulsor. The optimization is, thus, approached from both the propulsion system and individual propulsor perspectives. An optimal distribution of the thrust and power output is determined by how the system utilizes each passage's propulsive characteristics and its interaction with the airframe. These system level analysis and optimization are conducted using an actuator disk model to account for the propulsion-airframe integration numerically. With respect to the propulsor level, aerodynamic shape optimizations of the fan blades are performed in a sequential multi-objective optimization process for various design objectives, such as mass flow rate condition, fan pressure ratio, efficiency and the exit flow angle of the fan stage by using a genetic algorithm, NSGA-II. The radial chord distribution, and meanline distribution of the rotors are designed on the circumferentially averaged axi-symmetric inlet profiles and tested on the six inlet profiles from six divided sectors to reckon flow distortion. The performances of the counter rotating fans are, thus, evaluated accordingly for obtaining distortion tolerant fan. The performance of the distributed propulsion system is evaluated by two CFD tools, i.e., a multi-stage turbo-machinery CFD code and one propulsion-airframe integration flow solver coupled with a body-force model. The optimized boundary layer ingestion propulsion system of 16 distributed slots not only reaches the system target thrust, but also delivers a close to 20% fuel saving benefit against its counterpart 12 distributed clean inlet propulsion system.

Boundary-Layer-Ingestion Propulsion↗

Graph-learning approach to combine multiresolution seismic velocity models

SUMMARY The resolution of velocity models obtained by tomography varies due to multiple factors and variables, such as the inversion approach, ray coverage, data quality, etc. Combining velocity models with different resolutions can enable more accurate ground motion simulations. Toward this goal, we present a novel methodology to fuse multiresolution seismic velocity maps with probabilistic graphical models (PGMs). The PGMs provide segmentation results, corresponding to various velocity intervals, in seismic velocity models with different resolutions. Further, by considering physical information (such as ray path density), we introduce physics-informed probabilistic graphical models (PIPGMs). These models provide data-driven relations between subdomains with low (LR) and high (HR) resolutions. Transferring (segmented) distribution information from the HR regions enhances the details in the LR regions by solving a maximum likelihood problem with prior knowledge from HR models. When updating areas bordering HR and LR regions, a patch-scanning policy is adopted to consider local patterns and avoid sharp boundaries. To evaluate the efficacy of the proposed PGM fusion method, we tested the fusion approach on both a synthetic checkerboard model and a fault zone structure imaged from the 2019 Ridgecrest, CA, earthquake sequence. The Ridgecrest fault zone image consists of a shallow (top 1 km) high-resolution shear-wave velocity model obtained from ambient noise tomography, which is embedded into the coarser Statewide California Earthquake Center Community Velocity Model version S4.26-M01. The model efficacy is underscored by the deviation between observed and calculated traveltimes along the boundaries between HR and LR regions, 38 per cent less than obtained by conventional Gaussian interpolation. The proposed PGM fusion method can merge any gridded multiresolution velocity model, a valuable tool for computational seismology and ground motion estimation.

Geochemistry & Geophysics↗

Modeling instrumented test articles for Plasmatron X aerothermal environment

This work discusses the ablative response of thermal protection system materials when exposed to aerothermal environments. It focuses on modeling instrumented test article assemblies tested in the Plasmatron X inductively coupled plasma wind tunnel. For this purpose, as the samples are mounted to a water-cooled arm using numerous assembly components and are embedded with thermocouples for in-depth temperature measurements, the variation in ablative response of the samples as a result of these components is studied. The aerothermal boundary conditions for Plasmatron X are obtained for hypersonic flight-relevant conditions using the multi-physics aerothermal framework built at the University of Illinois at Urbana-Champaign. For simulating the material response, a material response code, Porous Material Analysis Toolbox based on OpenFOAM (PATO) is used. In PATO, numerous configurations and materials of thermocouples were studied. The optimal configuration for obtaining temperature measurements with minimal error using embedded thermocouples was found to be the U-shaped Type S, 30 American wire gauge thermocouple configuration.

Thermal Protection Systems↗

Modeling instrumented test articles for Plasmatron X aerothermal environment

This work discusses the ablative response of thermal protection system materials when exposed to aerothermal environments. It focuses on modeling instrumented test article assemblies tested in the Plasmatron X inductively coupled plasma wind tunnel. For this purpose, as the samples are mounted to a water-cooled arm using numerous assembly components and are embedded with thermocouples for in-depth temperature measurements, the variation in ablative response of the samples as a result of these components is studied. The aerothermal boundary conditions for Plasmatron X are obtained for hypersonic flight-relevant conditions using the multi-physics aerothermal framework built at the University of Illinois at Urbana-Champaign. For simulating the material response, a material response code, Porous Material Analysis Toolbox based on OpenFOAM (PATO) is used. In PATO, numerous configurations and materials of thermocouples were studied. The optimal configuration for obtaining temperature measurements with minimal error using embedded thermocouples was found to be the U-shaped Type S, 30 American wire gauge thermocouple configuration.

Thermal Protection Systems↗

A full potential flow analysis with realistic wake influence for helicopter rotor airload prediction

A 3-D, quasi-steady, full potential flow solver was adapted to include realistic wake influence for the aerodynamic analysis of helicopter rotors. The method is based on a finite difference solution of the full potential equation, using an inner and outer domain procedure for the blade flowfield to accommodate wake effects. The nonlinear flow is computed in the inner domain region using a finite difference solution method. The wake is modeled by a vortex lattice using prescribed geometry techniques to allow for the inclusion of realistic rotor wakes. The key feature of the analysis is that vortices contained within the finite difference mesh (inner domain) were treated with a vortex embedding technique while the influence of the remaining portion of the wake (in the outer domain) is impressed as a boundary condition on the outer surface of the finite difference mesh. The solution procedure couples the wake influence with the inner domain solution in a consistent and efficient solution process. The method has been applied to both hover and forward flight conditions. Correlation with subsonic and transonic hover airload data is shown which demonstrates the merits of the approach.

Egolf, T. Alan↗

Simulating Sources of Superstorm Plasmas

We evaluated the contributions to magnetospheric pressure (ring current) of the solar wind, polar wind, auroral wind, and plasmaspheric wind, with the surprising result that the main phase pressure is dominated by plasmaspheric protons. We used global simulation fields from the LFM single fluid ideal MHD model. We embedded the Comprehensive Ring Current Model within it, driven by the LFM transpolar potential, and supplied with plasmas at its boundary including solar wind protons, polar wind protons, auroral wind O+, and plasmaspheric protons. We included auroral outflows and acceleration driven by the LFM ionospheric boundary condition, including parallel ion acceleration driven by upward currents. Our plasmasphere model runs within the CRCM and is driven by it. Ionospheric sources were treated using our Global Ion Kinetics code based on full equations of motion. This treatment neglects inertial loading and pressure exerted by the ionospheric plasmas, and will be superceded by multifluid simulations that include those effects. However, these simulations provide new insights into the respective role of ionospheric sources in storm-time magnetospheric dynamics.

Fok, Mei-Ching↗

Transient, Small-Scale Field-Aligned Currents in the Plasma Sheet Boundary Layer During Storm Time Substorms

We report on field-aligned current observations by the four Magnetospheric Multiscale (MMS) spacecraft near the plasma sheet boundary layer (PSBL) during two major substorms on 23 June 2015. Small-scale field-aligned currents were found embedded in fluctuating PSBL flux tubes near the Separatrix region. We resolve, for the first time, short-lived earthward (downward) intense field-aligned current sheets with thicknesses of a few tens of kilometers, which are well below the ion scale, on flux tubes moving equatorward earth ward during outward plasma sheet expansion. They coincide with upward field-aligned electron beams with energies of a few hundred eV. These electrons are most likely due to acceleration associated with a reconnection jet or high-energy ion beam-produced disturbances. The observations highlight coupling of multiscale processes in PSBL as a consequence of magnetotail reconnection.

Nakamura, R.↗

Integrated Data Modeling and Simulation on the Joint Polar Satellite System Program

The Joint Polar Satellite System is a modern, large-scale, complex, multi-mission aerospace program, and presents a variety of design, testing and operational challenges due to: (1) System Scope: multi-mission coordination, role, responsibility and accountability challenges stemming from porous/ill-defined system and organizational boundaries (including foreign policy interactions) (2) Degree of Concurrency: design, implementation, integration, verification and operation occurring simultaneously, at multiple scales in the system hierarchy (3) Multi-Decadal Lifecycle: technical obsolesce, reliability and sustainment concerns, including those related to organizational and industrial base. Additionally, these systems tend to become embedded in the broader societal infrastructure, resulting in new system stakeholders with perhaps different preferences (4) Barriers to Effective Communications: process and cultural issues that emerge due to geographic dispersion and as one spans boundaries including gov./contractor, NASA/Other USG, and international relationships.

Roberts, Christopher J.↗

Global Sensitivity Analysis of Simulated Remote Sensing Polarimetric Observations Over Snow

This study presents a detailed theoretical assessment of the information content of passive polarimetric observations over snow scenes, using a global sensitivity analysis (GSA) method. Conventional sensitivity studies focus on varying a single parameter while keeping all other parameters fixed. In contrast, the GSA correctly addresses the covariance of state parameters across their entire parameter space, hence favoring a more correct interpretation of inversion algorithms and the optimal design of their state vectors. The forward simulations exploit a vector radiative transfer model to obtain the Stokes vector emerging at the top of the atmosphere for different solar zenith angles, when the bottom boundary consists of a vertically resolved snowpack of non-spherical grains. The presence of light-absorbing particulates (LAPs), either embedded in the snow or aloft in the atmosphere above in the form of aerosols, is also considered. The results are presented for a set of wavelengths spanning the visible (VIS), near-infrared (NIR), and shortwave infrared (SWIR) region of the spectrum. The GSA correctly captures the expected, high sensitivity of the reflectance to LAPs in the VIS–NIR and to grain size at different depths in the snowpack in the NIR–SWIR. With adequate viewing geometries, mono-angle measurements of total reflectance in the VIS–SWIR (akin to those of the Moderate Resolution Imaging Spectroradiometer, MODIS) resolve grain size in the top layer of the snowpack sufficiently well. The addition of multi-angle polarimetric observations in the VIS–NIR provides information on grain shape and microscale roughness. The simultaneous sensitivity in the VIS–NIR to both aerosols and snow-embedded impurities can be disentangled by extending the spectral range to the SWIR, which contains information on aerosol optical depth while remaining essentially unaffected when the same particulates are mixed with the snow. Multi-angle polarimetric observations can therefore (i) effectively partition LAPs between the atmosphere and the surface, which represents a notorious challenge for snow remote sensing based on measurements of total reflectance only and (ii) lead to better estimates of grain shape and roughness and, in turn, the asymmetry parameter, which is critical for the determination of albedo. The retrieval uncertainties are minimized when the degree of linear polarization is used in place of the polarized reflectance. The Sobol indices, which are the main metric for the GSA, were used to select the state parameters in retrievals performed on data simulated for multiple instrument configurations. Improvements in retrieval quality with the addition of measurements of polarization, multi-angle views, and different spectral channels reflect the information content, identified by the Sobol indices, relative to each configuration. The results encourage the development of new remote sensing algorithms that fully leverage multi-angle and polarimetric capabilities of modern remote sensors. They can also aid flight planning activities, since the optimal exploitation of the information content of multi-angle measurements depends on the viewing geometry. The better characterization of surface and atmospheric parameters in snow-covered regions advances research opportunities for scientists of the cryosphere and ultimately benefits albedo estimates in climate models.

remote sensing↗

An Embedded Statistical Method for Coupling Molecular Dynamics and Finite Element Analyses

The coupling of molecular dynamics (MD) simulations with finite element methods (FEM) yields computationally efficient models that link fundamental material processes at the atomistic level with continuum field responses at higher length scales. The theoretical challenge involves developing a seamless connection along an interface between two inherently different simulation frameworks. Various specialized methods have been developed to solve particular classes of problems. Many of these methods link the kinematics of individual MD atoms with FEM nodes at their common interface, necessarily requiring that the finite element mesh be refined to atomic resolution. Some of these coupling approaches also require simulations to be carried out at 0 K and restrict modeling to two-dimensional material domains due to difficulties in simulating full three-dimensional material processes. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the standard boundary value problem used to define a coupled domain. The method replaces a direct linkage of individual MD atoms and finite element (FE) nodes with a statistical averaging of atomistic displacements in local atomic volumes associated with each FE node in an interface region. The FEM and MD computational systems are effectively independent and communicate only through an iterative update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM). ESCM provides an enhanced coupling methodology that is inherently applicable to three-dimensional domains, avoids discretization of the continuum model to atomic scale resolution, and permits finite temperature states to be applied.

Saether, E.↗

Singular Perturbation-Based Large-Signal Order Reduction of Microgrids for Stability and Accuracy Synthesis With Control

The increasing penetration of distributed energy resources (DERs) highlights the growing importance of microgrids (MGs) in enhancing power system reliability. Employing electromagnetic transient (EMT) analysis in MGs becomes crucial for controlling the rapid transients. However, this requires an accurate but high-order model of power electronics and their underlying control loops, complexifying the stability analysis from the viewpoint of a higher control level. To overcome these challenges, this paper proposes a large-signal order reduction (LSOR) method for MGs with considerations of external control inputs and the detailed dynamics of underlying control levels based on singular perturbation theory (SPT). Specially, we innovatively proposed and strictly proved a general stability and accuracy assessment theorem that allows us to analyze the dynamic stability of a full-order nonlinear system by only leveraging our derived reduced-order model (ROM) and boundary layer model (BLM). Furthermore, this theorem furnishes a set of conditions that determine the accuracy of the developed ROM. Lastly, by embedding such a theorem into the SPT, we propose a novel LSOR approach with guaranteed accuracy and stability analysis equivalence. Case studies are conducted on MG systems to show the effectiveness of the proposed approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗

NBODY - a multipurpose trajectory optimization computer program

Documentation of the NBODY trajectory optimization program is presented in the form of a mathematical development plus a user's manual. Optimal multistage-launch ascent trajectories may be determined by variational thrust steering during the upper phase. Optimal low-thrust interplanetary spacecraft trajectories may also be calculated with solar power or constant power, all-propulsion or embedded coast arcs, fixed or optimal thrust angles, and a variety of terminal end conditions. A hybrid iteration scheme solves the boundary-value problem, while either transversality conditions or a univariate search scheme optimize vehicle or trajectory parameters.

Strack, W. C.↗

The free-wake prediction of rotor hover performance using a vortex embedding method

A method is developed to predict the rotor hover performance. This method solves the compressible mass conservation equation much like current full potential codes and can therefore predict the transonic flows on a rotor. However, the newly developed approach also allows for the free convection of shed vorticity and permits the computation of the entire hover wake system. The method uses a vortex embedding scheme in potential flow and has been implemented in a computer code, HELIX -I. To predict power we implement a simple boundary layer and two different induced-drag integration schemes. The induced-drag is obtained from surface pressure integration and an energy flux integral. Comparisons between computations and experiment show good agreement for the prediction of power polars, surface pressure distribution, and tip vortex geometry.

Ramachandran, K.↗