Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic 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

Control of unsteady separated flow associated with the dynamic stall of airfoils

A unique active flow-control device is proposed for the control of unsteady separated flow associated with the dynamic stall of airfoils. The device is an adaptive-geometry leading-edge which will allow controlled, dynamic modification of the leading-edge profile of an airfoil while the airfoil is executing an angle-of-attack pitch-up maneuver. A carbon-fiber composite skin has been bench tested, and a wind tunnel model is under construction. A baseline parameter study of compressible dynamic stall was performed for flow over an NACA 0012 airfoil. Parameters included Mach number, pitch rate, pitch history, and boundary layer tripping. Dynamic stall data were recorded via point-diffraction interferometry and the interferograms were analyzed with in-house developed image processing software. A new high-speed phase-locked photographic image recording system was developed for real-time documentation of dynamic stall.

Wilder, M. C.↗

From boundary data to bound states. Part II. Scattering angle to dynamical invariants (with twist)

We recently introduced in [9] a boundary-to-bound dictionary between gravitational scattering data and observables for bound states of non-spinning bodies. In this paper, we elaborate further on this holographic map. We start by deriving the following — remarkably simple — formula relating the periastron advance to the scattering angle: $ΔΦ(J,ε)=χ(J,ε)+χ(-J,ε)$, via analytic continuation in angular momentum and binding energy. Using explicit expressions from [9], we confirm its validity to all orders in the Post-Minkowskian (PM) expansion. Furthermore, we reconstruct the radial action for the bound state directly from the knowledge of the scattering angle. The radial action enables us to write compact expressions for dynamical invariants in terms of the deflection angle to all PM orders, which can also be written as a function of the PM-expanded amplitude. As an example, we reproduce our result in [9] for the periastron advance, and compute the radial and azimuthal frequencies and redshift variable to two-loops. Agreement is found in the overlap between PM and Post-Newtonian (PN) schemes. Last but not least, we initiate the study of our dictionary including spin. We demonstrate that the same relation between deflection angle and periastron advance applies for aligned-spin contributions, with J the (canonical) total angular momentum. Explicit checks are performed to display perfect agreement using state-of-the-art PN results in the literature. Using the map between test- and two-body dynamics, we also compute the periastron advance up to quadratic order in spin, to one-loop and to all orders in velocity. We conclude with a discussion on the generalized ‘impetus formula’ for spinning bodies and black holes as ‘elementary particles’. Our findings here and in [9] imply that the deflection angle already encodes vast amount of physical information for bound orbits, encouraging independent derivations using numerical and/or self-force methodologies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dynamic Adjustment of Model Parameters Using Assimilated Boundary Layer Height

The NASA Global Modeling and Assimilation Office is exploring the assimilation of boundary layer height (PBLH) derived from a variety of observations, including radiosondes, GNSS-RO, space- and ground-based lidar, and radar wind profilers. Previous work using the Goddard Earth Observing System (GEOS) has shown that within the lower troposphere, the representation of thermodynamic structure is strongly affected by model parameterizations. In assimilation experiments, parameterizations can compensate for analysis thermodynamic increments and reduce the observation information carried forward in time. To mitigate this effect, the present study explores a dynamic parameter adjustment based on the PBLH analysis field. The PBLH analysis increment is used to adjust model parameters that influence the depth of parameterized boundary layer mixing in both stable and convective regimes. Replay (nudging) experiments and cycled data assimilation experiments are used to evaluate the approach. The parameter adjustment is shown to reduce model compensation of the analysis thermodynamic tendencies and bring model parameterized mixing depths into closer agreement with observed PBLH.

Nathan Arnold↗

Structure and dynamics of Saturn's outer magnetosphere and boundary regions

