Search NASA⌕ Search

SEARCH · Search NASA

Results for “force matrix”

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 487 records · Page 27

SPARROW: A Steam Propelled Autonomous Retrieval Robot for Ocean Worlds

This paper presents the results of a NIAC Phase I study into the use of a propulsively hopping robot for the exploration of Europa’s rugged, icy surface. Named the“Steam Propelled Retrieval Robot for Ocean Worlds,” SPARROW is a multi-thruster robot passively gimballed within a protective, spherical shell, which enables it to freely rotate, self-right, and tumble over chaotic terrains. SPARROW is envisioned as a soccerball-sized payload to a primary lander mission. Europa’s abundant surface ice would be harvested as an in situ propellant source. The principal objective of SPARROW would be to increase the science return of a Europa landed asset by enabling access to distal, spatially distributed geologic units. The design of mobility systems for Europa is challenging, due in part to its almost entirely unconstrained surface topography and strength. Images returned by Voyager and Galileo yielded resolutions on the order of hundreds of meters per pixel, with localized regions reaching 6 meters per pixel—still far larger than a typical rover. A key benefit of SPARROW’s hopping, impact-tolerant design, is that it eliminates the need for a piori information on the terrain topography and surface strength; no surface reaction forces are required for motion. In this context, SPARROW is entirely terrain agnostic. In this paper we detail the results of three study objectives: i) to quantify the energy required to collect surface ice, change its phase, and maintain propellant temperature, ii) to identify control and estimation strategies that enable SPARROW to successfully reach, and return from, regions of scientific interest, and iii) to characterize the impact of SPARROW’s range on likely science return. Five water-based propellant architectures are presented alongside their mass, power, and volume requirements. Monte Carlo simulations of SPARROW hopping and tumbling over 1 km of glacial ice are summarized, characterizing SPARROW’s sensitivity to uncertainty in: initial conditions, thrust control, and cage-terrain interaction. Finally, a science traceability matrix is presented, which details the effect of sortie range on three science goals: constraining Europa’s evolutionary morphology, assessing sub-surface ocean habitability, and searching for life and/or biosignatures.

Phillips, Cynthia↗

SPARROW: Steam Propelled Autonomous Retrieval Robot for Ocean Worlds

The Steam Propelled Autonomous Retrieval Robot (SPARROW) for Ocean Worlds was a Phase I mission concept study funded under the NASA NIAC program. This report represents the findings of that study and recommendations for future work. SPARROW, envisioned as a soccer ball-sized payload to a primary lander mission, is a propulsively hopping robot for the exploration of Europa's rugged, icy surface. A multi-thruster, passively gimballed robot within a protective, spherical shell, SPARROW is able to freely rotate, self-right, and tumble over chaotic terrains. Europa's abundant surface ice would be harvested as an in situ propellant source. The principal objective of SPARROW is to increase the science return of a Europa landed asset by enabling access to distal, spatially distributed geologic units. The design of mobility systems for Europa is challenging, due in part to its almost entirely unconstrained surface topography and strength. Images returned by Voyager and Galileo yielded resolutions on the order of hundreds of meters per pixel, with localized regions reaching 6 meters per pixel—still far larger than a typical rover. A key benefit of SPARROW's hopping, impact-tolerant design, is that it eliminates the need for a priori information regarding terrain topography and surface strength; no surface reaction forces are required for motion. In this context, SPARROW is believed to be entirely terrain agnostic. In this report we detail the results of three study objectives: i) to quantify the energy required to collect surface ice, change its phase, and maintain propellant temperature, ii) to identify control and estimation strategies that enable SPARROW to successfully reach, and return from, regions of scientific interest, and iii) to characterize the impact of SPARROW's range on likely science return. Five water-based propellant architectures are presented alongside their mass, power, and volume requirements. Monte Carlo simulations of SPARROW hopping and tumbling over 1 km of glacial ice are summarized, characterizing SPARROW's sensitivity to uncertainty in: initial pose, thrust profile, and vehicle-terrain interaction. A science traceability matrix is presented, which details the effect of sortie range on three science goals: constraining Europa's evolutionary morphology, assessing sub-surface ocean habitability, and searching for life and/or biosignatures.

Autonomous↗

Evaluation of control laws and actuator locations for control systems applicable to deformable astronomical telescope mirrors

