Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “local motion stabilization”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Size limits of magnetic-domain engineering in continuous in-plane exchange-bias prototype films

Background: The application of superparamagnetic particles as biomolecular transporters in microfluidic systems for lab-on-a-chip applications crucially depends on the ability to control their motion. One approach for magnetic-particle motion control is the superposition of static magnetic stray field landscapes (MFLs) with dynamically varying external fields. These MFLs may emerge from magnetic domains engineered both in shape and in their local anisotropies. Motion control of smaller beads does necessarily need smaller magnetic patterns, i.e., MFLs varying on smaller lateral scales. The achievable size limit of engineered magnetic domains depends on the magnetic patterning method and on the magnetic anisotropies of the material system. Smallest patterns are expected to be in the range of the domain wall width of the particular material system. To explore these limits a patterning technology is needed with a spatial resolution significantly smaller than the domain wall width. Results: We demonstrate the application of a helium ion microscope with a beam diameter of 8 nm as a mask-less method for local domain patterning of magnetic thin-film systems. For a prototypical in-plane exchange-bias system the domain wall width has been investigated as a function of the angle between unidirectional anisotropy and domain wall. By shrinking the domain size of periodic domain stripes, we analyzed the influence of domain wall overlap on the domain stability. Finally, by changing the geometry of artificial two-dimensional domains, the influence of domain wall overlap and domain wall geometry on the ultimate domain size in the chosen system was analyzed. Conclusion: The application of a helium ion microscope for magnetic patterning has been shown. It allowed for exploring the fundamental limits of domain engineering in an in-plane exchange-bias thin film as a prototypical system. For two-dimensional domains the limit depends on the domain geometry. The relative orientation between domain wall and anisotropy axes is a crucial parameter and therefore influences the achievable minimum domain size dramatically.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Coastal Ohio Wind Project

