Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complete linearization method”

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 451 records · Page 25

A Massively Parallel Implementation of the CCSD(T) Method Using the Resolution-of-the-Identity Approximation and a Hybrid Distributed/Shared Memory Parallelization Model

In this work, a parallel algorithm is described for the coupled-cluster singles and doubles method augmented with a perturbative correction for triple excitations [CCSD(T)] using the resolution-of-the-identity (RI) approximation for two-electron repulsion integrals (ERIs). The algorithm bypasses the storage of four-center ERIs by adopting an integral-direct strategy. The CCSD amplitude equations are given in a compact quasi-linear form by factorizing them in terms of amplitude-dressed three-center intermediates. A hybrid MPI/OpenMP parallelization scheme is employed, which uses the OpenMP-based shared memory model for intranode parallelization and the MPI-based distributed memory model for internode parallelization. Parallel efficiency has been optimized for all terms in the CCSD amplitude equations. Two different algorithms have been implemented for the rate-limiting terms in the CCSD amplitude equations that entail and -scaling computational costs, where N O and N V denote the number of correlated occupied and virtual orbitals, respectively. One of the algorithms assembles the four-center ERIs requiring N V 4 and N O 2 N V 2 -scaling memory costs in a distributed manner on a number of MPI ranks, while the other algorithm completely bypasses the assembling of quartic memory-scaling ERIs and thus largely reduces the memory demand. It is demonstrated that the former memory-expensive algorithm is faster on a few hundred cores, while the latter memory-economic algorithm shows a better strong scaling in the limit of a few thousand cores. The program is shown to exhibit a near-linear scaling, in particular for the compute-intensive triples correction step, on up to 8000 cores. The performance of the program is demonstrated via calculations involving molecules with 24–51 atoms and up to 1624 atomic basis functions. As the first application, the complete basis set (CBS) limit for the interaction energy of the π-stacked uracil dimer from the S66 data set has been investigated. This work reports the first calculation of the interaction energy at the CCSD(T)/aug-cc-pVQZ level without local orbital approximation. The CBS limit for the CCSD correlation contribution to the interaction energy was found to be -8.01 kcal/mol, which agrees very well with the value -7.99 kcal/mol reported by Schmitz, Hättig, and Tew [ Phys. Chem. Chem. Phys. 2014 , 16 , 22167-22178]. The CBS limit for the total interaction energy was estimated to be -9.64 kcal/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Methodology issues concerning the accuracy of kinematic data collection and analysis using the ariel performance analysis system

Kinematics, the study of motion exclusive of the influences of mass and force, is one of the primary methods used for the analysis of human biomechanical systems as well as other types of mechanical systems. The Anthropometry and Biomechanics Laboratory (ABL) in the Crew Interface Analysis section of the Man-Systems Division performs both human body kinematics as well as mechanical system kinematics using the Ariel Performance Analysis System (APAS). The APAS supports both analysis of analog signals (e.g. force plate data collection) as well as digitization and analysis of video data. The current evaluations address several methodology issues concerning the accuracy of the kinematic data collection and analysis used in the ABL. This document describes a series of evaluations performed to gain quantitative data pertaining to position and constant angular velocity movements under several operating conditions. Two-dimensional as well as three-dimensional data collection and analyses were completed in a controlled laboratory environment using typical hardware setups. In addition, an evaluation was performed to evaluate the accuracy impact due to a single axis camera offset. Segment length and positional data exhibited errors within 3 percent when using three-dimensional analysis and yielded errors within 8 percent through two-dimensional analysis (Direct Linear Software). Peak angular velocities displayed errors within 6 percent through three-dimensional analyses and exhibited errors of 12 percent when using two-dimensional analysis (Direct Linear Software). The specific results from this series of evaluations and their impacts on the methodology issues of kinematic data collection and analyses are presented in detail. The accuracy levels observed in these evaluations are also presented.

Wilmington, R. P.↗

Creep fatigue life prediction for engine hot section materials (isotropic)

