Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic data collection”

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 1,495 records · Page 83

Models of Small-Scale Patchiness

Patchiness is perhaps the most salient characteristic of plankton populations in the ocean. The scale of this heterogeneity spans many orders of magnitude in its spatial extent, ranging from planetary down to microscale. It has been argued that patchiness plays a fundamental role in the functioning of marine ecosystems, insofar as the mean conditions may not reflect the environment to which organisms are adapted. Understanding the nature of this patchiness is thus one of the major challenges of oceanographic ecology. The patchiness problem is fundamentally one of physical-biological-chemical interactions. This interconnection arises from three basic sources: (1) ocean currents continually redistribute dissolved and suspended constituents by advection; (2) space-time fluctuations in the flows themselves impact biological and chemical processes, and (3) organisms are capable of directed motion through the water. This tripartite linkage poses a difficult challenge to understanding oceanic ecosystems: differentiation between the three sources of variability requires accurate assessment of property distributions in space and time, in addition to detailed knowledge of organismal repertoires and the processes by which ambient conditions control the rates of biological and chemical reactions. Various methods of observing the ocean tend to lie parallel to the axes of the space/time domain in which these physical-biological-chemical interactions take place. Given that a purely observational approach to the patchiness problem is not tractable with finite resources, the coupling of models with observations offers an alternative which provides a context for synthesis of sparse data with articulations of fundamental principles assumed to govern functionality of the system. In a sense, models can be used to fill the gaps in the space/time domain, yielding a framework for exploring the controls on spatially and temporally intermittent processes. The following discussion highlights only a few of the multitude of models which have yielded insight into the dynamics of plankton patchiness. In addition, this particular collection of examples is intended to furnish some exposure to the diversity of modeling approaches which can be brought to bear on the problem. These approaches range from abstract theoretical models intended to elucidate specific processes, to complex numerical formulations which can be used to actually simulate observed distributions in detail.

McGillicuddy, D. J.↗

Assessment of the Reconstructed Aerodynamics of the Mars Science Laboratory Entry Vehicle

On August 5, 2012, the Mars Science Laboratory entry vehicle successfully entered Mars atmosphere, flying a guided entry until parachute deploy. The Curiosity rover landed safely in Gale crater upon completion of the Entry Descent and Landing sequence. This paper compares the aerodynamics of the entry capsule extracted from onboard flight data, including Inertial Measurement Unit (IMU) accelerometer and rate gyro information, and heatshield surface pressure measurements. From the onboard data, static force and moment data has been extracted. This data is compared to preflight predictions. The information collected by MSL represents the most complete set of information collected during Mars entry to date. It allows the separation of aerodynamic performance from atmospheric conditions. The comparisons show the MSL aerodynamic characteristics have been identified and resolved to an accuracy better than the aerodynamic database uncertainties used in preflight simulations. A number of small anomalies have been identified and are discussed. This data will help revise aerodynamic databases for future missions and will guide computational fluid dynamics (CFD) development to improved prediction codes.

Schoenenberger, Mark↗

Design Concepts for an Outage Control Center Information Dashboard