The Coastal Ohio Wind Project intends to address problems that impede deployment of wind turbines in the coastal and offshore regions of Northern Ohio. The project evaluates different wind turbine designs and the potential impact of offshore turbines on migratory and resident birds by developing multidisciplinary research, which involves wildlife biology, electrical and mechanical engineering, and geospatial science. Firstly, the project conducts cost and performance studies of two- and three-blade wind turbines using a turbine design suited for the Great Lakes. The numerical studies comprised an analysis and evaluation of the annual energy production of two- and three-blade wind turbines to determine the levelized cost of energy. This task also involved wind tunnel studies of model wind turbines to quantify the wake flow field of upwind and downwind wind turbine-tower arrangements. The experimental work included a study of a scaled model of an offshore wind turbine platform in a water tunnel. The levelized cost of energy work consisted of the development and application of a cost model to predict the cost of energy produced by a wind turbine system placed offshore. The analysis found that a floating two-blade wind turbine presents the most cost effective alternative for the Great Lakes. The load effects studies showed that the two-blade wind turbine model experiences less torque under all IEC Standard design load cases considered. Other load effects did not show this trend and depending on the design load cases, the two-bladed wind turbine showed higher or lower load effects. The experimental studies of the wake were conducted using smoke flow visualization and hot wire anemometry. Flow visualization studies showed that in the downwind turbine configuration the wake flow was insensitive to the presence of the blade and was very similar to that of the tower alone. On the other hand, in the upwind turbine configuration, increasing the rotor blade angle of attack reduced the wake size and enhanced the vortices in the flow downstream of the turbine-tower compared with the tower alone case. Mean and rms velocity distributions from hot wire anemometer data confirmed that in a downwind configuration, the wake of the tower dominates the flow, thus the flow fields of a tower alone and tower-turbine combinations are nearly the same. For the upwind configuration, the mean velocity shows a narrowing of the wake compared with the tower alone case. The downwind configuration wake persisted longer than that of an upwind configuration; however, it was not possible to quantify this difference because of the size limitation of the wind tunnel downstream of the test section. The water tunnel studies demonstrated that the scale model studies could be used to adequately produce accurate motions to model the motions of a wind turbine platform subject to large waves. It was found that the important factors that affect the platform is whether the platform is submerged or surface piercing. In the former, the loads on the platform will be relatively reduced whereas in the latter case, the structure pierces the wave free surface and gains stiffness and stability. The other important element that affects the movement of the platform is depth of the sea in which the wind turbine will be installed. Furthermore, the wildlife biology component evaluated migratory patterns by different monitoring systems consisting of marine radar, thermal IR camera and acoustic recorders. The types of radar used in the project are weather surveillance radar and marine radar. The weather surveillance radar (1988 Doppler), also known as Next Generation Radar (NEXRAD), provides a network of weather stations in the US. Data generated from this network were used to understand general migratory patterns, migratory stopover habitats, and other patterns caused by the effects of weather conditions. At a local scale our marine radar was used to complement the datasets from NEXRAD and to collect additional monitoring parameters such as passage rates, flight paths, flight directions, and flight altitudes of nocturnal migrating species. Our work focused on the design and development of custom built marine radar that used t-bar and parabolic dish antennas. The marine radar used in the project was Furuno (XANK250) which was coupled with a XIR3000B digitizing card from Russell Technologies for collection of the radar data. The radar data was processed by open source radR processing software using different computational techniques and methods. Additional data from thermal IR imaging cameras were collected to detect heat emitted from objects and provide information on movements of birds and bats, data which we used for different animal flight behavior analysis. Lastly, the data from the acoustic recorders were used to provide the number of bird calls for assessing patterns and peak passage rates during migration. The development of the geospatial database included collection of different data sources that are used to support offshore wind turbine development. Many different data sets were collected and organized using initial version of web-based repository software tools that can accommodate distribution of rectified pertinent data sets such as the lake depth, lake bottom engineering parameters, extent of ice, navigation pathways, wind speed, important bird habitats, fish efforts and other layers that are relevant for supporting robust offshore wind turbine developments. Additional geospatial products developed during the project included few different prototypes for offshore wind farm suitability which can involve different stakeholders and participants for solving complex planning problems and building consensus. Some of the prototypes include spatial decision support system (SDSS) for collaborative decision making, a web-based Participatory Geographic Information System (PGIS) framework for evaluating importance of different decision alternatives using different evaluation criteria, and an Android application for collection of field data using mobile and tablet devices . In summary, the simulations of two- and three-blade wind turbines suggested that two-bladed machines could produce comparable annual energy as the three-blade wind turbines but have a lighter tower top weight, which leads to lower cost of energy. In addition, the two-blade rotor configuration potentially costs 20% less than a three blade configuration that produces the same power at the same site. The cost model analysis predicted a potential cost savings of approximately 15% for offshore two-blade wind turbines. The foundation design for a wind turbine in Lake Erie is likely to be driven by ice loads based on the currently available ice data and ice mechanics models. Hence, for Lake Eire, the cost savings will be somewhat smaller than the other lakes in the Great Lakes. Considering the size of cranes and vessels currently available in the Great Lakes, the cost optimal wind turbine size should be 3 MW, not larger. The surveillance data from different monitoring systems suggested that bird and bat passage rates per hour were comparable during heavy migrations in both spring and fall seasons while passage rates were significantly correlated to wind directions and wind speeds. The altitude of migration was higher during heavy migrations and higher over water relative to over land. Notable portions of migration on some spring nights occurred parallel the shoreline, often moving perpendicular to southern winds. The birds approaching the Western basin have a higher propensity to cross than birds approaching the Central basin of Lake Erie and as such offshore turbine development might be a better option further east towards Cleveland than in the Western basin. The high stopover density was more strongly associated with migration volume the following night rather than the preceding night. The processed mean scalar wind speeds with temporal resolutions as fine as 10-minute intervals near turbine height showed that August is the month with the weakest winds while December is the month, which typically has the strongest winds. The ice data suggests that shallow western basin of Lake Erie has higher ice cover duration many times exceeding 90 days during some winters.

17 WIND ENERGY↗
NASA NTRS1999

Eigenmodes of Ducted Flows With Radially-Dependent Axial and Swirl Velocity Components

This report characterizes the sets of small disturbances possible in cylindrical and annular ducts with mean flow whose axial and tangential components vary arbitrarily with radius. The linearized equations of motion are presented and discussed, and then exponential forms for the axial, circumferential, and time dependencies of any unsteady disturbances are assumed. The resultant equations form a generalized eigenvalue problem, the solution of which yields the axial wavenumbers and radial mode shapes of the unsteady disturbances. Two numerical discretizations are applied to the system of equations: (1) a spectral collocation technique based on Chebyshev polynomial expansions on the Gauss-Lobatto points, and (2) second and fourth order finite differences on uniform grids. The discretized equations are solved using a standard eigensystem package employing the QR algorithm. The eigenvalues fall into two primary categories: a discrete set (analogous to the acoustic modes found in uniform mean flows) and a continuous band (analogous to convected disturbances in uniform mean flows) where the phase velocities of the disturbances correspond to the local mean flow velocities. Sample mode shapes and eigensystem distributions are presented for both sheared axial and swirling flows. The physics of swirling flows is examined with reference to hydrodynamic stability and completeness of the eigensystem expansions. The effect of assuming exponential dependence in the axial direction is discussed.