Some of the major difficulties associated with large orbiting astronomical telescopes are the cost of manufacturing the primary mirror to precise tolerances and the maintaining of diffraction-limited tolerances while in orbit. One successfully demonstrated approach for minimizing these problem areas is the technique of actively deforming the primary mirror by applying discrete forces to the rear of the mirror. A modal control technique, as applied to active optics, has previously been developed and analyzed. The modal control technique represents the plant to be controlled in terms of its eigenvalues and eigenfunctions which are estimated via numerical approximation techniques. The report includes an extension of previous work using the modal control technique and also describes an optimal feedback controller. The equations for both control laws are developed in state-space differential form and include such considerations as stability, controllability, and observability. These equations are general and allow the incorporation of various mode-analyzer designs; two design approaches are presented. The report also includes a technique for placing actuator and sensor locations at points on the mirror based upon the flexibility matrix of the uncontrolled or unobserved modes of the structure. The locations selected by this technique are used in the computer runs which are described. The results are based upon three different initial error distributions, two mode-analyzer designs, and both the modal and optimal control laws.

Ostroff, A. J.↗

Exploring the Connection Between Sampling Problems in Bayesian Inference and Statistical Mechanics

The Bayesian and statistical mechanical communities often share the same objective in their work - estimating and integrating probability distribution functions (pdfs) describing stochastic systems, models or processes. Frequently, these pdfs are complex functions of random variables exhibiting multiple, well separated local minima. Conventional strategies for sampling such pdfs are inefficient, sometimes leading to an apparent non-ergodic behavior. Several recently developed techniques for handling this problem have been successfully applied in statistical mechanics. In the multicanonical and Wang-Landau Monte Carlo (MC) methods, the correct pdfs are recovered from uniform sampling of the parameter space by iteratively establishing proper weighting factors connecting these distributions. Trivial generalizations allow for sampling from any chosen pdf. The closely related transition matrix method relies on estimating transition probabilities between different states. All these methods proved to generate estimates of pdfs with high statistical accuracy. In another MC technique, parallel tempering, several random walks, each corresponding to a different value of a parameter (e.g. "temperature"), are generated and occasionally exchanged using the Metropolis criterion. This method can be considered as a statistically correct version of simulated annealing. An alternative approach is to represent the set of independent variables as a Hamiltonian system. Considerab!e progress has been made in understanding how to ensure that the system obeys the equipartition theorem or, equivalently, that coupling between the variables is correctly described. Then a host of techniques developed for dynamical systems can be used. Among them, probably the most powerful is the Adaptive Biasing Force method, in which thermodynamic integration and biased sampling are combined to yield very efficient estimates of pdfs. The third class of methods deals with transitions between states described by rate constants. These problems are isomorphic with chemical kinetics problems. Recently, several efficient techniques for this purpose have been developed based on the approach originally proposed by Gillespie. Although the utility of the techniques mentioned above for Bayesian problems has not been determined, further research along these lines is warranted

Pohorille, Andrew↗

Precipitation Strengthenable NiTiPd High Temperature Shape Memory Alloys

In binary NiTi alloys, it has long been known that Ni-rich alloys can be heat treated to produce precipitates which both strengthen the matrix against dislocations and improve the behavior of the material under thermal and mechanical cycling. Within recent years, the same effect has been observed in Ni-rich NiTiHf high temperature shape memory alloys and heat treatment regimens have been defined which will reliably produce improved properties. In NiTiPd alloys, precipitation has also been observed, but studies are still underway to define reliable heat treatments and compositions which will provide a balance of strengthening and good thermomechanical properties. For this study, a series of NiTi-32 at.Pd alloys was produced to determine the effect of changing nickeltitanium content on the transformation behavior and heat treatability of the material. Samples were aged at temperatures between 350C and 450C for times up to 100 hours. Actuation type behavior was evaluated using uniaxial constant force thermal cycling (UCFTC) to determine the effect of composition and aging on the material behavior. TEMSEM was used to evaluate the microstructure and determine the types of precipitates formed. The correlation between composition, heat treat, microstructure, and thermomechanical behavior will be addressed and discussed.

Transformation strain↗

Thermal Processing Effects on the Adhesive Strength of PS304 High Temperature Solid Lubricant Coatings

