Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalized predictive control”

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,117 records · Page 62

Towards fully predictive gyrokinetic full- f simulations: validation and triangularity studies in TCV

Designing economical magnetic confinement fusion power plants motivates computational tools that can estimate plasma behavior from engineering parameters without direct reliance on experimental measurement of the plasma profiles. In this work, we present full-f global long-wavelength gyrokinetic simulations of edge and scrape-off layer turbulence in tokamaks that use only magnetic geometry, heating power, and particle inventory as inputs. Unlike many modeling approaches that employ free parameters fitted to experimental data, raising uncertainties when extrapolating to reactor scales. This approach directly simulates turbulence and resulting profiles through gyrokinetics without such empirical adjustments. This is achieved via an adaptive sourcing algorithm in Gkeyll that strictly controls energy injection and emulates particle sourcing due to neutral recycling. We show that the simulated kinetic profiles compare reasonably well with Thomson scattering and Langmuir probe data for Tokamak á Configuration Variable (TCV) discharge #65125, and that the simulations reproduce characteristic features such as blob transport and self-organized electric fields. Applying the same framework to study triangularity effects suggests mechanisms contributing to the improved confinement reported for negative triangularity (NT). Simulations of TCV discharges #65125 and #65130 indicate that NT increases the E x B flow shear (by about 20% in these cases), which correlates with reduced turbulent losses and a modest change in the distribution of power exhaust to the vessel wall. While the physical models contain approximations that can be refined in future work, the predictive capability demonstrated here, evolving multiple profile relaxation times with kinetic electron and ion models in hundreds of GPU hours, indicates the feasibility of using Gkeyll to support design studies of fusion devices.