Kousen, Kenneth A.↗
DOE OSTI2019

Biochemical and physiological flexibility accompanies reduced cellulose biosynthesis in Brachypodium cesa1 S830N

Here, we present a study into the mechanisms of primary cell wall cellulose formation in grasses, using the model cereal grass Brachypodium distachyon. The exon found adjacent to the BdCESA1 glycosyltransferase QXXRW motif was targeted using Targeting Induced Local Lesions in Genomes (TILLING) and sequencing candidate amplicons in multiple parallel reactions (SCAMPRing) leading to the identification of the Bdcesa1 S830N allele. Plants carrying this missense mutation exhibited a significant reduction in crystalline cellulose content in tissues that rely on the primary cell wall for biomechanical support. However, Bdcesa1 S830N plants failed to exhibit the predicted reduction in plant height. In a mechanism unavailable to eudicotyledons, B. distachyon plants homozygous for the Bdcesa1 S830N allele appear to overcome the loss of internode expansion anatomically by increasing the number of nodes along the stem. Stem biomechanics were resultantly compromised in Bdcesa1 S830N . The Bdcesa1 S830N missense mutation did not interfere with BdCESA1 gene expression. However, molecular dynamic simulations of the CELLULOSE SYNTHASE A (CESA) structure with modelled membrane interactions illustrated that Bdcesa1 S830N exhibited structural changes in the translated gene product responsible for reduced cellulose biosynthesis. Molecular dynamic simulations showed that substituting S830N resulted in a stabilizing shift in the flexibility of the class specific region arm of the core catalytic domain of CESA, revealing the importance of this motion to protein function.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2000

A Dynamical Model for The Jupiter Thermocline

A new conceptual model for the observed banded wind structure of the giant outer planets is proposed, assuming a smoothly distributed potential vorticity (PV) atop a layer of deeply seated, latitudinally variable stratification. With a vanishing horizontal entropy contrast presumably enforced at the bottom of the flow layer by the underlying convection, the thermal-wind balance of observed cloud-top motions implies a mapping of constant potential temperature surfaces, nearly vertical at upper tropospheric levels, down to a deeper, flatter, but variably sloped 'thermocline.' For a generally colder-poleward distribution of isentropes consistent with the strong equatorial jets of both Jupiter and Saturn, this kind of mapping would also imply a poleward decreasing static stability. The proposed temperature-stability distribution is just the reverse of the Earth's troposphere, where the strongest latitudinal potential temperature gradients are at the bottom instead of the top, with static stability generally increasing toward the pole. The warmer-stabler correlation for Jupiter would be consistent, however, with a nearly monotonic distribution of PV from low to high latitudes, over planetary scales for which the planetary vorticity dominates the relative vorticity of the jets. In this way efficient PV mixing for the isentropically bounded thermocline can account for the dynamical maintenance of the cyclonic flanks of the equatorial jet. Over latitudinal intervals comparable to the internal deformation scale, however, the gradient of the absolute vorticity is dominated by the flow curvature, with the local latitudinal gradient of the static stability proportional to the local potential temperature gradient and therefore the local variation in geostrophic velocity. The resulting correlation of velocities and vorticity gradients implies an alternation of the flow over the internal deformation radius set by the vertical entropy contrast.

Allison, Michael↗
NASA NTRS1979

The dynamics and control of large flexible space structures, 2. Part A: Shape and orientation control using point actuators

The equations of planar motion for a flexible beam in orbit which includes the effects of gravity gradient torques and control torques from point actuators located along the beam was developed. Two classes of theorems are applied to the linearized form of these equations to establish necessary conditions for controlability for preselected actuator configurations. The feedback gains are selected: (1) based on the decoupling of the original coordinates and to obtain proper damping, and (2) by applying the linear regulator problem to the individual model coordinates separately. The linear control laws obtained using both techniques were evaluated by numerical integration of the nonlinear system equations. Numerical examples considering pitch and various number of modes with different combination of actuator numbers and locations are presented. The independent model control concept used earlier with a discretized model of the thin beam in orbit was reviewed for the case where the number of actuators is less than the number of modes. Results indicate that although the system is controllable it is not stable about the nominal (local vertical) orientation when the control is based on modal decoupling. An alternate control law not based on modal decoupling ensures stability of all the modes.