In this paper the effects of post deposition heat treatments on the cohesive and adhesive strength properties of PS304, a plasma sprayed nickel-chrome based, high temperature solid lubricant coating deposited on stainless steel, are studied. Plasma spray deposited coating samples were exposed in air at temperatures from 432 to 650 C for up to 500 hr to promote residual stress relief, enhance particle to particle bonding and increase coating to substrate bond strength. Coating pull-off strength was measured using a commercial adhesion tester that utilizes 13 mm diameter aluminum pull studs attached to the coating surface with epoxy. Pull off force was automatically recorded and converted to coating pull off strength. As deposited coating samples were also tested as a baseline. The as-deposited (untreated) samples either delaminated at the coating-substrate interface or failed internally (cohesive failure) at about 17 MPa. Samples heat treated at temperatures above 540 C for 100 hr or at 600 C or above for more than 24 hr exhibited strengths above 31 MPa, nearly a two fold increase. Coating failure occurred inside the body of the coating (cohesive failure) for nearly all of the heat-treated samples and only occasionally at the coating substrate interface (adhesive failure). Metallographic analyses of heat-treated coatings indicate that the Nickel-Chromium binder in the PS304 appears to have segregated into two phases, a high nickel matrix phase and a high chromium precipitated phase. Analysis of the precipitates indicates the presence of silicon, a constituent of a flow enhancing additive in the commercial NiCr powder. The exact nature and structure of the precipitate phase is not known. This microstructural change is believed to be partially responsible for the coating strength increase. Diffusion bonding between particles may also be playing a role. Increasing the heat treatment temperature, exposure time or both accelerate the heat treatment process. Preliminary measurements indicate that the heat treatment also results in a one time, permanent coating thickness increase of about 3%. Based upon these results, the incorporation of a heat treatment prior to final finishing has been incorporated in the application process of this coating technology.

DellaCorte, Christopher↗

Receptor-based differences in human aortic smooth muscle cell membrane stiffness

Cells respond to mechanical stimuli with diverse molecular responses. The nature of the sensory mechanism involved in mechanotransduction is not known, but integrins may play an important role. The integrins are linked to both the cytoskeleton and extracellular matrix, suggesting that probing cells via integrins should yield different mechanical properties than probing cells via non-cytoskeleton-associated receptors. To test the hypothesis that the mechanical properties of a cell are dependent on the receptor on which the stress is applied, human aortic smooth muscle cells were plated, and magnetic beads, targeted either to the integrins via fibronectin or to the transferrin receptor by use of an IgG antibody, were attached to the cell surface. The resistance of the cell to deformation ("stiffness") was estimated by oscillating the magnetic beads at 1 Hz by use of single-pole magnetic tweezers at 2 different magnitudes. The ratio of bead displacements at different magnitudes was used to explore the mechanical properties of the cells. Cells stressed via the integrins required approximately 10-fold more force to obtain the same bead displacements as the cells stressed via the transferrin receptors. Cells stressed via integrins showed stiffening behavior as the force was increased, whereas this stiffening was significantly less for cells stressed via the transferrin receptor (P<0.001). Mechanical characteristics of vascular smooth muscle cells depend on the receptor by which the stress is applied, with integrin-based linkages demonstrating cell-stiffening behavior.

Non-NASA Center↗

Molecular Dynamics Modeling of Carbon Nanotube Composite Fracture Using ReaxFF

Carbon nanotube (CNT) fiber reinforced composites with specific tensile strengths and moduli approaching those of aerospace grade carbon fiber composites have recently been reported. This achievement was enabled by the emerging availability of high N/tex yarns in kilometer-scale quantities. While the production of this yarn is an impressive advance, its strength is still much lower than that of the individual CNTs comprising the yarn. Closing this gap requires understanding load transfer between CNTs at the nanometer dimensional scale. This work uses reactive molecular dynamics simulations to gain an understanding at the nanometer scale of the key factors that determine CNT nanocomposite mechanical performance, and to place more realistic upper bounds on the target properties. While molecular dynamics simulations using conventional force fields can predict elastic properties, the ReaxFF reactive forcefield can also model fracture behavior because of its ability to accurately describe bond breaking and formation during a simulation. The upper and lower bounds of CNT composite properties are investigated by comparing systems composed of CNTs continuously connected across the periodic boundary with systems composed of finite length CNTs. These lengths, effectively infinite for the continuous tubes and an aspect ratio of 13 for the finite length case, result from practical limitations on the number of atoms that can be included in a simulation. Experimentally measured aspect ratios are typically on the order of 100,000, so the calculated results should represent upper and lower limits on experimental mechanical properties. Finally, the effect of various degrees of covalent crosslinking between the CNTs and amorphous carbon matrix is considered to identify the amount of CNT-matrix covalent bonding that maximizes overall composite properties.