In 1979-1981, the three USA spacecraft Pioneer 11 and Voyagers 1 and 2 discovered and explored the magnetosphere of Saturn to the limited extent possible on flyby trajectories. Considerable variation in the locations of the bow shock (BS) and magnetopause (MP) surfaces were observed in association with variable solar wind conditions and, during the Voyager 2 encounter, possible immersion in Jupiter's distant magnetic tail. The limited number of BS and MP crossings were concentrated near the subsolar region and the dawn terminator, and that fact, together with the temporal variability, makes it difficult to assess the three dimensional shape of the sunward magnetospheric boundary. The combined BS and MP crossing positions from the three spacecraft yield an average BS-to-MP stagnation point distance ratio of 1.29 +/- 0.10. This is near the 1.33 value for the Earth's magnetosphere, implying a similar sunward shape at Saturn. Study of the structure and dynamical behavior of the outer magnetosphere, both in the sunward hemisphere and the magnetotail region using combined plasma and magnetic field data, suggest that Saturn's magnetosphere is more similar to that of Earth than that of Jupiter.

Behannon, K. W.↗

Structure and dynamics of Saturn's outer magnetosphere and boundary regions

In 1979-1981, the three USA spacecraft Pioneer 11 and Voyagers 1 and 2 discovered and explored the magnetosphere of Saturn to the limited extent possible on flyby trajectories. Considerable variation in the locations of the bow shock (BS) and magnetopause (MP) surfaces were observed in association with variable solar wind conditions and, during the Voyager 2 encounter, possible immersion in Jupiter's distant magnetic tail. The limited number of BS and MP crossings were concentrated near the subsolar region and the dawn terminator, and that fact, together with the temporal variability, makes it difficult to assess the three dimensional shape of the sunward magnetospheric boundary. The combined BS and MP crossing positions from the three spacecraft yield an average BS-to-MP stagnation point distance ratio of 1.29 +/- 0.10. This is near the 1.33 value for the Earth's magnetosphere, implying a similar sunward shape at Saturn. Study of the structure and dynamical behavior of the outer magnetosphere, both in the sunward hemisphere and the magnetotail region using combined plasma and magnetic field data, suggest that Saturn's magnetosphere is more similar to that of Earth than that of Jupiter. Previously announced in STAR as N83-30346

Behannon, K. W.↗

Evaluation of Atmospheric Boundary Layer Height From Wind Profiling Radar and Slab Models and Its Responses to Seasonality of Land Cover, Subsidence, and Advection

In this study, we evaluated the effect of land cover, atmospheric subsidence, and advection on the annual dynamics of atmospheric boundary layer (ABL) height from two contrasting sites. The first site is the Walker Branch forest, a deciduous forest of temperate climate, complex topography, and cloudy summers. The second site is the Sacramento-San Joaquin River Delta, a site of Mediterranean climate, flat terrain on a local scale, and clear summers. After testing a new algorithm to calculate ABL heights from 915 MHz radar wind profilers, we evaluated a hierarchy of three slab models to recreate the diurnal and annual patterns of ABL growth. We found that the lower ABL heights in the Delta, particularly during late summer, are driven by the combined effects of increased atmospheric subsidence and marine air advection. In both sites, the annual pattern of ABL height was strongly correlated to total daily incoming radiation, and in the Delta, the annual pattern of ABL height closely followed the seasonal patterns of sensible heat flux from a mosaic of different land covers. A land composite of latent and sensible heat fluxes obtained through a meso-network of eddy covariance measurements and the ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) mission resulted in higher model skill, thus showing that land cover heterogeneity is an important driver of ABL growth. To conclude, model simulations show that in the Delta, restoring agricultural land to wetlands with large open water areas could result in a reduction of ABL height during those months with low subsidence and advection.

54 ENVIRONMENTAL SCIENCES↗

Advanced development of BEM for elastic and inelastic dynamic analysis of solids

Direct Boundary Element formulations and their numerical implementation for periodic and transient elastic as well as inelastic transient dynamic analyses of two-dimensional, axisymmetric and three-dimensional solids are presented. The inelastic formulation is based on an initial stress approach and is the first of its kind in the field of Boundary Element Methods. This formulation employs the Navier-Cauchy equation of motion, Graffi's dynamic reciprocal theorem, Stokes' fundamental solution, and the divergence theorem, together with kinematical and constitutive equations to obtain the pertinent integral equations of the problem in the time domain within the context of the small displacement theory of elastoplasticity. The dynamic (periodic, transient as well as nonlinear transient) formulations have been applied to a range of problems. The numerical formulations presented here are included in the BEST3D and GPBEST systems.