Bainum, P. M.↗
NASA NTRS2003

The Cytoskeleton and Force Response Mechanisms

The long term aim of this project was to define the mechanisms by which cells sense and respond to the physical forces experienced at 1g and missing in microgravity. Identification and characterization of the elements of the cells force response mechanism could provide pathways and molecules to serve as targets for pharmacological intervention to mitigate the pathologic effects of microgravity. Mechanical forces experienced by the organism can be transmitted to cells through molecules that allow cells to bind to the extracellular matrix and through other types of molecules which bind cells to each other. These molecules are coupled in large complexes of proteins to structural elements such as the actin cytoskeleton that give the cell the ability to sense, resist and respond to force. Application of small forces to tissue culture cells causes local elevation of intracellular calcium through stretch activated ion channels, increased tyrosine phosphorylation and a restructuring of the actin cytoskeleton. Using collagen coated iron oxide beads and strong magnets, we can apply different levels of force to cells in culture. We have found that force application causes the cells to polymerize actin at the site of mechanical deformation and unexpectedly, to depolymerize actin across the rest of the cell. Observations of GFP- actin expressing cells demonstrate that actin accumulates at the site of deformation within the first five minutes of force application and is maintained for many tens of minutes after force is removed. Consistent with the reinforcement of the cytoskeletal structures underlying the integrin-bead interaction, force also alters the motion of bound magnetic beads. This effect is seen following the removal of the magnetic field, and is only partially ablated by actin disruption with cytochalsin B. While actin is polymerizing locally at the site of force application, force also stimulates a global reduction in actin filament content within the cells. We have examined the roles of several actin filament disassembly factors in the global reduction of cellular actin filaments. The calcium regulated actin filament severing protein gelsolin is not necessary for the increased actin turnover, as cells derived from gelsolin null and wildtype mice still show a reduction in total actin filament content. Instead, our work suggests that the actin binding protein cofilin may be important for these changes in actin dynamics. Cofilin binds to and enhances the disassembly of actin filaments. Using immunological methods, we observe transient changes in the phosphorylation state of cofilin upon force application that suggests that cofilin may mediate actin filament turnover. Early after force application, cofilin is transiently dephosphorylated, activating its actin disassembly activity. Subsequently, we find a hyper-phosphorylation of cofilin, rendering it inactive. This reduction in cofilin activity may explain the stability of the force induced actin structuttes. In testing this hypothesis, we aimed to generate cells that express the constituitively active kinase (LIM-kinase) that phosphorylates cofilin. lnidial attempts in the cell lines used for the our previous studies proved unsuccessful. While we prepare this work for pubication, we are continuing to study other cell lines and tissue sources to determine whether they show a reduction in F-actin content after force application.

Allen, Philip Goodwin↗
DOE OSTI2017

Non-steady wind turbine response to daytime atmospheric turbulence

Relevant to drivetrain bearing fatigue failures, we analyse non-steady wind turbine responses from interactions between energy-dominant daytime atmospheric turbulence eddies and the rotating blades of a GE 1.5 MW wind turbine using a unique dataset from a GE field experiment and computer simulation. Time-resolved local velocity data were collected at the leading and trailing edges of an instrumented blade together with generator power, revolutions per minute, pitch and yaw. Wind velocity and temperature were measured upwind on a meteorological tower. The stability state and other atmospheric conditions during the field experiment were replicated with a large-eddy simulation in which was embedded a GE 1.5 MW wind turbine rotor modelled with an advanced actuator line method. Both datasets identify three important response time scales: advective passage of energy-dominant eddies (≈25–50 s), blade rotation (once per revolution (1P), ≈3 s) and sub-1P scale (<1 s) response to internal eddy structure. Large-amplitude short-time ramp-like and oscillatory load fluctuations result in response to temporal changes in velocity vector inclination in the aerofoil plane, modulated by eddy passage at longer time scales. Generator power responds strongly to large-eddy wind modulations. Lastly, we show that internal dynamics of the blade boundary layer near the trailing edge is temporally modulated by the non-steady external flow that was measured at the leading edge, as well as blade-generated turbulence motions.

17 WIND ENERGY↗
DOE OSTI2024

Tailoring Ion Transport in Li 3‐3y Ho 1+y Cl 6‐x Br x via Transition‐Metal Free Structural Planes and Charge Carrier Distribution