Jensen, Benjamin D.↗

The Gravity Field, Orientation, and Ephemeris of Mercury from MESSENGER Observations After Three Years in Orbit

We have analyzed three years of radio tracking data from the MESSENGER spacecraft in orbit around Mercury and determined the gravity field, planetary orientation, and ephemeris of the innermost planet. With improvements in spatial coverage, force modeling, and data weighting, we refined an earlier global gravity field both in quality and resolution, and we present here a spherical harmonic solution to degree and order 50. In this field, termed HgM005, uncertainties in low-degree coefficients are reduced by an order of magnitude relative to the earlier global field, and we obtained a preliminary value of the tidal Love number k(sub 2) of 0.451+/-0.014. We also estimated Mercury's pole position, and we obtained an obliquity value of 2.06 +/- 0.16 arcmin, in good agreement with analysis of Earth-based radar observations. From our updated rotation period (58.646146 +/- 0.000011 days) and Mercury ephemeris, we verified experimentally the planet's 3: 2 spin-orbit resonance to greater accuracy than previously possible. We present a detailed analysis of the HgM005 covariance matrix, and we describe some near-circular frozen orbits around Mercury that could be advantageous for future exploration.

MESSENGER↗

Optimal control of systems with capacity: Related noises

In the ordinary theory of optimal control (LQR and Kalman filter), the variances of the actuators and the sensors are assumed to be known (not related to the capacities of the devices). This assumption is not true in practice. Generally, a device with greater capacity to exert actuating forces and a sensor capable of sensing greater sensing range will generate noise of greater power spectral density. When the ordinary theory of optimal control is used to estimate the errors of the outputs in such cases it will lead to faulty results, because the capacities of such devices are unknown before the system is designed. The performance of the system designed by the ordinary theory will not be optimal as the variances of the sensors and the actuators are neither known nor constant. The interaction between the control system and structure could be serious because the ordinary method will lead to greater feedback (Kalman gain) matrices. Methods which can optimize the performance of systems when noises of the actuators and the sensors are related to their capacities are developed. These methods will result in smaller feedback (Kalman gain) matrix.

Ruan, Milfang↗

Design and Development of a Model to Simulate 0-G Treadmill Running Using the European Space Agency's Subject Loading System

Develop a model that simulates a human running in 0 G using the European Space Agency s (ESA) Subject Loading System (SLS). The model provides ground reaction forces (GRF) based on speed and pull-down forces (PDF). DESIGN The theoretical basis for the Running Model was based on a simple spring-mass model. The dynamic properties of the spring-mass model express theoretical vertical GRF (GRFv) and shear GRF in the posterior-anterior direction (GRFsh) during running gait. ADAMs VIEW software was used to build the model, which has a pelvis, thigh segment, shank segment, and a spring foot (see Figure 1).the model s movement simulates the joint kinematics of a human running at Earth gravity with the aim of generating GRF data. DEVELOPMENT & VERIFICATION ESA provided parabolic flight data of subjects running while using the SLS, for further characterization of the model s GRF. Peak GRF data were fit to a linear regression line dependent on PDF and speed. Interpolation and extrapolation of the regression equation provided a theoretical data matrix, which is used to drive the model s motion equations. Verification of the model was conducted by running the model at 4 different speeds, with each speed accounting for 3 different PDF. The model s GRF data fell within a 1-standard-deviation boundary derived from the empirical ESA data. CONCLUSION The Running Model aids in conducting various simulations (potential scenarios include a fatigued runner or a powerful runner generating high loads at a fast cadence) to determine limitations for the T2 vibration isolation system (VIS) aboard the International Space Station. This model can predict how running with the ESA SLS affects the T2 VIS and may be used for other exercise analyses in the future.

Caldwell, E. C.↗

