Search NASA⌕ Search

SEARCH · Search NASA

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 181 records · Page 10

A Model for Energy Buildup and Eruption Onset in Coronal Mass Ejections

Coronal mass ejections (CMEs) and eruptive flares (EFs) are the most energetic explosions in the solar system. Their underlying origin is the free energy that builds up slowly in the sheared magnetic field of a filament channel. We report the first end-to-end numerical simulation of a CME/EF, from zero-free-energy initial state through filament channel formation to violent eruption, driven solely by the magnetic-helicity condensation process. Helicity is the topological measure of linkages between magnetic flux systems, and is conserved in the corona, building up inexorably until it is ejected into interplanetary space. Numerous investigations have demonstrated that helicity injected by small-scale vortical motions, such as those observed in the photosphere, undergoes an inverse cascade from small scales to large, “condensing” at magnetic-polarity boundaries. Our new results verify that this process forms a filament channel within a compact bipolar region embedded in a background dipole field, and show for the first time that a fast CME eventually occurs via the magnetic-breakout mechanism. We further show that the trigger for explosive eruption is reconnection onset in the flare current sheet that develops above the polarity inversion line: this reconnection forms flare loops below the sheet and a CME flux rope above, and initiates high-speed outward flow of the CME. Our findings have important implications for magnetic self organization and explosive behavior in solar and other astrophysical plasmas, as well as for understanding and predicting explosive solar activity.

filaments↗

Evaluating Precipitation Features and Rainfall Characteristics in a Multi-scale Modeling Framework

Cloud and precipitation systems over the tropics and subtropics are simulated with a multi-scale modeling framework (MMF) and compared against the TRMM radar precipitation features (RPFs) product. A methodology, in close analogy to the TRMM RPFs, is developed to analyze simulated cloud precipitating structures from the embedded two-dimensional cloud-resolving models (CRMs) within an MMF. Despite the two-dimensionality of the CRMs, the simulated RPFs population distribution, and horizontal and vertical structure are in good agreement with TRMM observations. However, some deficits are also found in the model simulations. The model tends to overestimate mean convective precipitation rates for RPFs with a size less than 100 km, contributing to the excessive precipitation biases in the warm pool and western Pacific, western and northern India Ocean, and eastern Pacific commonly found in most MMFs. For large features with a size greater than 150 km, both convective and stratiform rain rates are underestimated. The distribution of maximum radar echo top heights as a function of RPF size is well simulated except the model tends to underestimate the occurrence frequency of maximum heights greater than 15 km. The maximum echo top heights for convective cells embedded within large RPFs with a size greater than 150 km are also underestimated. The cyclic lateral boundary with a limited model domain generates artificial occurrences for RPFs with a size close to the model domain size, producing a significant contribution to the total rainfall due to their sizes. This cyclic lateral boundary effect can be easily identified and quantified in both probability and cumulative distribution functions of RPFs. The geophysical distribution of the population of the largest RPFs in the control experiment shows they are mainly located in the Subtropics but also partially contribute to the common MMF biases of excessive precipitation in the Tropics. Sensitivity experiments using CRMs with different domain sizes and different grid spacings show larger domains (higher resolution) tend to shift the RPFs distribution to large (small) sizes. The cyclic lateral boundary biases increase as CRM domain size decreases. The impacts of model horizontal and vertical resolution on simulated convective systems are also investigated.

Chern, Jiun-Dar↗

Computation of unsteady flows over airfoils

Two methods are described for calculating unsteady flows over rapidly pitching airfoils. The first method is based on an interactive scheme in which the inviscid flow is obtained by a panel method. The boundary layer flow is computed by an interactive method that makes use of the Hilbert integral to couple the solutions of the inviscid and viscous flow equations. The second method is based on the solution of the compressible Navier-Stokes equations. The solution of these equations is obtained with an approximately factorized numerical algorithm, and with single block or multiple grids which enable grid embedding to enhance the resolution at isolated flow regions. In addition, the attached flow region can be computed by the numerical solution of compressible boundary layer equations. Unsteady pressure distributions obtained with both methods are compared with available experimental data.

Ekaterinaris, J. A.↗

Methods for analysis of cracks in three-dimensional solids

Various analytical and numerical methods used to evaluate the stress intensity factors for cracks in three-dimensional (3-D) solids are reviewed. Classical exact solutions and many of the approximate methods used in 3-D analyses of cracks are reviewed. The exact solutions for embedded elliptic cracks in infinite solids are discussed. The approximate methods reviewed are the finite element methods, the boundary integral equation (BIE) method, the mixed methods (superposition of analytical and finite element method, stress difference method, discretization-error method, alternating method, finite element-alternating method), and the line-spring model. The finite element method with singularity elements is the most widely used method. The BIE method only needs modeling of the surfaces of the solid and so is gaining popularity. The line-spring model appears to be the quickest way to obtain good estimates of the stress intensity factors. The finite element-alternating method appears to yield the most accurate solution at the minimum cost.