A series of high temperature strain controlled fatigue tests have been completed to study the effects of thermomechanical fatigue, multiaxial loading, reactive environments, and imposed mean stresses. The baseline alloy used in these tests was cast B1900+Hf (with and without coatings); a small number of tests of wrought INCO 718 are also included. A strong path dependence was demonstrated during the thermomechanical fatigue testing, using in-phase, out-phase, and non-proportional (elliptical and 'dogleg') strain-temperature cycles. The multiaxial tests also demonstrated cycle path to be a significant variable, using both proportional and non-proportional tension-torsion loading. Environmental screening tests were conducted in moderate pressure oxygen and purified argon; the oxygen reduced the specimen lives by two, while the argon testing produced ambiguous data. Both NiCoCrAlY overlay and diffusion aluminide coatings were evaluated under isothermal and TMF conditions; in general, the lives of the coated specimens were higher that those of uncoated specimens. Controlled mean stress TMF tests showed that small mean stress changes could change initiation lives by orders of magnitude; these results are not conservatively predicted using traditional linear damage summation rules. Microstructures were evaluated using optical, SEM and TEM methods.

Nelson, R. S.↗

Correction for Neutral Pressure-Driven Signal in Radiated Power Measurements on Proto-MPEX

A pinhole camera using four-channel resistive bolometers was designed and implemented on the Prototype Material Plasma Exposure eXperiment (Proto-MPEX) to measure the plasma radiated power. Quantification of the radiated power is required for comprehensive power balance assessments but faces unique challenges on linear plasma devices due to weak absolute signal, compared to typical use in toroidal fusion plasmas, and high neutral pressure changes. Around 1 mTorr, the pressure-driven signal contribution was shown to overwhelm the radiated power signal in Proto-MPEX, found to be ~3 μW . Since the pressure-driven signals could not be predicted from first principle arguments, the signals were corrected using empirical models that take two experimental pressure measurements and the raw bolometer signal to reveal the radiated power signals. Four “gas-only” Proto-MPEX shots were taken, without the radio frequency (rf) plasma discharge, to train, optimize, and validate empirical models. The data included two pressure measurements at different locations along with the bolometer signal. Multiple models were tested; partial least-squares regression with a cross correlation method to account for the time offset in the pressure measurements provided the best combination of accuracy and stability. The model was validated to have a modeling error of 0.0036 V, which was lower than 10% of the expected radiated power signal. The model was tested on a different “gas-only” shot with a complete time-series data. For plasma operation shots, the corrected bolometer measurements were compared with an absolute extreme ultraviolet photodiode measurement and similarities and differences are discussed.

Lee, Seungsup↗

Reconstructing large-scale structure with neutral hydrogen surveys

Upcoming 21-cm intensity surveys will use the hyperfine transition in emission to map out neutral hydrogen in large volumes of the universe. Unfortunately, large spatial scales are completely contaminated with spectrally smooth astrophysical foregrounds which are orders of magnitude brighter than the signal. This contamination also leaks into smaller radial and angular modes to form a foreground wedge, further limiting the usefulness of 21-cm observations for different science cases, especially cross-correlations with tracers that have wide kernels in the radial direction. Here in this paper, we investigate reconstructing these modes within a forward modeling framework. Starting with an initial density field, a suitable bias parameterization and non-linear dynamics to model the observed 21-cm field, our reconstruction proceeds by {combining} the likelihood of a forward simulation to match the observations (under given modeling error and a data noise model) {with the Gaussian prior on initial conditions and maximizing the obtained posterior}. For redshifts z=2 and 4, we are able to reconstruct 21cm field with cross correlation, r c > 0.8 on all scales for both our optimistic and pessimistic assumptions about foreground contamination and for different levels of thermal noise. The performance deteriorates slightly at z=6. The large-scale line-of-sight modes are reconstructed almost perfectly. We demonstrate how our method also provides a technique for density field reconstruction for baryon acoustic oscillations, outperforming standard methods on all scales. We also describe how our reconstructed field can provide superb clustering redshift estimation at high redshifts, where it is otherwise extremely difficult to obtain dense spectroscopic samples, as well as open up a wealth of cross-correlation opportunities with projected fields (e.g. lensing) which are restricted to modes transverse to the line of sight.

21 cm↗

Wildfire management decisions outweigh mechanical treatment as the keystone to forest landscape adaptation