Abstract Localized atomistic disorder in halide‐based solid electrolytes (SEs) can be leveraged to boost Li + mobility. In this study, Li + transport in structurally modified Li 3 HoCl 6 , via Br − introduction and Li + deficiency, is explored. The optimized Li 3‐3 y Ho 1+ y Cl 6‐ x Br x achieves an ionic conductivity of 3.8 mS cm −1 at 25 °C, the highest reported for holmium halide materials. 6,7 Li nuclear magnetic resonance and relaxometry investigations unveil enhanced ion dynamics with bromination, attaining a Li + motional rate neighboring 116 MHz. X‐ray diffraction analyses reveal mixed‐anion‐induced phase transitions with disproportionate octahedral expansions and distortions, creating Ho‐free planes with favorable energetics for Li + migration. Bond valence site energy analysis highlights preferred Li + transport pathways, particularly in structural planes devoid of Ho 3+ blocking effects. Molecular dynamics simulations corroborate enhanced Li + diffusion with Br − introduction into Li 3 HoCl 6 . Li‐Ho electrostatic repulsions in the (001) plane presumably drive Li + diffusion into the Ho‐free (002) layer, enabling rapid intraplanar Li + motion and exchange between the 2d and 4h sites. Li 3‐3 y Ho 1+ y Cl 6‐ x Br x also demonstrates good battery cycling stability. These findings offer valuable insights into the intricate correlations between structure and ion transport and will help guide the design of high‐performance fast ion conductors for all‐solid‐state batteries.

Chemistry↗
DOE OSTI2016

Aging kinetics of levoglucosan orientational glass as a rate dispersion process and consequences for the heterogeneous dynamics view

Aging kinetics of a glass is currently modeled in terms of slowing of its α-relaxation dynamics, whose features are interpreted in terms of dynamic heterogeneity, i.e., formation and decay of spatially and temporally distinct nm-size regions. To test the merits of this view, we studied the calorimetric effects of aging an orientational glass of levoglucosan crystal in which such regions would not form in the same way as they form in liquids, and persist in structural glasses, because there is no liquid-like molecular diffusion in the crystal. By measuring the heat capacity, C{sub p}, we determined the change in the enthalpy, H, and the entropy, S, during two aging-protocols: (a) keeping the samples isothermally at temperature, T{sub a}, and measuring the changes after different aging times, t{sub a}, and (b) keeping the samples at different T{sub a}s and measuring the changes after the same t{sub a}. A model-free analysis of the data shows that as t{sub a} is increased (procedure (a)), H and S decrease according to a dispersive rate kinetics, and as T{sub a} is increased (procedure (b)), H and S first increase, reach a local maximum at a certain T{sub a}, and then decrease. Even though there is no translational diffusion to produce (liquid-like) free volume, and no translational-rotational decoupling, the aging features are indistinguishable from those of structural glasses. We also find that the Kohlrausch parameter, originally fitted to the glass-aging data, decreases with decrease in T{sub a}, which is incompatible with the current use of the aging data for estimating the α-relaxation time. We argue that the vibrational state of a glass is naturally incompatible with its configurational state, and both change on aging until they are compatible, in the equilibrium liquid. So, dipolar fluctuations seen as the α-relaxation would not be the same motions that cause aging. We suggest that aging kinetics is intrinsically dispersive with its own characteristic rate constant and it does not yield the α-relaxation rate. In this view, thermodynamic and other properties define the fictive temperature; the real or imaginary components of a dynamic property do not define it. While particles’ overall motions may still play a crucial role in (structural) glass physics, we conclude that translational diffusion alone is not a requirement for structure stabilization on aging of a kinetically frozen state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

A method of surface marker location optimization for tumor motion estimation in lung stereotactic body radiation therapy