Banerjee, P. K.↗

Evaluating fine-resolution, regional outputs of a variable resolution global climate model

Climate models have been used to study water resources and regional hydrologic responses to climate change, but climate model outputs must be downscaled to provide relevant regional data. However, the accuracy of this regional data is limited by uncertainties across and within downscaling methods, uncertainty across global outputs, and discontinuities at downscaled boundaries. A new alternative to traditional downscaling is a variable resolution model that incorporates fine-resolution regions directly into a coarse-resolution, global climate simulation in order to capture contiguous dynamics across resolution boundaries. In this study, we used the Variable-Resolution Community Earth System Model (VR-CESM) to generate one-eighth degree (14 km) fine-resolution outputs for the western U.S. and eastern China from 1970-2006. We focus our evaluation on precipitaiton, temperature, snow pack, solar radiation, and wind. We compare the model outputs with remote-sensing-based precipitation data, and both reanalysis and gridded weather station data for precipitation and temperature. VR-CESM generally has a cold bias in winter and a warm bias in summer in the western U.S., which compensate each other to reduce the annual bias. In eastern China, however, the sign of temperature biases are more consistent throughout the year with cold biases in the higher mountains and warm biases throughout most of the rest of the region. Precipitation biases are dependent upon reference data, and show slight overestimation in high mountain regions in both the U.S. and China with respect to gridded weather station data. Simulated snow cover in the western U.S. is reasonable compared to remote sensing data, but snow cover and snow water equivalent have larger biases when compared to reanalysis data. In eastern China there are widespread snow cover biases compared to remote sensing data. VR-CESM underestimates downward shortwave radiation to a greater degree in summer than in winter, and underestimates surface layer windspeed over mountains to a greater degree than in other areas. Comparison between VR-CESM and a coarser simulation (1-degree Beijing Climate Center model) shows reduced precipitation biases in the mountainous regions with finer resolution, indicating the value of variable-resolution modeling for reigonal studies.

Di Vittorio, Alan↗

Transonic flutter study of a wind-tunnel model of an arrow-wing supersonic transport

A 1/20-size, low-speed flutter model of the SCAT-15F complete airplane was tested on cables to simulate a near free-flying condition. Only the model wing and fuselage were flexible. Flutter boundaries were measured for a nominal configuration and a configuration with wing fins removed at Mach numbers M from 0.76 to 1.2. For both configurations, the transonic dip in the wing flutter dynamic pressure q boundary was relatively small and the minimum flutter q occurred near M = 0.92. Removing the wing fins increased the flutter q about 14 percent and changed the flutter mode from symmetric to antisymmetric. Vibration and flutter analyses were made using a finite-element structural representation and subsonic kernel-function aerodynamics. For the nominal configuration, the analysis (using calculated modal data) predicted the experimental flutter q levels within 10 percent but did not predict the correct flutter mode at the higher M. For the configuration without wing fins, the analysis predicted 16 to 36 percent unconservative (higher than experimental) flutter q levels and showed extreme sensitivity to mass representation details that affected wing tip mode shapes. For high subsonic M, empennage aerodynamics had a significant effect on the predicted flutter boundaries of several symmetric modes.

Ruhlin, C. L.↗

Transonic flutter study of a wind-tunnel model of an arrow-wing supersonic transport

A 1/20-size, low-speed flutter model of the SCAT-15F complete airplane was tested in a wind-tunnel on cables to simulate a near free-flying condition. Only the model wing and fuselage were flexible. Flutter boundaries were measured for a nominal configuration and a configuration with wing fins removed at Mach numbers M from 0.76 to 1.2. For both configurations, the transonic dip in the wing flutter dynamic pressure q boundary was relatively small and the minimum flutter q occurred near M = 0.92. Removing the wing fins increased the flutter q about 14 percent and changed the flutter mode from symmetric to antisymmetric. Vibration and flutter analyses were made using a finite-element structural representation and subsonic kernel-function aerodynamics. For the nominal configuration, the analysis (using calculated modal data) predicted the experimental flutter q levels within 10 percent but did not predict the correct flutter mode at the higher M. For the configuration without wing fins, the analysis predicted 16 to 36 percent unconservative (higher than experimental) flutter q levels and showed extreme sensitivity to mass shapes.