Modern land management faces unprecedented uncertainty regarding future climates, novel disturbance regimes, and unanticipated ecological feedbacks. Mitigating this uncertainty requires a cohesive landscape management strategy that utilizes multiple methods to optimize benefits while hedging risks amidst uncertain futures. We used a process-based landscape simulation model (LANDIS-II) to forecast forest management, growth, climate effects, and future wildfire dynamics, and we distilled results using a decision support tool allowing us to examine tradeoffs between alternative management strategies. We developed plausible future management scenarios based on factorial combinations of restoration-oriented thinning prescriptions, prescribed fire, and wildland fire use. Results were assessed continuously for a 100-year simulation period, which provided a unique assessment of tradeoffs and benefits among seven primary topics representing social, ecological, and economic aspects of resilience. Projected climatic changes had a substantial impact on modeled wildfire activity. In the Wildfire Only scenario (no treatments, but including active wildfire and climate change), we observed an upwards inflection point in area burned around mid-century (2060) that had detrimental impacts on total landscape carbon storage. While simulated mechanical treatments (~ 3% area per year) reduced the incidence of high-severity fire, it did not eliminate this inflection completely. Scenarios involving wildland fire use resulted in greater reductions in high-severity fire and a more linear trend in cumulative area burned. Mechanical treatments were beneficial for subtopics under the economic topic given their positive financial return on investment, while wildland fire use scenarios were better for ecological subtopics, primarily due to a greater reduction in high-severity fire. Benefits among the social subtopics were mixed, reflecting the inevitability of tradeoffs in landscapes that we rely on for diverse and countervailing ecosystem services. This study provides evidence that optimal future scenarios will involve a mix of active and passive management strategies, allowing different management tactics to coexist within and among ownerships classes. Our results also emphasize the importance of wildfire management decisions as central to building more robust and resilient future landscapes.

54 ENVIRONMENTAL SCIENCES↗

Elucidating the Mechanical Energy for Cyclization of a DNA Origami Tile

DNA origami has emerged as a versatile method to synthesize nanostructures with high precision. This bottom-up self-assembly approach can produce not only complex static architectures, but also dynamic reconfigurable structures with tunable properties. While DNA origami has been explored increasingly for diverse applications, such as biomedical and biophysical tools, related mechanics are also under active investigation. Here we studied the structural properties of DNA origami and investigated the energy needed to deform the DNA structures. We used a single-layer rectangular DNA origami tile as a model system and studied its cyclization process. This origami tile was designed with an inherent twist by placing crossovers every 16 base-pairs (bp), corresponding to a helical pitch of 10.67 bp/turn, which is slightly different from that of native B-form DNA (~10.5 bp/turn). We used molecular dynamics (MD) simulations based on a coarse-grained model on an open-source computational platform, oxDNA. We calculated the energies needed to overcome the initial curvature and induce mechanical deformation by applying linear spring forces. We found that the initial curvature may be overcome gradually during cyclization and a total of ~33.1 kcal/mol is required to complete the deformation. These results provide insights into the DNA origami mechanics and should be useful for diverse applications such as adaptive reconfiguration and energy absorption.

59 BASIC BIOLOGICAL SCIENCES↗

Validation & Verification of Electrical Components B2 Test Facility

The research focus of this project is to assist in the closure of Measurement, Monitoring and Control System (MMCS) and other electrical requirements in support of the B2 Space Launch System (SLS) Core Stage Green Run Test Project. Alongside this goal I am to understand project management tools to analyze and control activities/tasks associated with those items. The project I am working on will assist electrical engineers in developing and or identifying closure rationale for the MMCS and electrical requirements. Also as serving the role of a training project manager, I am to develop and maintain a method for tracking progress with estimate completion dates specifically identifying those items during this performance window. The methods I used in order to perform the research consisted of worddocuments, excel sheets, and pdf documents needed for review. For example, documents consisted of: a SLS Core Stage Green Run Facility Requirement Document (FRD) identifying the SSC MMCS and electrical requirements, a Stage Controller to SSC Interface Control Document (ICD) identifying the SSC requirements, and the B2 SLS Core Stage Green Run Test project System Requirements Document (SRD) which identifies the MMCS and electrical requirements. The files were sent to me by both mentors, Mr. Barry Robinson and Ms. Dawn Davis. Overall, the review process for the word documents and excel spreadsheets proved to be successful. The data for the electrical components were accurate and were consistent with the original recorded data. However, the consistency with the pdf documents did not follow up all the way. Attention to this mistake was made and further revisions were done in order for the data to agree with each other. Furthermore, an itinerary was designed using Microsoft Outlook in order to track progress with estimated completion dates. This project contributes to NASA/Center Missions and Goals through the Waterfall Model. The Waterfall Model is a linear system used for engineering design. In the Waterfall Model, the fourth step is Verification and Validation. This involves installation, testing, and debugging of the B2 SLS Core Stage electrical components. For the SLS Core Stage Green Run Test Project, each requirement contained in this approved requirements set will have at least one Verification Item (VI) assigned to it. During the design phase, VIs will mostly consist of Analysis or Inspection types. Design phase VIs and closure information will be documented in Dynamics Data Management System (DDMS) in the form of analysis reports and design documentation.