Purpose: Accurately localizing lung tumor localization is essential for high-precision radiation therapy techniques such as stereotactic body radiation therapy (SBRT). Since direct monitoring of tumor motion is not always achievable due to the limitation of imaging modalities for treatment guidance, placement of fiducial markers on the patient’s body surface to act as a surrogate for tumor position prediction is a practical alternative for tracking lung tumor motion during SBRT treatments. In this work, the authors propose an innovative and robust model to solve the multimarker position optimization problem. The model is able to overcome the major drawbacks of the sparse optimization approach (SOA) model. Methods: The principle-component-analysis (PCA) method was employed as the framework to build the authors’ statistical prediction model. The method can be divided into two stages. The first stage is to build the surrogate tumor matrix and calculate its eigenvalues and associated eigenvectors. The second stage is to determine the “best represented” columns of the eigenvector matrix obtained from stage one and subsequently acquire the optimal marker positions as well as numbers. Using 4-dimensional CT (4DCT) and breath hold CT imaging data, the PCA method was compared to the SOA method with respect to calculation time, average prediction accuracy, prediction stability, noise resistance, marker position consistency, and marker distribution. Results: The PCA and SOA methods which were both tested were on all 11 patients for a total of 130 cases including 4DCT and breath-hold CT scenarios. The maximum calculation time for the PCA method was less than 1 s with 64 752 surface points, whereas the average calculation time for the SOA method was over 12 min with 400 surface points. Overall, the tumor center position prediction errors were comparable between the two methods, and all were less than 1.5 mm. However, for the extreme scenarios (breath hold), the prediction errors for the PCA method were not only smaller, but were also more stable than for the SOA method. Results obtained by imposing a series of random noises to the surrogates indicated that the PCA method was much more noise resistant than the SOA method. The marker position consistency tests using various combinations of 4DCT phases to construct the surrogates suggested that the marker position predictions of the PCA method were more consistent than those of the SOA method, in spite of surrogate construction. Marker distribution tests indicated that greater than 80% of the calculated marker positions fell into the high cross correlation and high motion magnitude regions for both of the algorithms. Conclusions: The PCA model is an accurate, efficient, robust, and practical model for solving the multimarker position optimization problem to predict lung tumor motion during SBRT treatments. Due to its generality, PCA model can also be applied to other imaging guidance system whichever using surface motion as the surrogates.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Integrated modeling applications for tokamak experiments with OMFIT

We developed one modeling framework for integrated tasks (OMFIT) is a comprehensive integrated modeling framework which will enable physics codes to interact in complicated workflows, and support scientists at all stages of the modeling cycle. The OMFIT development follows a unique bottom-up approach, where the framework design and capabilities organically evolve to support progressive integration of the components that are required to accomplish physics goals of increasing complexity. OMFIT provides a workflow for easily generating full kinetic equilibrium reconstructions that are constrained by magnetic and motional Stark effect measurements, and kinetic profile information that includes fast-ion pressure modeled by a transport code. It was found that magnetic measurements can be used to quantify the amount of anomalous fast-ion diffusion that is present in DIII-D discharges, and provide an estimate that is consistent with what would be needed for transport simulations to match the measured neutron rates. OMFIT was used to streamline edge-stability analyses, and evaluate the effect of resonant magnetic perturbation (RMP) on the pedestal stability, which have been found to be consistent with the experimental observations. Moreover, the development of a five-dimensional numerical fluid model for estimating the effects of the interaction between magnetohydrodynamic (MHD) and microturbulence, and its systematic verification against analytic models was also supported by the framework. OMFIT was used for optimizing an innovative high-harmonic fast wave system proposed for DIII-D. For a parallel refractive index n(parallel to) > 3, the conditions for strong electron-Landau damping were found to be independent of launched n(parallel to) and poloidal angle. OMFIT has been the platform of choice for developing a neural-network based approach to efficiently perform a non-linear multivariate regression of local transport fluxes as a function of local dimensionless parameters. Transport predictions for thousands of DIII-D discharges showed excellent agreement with the power balance calculations across the whole plasma radius and over a broad range of operating regimes. Concerning predictive transport simulations, the framework made possible the design and automation of a workflow that enables self-consistent predictions of kinetic profiles and the plasma equilibrium. It is found that the feedback between the transport fluxes and plasma equilibrium can significantly affect the kinetic profiles predictions. Finally, such a rich set of results provide tangible evidence of how bottom-up approaches can potentially provide a fast track to integrated modeling solutions that are functional, cost-effective, and in sync with the research effort of the community.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Flow visualizations of perpendicular blade vortex interactions

Helium bubble flow visualizations have been performed to study perpendicular interaction of a turbulent trailing vortex and a rectangular wing in the Virginia Tech Stability Tunnel. Many combinations of vortex strength, vortex-blade separation (Z(sub s)) and blade angle of attack were studied. Photographs of representative cases are presented. A range of phenomena were observed. For Z(sub s) greater than a few percent chord the vortex is deflected as it passes the blade under the influence of the local streamline curvature and its image in the blade. Initially the interaction appears to have no influence on the core. Downstream, however, the vortex core begins to diffuse and grow, presumably as a consequence of its interaction with the blade wake. The magnitude of these effects increases with reduction in Z(sub s). For Z(sub s) near zero the form of the interaction changes and becomes dependent on the vortex strength. For lower strengths the vortex appears to split into two filaments on the leading edge of the blade, one passing on the pressure and one passing on the suction side. At higher strengths the vortex bursts in the vicinity of the leading edge. In either case the core of its remnants then rapidly diffuse with distance downstream. Increase in Reynolds number did not qualitatively affect the flow apart from decreasing the amplitude of the small low-frequency wandering motions of the vortex. Changes in wing tip geometry and boundary layer trip had very little effect.