The nuclear industry, and the business world in general, is facing a rapidly increasing amount of data to be dealt with on a daily basis. In the last two decades, the steady improvement of data storage devices and means to create and collect data along the way influenced the manner in which we deal with information. Most data is still stored without filtering and refinement for later use. Many functions at a nuclear power plant generate vast amounts of data, with scheduled and unscheduled outages being a prime example of a source of some of the most complex data sets at the plant. To make matters worse, modern information and communications technology is making it possible to collect and store data faster than our ability to use it for making decisions. However, in most applications, especially outages, raw data has no value in itself; instead, managers, engineers and other specialists want to extract the information contained in it. The complexity and sheer volume of data could lead to information overload, resulting in getting lost in data that may be irrelevant to the task at hand, processed in an inappropriate way, or presented in an ineffective way. To prevent information overload, many data sources are ignored so production opportunities are lost because utilities lack the ability to deal with the enormous data volumes properly. Decision-makers are often confronted with large amounts of disparate, conflicting and dynamic information, which are available from multiple heterogeneous sources. Information and communication technologies alone will not solve this problem. Utilities need effective methods to exploit and use the hidden opportunities and knowledge residing in unexplored data resources. Superior performance before, during and after outages depends upon the right information being available at the right time to the right people. Acquisition of raw data is the easy part; instead, it is the ability to use advanced analytical, data processing and data visualization methods to turn the data into reliable information and comprehensible, actionable information. Techniques like data mining, filtering and analysis only work reliably for well-defined and well-understood problems. The path from data to decision is more complex. The ability to communicate knowledge during outages and emergent issues is crucial. This paper presents an approach to turn the unused data into an opportunity: applying principles from semiotics, human factors and visual analytics to transform the traditional way of processing outage data into media that will improve the collective situation awareness, knowledge, decisions, actions and overall performance of the entire outage team, and also support the reliability, quality and overall effectiveness of maintenance work. The application of the proposed visualization methods will become the medium of a semi-automated analytical process where humans and machines cooperate using their respective, distinct capabilities for the most effective results.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The Global Water Monitor: Lake, Wetland, and River Reach Monitoring for Resource Management and Hazard Observation

The Global Water Monitor is a monitoring program offering surface water-related products for lakes, reservoirs, river reaches, and wetlands, https://blueice.gsfc.nasa.gov/gwm. These products are derived from multiple satellite-based altimetry and multispectral imaging platforms. The primary measurements are water level, slope, and extent, derived from a series of NASA and other agency instruments (Sentinel-3,-6, ICESat-2, SWOT, MODIS). The system serves stakeholders, the US Department of Agriculture Foreign Agricultural Service and the US Geological Survey, as well as various international, national, state, and intelligence agencies, engineers, ecologists, and hydrological researchers. Considering the sensitivity of ground-based data and the remote accessibility of many basins, the satellites offer global coverage, accuracy, and continuity of measurements. Observation of storage fluctuations in lakes and reservoirs provides users with knowledge of short and long-term drought conditions that can affect both water and energy resources. Emphasis is on changing climatic conditions, water sharing between nations, and regional stability. Measurements across a wetland complex can assist with efforts to assess hydrological dynamics with a focus on conservation efforts that aid both natural ecology and the sharing of precious water resources – the latter potentially being required in a basin across the municipal, crop irrigation, aquaculture, livestock, and power station spectrum. Observations of river reaches assists in the collection of basic hydraulic information, particularly in gauge-poor regions, for modelling efforts and the estimation of discharge, and at high-latitudes satellite-based river observations serve projects that warn of spring melts and flood hazards. We look at how recent improvements in instrument technology and mission operations are assisting with acquiring true global coverage and cross-validation efforts, and how stakeholders interact with the project in terms of setting requirements, product format, and accessibility methods. Examples show the application of the products across a range of programs, and in addition to meeting operational requirements, how the multi-decadal surface water levels are creating a high-quality timeline of Earth Data Records.

Lakes↗

Parked aeroelastic field rotor response for a 20% scaled demonstrator of a 13‐MW downwind turbine

Abstract Aeroelastic parked testing of a unique downwind two‐bladed subscale rotor was completed to characterize the response of an extreme‐scale 13‐MW turbine in high‐wind parked conditions. A 20% geometric scaling was used resulting in scaled 20‐m‐long blades, whose structural and stiffness properties were designed using aeroelastic scaling to replicate the nondimensional structural aeroelastic deflections and dynamics that would occur for a lightweight, downwind 13‐MW rotor. The subscale rotor was mounted and field tested on the two‐bladed Controls Advanced Research Turbine (CART2) at the National Renewable Energy Laboratory's Flatiron Campus (NREL FC). The parked testing of these highly flexible blades included both pitch‐to‐run and pitch‐to‐feather configurations with the blades in the horizontal braked orientation. The collected experimental data includes the unsteady flapwise root bending moments and tip deflections as a function of inflow wind conditions. The bending moments are based on strain gauges located in the root section, whereas the tip deflections are captured by a video camera on the hub of the turbine pointed toward the tip of the blade. The experimental results are compared against computational predictions generated by FAST, a wind turbine simulation software, for the subscale and full‐scale models with consistent unsteady wind fields. FAST reasonably predicted the bending moments and deflections of the experimental data in terms of both the mean and standard deviations. These results demonstrate the efficacy of the first such aeroelastically scaled turbine test and demonstrate that a highly flexible lightweight downwind coned rotor can be designed to withstand extreme loads in parked conditions.