Bastian, Tyler↗

Epoxy Adhesives for Stator Magnet Assembly in Stirling Radioisotope Generators (SRG)

As NASA seeks to fulfill its goals of exploration and understanding through missions planned to visit the moons of Saturn and beyond, a number of challenges arise from the idea of deep space flight. One of the first problems associated with deep space travel is electrical power production for systems on the spacecraft. Conventional methods such as solar power are not practical because efficiency decreases substantially as the craft moves away from the Sun. The criterion for power generation during deep space missions are very specific, the main points requiring high reliability, low mass, minimal vibration and a long lifespan. A Stirling generator, although fairly old in concept, is considered to be a potential solution for electrical power generation for deep space flight. A Stirling generator works on the same electromagnetic principles of a standard generator, using the linear motion of the alternator through the stationary stator which produces electric induction. The motion of the alternator, however, is produced by the heating and cooling dynamics of pressurized gases. Essentially heating one end and cooling another of a contained gas will cause a periodic expansion and compression of the gas from one side to the other, which a displacer translates into linear mechanical motion. NASA needs to confirm that the materials used in the generator will be able to withstand the rigors of space and the life expectancy of the mission. I am working on the verification of the epoxy adhesives used to bond magnets to the steel lamination stack to complete the stator; in terms of in-service performance and durability under various space environments. Understanding the proper curing conditions, high temperature properties, and degassing problems as well as production difficulties are crucial to the long term success of the generator. system and steel substrate used in the stator. To optimize the curing conditions of the epoxies, modulated differential scanning calorimetry analysis was done as a function of cure time and temperatures. Adhesion bond strength was tested at various temperatures with lap shear samples using Hiperco 50 substrate to ensure that the proper adhesive is being used. To try and solve the problem of bondline thickness, micro glass beads of 0.0017" in diameter were investigated to see if any other physical properties of the epoxy were affected. Efforts will be made to develop a standard, optimized, fabrication process/procedure of sub-scale magnet-stator assemblies for various adhesive performance evaluation studies under simulated generator conditions. Also, accelerated aging testing will be done in a pressurized canister with stator assembly samples for three years to verify if any degassing or thermal degradation of the epoxy occurs. The necessity of verifying the correct epoxy adhesive system for the stator magnet in the SRG is crucial because failure of the stator assembly would jeopardize the electrical system, and thereby the entire mission itself. My work involves specimen fabrications, testing, and data analyses of the epoxy adhesive system for the stator magnet in the SRG is crucial because failure of the stator assembly would jeopardize the electrical system, and thereby the entire mission itself.

Cater, George M.↗

Upper bounds for convergence rates of vector extrapolation methods on linear systems with initial iterations

The application of the minimal polynomial extrapolation (MPE) and the reduced rank extrapolation (RRE) to a vector sequence obtained by the linear iterative technique x(sub j) + 1 = Ax(sub j) = b,j = 1,2,..., is considered. Both methods produce a two dimensional array of approximations s(sub n,k) to the solution of the system (I - A)x = b. Here, s(sub n,k) is obtained from the vectors x(sub j), n is less than or equal to j is less than or equal to n + k + 1. It was observed in an earlier publication by the first author that the sequence s(sub n,k), k = 1,2,..., for n greater than 0, but fixed, possesses better convergence properties than the sequence s(sub 0,k), k = 1,2,.... A detailed theoretical explanation for this phenomenon is provided in the present work. This explanation is heavily based on approximations by incomplete polynomials. It is demonstrated by numerical examples when the matrix A is sparse that cycling with s(sub n,k) for n greater than 0, but fixed, produces better convergence rates and costs less computationally than cycling with s(sub 0,k). It is also illustrated numerically with a convection-diffusion problem that the former may produce excellent results where the latter may fail completely. As has been shown in an earlier publication, the results produced by s(sub 0,k) are identical to the corresponding results obtained by applying the Arnoldi method or generalized minimal residual scheme (GMRES) to the system (I - A)x = b.

Sidi, Avram↗

Use of Minute-by-Minute Cardiovascular Measurements During Tilt Tests to Strengthen Inference on the Effect of Long-Duration Space Flight on Orthostatic Hypotension