Systematic Relationships Between Lidar Observables and Sizes And Mineral Composition Of Dust Aerosols

The physical and chemical properties of soil dust aerosol particles fundamentally affect their interaction with climate, including shortwave absorption and radiative forcing, nucleation of cloud droplets and ice crystals, heterogeneous formation of sulfates and nitrates on the surface of dust particles, and atmospheric processing of iron into bioavailable forms that increase the productivity of marine phytoplankton. Lidar measurements, such as extinction-to-backscatter, color and depolarization ratios, are frequently used to distinguish between aerosol types with different physical and chemical properties. The chemical composition of aerosol particles determines their complex refractive index, hence affecting their backscattering properties. Here we present a study on how dust aerosol backscattering and depolarization properties at wavelengths of 355, 532 and 1064 nm are related to size and complex refractive index, which varies with the mineral composition of the dust. Dust aerosols are represented by collections of spheroids with a range of prolate and oblate aspect ratios and their optical properties are obtained using T-matrix calculations. We find simple, systematic relationships between lidar observables and the dust size and complex refractive index that may aid the use of space-based or airborne lidars for direct retrieval of dust properties or for the evaluation of chemical transport models using forward simulated lidar variables. In addition, we present first results on the spatial variation of forward-simulated lidar variables based on a dust model that accounts for the atmospheric cycle of eight different mineral types plus internal mixtures of seven mineral types with iron oxides, which was recently implemented in the NASA GISS Earth System ModelE2.

Dust↗

The Omics of Stem Cell Mediated Regeneration: A Pilot Single Cell RNA-Seq Study of Mechanotransduction

Mechanical forces are potent modulators of stem cell based tissue regenerative mechanisms, inducing cell fate decisions and tissue specific commitment. A unique platform for investigating mechanotransduction is spaceflight, where microgravity and altered fluid mechanics provide a loading-null experimental condition. Seminal investigations of regenerative capacity in a wholly regenerative species, the newt model, and in a variety of totipotent and adult stem cell populations have demonstrated the detrimental effects of unloading on maintenance of stem cell based regeneration. Of particular interest is the observation that unloading interferes with the transition of stem cell pools from proliferative state to differentiation commitment. In this work we sought to test the hypothesis that gravity mechanotransduction regulates stem cell tissue regenerative processes by modulating stem cell proliferation and differentiation fates at specific cell cycle stages. To do this, clonally-derived ESCs were plated on a collagen matrix and expanded for 36 hours before re-plating on a non-adherent culture dish in the absence of leukemia inhibitory factor (LIF) to form spheroid aggregate EBs. After formation, the EBs were transferred to a collagen matrix coated culture dishes and given 4 days to allow implantation and outgrowth. In parallel, totipotent ESCs were plated 24 hours before mechanical stimulation on collagen matrix culture dishes in the presence of LIF to maintain totipotency and serve as un-differentiation committed controls. The EBs and ESCs were then subjected to either a 60 minute pulse of gravity (static loading) or 60 minutes of cyclic stretch (dynamic loading) mechanotransduction. Six hours post-stimulation, we used a 10X Genomics Single Cell controller to generate bar-coded single cell Illumina libraries and sequenced expressomes for 5,000 static loaded cells, representative of a change in gravity mechanotransduction, 5,000 dynamic loaded cells, representative of tissue loading associate with physiologic function, and 5,000 unstimulated 1g control cells. The comparison of these 3 libraries by cluster assignment based on like gene expression patterns show substantial alteration in cluster geometry due to mechanical loading. Specifically the mechanically loaded EB outgrowth cells to retain potency markers (PAX6, SOX2, CD34) and suppress early commitment markers (Dhh, VCAN, Igf1). Whereas the EBs cultured under the non-stimulated conditions display clear departure from the ESC expressome with lineage commitment markers upregulated and several tissue specific markers being expressed (BMP "early musculoskeletal development, Mesp1" early cardiovascular cell lineage). These markers are not seen in the mechano-stimulated cultures or the totipotent ESC cultures. Comparison of like clusters between our experimental conditions revealed an array of regenerative and stem cell genes are significantly mechano-regulated. Of particular importance CDKN1a/p21, a gene shown by previous investigation of our research team to be significantly upregulated in unloading, was suppressed in the static and dynamic loaded EBS. In addition to CDKN1a/p21 many genes related to cell cycle and transitory differentiation markers had elevated expression in the mechano-stimulated EBs, but surprisingly these trends were not observed in the ESC cultures. This study is the first of its kind investigating for mechano-signaling and mechano-regulated pathways, and has alre