Hoffmann, Antoine Cyril David [Princeton Plasma Ph↗

Partially averaged Navier-Stokes closure modeling for variable-density turbulent flow

We report this work extends the framework of the partially averaged Navier-Stokes (PANS) equations to variable-density flow, i.e., multimaterial and/or compressible mixing problems with density variations and production of turbulence kinetic energy by both shear and buoyancy mechanisms. The proposed methodology is utilized to derive the PANS BHR-LEVM closure. This includes a priori testing to analyze and develop guidelines toward the efficient selection of the parameters controlling the physical resolution and, consequently, the range of resolved scales of PANS. Two archetypal test-cases involving transient turbulence, hydrodynamic instabilities, and coherent structures are used to illustrate the accuracy and potential of the method: the Taylor-Green vortex at Reynolds number Re = 3000, and the Rayleigh-Taylor flow at Atwood number 0.5 and (Re) max ≈ 500. These representative problems, for which turbulence is generated by shear and buoyancy processes, constitute the initial validation space of the new model, and their results are comprehensively discussed in two subsequent studies. The computations indicate that PANS can accurately predict the selected flow problems, resolving only a fraction of the scales of large-eddy simulation and direct numerical simulation strategies. The results also reiterate that the physical resolution of the PANS model must guarantee that the key instabilities and coherent structures of the flow are resolved. The remaining scales can be modeled through an adequate turbulence scale-dependent closure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Socioeconomically-inspired modeling to justify use of fine-grain mobility data

When designing measures to control infectious disease spread, it is crucial to understand the structure of the population for which interventions are being implemented. Recent work has highlighted the need for models that incorporate demographic heterogeneity not just in age structure but also by socioeconomic status (SES). Appropriately capturing additional sources of population heterogeneity requires considerable data and model development. To understand the potential disagreement between SES-explicit or SES-agnostic disease models, we adapted Sandia’s Adaptive Recovery Model (ARM) model to consider differences in contact structure and mortality by Social Vulnerability Index (SVI) on a theoretical network. We also incorporated an Average network that did not consider SVI. By exploring disparities in vaccine and PPE uptake by SES and comparing to Average networks, as well as analyzing the influence of global vs. local contact, we found that the two model constructions often predicted different outcomes. Whether these differences are truly reflective of incorporating SES, and which model most closely represents reality, merits further investigation.

99 GENERAL AND MISCELLANEOUS↗

Processes Controlling Water Vapor in the Winter Arctic Stratospheric Middleworld

Water vapor in the winter arctic stratospheric middleworld is import-an: for two reasons: (1) the arctic middleworld is a source of air for the upper Troposphere because of the generally downward motion, and thus its water vapor content helps determine upper tropospheric water, a critical part of the earth's radiation budget; and (2) under appropriate conditions, relative humidities will be large, even to the point of stratospheric cirrus cloud formation, leading to the production of active chlorine species that could destroy ozone. On a number of occasions during SOLVE, clouds were observed in the stratospheric middleworld by the DC-8 aircraft. These tended to coincide with regions of low temperatures, though some cases suggest water vapor enhancements due to troposphere-to-stratosphere transport. The goal of this work is to understand the importance of processes in and at the edge of the arctic stratospheric middleworld in determining water vapor at these levels. Specifically, is water vapor at these levels determined largely by the descent of air from above, or are clouds both within and at the edge of the stratospheric middleworld potentially important? How important is troposphere-to-stratosphere transport of air in determining stratospheric middleworld water vapor content? To this end, we will first examine the minimum saturation mixing ratios along theta/EPV tubes during the SOLVE winter and compare these with DC-8 water vapor observations. This will be a rough indicator of how high relative humidities can get, and the likelihood of cirrus cloud formation in various parts of the stratospheric middleworld. We will then examine saturation mixing ratios along both diabatic and adiabatic trajectories, comparing these values with actual aircraft water vapor observations, both in situ and remote. Finally, we will attempt to actually predict water vapor using minimum saturation mixing ratios along trajectories, cloud injection (derived from satellite imagery) along trajectories, and suitable initial conditions.

Pfister, Leonhard↗

Interplay of coulomb and exciton–phonon coupling controls singlet fission dynamics in two pentacene polymorphs

Pentacene is an important model organic semiconductor in both the singlet exciton fission (SF) and organic electronics communities. We have investigated the effect of changing crystal structure on the SF process, generating multiple triplet excitons from an initial singlet exciton, and subsequent triplet recombination. Unlike for similar organic semiconductors that have strong SF sensitive to polymorphism, we find almost no quantitative difference between the kinetics of triplet pair (TT) formation in the two dominant polymorphs of pentacene. Both pairwise dimer coupling and momentum-space crystal models predict much faster TT formation from the bright singlet excited state of the Bulk vs ThinFilm polymorph, contrasting with the experiment. GW and Bethe–Salpeter equation calculations, including exciton–phonon coupling, reveal that ultrafast phonon-driven transitions in the ThinFilm polymorph compensate the intrinsically slower purely Coulomb-mediated TT formation channel, rationalizing the similarity in observed rates. Taking into account the influence of subtle structural distinctions on both the direct and phonon-mediated SF pathways reveals a predictive capability to these methods, expected to be applicable to a wide variety of molecular crystals.

14 SOLAR ENERGY↗

A Physics-Constrained Deep Learning Model for Simulating Multiphase Flow in 3D Heterogeneous Porous Media

Physics-based simulators for multiphase flow in porous media emulate nonlinear processes with coupled physics, and usually require extensive computational resources for software development, maintenance and simulation execution. As a result, a huge demand exists for fast modeling of coupled processes in a wide range of subsurface applications including geological sequestration, hydrocarbon recovery and geothermal energy extraction. In this work, an efficient physics-constrained deep learning model is developed for solving multiphase flow in 3-Dimensional (3D) heterogeneous porous media. The model fully leverages the spatial topology predictive capability of convolutional neural networks, specifically U-Net with successive contracting and expansive steps, and is coupled with an efficient continuity-based smoother to predict flow responses that need spatial continuity. Furthermore, the transient regions are penalized to steer the training process such that the model can accurately capture flow in these regions. The model takes inputs including properties of porous media, fluid properties and well controls, and predicts the temporal-spatial evolution of the state variables (pressure and saturation). While maintaining the continuity of fluid flow, the 3D spatial domain is decomposed into 2D images for reducing training cost, and the decomposition results in an increased number of training data samples and better training efficiency. Additionally, a surrogate model is separately constructed as a postprocessor to calculate well flow rate based on the predictions of state variables from the deep learning model. We use the example of CO 2 injection into saline aquifers, and apply the physics-constrained deep learning model that is trained from physics-based simulation data and emulates the physics process. The model performs prediction with a speedup of ~ 1400 times compared to physics-based simulations, and the average temporal errors of predicted pressure and saturation plumes are 0.27% and 0.099% respectively. Furthermore, water production rate is efficiently predicted by a surrogate model for well flow rate, with a mean error less than 5%. Therefore, with its unique scheme to cope with the fidelity in fluid flow in porous media, the physics-constrained deep learning model can become an efficient predictive model for computationally demanding inverse problems or other coupled processes.

58 GEOSCIENCES↗

Investigation of the Influence of Microgravity on Transport Mechanism in a Virtual Spaceflight Chamber: A Flight Definition Program

A need exists for understanding precisely how particles move and interact in a fluid in the absence of gravity. Such understanding is required, for example, for modeling and predicting crystal growth in space where crystals grow from solution around nucleation sites as well as for any study of particles or bubbles in liquids or in experiments where particles are used as tracers for mapping microconvection. We have produced an exact solution to the general equation of motion of particles at extremely low Reynolds number in microgravity that covers a wide range of interesting conditions. We have also developed diagnostic tools and experimental techniques to test the validity of the general equation . This program, which started in May, 1998, will produce the flight definition for an experiment in a microgravity environment of space to validate the theoretical model. We will design an experiment with the help of the theoretical model that is optimized for testing the model, measuring g, g-jitter, and other microgravity phenomena. This paper describes the goals, rational, and approach for the flight definition program. The first objective of this research is to understand the physics of particle interactions with fluids and other particles in low Reynolds number flows in microgravity. Secondary objectives are to (1) observe and quantify g-jitter effects and microconvection on particles in fluids, (2) validate an exact solution to the general equation of motion of a particle in a fluid, and (3) to characterize the ability of isolation tables to isolate experiments containing particle in liquids. The objectives will be achieved by recording a large number of holograms of particle fields in microgravity under controlled conditions, extracting the precise three-dimensional position of all of the particles as a function of time and examining the effects of all parameters on the motion of the particles. The feasibility for achieving these results has already been established in the ongoing ground-based NRA, which led to the "virtual spaceflight chamber" concept.

Trolinger, James D.↗

An Advanced Trajectory-Based Operations Prototype Tool and Focus Group Evaluation

Trajectory-based operations (TBO) is a key concept in the Next Generation Air Transportation System transformation of the National Airspace System (NAS) that will increase the predictability and stability of traffic flows, support a common operational picture through the use of digital data sharing, facilitate more effective collaborative decision making between airspace users and air navigation service providers, and enable increased levels of integrated automation across the NAS. The National Aeronautics and Space Administration (NASA) has been developing trajectory-based systems to improve the efficiency of the NAS during specific phases of flight and is now also exploring Advanced 4-Dimensional Trajectory (4DT) operational concepts that will integrate these technologies and incorporate new technology where needed to create both automation and procedures to support gate-to-gate TBO. A TBO Prototype simulation toolkit has been developed that demonstrates initial functionality that may reside in an Advanced 4DT TBO concept. Pilot and controller subject matter experts (SMEs) were brought to the Air Traffic Operations Laboratory at NASA Langley Research Center for discussions on an Advanced 4DT operational concept and were provided an interactive demonstration of the TBO Prototype using four example scenarios. The SMEs provided feedback on potential operational, technological, and procedural opportunities and concerns. After viewing the interactive demonstration scenarios, the SMEs felt the operational capabilities demonstrated would be useful for performing TBO while maintaining situation awareness and low mental workload. The TBO concept demonstrated produced defined routings around weather which resulted in a more organized, consistent flow of traffic where it was clear to both the controller and pilot what route the aircraft was to follow. In general, the controller SMEs felt that traffic flow management should be responsible for generating and negotiating the operational constraints demonstrated, in cooperation with the Air Traffic Control System Command Center, while air traffic control should be responsible for the implementation of those constraints. The SMEs also indicated that digital data communications would be very beneficial for TBO operations and would result in less workload due to reduced communications, would eliminate issues due to language barriers and frequency problems, and would make receiving, loading, accepting, and executing clearances easier, less ambiguous, and more expeditious. This paper describes an Advanced 4DT operational concept, the TBO Prototype, the demonstration scenarios and methods used, and the feedback obtained from the pilot and controller SMEs in this focus group evaluation.

Guerreiro, Nelson M.↗

Investigating Material Properties of Subsurface Rock Formations Modified by Engineering Mineral Precipitation (Final Scientific and Technical Report)

Montana State University’s (MSU) Energy Research Institute (ERI), in collaboration with the Center for Biofilm Engineering (CBE) and the Department of Civil Engineering (CE), has conducted a long‐term research program aimed at developing a novel cementing agent to address wellbore integrity and reduce the unwanted upward migration of fluids and greenhouse gases from the subsurface. The primary technology developed through this research program is known as ureolysis‐induced calcite precipitation (UICP), which harnesses bio‐chemical processes to precipitate calcium carbonate (CaCO 3 ). The same general process can also be called microbially-induced calcium carbonate precipitation (MICP) when microbes provide the process-catalyzing urease enzyme. Both terms are used in this report. Results have conclusively demonstrated that, if properly controlled, UICP can successfully seal fractures, high permeability zones, and compromised cement in the vicinity of wellbores and in nearby caprock. This technology has been successfully deployed to mitigate annular leakage in two test wells and over sixty commercial wells with a 100% success rate. This success in downhole deployment generates consideration of other subsurface applications where UICP could provide benefit to the energy sector, such as shale property modification for unconventional oil and gas recovery. The focus of this research project was to investigate fundamental material and mechanical properties of select shale cores and analyze how these properties change due to engineered mineral precipitation with the intent to control these properties to achieve a range of engineering objectives. Ultimately, the project aim was to identify valuable new areas where application of UICP might contribute to national energy security and environmental protection. The research workplan coupled UICP treatment of core samples, nuclear magnetic resonance (NMR) characterization, and mechanical strength testing at MSU with advanced X‐Ray micro-computed tomography (μCT) imaging and numerical modeling performed by collaborators at two national laboratories, the National Energy Technology Laboratory (NETL) and Lawrence Berkeley National Laboratory (LBNL). Experimental results are useful to inform geo-mechanical models which could be applied to predict mineralized rock formation behavior at field scale. Our findings suggest that NMR and μCT methods to detect and quantify biomineral formation in shale fractures are complementary and consistent with each other. Either could be used to estimate the volume of new mineral formed by UICP in shale fractures. The use of surfactants and guar gum to enhance biomineral precipitation in shale fractures merits further research. UICP can, under some conditions, increase the tensile strength of sealed shale fractures beyond that of the intact shale. These findings demonstrate that continued research in this area may be valuable to understanding and improving shale resource recovery techniques.

58 GEOSCIENCES↗

Multi-harmonic Rutherford island theory

Rutherford island theory, which governs the nonlinear evolution of tearing modes in tokamak plasmas, is generalized to take into account situations in which the conventional one-harmonic approximation is not valid. The analysis incorporates non-inductive currents driven by radio frequency (RF) electromagnetic waves injected into the plasma. A multi-harmonic tearing mode dispersion relation is derived that takes the form of a nonlinear inhomogeneous matrix eigenvalue problem. The dispersion relation is solved in the so-called two-harmonic approximation, in which only the principal Fourier harmonic of the perturbed magnetic flux and its first overtone are included in the calculation. In the absence of RF current drive, the nonlinear behavior of a tearing mode predicted in the two-harmonic approximation does not differ substantially from that predicted in the one-harmonic approximation. On the other hand, RF current drive that is sufficiently localized in the vicinity of the O-points of the mode's magnetic island chain is capable of triggering bifurcations of the O-points (which is impossible in the one-harmonic approximation). However, the current drive is incapable of triggering bifurcations of the island X-points. This finding is significant because Bardóczi and Evans recently observed bifurcations of magnetic island chain O-points in the presence of RF current drive in the DIII-D tokamak but did not observe bifurcations of the X-points. Finally, the changes in the topology of the magnetic island flux-surfaces induced by RF current drive are found to facilitate the stabilization of the tearing mode.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

High-Efficiency High-Resolution Global Model Developments at the NASA Goddard Data Assimilation Office

The Data Assimilation Office (DAO) has been developing a new generation of ultra-high resolution General Circulation Model (GCM) that is suitable for 4-D data assimilation, numerical weather predictions, and climate simulations. These three applications have conflicting requirements. For 4-D data assimilation and weather predictions, it is highly desirable to run the model at the highest possible spatial resolution (e.g., 55 km or finer) so as to be able to resolve and predict socially and economically important weather phenomena such as tropical cyclones, hurricanes, and severe winter storms. For climate change applications, the model simulations need to be carried out for decades, if not centuries. To reduce uncertainty in climate change assessments, the next generation model would also need to be run at a fine enough spatial resolution that can at least marginally simulate the effects of intense tropical cyclones. Scientific problems (e.g., parameterization of subgrid scale moist processes) aside, all three areas of application require the model's computational performance to be dramatically improved as compared to the previous generation. In this talk, I will present the current and future developments of the "finite-volume dynamical core" at the Data Assimilation Office. This dynamical core applies modem monotonicity preserving algorithms and is genuinely conservative by construction, not by an ad hoc fixer. The "discretization" of the conservation laws is purely local, which is clearly advantageous for resolving sharp gradient flow features. In addition, the local nature of the finite-volume discretization also has a significant advantage on distributed memory parallel computers. Together with a unique vertically Lagrangian control volume discretization that essentially reduces the dimension of the computational problem from three to two, the finite-volume dynamical core is very efficient, particularly at high resolutions. I will also present the computational design of the dynamical core using a hybrid distributed-shared memory programming paradigm that is portable to virtually any of today's high-end parallel super-computing clusters.

Lin, Shian-Jiann↗

High-Efficiency High-Resolution Global Model Developments at the NASA Goddard Data Assimilation Office

The Data Assimilation Office (DAO) has been developing a new generation of ultra-high resolution General Circulation Model (GCM) that is suitable for 4-D data assimilation, numerical weather predictions, and climate simulations. These three applications have conflicting requirements. For 4-D data assimilation and weather predictions, it is highly desirable to run the model at the highest possible spatial resolution (e.g., 55 kin or finer) so as to be able to resolve and predict socially and economically important weather phenomena such as tropical cyclones, hurricanes, and severe winter storms. For climate change applications, the model simulations need to be carried out for decades, if not centuries. To reduce uncertainty in climate change assessments, the next generation model would also need to be run at a fine enough spatial resolution that can at least marginally simulate the effects of intense tropical cyclones. Scientific problems (e.g., parameterization of subgrid scale moist processes) aside, all three areas of application require the model's computational performance to be dramatically improved as compared to the previous generation. In this talk, I will present the current and future developments of the "finite-volume dynamical core" at the Data Assimilation Office. This dynamical core applies modem monotonicity preserving algorithms and is genuinely conservative by construction, not by an ad hoc fixer. The "discretization" of the conservation laws is purely local, which is clearly advantageous for resolving sharp gradient flow features. In addition, the local nature of the finite-volume discretization also has a significant advantage on distributed memory parallel computers. Together with a unique vertically Lagrangian control volume discretization that essentially reduces the dimension of the computational problem from three to two, the finite-volume dynamical core is very efficient, particularly at high resolutions. I will also present the computational design of the dynamical core using a hybrid distributed- shared memory programming paradigm that is portable to virtually any of today's high-end parallel super-computing clusters.

Lin, Shian-Jiann↗

One-Carbon Metabolism and SANS 2022 Update

Spaceflight Associated Neuro-ocular Syndrome, or SANS, affects a subset of astronauts (1, 2), and biochemical evidence has documented differences in those astronauts (3). Specifically, they had higher circulating concentrations of metabolites of the one-carbon metabolic pathway (1C), including homocysteine, and these concentrations were higher before flight (3). After ruling out many potential confounding factors in these otherwise healthy individuals (e.g., sex, kidney function, vitamin status, coffee consumption), a study of genetics was warranted. In an initial pilot effort, we documented a genetic predisposition to develop ophthalmic changes after long-duration space flight (4). That is, from a limited study of 5 single-nucleotide polymorphisms (SNPs), we found that the G allele for the MTRR A66G SNP was associated with a greater risk of choroidal folds and cotton-wool spots after flight, and the C allele for SHMT1 C1420T was protective against optic disc edema (4). These data provide a potential pathway for understanding why some individuals develop SANS, while others do not. The initial pilot study of 5 SNPs yielded striking findings, but the 1C pathway is far more complex. An effort was undertaken to examine more than 500 1C SNPs to see if a broader examination could help illuminate this association. That work is ongoing. The astronaut findings led us to advocate for the inclusion of 1C pathway genetic and biochemistry testing on other SANS-related projects, noting that genetics might help identify responders, non-responders, or outliers. The first such effort yielded evidence of an association of specific forms of the MTRR and SHMT-1 SNPs and vitamin B12 status with end-tidal CO2 after acute carbon dioxide exposure (5). The second such effort led to the identification that individuals exposed to strict head-down tilt and CO2 for 30-d who developed optic disc edema also had risk alleles for the two SNPs described above (6). Additionally, we identified a clinical population with many characteristics either attributed or purported to be involved in the ocular changes seen in affected astronauts: women with polycystic ovary syndrome (PCOS). PCOS is a condition of androgen excess and anovulatory menstrual cycles. The shared characteristics and clinical findings between SANS and PCOS generally include higher circulating homocysteine concentrations, increased retinal nerve fiber layer thickness, increased androgen concentrations (or responses), and altered carbohydrate metabolism. To our knowledge, no study has examined whether women with PCOS have asymptomatic ophthalmic anomalies observed in astronauts with SANS. While researchers have evaluated the one-carbon metabolism pathway polymorphisms of PCOS patients, and initial studies show an association with certain one-carbon polymorphisms, none have looked at the set of SNPs identified in our studies that are associated with ophthalmic changes in astronauts. Accordingly, we designed a study to evaluate the association of one-carbon pathway SNPs and ophthalmic findings in patients with PCOS and/or IIH compared to controls. Subjects provided blood samples for vitamin and one carbon biochemistry analyses, an extensive analysis of >500 SNPs associated with one carbon metabolism and had eye examinations and ocular imaging. Data analysis are underway. The data collected to date have shown associations between one carbon pathway biochemistry and genetics and incidence of SANS. The mechanisms for SANS has yet to be identified, although many hypotheses exist. Based on our data, we have developed (8, 9) and expanded (6) a multi-hit hypothesis for how these seemingly disparate findings could be linked. While intriguing, the hypothesis represents the starting point for further research. We aim to clarify the relationship between B-vitamin status and genetics with regard to the risk of SANS. Ultimately, understanding the mechanism(s) behind this will provide a means to predict, prevent, or treat these ophthalmologic pathologies in astronauts, and terrestrial populations.

S M Smith↗

Monitoring Airspace Complexity and Determining Contributing Factors

The national airspace has evolved over many years to accommodate increased traffic demand [1] while simultaneously maintaining one of the safest forms of transportation [2], [3]. One of the reasons for this success is the ability of the system and the operators to adapt and accommodate to situations that routinely disrupt optimal operations. These situations may include: adverse weather, delays, early arrivals, equipment outages, and other factors that are outside the operators ability to control. These factors can lead to states where automation is unable to properly handle these issues and therefore air traffic controllers and pilots have to intervene, ultimately increasing communication between operators resulting in higher workload. As controller workload increases to handle sub-optimal operating conditions this can be viewed as an increase in complexity. The reasoning for this is because humans are now required to make tactical decisions in response to external factors, resulting in a departure from the strategic plan where operations would be more efficiently managed. Human operators control airspace complexity under rigid regulations that are constantly changing. The airspace is divided into sectors and the number of aircraft assigned to each controller is limited for safe handling. There has been past work that devised airspace complexity metrics in commercial aviation and related these metrics to controller workload (e.g., [4],[5]). The upper bounds on the system load are pre-determined. Such bounds on complexity make for a safe system, but the system cannot scale and adapt to autonomous, dense, and heterogeneous traffic, including the many types of Unmanned Aerial Vehicles (UAVs) envisioned to be added to the operations. We hypothesize that, as traffic density and heterogeneity grow, and other key metrics change, there will be phase transitions at which the way traffic should be managed changes significantly [6]. We offer a method for in-time detection of contributing factors that lead to phase transitions, characterized by increased complexity. To the best of our knowledge, there is no tool similar to our proposed effort that identifies such contributing factors or precursor patterns. To define the scope we are proposing to measure complexity from the viewpoint of the Terminal Radar Approach Control Facilities (TRACON) controller’s perspective. In particular we are analyzing arrivals into KSFO. With safety as the top concern for airspace operators, it is important to recognize that as density and heterogeneity grow, the focus of the system will change. Times of the day when the airspace has low density and heterogeneity, the flights will follow more efficient paths where the aircraft move on established routes that are more or less directly to the destination. However, when density and heterogeneity increases, the system will begin changing focus to avoiding conflicts and collisions and route the flights in a more flexible way. Higher flexibility requires more communication and coordination between controllers and pilots which the current automation is unable to handle. This paper proposes a novel approach that monitors airspace complexity at multiple scales, uses a Machine Learning-based tool that predicts when operations will transition to a regime of greater complexity, and identifies actions that can reduce the complexity while still maintaining efficient and safe operations. We demonstrate our proposed approach using data from multiple complementary sources. This includes, but is not limited to: historical aircraft surveillance data from NASA’s Sherlock Data Warehouse [7], METAR weather data, and airport configuration data from Aviation System Performance Metrics (ASPM). The surveillance data flight paths are sampled at a variable sample rate — increasing as the aircraft approaches the airport. This is due to how Sherlock manages flight track stitching between different radar facilities which have different sampling rates. The weather and performance data are logged at defined intervals throughout the day at a courser refresh rate. In addition to the logged data and metrics, we leverage pre-defined Standard Terminal Arrival Routes (STARs) procedures to characterize the path of each flight. Each flight files for one of these routes in the flight plan well before entering the terminal airspace, and approximately follows the route until it leaves the STAR, typically on the final fix of a runway transition. However, most flights do not always fly the full STAR procedure to completion [8], but the majority do adhere to the fixes within the common route of the procedure. Our approach leverages fixes in the common route of each of the STARs to build a reference path to the airport. This allows us to characterize the flight paths in what we are defining as the “maneuvering area” (the airspace between the STAR and before the flight is lined up on the runway’s final approach) to determine how off nominal the flights are to calculate its complexity score. Determining airspace complexity is a concept that does not have a concrete answer. In designing this metric, we consider what increases the workload for the air traffic controllers. Consequently more specialized vectoring maneuvers results in higher workload. Accordingly, we start with a theory: each flight has a direct path it takes from the STAR’s common route to the final approach’s outer marker fix for the flight’s landing runway. It is important to note that the direct path is only used as a reference. If the majority of the flights have a large consistent offset as compared to other routes it does not necessarily mean that those flights have higher complexity. We are merely building a distribution based on this direct path for that particular STAR and runway pair to determine the normal mode of operations for that route. Flights that are in the upper tail of these distributions will result in higher complexity scores and flights that fly in the median will represent the normal mode of operations and therefore will have lower complexity scores. Since flights following each STAR route take different paths to the airport, we have a different distribution for each STAR route and therefore can model these distributions to compute a complexity score from their respective normalized distributions. To evaluate the effectiveness of our proposed airspace complexity metric we will compare against an established approach based on trajectory clustering [9]. This unsupervised learning technique consists of the following steps: (1) identify the general maneuvering areas (waypoints) by performing $\kappa$-means or DBSCAN clustering on locations where aircraft frequently turn based on the surveillance radar track data, (2) map flight trajectories onto sequences of waypoints, and (3) cluster the sequences based on their common subsequences. From a high-level perspective, this baseline model learns nominal operations in the airspace through the sequence of waypoints that are representative of where aircraft change direction and defines deviations from the nominal operations as “complex.” Therefore, more deviations from the nominal operations correspond to higher complexity values. For our validation, we re-implemented this technique and tune model hyper-parameters to correctly detect waypoints for the arrival traffic into the San Francisco bay area. We will compute the complexity measure over a one-year period using our proposed technique as well as the baseline. Our validation will be based on each technique’s ability to detect a set of undesirable outcomes (e.g., go-arounds, holding patterns, average time in the airspace, etc.). Since our current complexity metric is derived from the offset from the direct reference path, it’s important to understand what causes these offsets. In many of the flights with high offset distance, flights performing holding patterns and S turns can be observed. These maneuvering tactics are utilized to add distance between the aircraft and the destination runway to prevent multiple flights from having conflicting arrival times. In order to predict a rise in complexity (or the precursor to complexity), it’s necessary to be able to identify these potential conflicts (which in turn, result in higher offsets). To do this, we define a “representative flight” for each STAR route and runway pair. This flight is approximately the path the flight would take if there was a clear path with no other flights in the airspace — including the time remaining to the airport. We first identify the flights for a given STAR runway pair using the offset to the reference path distributions that fall between the 44-55 percentiles. This yields the flights that conform to the most normal mode of operation. Each of these flights is partitioned based on the percent complete from the entry point into the maneuvering areas from 0\% – 100\% complete. Then for each percent “bin”, we take the median value of the flight’s latitude/longitude coordinates, airspeed, and (non causal) time remaining to the airport to construct a lookup table for each percent complete bin on a given route. As a flight enters the maneuvering area, we can find the estimated arrival time of a flight to the airport by finding the closest point to the representative path’s percent complete bin (relative to the flight’s current position at any snapshot in the airspace) and therefore retrieve the corresponding remaining time left on the “representative path”. We assume that the flight will follow the representative path to completion when deriving these estimates. We can then compare these estimated arrival times against other flights for the same snapshot in time to identify potential conflicts. If more flights are estimated to arrive within a tolerance window than there are runways available, then we have a potential conflict. We can use this derived measure along with other factors expected to add disruption to the operation such as weather and runway configuration changes as an input to machine learning tools to detect precursors that increases in our complexity measure. This novel method will assist in uncovering insights into the contributing factors that lead to increased complexity that may allow for in-time responses to avoid reaching a high complexity state in the airspace.

complexity↗

The Effect of Habitual Smoking on VO2max

VO2max is associated with many factors, including age, gender, physical activity, and body composition. It is popularly believed that habitual smoking lowers aerobic fitness. PURPOSE: to determine the effect of habitual smoking on VO2max after controlling for age, gender, activity and BMI. METHODS: 2374 men and 375 women employed at the NASA/Johnson Space Center were measured for VO2max by indirect calorimetry (RER>=1.1), activity by the 11 point (0-10) NASA Physical Activity Status Scale (PASS), BMI and smoking pack-yrs (packs day*y of smoking). Age was recorded in years and gender was coded as M=1, W=0. Pack.y was made a categorical variable consisting of four levels as follows: Never Smoked (0), Light (1-10), Regular (11-20), Heavy (>20). Group differences were verified by ANOVA. A General Linear Models (GLM) was used to develop two models to examine the relationship of smoking behavior on VO2max. GLM #1(without smoking) determined the combined effects of age, gender, PASS and BMI on VO2max. GLM #2 (with smoking) determined the added effects of smoking (pack.y groupings) on VO2max after controlling for age, gender, PASS and BMI. Constant errors (CE) were calculated to compare the accuracy of the two models for estimating the VO2max of the smoking subgroups. RESULTS: ANOVA affirmed the mean VO2max of each pack.y grouping decreased significantly (p<0.01) as the level of smoking exposure increased. GLM #1 showed that age, gender, PASS and BMI were independently related with VO2max (R2 = 0.642, SEE = 4.90, p<0.001). The added pack.y variables in GLM #2 were statistically significant (R2 change = 0.7%, p<0.01). Post hoc analysis showed that compared to Never Smoked, the effects on VO2max from Light and Regular smoking habits were -0.83 and -0.85 ml.kg- 1.min-1 respectively (p<0.05). The effect of Heavy smoking on VO2max was -2.56 ml.kg- 1.min-1 (p<0.001). The CE s of each smoking group in GLM #2 was smaller than the CE s of the smoking group counterparts in GLM #1. CONCLUSIONS: After accounting for the effects of gender, age, PASS and BMI the effect of habitual smoking on reducing VO2max is minimal, about 0.85 ml/kg/min, until the habit exceeds 20 pack.y at which point an additional decrease of 1.71 ml/kg/min is noted. Adding pack.y data improves the accuracy of predicting the VO2max of smokers.

Wier, Larry T.↗

Credible practice of modeling and simulation in healthcare: ten rules from a multidisciplinary perspective

The complexities of modern biomedicine are rapidly increasing. Thus, modeling and simulation have become increasingly important as a strategy to understand and predict the trajectory of pathophysiology, disease genesis, and disease spread in support of clinical and policy decisions. In such cases, inappropriate or ill-placed trust in the model and simulation outcomes may result in negative outcomes, and hence illustrate the need to formalize the execution and communication of modeling and simulation practices. Although verification and validation have been generally accepted as significant components of a model’s credibility, they cannot be assumed to equate to a holistic credible practice, which includes activities that can impact comprehension and in-depth examination inherent in the devel-opment and reuse of the models. For the past several years, the Committee on Credible Practice of Modeling and Simulation in Healthcare, an interdisciplinary group seeded from a U.S. interagency initiative, has worked to codify best practices. Here, we provide Ten Rules for credible practice of modeling and simulation in healthcare developed from a comparative analysis by the Committee’s multidisciplinary membership, followed by a large stakeholder com-munity survey. These rules establish a unified conceptual framework for modeling and simulation design, implementation, evaluation, dissemination and usage across the modeling and simulation life-cycle. While biomedical science and clinical care domains have somewhat different requirements and expectations for credible practice, our study converged on rules that would be useful across a broad swath of model types. In brief, the rules are: (1) Define context clearly. (2) Use contextually appropriate data. (3) Evaluate within context. (4) List limitations explicitly. (5) Use version control. (6) Document appropriately. (7) Disseminate broadly. (8) Get independent reviews. (9) Test competing imple-mentations. (10) Conform to standards. Although some of these are common sense guidelines, we have found that many are often missed or misconstrued, even by seasoned practitioners. Computational models are already widely used in basic science to generate new biomedical knowledge. As they penetrate clinical care and healthcare policy, contributing to personalized and precision medicine, clinical safety will require established guidelines for the credible practice of modeling and simulation in healthcare.

Credibility↗

Lessons Learned from the Aquarius Radiometer Engineering Model

The Aquarius instrument is a combination microwave radiometer and scatterometer operating in the L-band for the measurement of sea surface salinity (SSS) over the open-ocean. "The Aquarius science goals are to observe and model the processes that relate salinity variations to climatic changes in the global cycling of water and to understand how these variations influence the general ocean circulation. By measuring salinity globally and synoptically every month for 3 years, Aquarius will provide an unprecedented new view of the ocean's role in climate." (http://aquarius.asfc.nasa.gov on Jan. 4, 2007). To launch in 2009, Aquarius will provide 8-day global coverage and monthly-averaged salinity maps at approx.120-km spatial resolution. Because the microwave emissivity of the ocean surface is weakly dependent upon salinity, the Aquarius radiometer remain sufficiently stable to provide accurate salinity measurements. An engineering model (EM) of one radiometer receiver was developed and tested for functionality and performance at NASA's Goddard Space Flight Center. It was tested for long-term stability starting July 27,2006. Two discrepancies where found at the start of testing and a four-month correction process ensued. This paper covers the testing and debugging of the EM radiometer through resolution of the problems. The radiometer is required to have calibration stability of less than or equal to 0.15 K over 7 days. The radiometer is in a three-beam push-broom configuration with each beam being an independent three-state Dicke radiometer (antenna, reference load, and additive noise). Internal noise diodes and reference loads provide periodic and rapid calibration of the receiver. Front-end loss contributions must be corrected as the lossy components are outside the calibration loop. A key design feature for obtaining stability is active thermal control of the front-end lossy components and the RF receiver. The thermal stability is predicted to be about 0.1 C over an orbit and no more than 1 C over a year. The radiometric stability requirement is 0.13 K rms 1 7 days without using temperature corrections. Temperature corrections and post-launch calibration are expected to extend the stability period through the life of the mission.

Piepmeier, Jeffrey↗

Quantifying uncertainty in uranium concentration measurements via K-edge densitometry

This study quantifies the uncertainty in uranium concentration predictions of fluoride and chloride-based salts within a steel pipe using K-edge densitometry. Modeling and simulation was conducted with the Monte Carlo N-Particle Transport (MCNP) code. The quality of of this technique’s prediction in a pipe requires proper characterization of the pipe’s thickness, which is dependent on the source size and axial offset from the pipe centerline. The thickness was determined as either the center-line thickness seen by the X-ray source or an average value determined through random sampling. Generally, the predicted concentrations were slightly better at lower offset with the random sampling thickness and using the center-line thickness for the highest offsets. For a line-beam source and varying axial offsets, the relative error of concentration was within 1% of the true value but uncertainty increased by 2 orders of magnitude. Similarly, for no axial offset, the relative error was significantly less than 1% while no trend for uncertainty was found. However, at the largest possible offset for a given source size, the concentrations become erroneous and greater than the allowable 1% relative error. Furthermore, high offsets tended to increase the variance of the transmission spectra by 3 orders of magnitude.

Characterization and Analytical Technique↗