17 WIND ENERGY↗

Amplitude analysis of ψ3686→γKS0KS0

Using (2712 ± 14) × 106ψ(3686) events collected with the BESIII detector, we perform the first amplitude analysis of the radiative decay ψ3686→γKS0KS0$$ \psi (3686)\to \gamma {K}_S^0{K}_S^0 $$ within the mass region MKS0KS0<2.8$$ {M}_{K_S^0{K}_S^0}<2.8 $$ GeV/c2. Employing a one-channel K-matrix approach for the description of the dynamics of the KS0KS0$$ {K}_S^0{K}_S^0 $$ system, the data sample is well described with four poles for the f0-wave and three poles for the f2-wave. The determined pole positions are consistent with those of well-established resonance states. The observed f0 and f2 states are found to be in agreement with those produced in radiative J/ψ decays. The production behaviors of f0 and f2 poles in ψ3686→γKS0KS0$$ \psi (3686)\to \gamma {K}_S^0{K}_S^0 $$ are qualified with their residues and the converted branching fractions. By comparing with J/ψ→γKS0KS0$$ J/\psi \to \gamma {K}_S^0{K}_S^0 $$ decay, the ratios Bψ3686→γf0,2BJ/ψ→γf0,2$$ \frac{\mathcal{B}\left(\psi (3686)\to \gamma {f}_{0,2}\right)}{\mathcal{B}\left(J/\psi \to \gamma {f}_{0,2}\right)} $$ are determined, which provides crucial experimental inputs on the internal structure of the f0,2 states, especially their potential mixing with glueball components.

Ablikim, M↗

Occupant comfort and behavior: High-resolution data from a 6-month field study of personal comfort systems with 37 real office workers

Personal Comfort Systems (PCS) provide individual occupants local heating and cooling to meet their comfort needs without affecting others in the same space. It saves energy by relaxing ambient temperature requirements for the HVAC system. Aside from these benefits, PCS offers a wealth of data that can describe how individuals interact with heating/cooling devices in their own environment. Recently developed Internet-connected PCS chairs have unlocked this opportunity by generating continuous streams of heating and cooling usage data, along with occupancy status and environmental measurements via embedded sensors. The data allow individuals' comfort and behavior to be learned, and can inform centralized systems to provide ‘just the right’ amount of conditioning to meet occupant needs. In summer 2016, we carried out a study with PCS chairs involving 37 occupants in an office building in California. During the study period, we collected >5 million chair usage data-points and 4500 occupant survey responses, as well as continuous measurements of environmental and HVAC system conditions. The data analysis shows that (1) local temperatures experienced by individual occupants vary quite widely across different parts of the building, even within the same thermal zone; (2) occupants often have different thermal preferences even under the same thermal conditions; (3) PCS control behavior can dynamically describe individuals' thermal preferences; (4) PCS chairs produce much higher comfort satisfaction (96%) than typically achieved in buildings. We conclude that PCS not only provide personalized comfort solutions but also offer individualized feedback that can improve comfort analytics and control decisions in buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Structural Evolution of Polylactide Molecular Bottlebrushes: Kinetics Study by Size Exclusion Chromatography, Small Angle Neutron Scattering, and Simulations