Rife, Michael C.↗
DOE OSTI2014

The onset of spiral structure in the universe

The onset of spiral structure in galaxies appears to occur between redshifts 1.4 and 1.8 when disks have developed a cool stellar component, rotation dominates over turbulent motions in the gas, and massive clumps become less frequent. During the transition from clumpy to spiral disks, two unusual types of spirals are found in the Hubble Ultra Deep Field that are massive, clumpy, and irregular like their predecessor clumpy disks, yet spiral-like or sheared like their descendants. One type is 'woolly' with massive clumpy arms all over the disk and is brighter than other disk galaxies at the same redshift, while another type has irregular multiple arms with high pitch angles, star formation knots, and no inner symmetry like today's multiple-arm galaxies. The common types of spirals seen locally are also present in a redshift range around z ∼ 1, namely grand design with two symmetric arms, multiple arm with symmetry in the inner parts and several long, thin arms in the outer parts, and flocculent, with short, irregular, and patchy arms that are mostly from star formation. Normal multiple-arm galaxies are found only closer than z ∼ 0.6 in the Ultra Deep Field. Grand design galaxies extend furthest to z ∼ 1.8, presumably because interactions can drive a two-arm spiral in a disk that would otherwise have a more irregular structure. The difference between these types is understandable in terms of the usual stability parameters for gas and stars, and the ratio of the velocity dispersion to rotation speed.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2016

Development of Semi-Empirical Damping Equation for Baffled Tank with Oblate Spheroidal Dome

Propellant slosh is a potential source of disturbance that can significantly impact the stability of space vehicles. The slosh dynamics are typically represented by a mechanical model of a spring-mass-damper. This mechanical model is then included in the equation of motion of the entire vehicle for Guidance, Navigation and Control analysis. The typical parameters required by the mechanical model include natural frequency of the slosh, slosh mass, slosh mass center location, and the critical damping ratio. A fundamental study has been undertaken at NASA MSFC to understand the fluid damping physics from a ring baffle in the barrel section of a propellant tank. An asymptotic damping equation and CFD blended equation have been derived by NASA MSFC team to complement the popularly used Miles equation at different flow regimes. The new development has found success in providing a nonlinear damping model for the Space Launch System. The purpose of this study is to further extend the semi-empirical damping equations into the oblate spheroidal dome section of the propellant tanks. First, previous experimental data from the spherical baffled tank are collected and analyzed. Several methods of taking the dome curvature effect, including a generalized Miles equation, area projection method, and equalized fill height method, are assessed. CFD simulation is used to shed light on the interaction of vorticity around the baffle with the locally curved wall and liquid-gas interface. The final damping equation will be validated by a recent subscale test with an oblate spheroidal dome conducted at NASA MSFC.

Yang, H. Q.↗
NASA NTRS1999

Hydrodynamic Stability and Magnetic Reconnection in Disks and Stars