Typical methodology for evaluating the effects of spaceflight on orthostatic hypotension (OH) has been survival analysis of tolerance times from 80 head-up tilt tests. However when scheduled test durations are short, there may not be enough failures to allow survival analysis to adequately estimate and compare the effects of flight phase (e.g. pre-flight, number of days post-flight), flight duration, and their interaction, as well as interactions with effects of interventions or countermeasures. The problem is exacerbated in the presence of a repeated measures design, in which subjects participate in tilt tests during various flight phases. Here we show how it is possible to dramatically improve the efficiency of statistical inference in this setting by making use of the additional information contained in minute-by-minute observations of cardiovascular parameters thought to be reflective of progression towards presyncope during tilt testing. Methods: We retrospectively examined operational tilt test (OTT; 10 -min 80 head-up tilt) data from 20 International Space Station (ISS) and 66 Shuttle astronauts 10 d before launch (L-10), on landing day (R+0) and during recovery (R+1, R+3, R+6-10) depending on the level of participation. Data from 5 ISS astronauts tested on R+0 or R+1 who used non-standard countermeasures were excluded. In addition to OTT survival time, 8 cardiovascular parameters (CP: heart rate, systolic, diastolic, and mean arterial blood pressure, pulse pressure, stroke volume, cardiac output, and total peripheral resistance) that might be predictive of progression towards presyncope were measured every minute of each OTT. Statistical analysis was predicated on a two ]stage model of causation. In the first stage, flight duration and time from landing affect the astronauts' degree of OH, which is manifested in the time trends and variation of the above CPs during OTTs. In the second stage, the behavior of these parameters directly affects OTT survival time. Actual analysis proceeded in the opposite direction. First we identified those CPs or linear combinations that best predicted OTT survival regardless of what spaceflight conditions led to OTT completion or presyncope. From these, we calculated a summary statistic (one per OTT) that best predicted survival. We then used mixed ]model regression analysis to relate changes in the summary statistic to flight phase and duration. Inference on the effects of phase, duration, and their interaction on OH follows directly from this second analysis. Results: A linear combination (W) of diastolic blood pressure (DBP) and stroke volume (SV) was found to be the best predictor of OTT survival using the complete data set of minute-by-minute observations of CPs for each OTT. Furthermore, the log-transformed standard deviation of W (Z = log SW) was found to be a strong predictor of survival in the reduced data set consisting of one observation per OTT. In other words, this measure of variability of W during an OTT was the best indicator of whether or not the subject could complete the 10-min test, with higher variability (i.e. higher values of Z) being associated with greater probability of failure. In the mixed-model regression analysis where Z was now treated as a outcome with flight phase and duration groups (ISS and STS) as predictors, we found that there was a significantly more variability in W (higher values of Z) for both groups at R+0, but with no evidence of an interaction until R+3, when the ISS group still had inflated variability, but not the STS group. Conclusions: Variability of the cardiovascular index W recovers more slowly after long-compared to short-duration spaceflight. Since high variability of W has also been shown to be predictive of OTT failure, a primary manifestation of OH, a logical conclusion is that recovery from OH also is slower after long-duration compared to short-duration spaceflights.

Feiveson, Alan H.↗

Dual-energy CT in the obese: a preliminary retrospective review to evaluate quality and feasibility of the single-source dual-detector implementation

Purpose: To determine the feasibility of performing dual-energy CT with a single-source spectral detector system in obese patients. Materials and methods: Retrospective, IRB-approved review of 28 patients weighing ≥ 270 lbs (122 kg) who underwent CT of the abdomen on a single-source spectral detector system was performed. Two blinded, independent radiologists rated relative preference between conventional CT images taken at 120 kVp (CCT120) and monoenergetic 70 keV equivalent (MonoE70) as well as iodine map image quality in the spleen, pancreas, kidneys, and liver. Signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) were compared between conventional CT and MonoE70 images and correlated with body habitus markers of weight, height, and abdominal diameter. Results: MonoE70 images were preferred by radiologists 100% of the time (1-sample t test, p < 0.0001) over conventional CCT120 images. Noise was significantly lower; SNR and CNR were significantly higher in MonoE70 images than in CCT120 images (paired t tests, p < 0.0001). Mean iodine map rating (scale 1–5) was 4.54 ± 0.58, denoting near homogenous and complete iodine mapping through the spleen, pancreas, kidneys, and liver for the majority of patients. Body habitus markers were not significantly correlated with image preference score; noise; MonoE70 SNR; MonoE70 CNR; change in noise, SNR, or CNR from CCT120 to MonoE70, or iodine map quality; ordinal and linear regression, p = 0.2547, p = 0.6837, p = 0.1888, p = 0.5489, p = 0.9830, p = 0.8849, p = 0.8741, p = 0.1522, respectively. Conclusion: The single-source spectral detector implementation of dual-energy CT provides viable, high-quality imaging for obese patients.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Predicted vibration responses of Apollo structure and effects of pressure correlation lengths on response