In this paper, structural evolution from poly(lactide) (PLA) macromonomer to resultant PLA molecular bottlebrush during ring opening metathesis polymerization (ROMP) was investigated for the first time by combining size exclusion chromatography (SEC), small-angle neutron scattering (SANS), and coarse-grained molecular dynamics (CG-MD) simulations. Multiple aliquots were collected at various reaction times during ROMP and subsequently analyzed by SEC and SANS. These complementary techniques enable the understanding of systematic changes in conversion, molecular weight and dispersity as well as structural details of PLA molecular bottlebrushes. CG-MD simulation not only predicts the experimental observations, but it also provides further insight into the analysis and interpretation of data obtained in SEC and SANS experiments. Finally, we find that PLA molecular bottlebrushes undergo three conformational transitions with increasing conversion (i.e., increasing the backbone length): (1) from an elongated to a globular shape due to longer side chain at low conversion, (2) from a globular to an elongated shape at intermediate conversion caused by excluded volume of PLA side chain, and (3) the saturation of contour length at high conversion due to chain transfer reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular quantum magnetism in LiZn2Mo3O8

Inelastic neutron scattering for temperatures below 30 K from a powder of LiZn2Mo3O8 demonstrates this triangular-lattice antiferromagnet hosts collective magnetic excitations from spin 1/2 Mo3O13 molecules. Apparently gapless ( \Delta < 0.2 meV) and extending at least up to 2.5 meV, the low energy magnetic scattering cross section is surprisingly broad in momentum space and involves one third of the spins present above 100 K. The data are compatible with the presence of valencebonds involving nearest-neighbor and next-nearest-neighbor spins forming a disordered or dynamic state.

Mourigal, Martin↗

FOSWEC v2: System identification and control testing at Hinsdale Wave Research Laboratory

This report summarizes the design and execution of a wave tank test of the floating oscillating surge wave energy converter (FOSWEC) in the O.H. Hinsdale Wave Research Laboratory Directional Wave Basin at Oregon State University. This device, which uses two "flaps" that pivot about a central platform when excited by waves, has a natural frequency within the range of the waves by which it is excited. The FOSWEC was originally considered to be a 1:33 scale device, however, for the current tests, no fixed relative scale is used (i.e., the WEC is considered to be scaled for the basin’s wave environment in which it operates). The primary goal of this test was to assess the degree to which previously developed modeling, experimentation, and control design methods could be applied to a broad range of wave energy converter designs. Testing was conducted to identify a dynamic model for the impedance and excitation behavior of the device. Using these models, a series of closed loop tests were conducted using a causal impedance matching controller. This report provides a brief description of the results, as well as a summary of the device and experimental design. The results show that the methods applied to this experimental device perform well and should be broadly applicable. The data collected during testing is compressed into FOSWEC.zip. Please refer to Appendix C (pages 61-63) of the test report for descriptions of each test ID corresponding to the compressed files.

16 TIDAL AND WAVE POWER↗

Global Carbon Budget 1800-2018 Loss of Additional Sink Capacity and Present Transient Difference

A dataset of two global carbon fluxes: the loss of additional sink capacity (LASC) and present versus transient difference (PTD). The loss of additional sink capacity is the loss of indirect anthropogenic sink capacity in land ecosystems caused by the loss of ecosystems that can increase carbon in response to changing environments, primarily forests. The present transient difference refers to the direct anthropogenic flux in land ecosystems (from land cover and land use change) and the difference in that flux when calculated assuming "present day" (actually more like 1980) ecosystem carbon stocks or when calculated assuming transient carbon stocks (i.e. caused by environmental change). In this dataset, the LASC and PTD fluxes were calculated using the TRENDY ensemble of Dynamic Global Vegetation Models (DGVMs; Sitch et al., 2024) following the methods described therein and in Obermeier et al. (2021). For further details see Obermeier et al. (2021), Friedlingstein et al. (2019), Sitch et al. (2024), and Walker et al. (2025). The dataset includes ensemble means and standard deviations for annual fluxes, fluxes cumulated annually over the whole time period, and cumulated since 1959. These data were originally collected to compare different approaches (DGVMs versus bookkeeping models) for estimating land use and land cover change (LULCC) emissions (Obermeier et al. 2021). More recently, and similarly, these data were used to harmonize the estimates of the direct and indirect anthropogenic fluxes in the global carbon cycle and to recalculate the global carbon budget (Walker et al., 2025).