Ruhlin, C. L.↗

A Poisson equation method for prescribing fully developed non-Newtonian inlet conditions for computational fluid dynamics simulations in models of arbitrary cross-section

Prescribing inlet boundary conditions for computational fluid dynamics (CFD) simulations of internal flow in complex geometries such as anatomical vascular models is challenging. In the absence of patient-specific inlet velocity data, a common approach for long blood vessels is to assume that the inlet flow is fully developed. In vessels of irregular cross section, however, prescribing fully developed conditions is complicated due to the lack of a general closed-form analytical solution. In this study, we develop a simple Poisson equation method for prescribing fully developed inlet conditions for the flow of either Newtonian or non-Newtonian fluids in CFD models of arbitrary cross-section. We first derive the generalized Poisson equation for fully developed flow of a non-Newtonian fluid and we then develop and verify a methodology for numerically computing the solution on any planar boundary domain. In addition, we develop a simple extension of the method for prescribing a non-orthogonal inlet velocity that represents fully developed flow from an upstream tube that is connected to the CFD inlet at a non-orthogonal angle. This may be used to investigate a common source of uncertainty in CFD simulations of internal flow that is due to a lack of information concerning the exact streamwise flow direction at the inlets. Comparison to several Newtonian and non-Newtonian benchmark verification solutions shows the method to be extremely accurate. As a practical demonstration case, we use the method to prescribe fully developed conditions on multiple non-circular inlets for the non-Newtonian flow of blood in a patient-specific model of the inferior vena cava (IVC). Finally, we further demonstrate the utility of the method by performing a sensitivity study using the patient-specific IVC model, wherein we investigate the influence of inlet velocity flow direction on the non-Newtonian IVC hemodynamics. Given its simplicity and computational efficiency, the method is shown to be far superior to alternative approaches for prescribing fully developed inlet conditions in such complicated geometries. In conclusion, to facilitate the adoption of our Poisson equation method, we have distributed our OpenFOAM source code and the associated test cases from this study as open-source software.

97 MATHEMATICS AND COMPUTING↗

Investigation of Inlet Concepts for Maneuver Improvement at Transonic Speeds

A 15 percent scale lightweight fighter type inlet forebody was tested in the Ames 14 foot transonic wind tunnel at Mach numbers of 0.7, 0.9, and 1.04. The inlet was a two dimensional horizontal ramp system designed for a Mach number of 2.2. Four inlet devices designed to prevent or delay cowl-lip boundary layer separation or to improve the inlet internal flow characteristics at high angles of attack were investigated. The devices used to control cowl-lip separation consisted of cowl leading edge flaps, slotted flaps, and tangential blowing. To improve the internal flow characteristics, discrete jet nozzle flows were directed downstream and parallel to the duct surface in the subsonic diffuser to energize the wall boundary layer. The discrete jets used in the subsonic diffuser were also tested in combination with each of the cowl leading edge devices. Test measurements included engine-face total pressure recovery, steady state distortion, dynamic distortion, duct boundary layer profiles, and duct-surface static pressures.

Latham, E.↗

Modeling of an Atmospheric-Boundary-Layer Profile in Support of Experiments in the NASA Langley Transonic Dynamics Tunnel

This paper presents the results of the Computational Fluid Dynamics (CFD) analyses to model two atmospheric-boundary-layer (ABL) profiles inside the NASA Langley Research Center Transonic Dynamics Tunnel (TDT). The CFD models include tunnel walls and all additional hardware needed to simulate the ABL profiles. This hardware includes Irwin spires and floor-mounted roughness elements. The numerical simulations show that the application of higher-fidelity numerical methods is necessary to compute boundary-layer and turbulenceintensity profiles that match experimental data. The ABL profiles are computed both inside the numerical model of the TDT and in a classical free-air model. Both Unsteady Reynolds Averaged Navier-Stokes (URANS) with Spalart-Allmaras (SA) turbulence model and Modified Delayed Detached Eddy (MDDES) simulation methods are used. The results show that the MDDES-simulation results match the experimental data very well while URANS-SA does not. The results also show that the wind-tunnel walls have a significant effect on ABL prediction.