Progressive wove and reverberant acoustic fields ore often used in the laboratory for qualification testing of aerospace structures that are subjected to unsteady aerodynamic environments which occur during flight through the atmosphere. Acoustic simulation of these in-flight environments i s achieved i f the structural vibration response to acoustic excitation is equivalent, over the frequency range of interest, to the response caused by aerodynamic turbulence. To determine the degree of response simulation that could be expected for such qualification testing of the Apollo Spacecraft, responses of a segment of the Apollo structure were estimated theoretically for both acoustic and flight environments. The vibration analysis presented in this report was performed for the Spacecraft Lunar Module Adaptor; and the three environments treated are turbulence at Mach 1 and Mach 2, a reverberant acoustic field, and a modified acoustic progressive wave field. The progressive wave field is assumed to be generated within a specially designed shroud of sixteen axially oriented ducts which is constructed around and contoured to the external skin of the vehicle. Each duct is open on the side adjacent to the skin and is driven acoustically at one end by an independent noise source, thus allowing for different pressure correlation patterns around the circumference of the vehicle. The purpose of the analysis was to determine the effects of surface pressure correlation lengths on the response characteristics of representative Apollo structure, and to select on this basis, the optimum acoustic environment for structural vibration qualification tests of the Apollo in the Spacecraft Acoustic Laboratory at MSC. The method of analysis employed consists of determining, independently, the space average, mean-square acceleration spectral density for each of the significant modes of vibration of an equivalent cylindrical shell, and then expressing the total structural response as a linear summation of these modal responses. Approximately 570 classical modes of a pinned-end cylinder were used in order to find the spectral response characteristics of the structure throughout the frequency range of 10 - 1,000 Hz. The effective forcing functions or joint acceptances of the various modes of vibration were computed and are presented graphically for each of the pressure excitation environments. A brief development of the response equations used is presented in the report. Space average acceleration response spectra were computed for twelve different cases including aerodynamic turbulence at Mach 1 and Mach 2, reverberant acoustic field, five duct correlation patterns, two cases of an axially damped progressive wave field, and two cases in which the structural damping constant was varied from the expected value. The response spectra were computed and plotted by a high speed digital computer far a flat excitation pressure spectrum; and the response spectra were numerically integrated to give one-third octave bond average responses. A complete set of digital computer programs has been developed for both flat plate and cylindrical shell structures for the above environments, and each requires run times of three minutes or less. Thus, the work presented herein, along with the results, constitutes a practical advancement to the existing state-of-the-art of vibration prediction.

R W White↗

Weather Impact on Airport Arrival Meter Fix Throughput

Time-based flow management provides arrival aircraft schedules based on arrival airport conditions, airport capacity, required spacing, and weather conditions. In order to meet a scheduled time at which arrival aircraft can cross an airport arrival meter fix prior to entering the airport terminal airspace, air traffic controllers make regulations on air traffic. Severe weather may create an airport arrival bottleneck if one or more of airport arrival meter fixes are partially or completely blocked by the weather and the arrival demand has not been reduced accordingly. Under these conditions, aircraft are frequently being put in holding patterns until they can be rerouted. A model that predicts the weather impacted meter fix throughput may help air traffic controllers direct arrival flows into the airport more efficiently, minimizing arrival meter fix congestion. This paper presents an analysis of air traffic flows across arrival meter fixes at the Newark Liberty International Airport (EWR). Several scenarios of weather impacted EWR arrival fix flows are described. Furthermore, multiple linear regression and regression tree ensemble learning approaches for translating multiple sector Weather Impacted Traffic Indexes (WITI) to EWR arrival meter fix throughputs are examined. These weather translation models are developed and validated using the EWR arrival flight and weather data for the period of April-September in 2014. This study also compares the performance of the regression tree ensemble with traditional multiple linear regression models for estimating the weather impacted throughputs at each of the EWR arrival meter fixes. For all meter fixes investigated, the results from the regression tree ensemble weather translation models show a stronger correlation between model outputs and observed meter fix throughputs than that produced from multiple linear regression method.