54 ENVIRONMENTAL SCIENCES↗

Achieving Tier 4 Emissions in Biomass Cookstoves

Previous literature on top-lit updraft (TLUD) gasifier cookstoves suggested that these stoves have the potential to be the lowest emitting biomass cookstove. However, the previous literature also demonstrated a high degree of variability in TLUD emissions and performance, and a lack of general understanding of the TLUD combustion process. The objective of this study was to improve understanding of the combustion process in TLUD cookstoves. In a TLUD, biomass is gasified and the resulting producer gas is burned in a secondary flame located just above the fuel bed. The goal of this project is to enable the design of a more robust TLUD that consistently meets Tier 4 performance targets through a better understanding of the underlying combustion physics. The project featured a combined modeling, experimental and product design/development effort comprised of four different activities: Development of a model of the gasification process in the biomass fuel bed; Development of a CFD model of the secondary combustion zone; Experiments with a modular TLUD test bed to provide information on how stove design, fuel properties, and operating mode influence performance and provide data needed to validate the fuel bed model; Planar laser-induced fluorescence (PLIF) experiments with a two-dimensional optical test bed to provide insight into the flame dynamics in the secondary combustion zone and data to validate the CFD model; Design, development and field testing of a market ready TLUD prototype. Over 180 tests of 40 different configurations of the modular TLUD test bed were performed to demonstrate how stove design, fuel properties and operating mode influences performance, and the conditions under which Tier 4 emissions are obtainable. Images of OH and acetone PLIF were collected at 10 kHz with the optical test bed. The modeling and experimental results informed the design of a TLUD prototype that met Tier 3 to Tier 4 specifications in emissions and Tier 2 in efficiency. The final prototype was field tested in India.

09 BIOMASS FUELS↗

Quantifying Subsurface Biogeochemical Variability in a High Altitude Watershed During Winter Isolation

Shallow subsurface microbial and geochemical processes in watersheds are dynamic and responsive to seasonal and long-term environmental change. At high elevation or latitude, such changes may occur during winter months when normal sampling is impossible, impractical, unsafe, or may actually alter the processes being studied. Yet accurate modeling of biogeochemical parameters in these environments requires sampling that can capture events and persistent cold climate trends. Gaps in current models exist where data have not been collected during extended snow- or ice-covered periods. The goals of our exploratory study were to examine the biogeochemical processes occurring in a well-studied high altitude watershed during winter when the system is largely inaccessible. We tested the hypothesis that during snow cover in the East River (ER) watershed, episodic excursions of microbial community structure and biogeochemical processes and concentrations fluctuated from values extrapolated from pre- and post-snow time periods. Through our research we identified previously unresolved subsurface-surface transport phenomena and a heretofore under-appreciated methane cycle occurring during wintertime months, thus providing valuable insight into watershed processes. Our research showed that autonomous sampling techniques, first used in marine systems, is extremely versatile and can be adapted for use in harsh continental settings to examine temporal changes in freshwater chemistry and microbial community assemblage. Further, we found that by sampling during periods when systems are hard to access, will lend insight and better accountability for global budgets for methane, an important greenhouse gas. By strengthening our knowledge of wintertime biogeochemical relationships in high altitude watersheds our research leads to more complete descriptions of cryptic processes and ultimately will inform and improve the reactive transport models that simulate processes in these systems.

54 ENVIRONMENTAL SCIENCES↗

Stereo Camera Deployment in Support of TRACER (Field Campaign Report)