Chwalowski, Pawel↗

Numerical approximation of boundary conditions with applications to inviscid equations of gas dynamics

A comprehensive overview of the state of the art of well-posedness and stability analysis of difference approximations for initial boundary value problems of the hyperbolic type is presented. The applicability of recent theoretical development to practical calculations for nonlinear gas dynamics is examined. The one dimensional inviscid gas dynamics equations in conservation law form are selected for numerical experiments. The class of implicit schemes developed from linear multistep methods in ordinary differential equations is chosen and the use of linear extrapolation as an explicit or implicit boundary scheme is emphasized. Specification of boundary data in the primitive variables and computation in terms of the conservative variables in the interior is discussed. Some numerical examples for the quasi-one-dimensional nozzle are given.

Yee, H. C.↗

Representation of Nucleation Mode Microphysics in a Global Aerosol Model with Sectional Microphysics

In models, nucleation mode (1 nm<Dp <10 nm) particle microphysics can be represented explicitly with aerosol microphysical processes or can be parameterized to obtain the growth and survival of nuclei to the model's lower size boundary. This study investigates how the representation of nucleation mode microphysics impacts aerosol number predictions in the TwO-Moment Aerosol Sectional (TOMAS) aerosol microphysics model running with the GISS GCM II-prime by varying its lowest diameter boundary: 1 nm, 3 nm, and 10 nm. The model with the 1 nm boundary simulates the nucleation mode particles with fully resolved microphysical processes, while the model with the 10 nm and 3 nm boundaries uses a nucleation mode dynamics parameterization to account for the growth of nucleated particles to 10 nm and 3 nm, respectively.We also investigate the impact of the time step for aerosol microphysical processes (a 10 min versus a 1 h time step) to aerosol number predictions in the TOMAS models with explicit dynamics for the nucleation mode particles (i.e., 3 nm and 1 nm boundary). The model with the explicit microphysics (i.e., 1 nm boundary) with the 10 min time step is used as a numerical benchmark simulation to estimate biases caused by varying the lower size cutoff and the time step. Different representations of the nucleation mode have a significant effect on the formation rate of particles larger than 10 nm from nucleated particles (J10) and the burdens and lifetimes of ultrafinemode (10 nm≤Dp ≤70 nm) particles but have less impact on the burdens and lifetimes of CCN-sized particles. The models using parameterized microphysics (i.e., 10 nm and 3 nm boundaries) result in higher J10 and shorter coagulation lifetimes of ultrafine-mode particles than the model with explicit dynamics (i.e., 1 nm boundary). The spatial distributions of CN10 (Dp ≥10 nm) and CCN(0.2 %) (i.e., CCN concentrations at 0.2%supersaturation) are moderately affected, especially CN10 predictions above ∼700 hPa where nucleation contributes most strongly to CN10 concentrations. The lowermost-layer CN10 is substantially improved with the 3 nm boundary (compared to 10 nm) in most areas. The overprediction in CN10 with the 3 nm and 10 nm boundaries can be explained by the overprediction of J10 or J3 with the parameterized microphysics, possibly due to the instantaneous growth rate assumption in the survival and growth parameterization. The errors in CN10 predictions are sensitive to the choice of the lower size boundary but not to the choice of the time step applied to the microphysical processes. The spatial distribution of CCN(0.2 %) with the 3 nm boundary is almost identical to that with the 1 nm boundary, but that with the 10 nm boundary can differ more than 10-40% in some areas. We found that the deviation in the 10 nm simulations is partly due to the longer time step (i.e., 1 h time step used in the 10 nm simulations compared to 10 min time step used in the benchmark simulations), but, even with the same time step, the 10 nm cutoff showed noticeably higher errors than the 3 nm cutoff. In conclusion, we generally recommend using a lower diameter boundary of 3 nm for studies focused on aerosol indirect effects but down to 1 nm boundary for studies focused on CN10 predictions or nucleation.