Machine Learning Model.↗

Enhancement of EM Pump Performance through Modeling and Testing

Liquid sodium has seen development as a nuclear reactor coolant in the United States for over 50 years. While sodium has many desirable characteristics as a heat transfer fluid, its reactivity poses numerous engineering challenges for even the most basic plant components like pumps. Despite the challenges, mechanical pumps have seen success in sodium service albeit at the expense of simplicity and flexibility. In fact, significant failures of these pumps have been recorded which can be attributed to their complexity. Notwithstanding the complexity, mechanical pumps are still widely proposed for the main circuit pump of sodium cooled reactors. However, an alternative to mechanical pumps are Electromagnetic Induction Pumps (EMIPs) which use the principle of the Lorentz Force. EMIPs avoid mechanically moving parts in contact with sodium which simplifies the design by avoiding mechanical wear and seals. Therefore, these pumps are generally safer and more reliable than their mechanical counterparts. That being said, EMIPs suffer from relatively low efficiencies that range from 5% to 45%. Some of the inefficiencies are inherent to the EMIPs pumps; namely, inefficiencies arise in the form of resistive heating losses in the fluid itself. However, other inefficiencies arise from Edge Effects which include Finite Length and Finite Width effects. Due to the significant safety advantages EMIPs have compared to mechanical pumps, there is a strong motivation to integrate these technologies into the design of Generation IV liquid-metal-cooled nuclear power plants. However, EMIPs have low efficiencies partly attributed to Edge Effects. Much of the literature focusing on Edge Effects is dated or simply does not exist in the public domain. More recently, numerical models of EMIPs have shown promising results. However, these models are often custom solvers or use complex software suits which limit their application. This report describes the work completed to understand Edge Effects in two EMIPs designs. These designs are the Annular Linear Induction Pump (ALIP) and Permanent Magnet Induction Pump (PMIP) also known as a Moving Magnet Pump (MMP). The methodology developed during this project aim at measuring key performance parameters (pressure-flowrate and efficiency-slip performance) using high temperature liquid-sodium. Baseline experimental measurements were then used to verify the accuracy of numerical models developed in the software suit ANSYS/FLUENT. On the experimental side, instrumentation was developed to characterize the pressure-flowrate and efficiency-slip performance of both the ALIP and PMIP in high-temperature liquid-sodium. Additionally, instrumentation was developed to measure potential output pressure pulsations and the magnetic field distributions. Moreover, a specialized ALIP was constructed to study the impacts magnetic field tapering, also known as coil grading, on the Finite Length Effect. Lastly, experimental performance curves of the PMIP were collected at 300 oC and 600 oC and compared to analytical theory which accounted for Finite Width Effects. These results showed that the Finite Width Effect and non-uniform magnetic field in a PMIP impact the performance. On the numerical side, models were developed using the software suit FLUENT/ANSYS. These methods were initially compared against a simple MHD flow in a Permanent Magnet Electromagnet Flowmeter (PMEMFM). The results were significant in that they showed agreeance between the numerical model and demonstrated that PMEMFM’s could be calibrated numerically. These modeling methods were then applied to a prototypic PMIP. Results showed that this simple model provided reasonably accurate results of the PMIP p-Q performance.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Practical Aerodynamic Design Optimization Based on the Navier-Stokes Equations and a Discrete Adjoint Method

Compressible and incompressible versions of a three-dimensional unstructured mesh Reynolds-averaged Navier-Stokes flow solver have been differentiated and resulting derivatives have been verified by comparisons with finite differences and a complex-variable approach. In this implementation, the turbulence model is fully coupled with the flow equations in order to achieve this consistency. The accuracy demonstrated in the current work represents the first time that such an approach has been successfully implemented. The accuracy of a number of simplifying approximations to the linearizations of the residual have been examined. A first-order approximation to the dependent variables in both the adjoint and design equations has been investigated. The effects of a "frozen" eddy viscosity and the ramifications of neglecting some mesh sensitivity terms were also examined. It has been found that none of the approximations yielded derivatives of acceptable accuracy and were often of incorrect sign. However, numerical experiments indicate that an incomplete convergence of the adjoint system often yield sufficiently accurate derivatives, thereby significantly lowering the time required for computing sensitivity information. The convergence rate of the adjoint solver relative to the flow solver has been examined. Inviscid adjoint solutions typically require one to four times the cost of a flow solution, while for turbulent adjoint computations, this ratio can reach as high as eight to ten. Numerical experiments have shown that the adjoint solver can stall before converging the solution to machine accuracy, particularly for viscous cases. A possible remedy for this phenomenon would be to include the complete higher-order linearization in the preconditioning step, or to employ a simple form of mesh sequencing to obtain better approximations to the solution through the use of coarser meshes. An efficient surface parameterization based on a free-form deformation technique has been utilized and the resulting codes have been integrated with an optimization package. Lastly, sample optimizations have been shown for inviscid and turbulent flow over an ONERA M6 wing. Drag reductions have been demonstrated by reducing shock strengths across the span of the wing. In order for large scale optimization to become routine, the benefits of parallel architectures should be exploited. Although the flow solver has been parallelized using compiler directives. The parallel efficiency is under 50 percent. Clearly, parallel versions of the codes will have an immediate impact on the ability to design realistic configurations on fine meshes, and this effort is currently underway.