An improved understanding of the salient environmental controls on cloud formation, evolution, and eventual dissipation is critical to address ongoing challenges with cloud process and parameterization representations in global climate and Earth system models. One goal for the recent U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s Tracking Aerosol Convection Interactions Experiment (TRACER) campaign was to collect a comprehensive data set that enabled such convective cloud process studies and key demonstrations for those controls that influenced cloud life cycle (i.e., aerosols, thermodynamics) in the Houston, Texas region. To help accomplish this, ARM instrumentation during TRACER was tasked with tracking a large number of individual, isolated convective cells – and to follow the evolution of these cells at high spatiotemporal resolution to characterize changes in cloud dynamic and microphysical properties. Since Houston experiences a range of convective clouds, it was known that the standard ARM Mobile Facility (AMF) instruments may not be sufficient to completely document initiating, transient, or dissipating low- or shallow-cloud behaviors that were also expected during this campaign (in terms of sensitivity, resolution, and/or operational availability). As one partial solution, a supplemental stereo camera deployment (this sub-campaign) was requested to augment the ARM AMF instrumentation to better address shallow and shallow-to-deep transitional types of cloud process drivers during TRACER (ARM stereo cameras for clouds [STEREOCAM]; Romps and Öktem 2018). The primary scientific focus was the relationships between cloud properties and the ambient conditions, which points to several key TRACER science questions including: ‘What is the relationship between cloud size or updraft intensity to the environmental wind shear and/or humidity?’ Overall, the ARM TRACER stereo camera deployment demonstrated unique effectiveness in observing a wide range of critical shallow, congestus, and transitioning or time-evolving cloud characteristics. The data sets from these cameras include information on the clouds' horizontal dimensions, elevations, and depths, while also enabling potential products for cloud initiation and dissipation rates, and vertical velocities. Stereo cameras simultaneously inform on cloud life cycle stage and spatial properties such as cloud fractional coverage, which should provide complementary information for ARM users when combined with TRACER cloud radars, lidar, and/or other profiling sensors. Moreover, camera products offer large-eddy simulation (LES),-scale-appropriate cloud coverage, depth, and spatial variability estimates, while opening additional avenues to challenge difficult process questions on cloud updrafts/entrainment and their covariability with environmental controls such as wind shear and humidity.

47 OTHER INSTRUMENTATION↗

MHD Waves and Turbulence in the Sun and Interplanetary Medium

The existence, nature, and dynamics of magnetohydrodynamic waves and turbulence in the solar atmosphere and interplanetary medium are addressed. Remote sensing observations of global oscillations of the Sun and their possible interpretation in terms of waves or turbulence are considered. The region of solar wind acceleration as described by coronal imaging, resonance-line spectrometry, and radio techniques is discussed. Fluctuations and discontinuities in the interplanetary medium are considered using primary data collected by spaceborne magnetometers and plasma analyzers.

Source record↗

Global Design Optimization for Aerodynamics and Rocket Propulsion Components

Modern computational and experimental tools for aerodynamics and propulsion applications have matured to a stage where they can provide substantial insight into engineering processes involving fluid flows, and can be fruitfully utilized to help improve the design of practical devices. In particular, rapid and continuous development in aerospace engineering demands that new design concepts be regularly proposed to meet goals for increased performance, robustness and safety while concurrently decreasing cost. To date, the majority of the effort in design optimization of fluid dynamics has relied on gradient-based search algorithms. Global optimization methods can utilize the information collected from various sources and by different tools. These methods offer multi-criterion optimization, handle the existence of multiple design points and trade-offs via insight into the entire design space, can easily perform tasks in parallel, and are often effective in filtering the noise intrinsic to numerical and experimental data. However, a successful application of the global optimization method needs to address issues related to data requirements with an increase in the number of design variables, and methods for predicting the model performance. In this article, we review recent progress made in establishing suitable global optimization techniques employing neural network and polynomial-based response surface methodologies. Issues addressed include techniques for construction of the response surface, design of experiment techniques for supplying information in an economical manner, optimization procedures and multi-level techniques, and assessment of relative performance between polynomials and neural networks. Examples drawn from wing aerodynamics, turbulent diffuser flows, gas-gas injectors, and supersonic turbines are employed to help demonstrate the issues involved in an engineering design context. Both the usefulness of the existing knowledge to aid current design practices and the need for future research are identified.