microphysics↗

Large-Scale Structure of Magnetospheric Plasma

Recent investigations of magnetospheric plasma structure are summarized under the broad categories of empirical models, transport across boundaries, formation, and dynamics of the plasma sheet. This report reviews work in these areas during the period 1991 to 1993. Fully three-dimensional empirical models and simulations have become important contributors to our understanding of the magnetospheric system. Some new structural concepts have appeared in the literature: the 'entry boundary' and 'geo-pause', the plasma sheet 'region 1 vortices', the 'low-energy layer', the 'adia-baticity boundary' or 'wall region', and a region in the tail to which we refer as the 'injection port'. Traditional structural concepts have also been the subject of recent study, notably the plasmapause, the magnetopause, and the plasma sheet. Significant progress has been made in understanding the nature of plasma sheet formation and dynamics, but the acceleration of electrons to high energy remains somewhat mysterious.

Moore, T. E.↗

Performance Analysis and Tank Test Validation of a Hybrid Wave-Current Energy Converter with a Single Power Takeoff

Marine and hydrokinetic (MHK) energy, including ocean waves, tidal current, ocean current and river current, has been recognized as a promising power source due to its full-day availability and high energy potential. At this stage, ocean current energy, tidal energy and ocean wave energy are currently the most competitive sourves among all the categories of MHK. The state of art MHK energy harvesting technology mainly focus on harvesting either ocean wave energy or current energy, but not both. However, a significant amount of ocean waves and tidal/ ocean current coexist in many sites and traditional devices that harvest from a single form of MHK energy, cannot make full use of the coexisting ocean energy. Furthermore, MHK energy harvesting devices need to advance to be cost-effective and competitive with other energy sources. This is difficult to achieve. Ocean wave excitation is irregular, which means that ocean wave height and wave periods are unpredictable and excitation forces on energy harvesting devices can have large variance in amplitude and frequency. Such problems/ restrictions can be possibly addressed by the concept of a hybrid energy converter. In this sense, a hybrid wave-current ocean energy conveter (HWCEC) that simutaneously harvests energy from current and wave with one single power takeoff (PTO) is designed.The wave energy is extracted through relative heaving motion between a floating buoy and a submerged second body, while the current energy is extracted using a marine current turbine (MCT). Energy from both sources are integrated by a hybrid PTO whose concept is based on a mechanical motion rectifier (MMR). In this study, different working modes are investigated together with switching criteria.Simulations were conducted with hydrodynamic coefficients obtained from computational fluid dynamics analysis and boundary element method. Tank tests were conducted for a HWCEC under coexisting wave and current inputs. For comparison, separate baseline tests of a turbine and a two-body point absorber, each acting in isolation, are conducted. Experimental results validate the dynamic modeling and show that a HWCEC can increase the output power with a range between 29-87 percent over either current turbine and wave energy converter acting individually, and it can reduce by up to 70 percent the peak-to-average power ratio compared with the wave energy converter on the tested conditions.Such results demonstrate the potential of the HWCEC as an efficient and cost-effective design.

16 TIDAL AND WAVE POWER↗

Measuring dynamics of differentially rotating, unmagnetized, free-boundary plasma produced by soft x-ray irradiation

We present results from a new experiment, fielded on the MAGPIE pulsed power generator, producing a rotating, boundary-free plasma with no pre-imposed magnetic field. Angular momentum is introduced to the rotating system by the oblique collision of multiple plasma flows, which are driven by x-ray ablation of solid targets using the soft x-rays emitted from the implosion of wire array z-pinches. This produces a hydrodynamically stable plasma which undergoes ~2-3 rotations over the duration of the experiment, significantly more than previous pulsed-power platforms. Estimating the angular frequency from the electron density profile in the inner part of the rotating plasma shows that it has a quasi-Keplerian rotation profile. The system also allows for the addition of a controllable magnetic field. This, combined with the sufficiently large Reynolds (~10 5 ) and magnetic Reynolds (~10) numbers, will enable investigation of the effect of magnetic field on the structure and stability of the rotating plasma.

Physics - Plasma physics↗