The purpose of this grant is to study parametric instability. The simplest example of parametric instability is a harmonic oscillator with a periodic modulation of the spring constant. If the modulation frequency is close to twice the natural frequency of the oscillator, the amplitude of oscillation tends to grow exponentially. The growth rate is proportional to the strength of the modulation, but it also depends upon the closeness to resonance of the two frequencies, and upon natural damping rate or "Q" of the oscillator. Parametric instabilities are very common in physics. A familiar example is a jogger's ponytail--normally a very strongly damped pendulum, it can be destabilized by the variation in effective gravity during the jogger's stride. Observation confirms that the period of the pendulum is half that of the jogger's vertical motion. In astrophysics, parametric instability may occur by external tidal forcing, or by interaction among eigenmodes. In the latter case, an energetic eigenmode may destabilize modes of half its frequency, provided some weak nonlinearity exists to couple them. Under a previous Astrophysical Theory grant (NAGW-2419), the PI discovered a parametric instability of tidally forced disks such as the accretion disks in cataclysmic variables and X ray binaries [2]. The destabilized modes are tightly-wound, incompressible, three-dimensional waves analogous to g-modes and r-modes in stars. Later work has confirmed our analysis [4]. It was hoped that these modes might provide a source of turbulence and angular momentum transport in accretion disks. However, a follow-up investigation of this instability by local numerical simulations, although confirming the analytically estimated growth rates, found negligible angular momentum flux [3]. Other work, partly supported by the ATP, now strongly indicates that the transport mechanism in such disks is magnetohydrodynamic turbulence [6]. Nevertheless, the parametric mechanism may truncate the outer edges of disks in close binaries [2], and it may be important in disks of very low ionization such as protostellar disks, or even cataclysmic-variable disks in quiescence where the MHD mechanism may be ineffective [5]. All analyses up to 1996 were done in a local approximation where the orbital frequency, shear rate, and tidal field were treated as constants. The locally computed growth rate turns out to depend strongly on radius, and it was unclear how to average these local rates to obtain the correct global rate. This is a critical issue for accretion disks in close binaries, because the local growth rate is comparable to the orbital frequency towards the outer edge of the disk but decreases rapidly inwards. Paper #1 examined this issue in a simplified global model where the destabilizing terms vary with position. We found that the global growth rate is essentially equal to the maximum local rate, provided that the latter is smoothed over a radial range equal to the distance that the destabilized wave propagates at its group speed in one growth time. Thus, in an accretion disk, waves would grow rapidly in the outer parts but would propagate both inwards and outwards at a maximum group speed of order the disk thickness divided by the orbital period.

Goodman, Jeremy↗
NASA NTRS2012

Aeroelastic Modeling of Elastically Shaped Aircraft Concept via Wing Shaping Control for Drag Reduction

Lightweight aircraft design has received a considerable attention in recent years as a means for improving cruise efficiency. Reducing aircraft weight results in lower lift requirements which directly translate into lower drag, hence reduced engine thrust requirements during cruise. The use of lightweight materials such as advanced composite materials has been adopted by airframe manufacturers in current and future aircraft. Modern lightweight materials can provide less structural rigidity while maintaining load-carrying capacity. As structural flexibility increases, aeroelastic interactions with aerodynamic forces and moments become an increasingly important consideration in aircraft design and aerodynamic performance. Furthermore, aeroelastic interactions with flight dynamics can result in issues with vehicle stability and control. Abstract This paper describes a recent aeroelastic modeling effort for an elastically shaped aircraft concept (ESAC). The aircraft model is based on the rigid-body generic transport model (GTM) originally developed at NASA Langley Research Center. The ESAC distinguishes itself from the GTM in that it is equipped with highly flexible wing structures as a weight reduction design feature. More significantly, the wings are outfitted with a novel control effector concept called variable camber continuous trailing edge (VCCTE) flap system for active control of wing aeroelastic deflections to optimize the local angle of attack of wing sections for improved aerodynamic efficiency through cruise drag reduction and lift enhancement during take-off and landing. The VCCTE flap is a multi-functional and aerodynamically efficient device capable of achieving high lift-to-drag ratios. The flap system is comprised of three chordwise segments that form the variable camber feature of the flap and multiple spanwise segments that form a piecewise continuous trailing edge. By configuring the flap camber and trailing edge shape, drag reduction could be achieved. Moreover, some parts of the flap system can be made to have a high frequency response for roll control, gust load alleviation, and aeroservoelastic (ASE) modal suppression control. Abstract The aeroelastic model of the ESAC is based on one-dimensional structural dynamic theory that captures the aeroelastic deformation of a wing structure in a combined motion that involves flapwise bending, chordwise bending, and torsion. The model includes the effect of aircraft propulsion due to wing flexibility which causes the propulsive forces and moments to couple with the wing elastic motion. Engine mass is also accounted in the model. A fuel management model is developed to describe the wing mass change due to fuel usage in the main tank and wing tanks during cruise. Abstract The model computes both static and dynamic responses of the wing structures. The static aeroelastic deflections are used to estimate the effect of wing flexibility on induced drag and the potential drag reduction by the VCCTE flap system. A flutter analysis is conducted to estimate the flutter speed boundary. Gust load alleviation via adaptive control has been recently investigated to address flexibility of aircraft structures. A multi-objective flight control approach is presented for drag reduction control. The approach is based on an optimal control framework using a multi-objective cost function. Future studies will demonstrate the potential benefits of the approach.

Aeroelasticity↗