Shyy, Wei↗

Airborne Tomographic Swath Ice Sounding Processing System

Glaciers and ice sheets modulate global sea level by storing water deposited as snow on the surface, and discharging water back into the ocean through melting. Their physical state can be characterized in terms of their mass balance and dynamics. To estimate the current ice mass balance, and to predict future changes in the motion of the Greenland and Antarctic ice sheets, it is necessary to know the ice sheet thickness and the physical conditions of the ice sheet surface and bed. This information is required at fine resolution and over extensive portions of the ice sheets. A tomographic algorithm has been developed to take raw data collected by a multiple-channel synthetic aperture sounding radar system over a polar ice sheet and convert those data into two-dimensional (2D) ice thickness measurements. Prior to this work, conventional processing techniques only provided one-dimensional ice thickness measurements along profiles.

Wu, Xiaoqing↗

Implicit Formulation of Muscle Dynamics in OpenSim

Astronauts lose bone and muscle mass during spaceflight. Exercise countermeasure is the primary method for counteracting bone and muscle mass loss in space. New spacecraft exercise device concepts are currently being developed for the NASAs new crew exploration vehicle. The NASA Digital Astronaut Project (DAP) uses computational modeling to help determine if the new exercise devices will be effective as countermeasures. The NASA Digital Astronaut Project is developing the ability to utilize predictive simulation to provide insight into the change in kinematics and kinetics with a change in device and gravitational environment (1-g versus 0-g). For example, in space exercise the subject's body weight is applied in addition to the loads prescribed for musculoskeletal maintenance. How and where these loads are applied obviously directly impacts bone and tissue loads. Additionally, due to space vehicle structural requirements, exercise devices are often placed on vibration isolation systems. This changes the apparent impedance or stiffness of the device as seen by the user. Data collection under these conditions is often impractical and limited. Predictive modeling provides a means to have a virtual subject to test hypotheses. Predictive simulation provides a virtual subject for which we are able to perform studies such as sensitivity to device loading and vibration isolation without the need for laboratory kinematic or kinetic test data.Direct Collocation optimization provides an efficient means to perform task based optimization and predictive modeling. It is relatively straight forward to structure a physical exercise task in a Direct Collocation mathematical formulation: perform a motion such that you start at an initial pose, achieve a given amount of deflection i.e a squat, return to the initial pose, and minimize muscle activation cost. Direct Collocation is advantageous in that it does not require numerical integration to evaluate the objective function. Instead, the system dynamics are transformed to discrete time and the optimizer is constrained such that the solution is not considered to be a valid unless the dynamic equations are satisfied at all time points. The simulation and optimization are effectively done simultaneously. Due to the implicit integration, time steps can be more coarse than in a differential equation solver. In a gait scenario this means that that the model constraints and cost function are evaluated at 100 nodes in the gait cycle versus 10,000 integration steps in a variable-step forward dynamic simulation. Furthermore, no time is wasted on accurate simulations of movements that are far from the optimum. Constrained optimization algorithms require a Jacobian matrix that contains the partial derivatives of each of the dynamic constraints with respect to of each of the state and control variables at all time points. This is a large but sparse matrix. An implicit dynamics formulation requires computation of the dynamic residuals f as a function of the states x and their derivatives, and controls u:f(x, dxdt, u) 0If the dynamics of musculoskeletal system are formulated implicitly, the Jacobian elements are often available analytically, eliminating the need for numerical differentiation; this is obviously computationally advantageous. Additionally, implicit formulation of musculoskeletal dynamics do not suffer from singularities from low mass bodies, zero muscle activation, or other stiff system or

physical exercise↗