Grossman, Bernard↗

Applications of Data Assimilation to Analysis of the Ocean on Large Scales

It is commonplace to begin talks on this topic by noting that oceanographic data are too scarce and sparse to provide complete initial and boundary conditions for large-scale ocean models. Even considering the availability of remotely-sensed data such as radar altimetry from the TOPEX and ERS-1 satellites, a glance at a map of available subsurface data should convince most observers that this is still the case. Data are still too sparse for comprehensive treatment of interannual to interdecadal climate change through the use of models, since the new data sets have not been around for very long. In view of the dearth of data, we must note that the overall picture is changing rapidly. Recently, there have been a number of large scale ocean analysis and prediction efforts, some of which now run on an operational or at least quasi-operational basis, most notably the model based analyses of the tropical oceans. These programs are modeled on numerical weather prediction. Aside from the success of the global tide models, assimilation of data in the tropics, in support of prediction and analysis of seasonal to interannual climate change, is probably the area of large scale ocean modeling and data assimilation in which the most progress has been made. Climate change is a problem which is particularly suited to advanced data assimilation methods. Linear models are useful, and the linear theory can be exploited. For the most part, the data are sufficiently sparse that implementation of advanced methods is worthwhile. As an example of a large scale data assimilation experiment with a recent extensive data set, we present results of a tropical ocean experiment in which the Kalman filter was used to assimilate three years of altimetric data from Geosat into a coarsely resolved linearized long wave shallow water model. Since nonlinear processes dominate the local dynamic signal outside the tropics, subsurface dynamical quantities cannot be reliably inferred from surface height anomalies. Because of its potential for large scale synoptic coverage of the deep ocean, acoustic travel time data should be a natural complement to satellite altimetry. Satellite data give us vertical integrals associated with thermodynamic and dynamic processes.

Miller, Robert N.↗

Gravity Sensors and the Role of 3-D Visualization and Simulation in Biomedical Research

Ross 3.0 software was developed in the Biocomputation Center for semi-automated reconstruction of objects from serial thin sections. Data are captured directly from a transmission electron microscope via a video camera to a graphics workstation where the sections are mosaicked and contours are traced, registered and displayed by semi-automated methods. For the first time, macular type II cells are described completely for their innervation patterns. The purposes are to learn more about the fundamental circuitry of the macula and to demonstrate whether the terminals are altered morphologically by space flight. Current examples, from the medial part of the macula, are from maculas collected in-flight on the Space Life Sciences-2 mission, 4.5 hrs post-flight, and from a ground control. Results show that the typical type 11 cell receives processes from up to six nearby calyces or afferents. Nearly all the processes are elongated; some have bouton-like swellings and numerous vesicles. Multiple (2 to 4) processes from a single calyx to a type II cell are common, and approx. 1/3 of the processes innervate 2 type II cells of a neighboring cluster of 3 cells. About 2% of type II cells resemble type I cells morphologically and are surrounded by demicalyces. Differences in size or shape of the terminals under flight conditions could not be determined because the sample size is still too small; but it is clear that reconstruction methods provide insights into macular circuitry not obtainable by other techniques. The results demonstrate a morphological basis for interactions between adjacent receptive fields, through feedback-feedforward connections, during preprocessing of linear acceleratory information by the maculas. While the methods are currently being tested using vestibular maculas as the model system, it is clear that the technology is applicable to any tissue that can be physically or optically sectioned. ROSS software has already been implemented for reconstructing objects from tissues studied by confocal and by transmitted light microscopy, and research into magnetic resonance imaging-computational tomography combined visualization are underway.

Ross, Muriel D.↗