Single cell RNA-sequecing↗

Validation of three-dimensional incompressible spatial direct numerical simulation code: A comparison with linear stability and parabolic stability equation theories for boundary-layer transition on a flat plate

Spatially evolving instabilities in a boundary layer on a flat plate are computed by direct numerical simulation (DNS) of the incompressible Navier-Stokes equations. In a truncated physical domain, a nonstaggered mesh is used for the grid. A Chebyshev-collocation method is used normal to the wall; finite difference and compact difference methods are used in the streamwise direction; and a Fourier series is used in the spanwise direction. For time stepping, implicit Crank-Nicolson and explicit Runge-Kutta schemes are used to the time-splitting method. The influence-matrix technique is used to solve the pressure equation. At the outflow boundary, the buffer-domain technique is used to prevent convective wave reflection or upstream propagation of information from the boundary. Results of the DNS are compared with those from both linear stability theory (LST) and parabolized stability equation (PSE) theory. Computed disturbance amplitudes and phases are in very good agreement with those of LST (for small inflow disturbance amplitudes). A measure of the sensitivity of the inflow condition is demonstrated with both LST and PSE theory used to approximate inflows. Although the DNS numerics are very different than those of PSE theory, the results are in good agreement. A small discrepancy in the results that does occur is likely a result of the variation in PSE boundary condition treatment in the far field. Finally, a small-amplitude wave triad is forced at the inflow, and simulation results are compared with those of LST. Again, very good agreement is found between DNS and LST results for the 3-D simulations, the implication being that the disturbance amplitudes are sufficiently small that nonlinear interactions are negligible.

Joslin, Ronald D.↗

Continuation of tailored composite structures of ordered staple thermoplastic material

The search for the cost effective composite structure has motivated the investigation of several approaches to develop composite structure from innovative material forms. Among the promising approaches is the conversion of a planar sheet to components of complex curvature through sheet forming or stretch forming. In both cases, the potential for material stretch in the fiber direction appears to offer a clear advantage in formability over continuous fiber systems. A framework was established which allows the simulation of the anisotropic mechanisms of deformation of long discontinuous fiber laminates wherein the matrix phase is a viscous fluid. Predictions for the effective viscosities of a hyper-anisotropic medium consisting of collimated, discontinuous fibers suspended in viscous matrix were extended to capture the characteristics of typical polymers including non-Newtonian behavior and temperature dependence. In addition, the influence of fiber misorientation was also modeled by compliance averaging to determine ensemble properties for a given orientation distribution. A design tool is presented for predicting the effect of material heterogeneity on the performance of curved composite beams such as those used in aircraft fuselage structures. Material heterogeneity can be induced during manufacturing processes such as sheet forming and stretch forming of thermoplastic composites. This heterogeneity can be introduced in the form of fiber realignment and spreading during the manufacturing process causing radial and tangential gradients in material properties. Two analysis procedures are used to solve the beam problems. The first method uses separate two-dimensional elasticity solutions for the stresses in the flange and web sections of the beam. The separate solutions are coupled by requiring that forces and displacements match section boundaries. The second method uses an approximate Rayleigh-Ritz technique to find the solutions for more complex beams. Analyses are performed for curved beams of various cross-sections loaded in pure bending and with a uniform distributed load. Preliminary results show that the geometry of the beam dictates the effect of heterogeneity on performance. The role of heterogeneity is larger in beams with a small average radius-to-depth ration, R/t, where R is the average radius of the beam and t is the difference between the inside and outside radii. Results of the anlysis are in the form of stresses and displacements and are compared to both mechanics of materials and numerical solutions obtained using finite element analysis.

Santare, Michael H.↗

3D Navier-Stokes Time Accurate Solutions Using Multipartitioning Parallel Computation Methodology