Raju, I. S.↗

Thermospheric dynamics during November 21-22, 1981 - Dynamics Explorer measurements and thermospheric general circulation model predictions

Time-dependent aurora and magnetospheric convection parameterizations have been derived from solar wind and aurora particle data for November 21-22, 1981, and are used to drive the auroral and magnetospheric convection models that are embedded in the National Center for Atmospheric Research thermospheric general circulation model (TGCM). Neutral wind speeds and transition boundaries between the midlatitude solar-driven circulation and the high-latitude magnetospheric convection-driven circulation are examined on an orbit-by-orbit basis. The results show that TGCM-calculated winds and reversal boundary locations are in generally good agreement with Dynamics Explorer 2 measurements for the orbits studied. This suggests that, at least for this particular period of relatively moderate geomagnetic activity, the TGCM parameterizations on the eveningside of the auroral oval and polar cap are adequate.

Roble, R. G.↗

The 5-6 December 1991 FIRE IFO II jet stream cirrus case study: Possible influences of Volcanic Aerosols

In presenting an overview of the cirus clouds comprehensively studied by ground-based and airborne sensors from Coffeyville, Kansas, during the 5-6 December 1992 Project First ISCCP Region Experiment (FIRE) Intensive Fields Observation (IFO) II case study period, evidence is provided that volcanic aerosols friom the June 1991 Pinatubo eruptions may have significantly influenced the formation and maintenance of the cirrus. Following the local appearance of a spur of stratospheric volcanic debris from the subtropics, a series of jet streaks subsequently conditioned the troposphere through tropopause foldings with sulfur-based particles that became effective cloud-forming nuclei in cirrus clouds. Aerosol and ozone measurements suggest a complicated history of stratospheric-tropospheric exchanges embedded within the upper-level flow, and cirrus cloud formation was noted to occur locally at the boundaries of stratospheric aerosol-enriched layers that became humidified through diffusion, precipitation, or advective processes. Apparent cirrus cloud alterations include abnormally high ice crystal concentrations (up to approximately 600/L), complex radial ice crystal types, and relatively large haze particles in cirrus uncinus cell heads at temperatures between -40 and -50 C. Implications for volcanic-cirrus cloud climate effects and usual (nonvolcanic aerosol) jet stream cirrus cloud formation are discussed.

Sassen, Kenneth↗

A Hybrid Boundary Element-Finite Volume Method for Unsteady Transonic Airfoil Flows

A hybrid boundary element finite volume method for unsteady transonic flow computation has been developed. In this method, the unsteady Euler equations in a moving frame of reference are solved in a small embedded domain (inner domain) around the airfoil using an implicit finite volume scheme. The unsteady full-potential equation, written in the same frame of reference and in the form of the Poisson equation. is solved in the outer domain using the integral equation boundary element method to provide the boundary conditions for the inner Euler domain. The solution procedure is a time-accurate stepping procedure, where the outer boundary conditions for the inner domain are updated using the integral equation -- boundary element solution over the outer domain. The method is applied to unsteady transonic flows around the NACA0012 airfoil undergoing pitching oscillation and ramp motion. The results are compared with those of an implicit Euler equation solver, which is used throughout a large computational domain, and experimental data.

Hu, Hong↗

Flow-Field Measurement of Device-Induced Embedded Streamwise Vortex on a Flat Plate

Detailed flow-field measurements were performed downstream of a single vortex generator (VG) using an advanced Stereo Digital Particle Image Velocimetry system. Thc passive flow-control devices examined consisted of a low-profile VG with a device height, h, approximately equal to 20 percent of the boundary-layer thickness, sigma, and a conventional VG with h is approximately sigma. Flow-field data were taken at twelve cross-flow planes downstream of the VG to document and quantify the evolution of embedded streamwise vortex. The effects of device angle of attack on vortex development downstream were compared between the low-profile VG and the conventional VG. Key parameters including vorticity, circulation, trajectory, and half-life radius - describing concentration, strength, path, and size, respectively--of the device-induced streamwise vortex were extracted from the flow-field data. The magnitude of maximum vorticity increases as angle of attack increases for the low-profile VG, but the trend is reversed for the conventional VG, probably due to flow stalling around the larger device at higher angles of attack. Peak vorticity and circulation for the low-profile VG decays exponentially and inversely proportional to the distance downstream from the device. The device-height normalized vortex trajectories for the low-profile VG, especially in the lateral direction, follow the general trends of the conventional VG. The experimental database was used to validate the predictive capability of computational fluid dynamics (CFD). CFD accurately predicts the vortex circulation and path; however, improvements are needed for predicting the vorticity strength and vortex size.

Yao, Chung-Sheng↗

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↗

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↗