A parallel CFD code solving 3D time accurate Navier-Stokes equations with multipartitioning parallel Methodology is being developed in collaboration with Ohio State University within the Air Vehicle Directorate, at Wright Patterson Air Force Base. The advantage of the multipartitioning parallel method is that the domain decomposition will not introduce domain boundaries for the implicit operators. A ring structure data communication is employed so that the implicit time accurate method can be implemented for multi-processors with the same accuracy as for the single processor. No sub-iteration is needed at the domain boundaries. The code has been validated for some typical unsteady flows, which include Coutte Flow, flow passing a cylinder. The code now is being employed for a large scale time accurate wall jet transient flow computation. 'ne preliminary results are promising. The mesh has been refined to capture more details of the flow field. The mesh refinement computation is in progress and would be difficult to successfully implement without the parallel computation techniques used. A modified version of the code with more efficient inversion of the diagonalized block matrix is currently being tested.

Zha, Ge-Cheng↗

Multi-Operator Multi-UAV (MOMU) Control: Exploring the Influence of Sensor Tools and Playbook Task Delegation

New concepts of operations for Unmanned Aerial Vehicles (UAVs) will require a change from the current 2:1 operator to vehicle crew configuration. One particular control paradigm, largely driven by logistics, manpower, and training burdens, as well as the desire to force multiply, involves a single operator simultaneously managing multiple UAVs. This mode of operations has shown to significantly increase cognitive workload and decrease situation awareness, as operators are required to simultaneously attend to multiple sources of information. One potential way to mitigate potential drawbacks of multi-vehicle control by a single operator is to migrate to a multi-operator multi-UAV (MOMU) crew configuration, whereby M operators control N (> M) vehicles. This type of crew configuration can be organized in several ways to dynamically manage cognitive workload, match operator qualifications and skills to mission requirements, increase utilization of available assets, and thereby achieve maximum force multiplication. The present experiment examined task performance in a simulated MOMU environment and evaluated the potential benefits of sensor management aids ("Tools") as well as integrated sensor and flight automation ("Plays") compared to a fully manual condition ("Manual"). Tools support the operator by facilitating rapid understanding and management of sensor information, while the Plays support the operator by offloading/automating subtasks. Six pairs of participants were recruited for this study and tasked with sharing a pool of UAVs in order to conduct reconnaissance, surveillance, and target acquisition (RSTA) missions in adjacent Areas of Operation (AOs). Participants were given four tasks to accomplish, in order of priority: 1) prosecute High Value Targets (HVTs); 2) identify / track targets (military vehicles); 3) identify / mark civilian vehicles; and 4) respond to chat messages. Performance on the mission tasks was measured in terms of accuracy and reaction time. A composite mission score was also calculated using a payoff matrix that weighted each task according to priority. The results indicate that Playbook demonstrated better performance overall with higher accuracy rates and the highest composite score compared to Tools and Manual. The implications of these results to supporting future MOMU concepts of operations is discussed.

Playbook↗

A General Interface Method for Aeroelastic Analysis of Aircraft

The aeroelastic analysis of an aircraft requires an accurate and efficient procedure to couple aerodynamics and structures. The procedure needs an interface method to bridge the gap between the aerodynamic and structural models in order to transform loads and displacements. Such an interface method is described in this report. This interface method transforms loads computed by any aerodynamic code to a structural finite element (FE) model and converts the displacements from the FE model to the aerodynamic model. The approach is based on FE technology in which virtual work is employed to transform the aerodynamic pressures into FE nodal forces. The displacements at the FE nodes are then converted back to aerodynamic grid points on the aircraft surface through the reciprocal theorem in structural engineering. The method allows both high and crude fidelities of both models and does not require an intermediate modeling. In addition, the method performs the conversion of loads and displacements directly between individual aerodynamic grid point and its corresponding structural finite element and, hence, is very efficient for large aircraft models. This report also describes the application of this aero-structure interface method to a simple wing and an MD-90 wing. The results show that the aeroelastic effect is very important. For the simple wing, both linear and nonlinear approaches are used. In the linear approach, the deformation of the structural model is considered small, and the loads from the deformed aerodynamic model are applied to the original geometry of the structure. In the nonlinear approach, the geometry of the structure and its stiffness matrix are updated in every iteration and the increments of loads from the previous iteration are applied to the new structural geometry in order to compute the displacement increments. Additional studies to apply the aero-structure interaction procedure to more complicated geometry will be conducted in the second phase of the present contract.

Tzong, T.↗