Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation, Modeling, Performance, and Visualization”

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.

748 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Pre-commercial Demonstration of Direct Non-iterative State Estimator for Operational Use of Synchrophasor Management Systems

The objective of this project was to demonstrate functionality and performance of a Direct Non-iterative State Estimator (DNSE) integrated with NYPA’s Energy Management System (EMS) and with enhanced Real Time Dynamics Monitoring System (RTDMS) synchrophasor platform from Electric Power Group (EPG). DNSE is designed to overcome a major obstacle to operational use of Synchro-Phasor Management Systems (SPMS) by providing to synchrophasor management systems (SPMS) applications a consistent and a complete synchrophasor data foundation in the same way that a traditional EMS State Estimator (SE) provides to EMS applications. Specifically, DNSE is designed to use synchrophasor measurements collected by a central PDC, Supervisory Control and Data Acquisition (SCADA) measurements, and Energy Management System (EMS) network model, to obtain the complete state of the utility’s operating model at rates that are close to the synchrophasor data rates. In this way, the system is comprehensive in that it does not only cover the part of the network that is “visible” via synchrophasors, but also the part that is only “visible” through the SCADA measurements. Visualization needs associated with the use of DNSE results are fulfilled through suitably enhanced Real Time Dynamics Monitoring System (RTDMS), with the enhancements implemented by EPG. This project had the following goals in mind: To advance the deployment of commercial grade DNSE software application that relies on synchrophasor and SCADA data ; Apply DNSE at other utilities, to address a generic and fundamental need for “clean” operational data for synchrophasor applications; Provide means for “live” estimated data access by control system operators; Enhance potential for situational awareness through full system operational model coverage; Sub-second execution rate of the Direct Non-iterative State Estimator, eventually at a near-phasor data rate execution speed, i.e. < 0.1 sec. Anticipated benefits from this projects are: Enhanced reliability and improvements in the economic efficiency of bulk power system planning and operations; Providing “clean” data to other synchrophasor applications; Enhancement of situational awareness by providing the full operational model updated at near synchrophasor rate; A production-grade software tool that incorporate synchrophasor and SCADA data; Provides a basis for development of next generation monitoring and control applications, based on both SCADA and PMU data. Quanta Technology (QT) team worked in collaboration with Electric Power Group (EPG) whose team has enhanced its commercial Real Time Dynamics Monitoring System (RTDMS) to accommodate the requirements posed by DNSE application. EPG also provided its ePDC and Model-less Data Conditioning (PDVC) software for integration with DNSE+. QT developed the system requirements for DNSE; developed system architecture and defined interfaces between internal DNSE components. The core DNSE algorithm with all surrounding interfaces was named DNSE+. Since the DNSE development was done in a simulated system environment, QT used its PMU simulator that was enhanced during this project for development and factory acceptance testing (FAT). SCADA data in this stage was simulated by commercial PSS/e software. The output of DNSE are estimates of System states in C37.118-2 format, sent to RTDMS for further processing and display. As the number of these states is large, it was necessary to expand the C37.111-2 standard to accommodate large data sets. This enhancement was implemented in RTDMS. The demonstration of pre-production DNSE technology was done at NYPA using streaming field data from NYPA PMUs and from its RTUs through their SCADA system. NYPA provided ICCP interface as well as Common Information Model (CIM). The relevance of the DNSE+ application is that it provides state estimation of the power system based on hybrid set of data, consisting of both available PMU data and SCADA measurements. As this is a direct, non-iterative method of calculation of the system states, if does not suffer from convergence issues which is potential problem for conventional state estimators. Also, it can take any available PMU measurements, so it does not need to have a high percentage of PMU coverage needed in the case of Linear State Estimator. As the DNSE calculates synchrophasors of the system states (both phase and absolute value) as sub-second rate, this application can provide a basis for development of next generation of applications based both on SCADA and PMU data.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

PipeSight: A High-Performance Computing Platform for Pipeline Integrity Management

The Phase I feasibility study completed as part of this project has led to a number of innovative technologies being developed and has laid the foundation for a successful Phase II effort to commercialize a platform for managing the integrity of pipelines for the damage mechanisms of the new, hybrid-energy based economy. To ground the development efforts and direction of the project, an extensive market research and customer discovery effort was undertaken early in Phase I. Through this effort, a number of pipeline owners and operators were interviewed, and the following key findings were discovered about the pipeline industry: • Small pipeline operators do not have the central engineering groups necessary to perform their own independent analysis of inspection data, but instead rely on summarized tally sheets provided to them by inspection service providers. • The time it takes to go from an inspection to a completed engineering assessment, even for small segments of pipeline, can take anywhere from 30-120 days. During this delay, critical threats can (and have been known to) cause failures. • Uncertainty is often not accounted for in the assessment of pipeline integrity. The tally sheets provided by third-party service providers are almost always deterministic in nature, identifying threats that present a concern only to the current (not the future) integrity of the pipeline. • It is uncommon to apply the latest technologies to perform advanced assessments of damaged pipelines. There is a desire to use more advanced analysis capabilities to assess threats. Many pipeline operators indicated that they would often excavate a pipeline to perform an inspection and find that the damage was not as bad as they anticipated, thus using limited resources unnecessarily. Companies are not consistent in their use of inspection data to determine corrosion rates, and those that do only calculate deterministic corrosion rates. • The industry has prominently relied on time-based inspections but has recently started to transition to risk-based inspections. However, there appears to be no uniform guidance on how to do so while properly accounting for all sources of uncertainty. • Companies are not storing inspection data in a manner that allows for the ready determination of temporal trends. • Predictive maintenance principles and practices are beginning to be used by early adopters • Some pipelines are being re-purposed to transport different process fluids than they were designed for, e.g., H 2 and CO 2 rich process streams to serve the new hybrid-energy based economy, which are presenting new integrity concerns for the existing pipeline network that crisscrosses the United States. As a result of these discoveries, we were able to target the development efforts in Phase I to best serve the needs of the industry. In Phase I, we developed a way to correlate multiple large-scale scans of the pipeline to determine a probabilistic corrosion rate that accounts for all sources of error and uncertainty in the inspection process. This probabilistic corrosion rate can be used to predict the future thickness distribution of the pipe wall. We demonstrate how this analysis may be performed in an analytical fashion and has been implemented in such a manner that it can be readily distributed using GPU computing through integration of the Kokkos programming model. We also make a very novel extension of the analytical corrosion rate model to Bayesian Networks (an explainable AI technique) that can account for non-parametric distributions of corrosion rates. With the predictions made above for the probabilistic corrosion rate and corresponding future distribution of the pipe wall thickness, we can assess the integrity of the pipeline through the use of a probabilistic engineering assessment. We developed a novel screening data analysis approach that can rapidly identify ‘hotspots’ (local thin areas) where the integrity of the pipeline is a concern. Once more, we implemented this screening approach in C++ to leverage GPU computing via the Kokkos programming model. After the critical hotspots are identified, we developed a program that can automatically generate an advanced finite element model of the damaged regions. Since the number of damaged regions that require advanced analysis can number in the thousands, we integrated an open-source container-native workflow engine for orchestrating parallel jobs on the cloud. Initially, these advanced numerical models were only designed to account for loading due to internal pressure. However, in a slight pivot from the initial Phase I proposal, we developed a complete pipe stress analysis program (called Simflex) which can simulate the complete pipeline and its response to thermal expansion, pressure, thermal bowing, weight, wind, earthquake, support displacement, support friction and external forces. This pipe stress analysis program was written generically, to handle any piping system, but contains the features needed to model long pipelines (i.e., it incorporates a model for soil mechanics and can account for the nonlinear boundary conditions necessary to simulate long underground pipelines). This pipe stress analysis program can simulate any segment of the pipeline (simple or complex) under any set of conditions and loads, to determine the supplemental loads (axial forces and bending moments) at the location of damage. This enables the most accurate state of stress to be accounted for in the pipeline, which can prove critical when evaluating the integrity of a damaged region. In the process of developing the technologies to perform the integrity assessment of the pipeline, we also extended one of the industry standard approaches for performing the assessment of local thin areas that extend more in the circumferential direction than the longitudinal direction of the pipeline. This approach was presented to the API 579-1/AS ME FFS-1 steering committee in November 2021 for consideration in the next edition of the industry standard for Fitness-For-Service (expected to be released in 2023). To help pipeline operators make decisions with the results on any integrity assessment, we developed a new approach to the life-cycle management of pipelines which uses a Bayesian Decision Network. The network is designed to help pipeline operators plan and prioritize inspection activities and ultimately make smarter, more cost-effective decisions. The Bayesian approach accounts for all sources of uncertainty and carries them through to the final optimal decisions, providing a probabilistic framework for optimizing inspection intervals. The proof-of-concept networks developed in the feasibility study are complete, verified, and are focused on a subset of the pipeline. To expand this novel approach to the scale necessary for an entire network of pipelines in Phase II, we will leverage the DOE-funded Bengi solver for industrial-scale decision making with Bayesian Networks [22]. Once implemented, we will be able to provide the pipeline industry with a much-needed tool for optimal inspection planning using truly explainable artificial intelligence (XAI). To handle all of these advanced capabilities into a cloud-based platform, the architecture of the Equity Engineering Cloud (EEC) was extended to include Argo Workflows, a framework capable of distributing and managing a massive number of jobs that consume their own resources, such that thousands of serial finite element simulations can be run in parallel. As part of this substantial undertaking, we also integrated Argo Continuous Delivery (CD) into the EEC, to aid with the rapid prototyping and iterations that will be imperative to the success of the PipeSight platform’s Agile development process in Phase II. As part of the pipe stress analysis program, we also developed a custom visualizer that leverages the DOE-funded VTK visualization library. We added custom contouring capabilities and a means for interacting visually with both the inputs and outputs of the pipe stress analysis program. We also developed routines for automating the post-processing of the finite element simulations to determine if any failure criteria are met and to visualize the deformations, stresses and strains in ParaView using the exodus II file format (a subset of netCDF).

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2019

First Annual Report on Pulsed Thermal Tomography Nondestructive Evaluation of Additively Manufactured Reactor Materials and Components: Pulsed Thermal Tomography Nondestructive Examination of Additively Manufactured Reactor Materials and Components

Additive manufacturing (AM, or 3D printing) for commercial nuclear energy applications is an emerging method for cost-efficient manufacturing aimed at replacing aging nuclear reactor parts and reducing costs for new construction. Because of the geometry of metallic structures of interest for nuclear applications, which consist of planar primitives with no symmetry of revolution, limited options are available for nondestructive evaluation (NDE), either during or post manufacturing. Known material flaws in AM include low-density regions consisting of non-sintered powder, which have to be detected to ensure the safety of long-term performance nuclear reactor components. As a solution to NDE of AM, we are developing pulsed thermal tomography (PTT) models and depth inversion algorithms for 3D imaging. PTT has many advantages because the method is non-contact and allows for in-service, NDE of AM nuclear reactor parts. By analyzing transients of surface temperature response due to internal thermal resistances, one can obtain 3D reconstructions of material effusivity using a unique inversion algorithm developed at Argonne. This report provides results of preliminary performance evaluation of PTT capability in detection of flaws in metallic structures, as well imaging complex geometry AM structures. Performance of PTT was investigated through modeling and experiments, Computer models of PTT were developed using COMSOL Heat Transfer module. Defects were modeled as cylindrical flat bottom holes (FBH), which is a common model of calibrated material flaws in thermal tomography experiments. Materials considered in this study include stainless steel 316 (SS316), stainless steel 304 (SS304), and Inconel 718. Theoretical analyses were conducted to validate inversion of simulated PTT data with COMSOL for a plate. Subsequently, 3D reconstructions were performed on COMSOL simulations for FBH, revealing a decrease in spatial resolution over depth due to thermal diffusion. The results of this study show that the performance of the inversion algorithm for detecting smaller defects depends strongly on the depth of the defect as well as the incident heat flux. The size of detectable defect was estimated by fitting a Gaussian function to surface temperature profile. The criteria for detectability was taken as 20mK noise equivalent temperature difference (NETD), which is currently the sensitivity limit of high-performance infrared cameras. It was determined through computer simulations that the smallest detectable FBH in SS316 has a 50µm diameter and is located 0.5mm below the plate surface. Preliminary experimental studies were performed to evaluate PTT performance in detecting calibrated flaws using metallic specimens with AM defects simulated as drilled-in FBH’s. High strength Stainless Steel 316 and Inconel 718 alloys were considered, as well as lower grade Stainless Steel 304, Nickel 200, and Hastelloy C276. Specimens investigated in this report consisted of approximately 1/4in-thick plates made out of these alloys using conventional manufacturing methods. The diameters of FBH’s varied from 1mm to 8mm, and their depths below the plate flat surface varied between 1mm and 6mm. The size of the smallest FBH was limited to 1mm because conventional mechanical drills were used for creating the holes. PTT imaging results have shown that 1mm-diameter FBH located 1mm and 2mm below the surface were detectable. Larger size FBH were detectable at greater depth. For example, 6mm-diameter FBH could be detected at 8mm depth. Image contrast varied slightly between the specimens, with the best reconstructions obtained in SS316 and C276 plates. Finally, PTT capability in imaging several AM structures fabricated from Inconel 718 (IN718) powder feedstock was investigated. The AM structures have complex geometry, but do not have calibrated internal defects. The objective is to determine PTT settings, such as the total integration time, for imaging of representative AM structures. A new FLIR X8500sc IR camera with higher spatial resolution (1280x1024 at 181Hz frame rate) was integrated into the PTT laboratory setup. Software tools were developed to ensure compatibility of new camera files with legacy effusivity reconstruction codes. The new camera was used to obtain 3-D reconstructions of AM IN718 nozzle plate and a particle filter plate. It was shown that PTT can scan through 2/3in-thick Inconel 718 plates in approximately 15s. Several modes of 3-D data visualization were explored, including using ImageJ and MATLAB software packages. The results demonstrate the capability of PTT in 3-D imaging of complex geometry AM structures. Intermediate results demonstrated in this report will be developed further to investigate detection of flaws in AM specimens.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Understanding sources of organic aerosol during CalNex-2010 using the CMAQ-VBS

Abstract. Community Multiscale Air Quality (CMAQ) model simulations utilizing the traditional organic aerosol (OA) treatment (CMAQ-AE6) and a volatility basis set (VBS) treatment for OA (CMAQ-VBS) were evaluated against measurements collected at routine monitoring networks (Chemical Speciation Network (CSN) and Interagency Monitoring of Protected Visual Environments (IMPROVE)) and those collected during the 2010 California at the Nexus of Air Quality and Climate Change (CalNex) field campaign to examine important sources of OA in southern California. Traditionally, CMAQ treats primary organic aerosol (POA) as nonvolatile and uses a two-product framework to represent secondary organic aerosol (SOA) formation. CMAQ-VBS instead treats POA as semivolatile and lumps OA using volatility bins spaced an order of magnitude apart. The CMAQ-VBS approach underpredicted organic carbon (OC) at IMPROVE and CSN sites to a greater degree than CMAQ-AE6 due to the semivolatile POA treatment. However, comparisons to aerosol mass spectrometer (AMS) measurements collected at Pasadena, CA, indicated that CMAQ-VBS better represented the diurnal profile and primary/secondary split of OA. CMAQ-VBS SOA underpredicted the average measured AMS oxygenated organic aerosol (OOA, a surrogate for SOA) concentration by a factor of 5.2, representing a considerable improvement to CMAQ-AE6 SOA predictions (factor of 24 lower than AMS). We use two new methods, one based on species ratios (SOA/ΔCO and SOA/Ox) and another on a simplified SOA parameterization, to apportion the SOA underprediction for CMAQ-VBS to slow photochemical oxidation (estimated as 1.5 × lower than observed at Pasadena using −log(NOx : NOy)), low intrinsic SOA formation efficiency (low by 1.6 to 2 × for Pasadena), and low emissions or excessive dispersion for the Pasadena site (estimated to be 1.6 to 2.3 × too low/excessive). The first and third factors are common to CMAQ-AE6, while the intrinsic SOA formation efficiency for that model is estimated to be too low by about 7 × . From source-apportioned model results, we found most of the CMAQ-VBS modeled POA at the Pasadena CalNex site was attributable to meat cooking emissions (48 %, consistent with a substantial fraction of cooking OA in the observations). This is compared to 18 % from gasoline vehicle emissions, 13 % from biomass burning (in the form of residential wood combustion), and 8 % from diesel vehicle emissions. All "other" inventoried emission sources (e.g., industrial, point, and area sources) comprised the final 13 %. The CMAQ-VBS semivolatile POA treatment underpredicted AMS hydrocarbon-like OA (HOA) + cooking-influenced OA (CIOA) at Pasadena by a factor of 1.8 compared to a factor of 1.4 overprediction of POA in CMAQ-AE6, but it did capture the AMS diurnal profile of HOA and CIOA well, with the exception of the midday peak. Overall, the CMAQ-VBS with its semivolatile treatment of POA, SOA from intermediate volatility organic compounds (IVOCs), and aging of SOA improves SOA model performance (though SOA formation efficiency is still 1.6–2 × too low). However, continued efforts are needed to better understand assumptions in the parameterization (e.g., SOA aging) and provide additional certainty to how best to apply existing emission inventories in a framework that treats POA as semivolatile, which currently degrades existing model performance at routine monitoring networks. The VBS and other approaches (e.g., AE6) require additional work to appropriately incorporate IVOC emissions and subsequent SOA formation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Leveraging 3D CAD and Virtual Reality in Design of the Calcine Disposition Project - 20512

Fluor Idaho's Calcine Disposition Project (CDP) is leveraging the power of three-dimensional (3D) visualization, light detection and ranging (Lidar), and virtual reality to improve designs and worker safety. With the use of computer-aided design (CAD) modeling software, the CDP has created 3D models of the facility and systems to aid in the design, development, and operation of its retrieval and transfer system. This by itself is nothing new and has been done in industry for years. The CDP project however is taking this a few steps farther with the use of Lidar and virtual reality software. The CRP is tasked with removing a radioactive granular material from stainless steel bins located in concrete vaults that were constructed in the late 1950's and early 1960's. While construction drawings are available, it is not certain that the drawings are as-built or how accurate they truly are. Part of the project requires precise placement of equipment on to bins around multiple pipe, electrical lines, and bin stiffeners. Relying solely on the existing drawings and models created from these drawings is highly risky and prone to failure. The first step to determine the precise location and orientation of all obstructions in the vault the project will use a Lidar system to scan the vault structure and precisely locate all vault and bin components. Once this information is obtained, the data will be integrated into the CAD model. The model will then be verified, and a true as-built model developed. Precise component location, size, and orientation will be used for final design and placement of retrieval components. The second step will be to take the verified 3D model and scan data and import it into virtual reality software creating a virtual world. While some work on the project requires all work to be completed remotely due to high radiation fields, some areas have low enough radiation levels that personnel can enter and perform work. With the use of a virtual reality headset, operators and technicians will be able to enter our virtual world and become familiar with the surroundings and perform work prior to ever entering the radiation environment. This training is expected to pay dividends by improving worker efficiency, reducing errors, and improving confidence that the work can be performed as expected. It will play a key aspect in as low as reasonably achievable (ALARA) principles. An additional aspect in using the Lidar scans and virtual environment, is that it will give engineers and managers an opportunity to revisit the facility, especially radiation areas, at any time to obtain design information and measurements that may be costly or impossible to obtain once the facility transitions into radiation operations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2014

Novel Supercritical Carbon Dioxide Power Cycle Utilizing Pressured Oxy-combustion in Conjunction with Cryogenic Compression

The team of Southwest Research Institute ® (SwRI) and Thar Energy LLC (Thar) applied technology engineering and economic analysis to evaluate two advanced oxy-combustion power cycles, the Cryogenic Pressurized Oxy-combustion Cycle (CPOC), and the Supercritical Oxy-combustion Cycle. This assessment evaluated the performance and economic cost of the two proposed cycles with carbon capture and included a technology gap analysis of the proposed technologies to determine the technology readiness level of the cycle and the cycle components. The results of the engineering and economic analysis and the technology gap analysis were used to identify the next steps along the technology development roadmap for the selected cycle. The project objectives, as outlined in the FOA, were 90% CO 2 removal at no more than a 35% increase in cost of electricity (COE) as compared to a Supercritical Pulverized Coal Plant without CO 2 capture. The supercritical oxy-combustion power cycle with 99% carbon capture achieves a COE of $\$$121/MWe. This revised COE represents a 21% reduction in cost as compared to supercritical steam with 90% carbon capture ($\$$137/MWe). However, this represents a 49% increase in the COE over supercritical steam without carbon capture ($\$$80.95/MWe), exceeding the 35% target. The supercritical oxy-combustion cycle with 99% carbon capture achieved a 37.9% HHV plant efficiency (39.3% LHV plant efficiency), when coupling a supercritical oxy-combustion thermal loop to an indirect supercritical CO 2 (sCO 2 ) power block. In this configuration, the power block achieved 48% thermal efficiency for turbine inlet conditions of 650°C and 290 atm. Power block efficiencies near 60% are feasible with higher turbine inlet temperatures, however a design tradeoff to limit firing temperature to 650°C was made in order to use austenitic stainless steels for the high temperature pressure vessels and piping and to minimize the need for advanced turbomachinery features such as blade cooling. The overall technical readiness of the supercritical oxy-combustion cycle is TRL 2, Technology Concept, due to the maturity level of the supercritical oxy-combustor for solid fuels, and several critical supporting components, as identified in the Technical Gap Analysis. The supercritical oxycombustor for solid fuels operating at pressures near 100 atm is a unique component of the supercritical oxy-combustion cycle. In addition to the low TRL supercritical oxy-combustor, secondary systems were identified that would require adaptation for use with the supercritical oxycombustion cycle. These secondary systems include the high pressure pulverized coal feed, high temperature cyclone, removal of post-combustion particulates from the high pressure cyclone underflow stream, and micro-channel heat exchangers tolerant of particulate loading. Bench scale testing was utilized to measure coal combustion properties at elevated pressures in a CO 2 environment. This testing included coal slurry preparation, visualization of coal injection into a high pressure fluid, and modification of existing test equipment to facilitate the combustion properties testing. Additional bench scale testing evaluated the effectiveness of a rotary atomizer for injecting a coal-water slurry into a fluid with similar densities, as opposed to the typical application where the high density fluid is injected into a low density fluid. The swirl type supercritical oxy-combustor was developed from initial concept to an advanced design stage through numerical simulation using FLUENT and Chemkin to model the flow through the combustor and provide initial assessment of the coal combustion reactions in the flow path. This effort enabled the initial combustor mechanical layout, initial pressure vessel design, and the conceptual layout of a pilot scale test loop. A pilot scale demonstration of the supercritical oxy-combustion cycle is proposed as the next step in the technology development. This demonstration would advance the supercritical oxy-combustion cycle and the supercritical oxy-combustor from a current TRL of 2, Technology Concept, to TRL 6, Pilot Scale System Demonstrated in a Relevant Environment, and enable the evaluation and continued refinement of the supercritical oxy-combustor and critical secondary systems.

42 ENGINEERING↗
NASA NTRS2015

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗
NASA NTRS2016

Forward Skirt Structural Testing on the Space Launch System (SLS) Program

Structural testing was performed to evaluate heritage forward skirts from the Space Shuttle program for use on the Space Launch System (SLS) program. One forward skirt is located in each solid rocket booster. Heritage forward skirts are aluminum 2219 welded structures. Loads are applied at the forward skirt thrust post and ball assembly. Testing was needed because SLS ascent loads are roughly 40% higher than Space Shuttle loads. Testing objectives were to determine margins of safety, demonstrate reliability, and validate analytical models. Two forward skirts were structurally tested using the test configuration. The test stand applied loads to the thrust post. Four hydraulic actuators were used to apply axial load and two hydraulic actuators were used to apply radial and tangential loads. The first test was referred to as FSTA-1 (Forward Skirt Structural Test Article) and was performed in April/May 2014. The purpose of FSTA-1 was to verify the ultimate capability of the forward skirt subjected to ascent ultimate loads. Testing consisted of two liftoff load cases taken to 100% limit load followed by an ascent load case taken to 110% limit load. The forward skirt was unloaded to no load after each test case. Lastly, the forward skirt was tested to 140% limit and then to failure using the ascent loads. The second test was referred to as FSTA-2 and performed in July/August of 2014. The purpose of FSTA-2 was to verify the ultimate capability of the forward skirt subjected to liftoff ultimate loads. Testing consisted of six liftoff load cases taken to 100% limit load followed by the six liftoff cases taken to 140% limit load. Two ascent load cases were then tested to 100% limit load. The forward skirt was unloaded to no load after each test case. Lastly, the forward skirt was tested to 140% limit and then to failure using the ascent loads. The forward skirts on FSTA-1 and FSTA-2 successfully carried all applied liftoff and ascent load cases. Both FSTA-1 and FSTA-2 were tested to failure by increasing the ascent loads. Failure occurred in the forward skirt thrust post radius. The forward skirts on FSTA-1 and FSTA-2 had nearly identical failure modes. FSTA-1 failed at 1.72 times limit load and FSTA-2 failed at 1.62 times limit load. This difference is primarily attributed to variation in material properties in the thrust post region. Test data were obtained from strain gages, deflection gages, ARAMIS digital strain measurement, acoustic emissions, and high-speed video. Strain gage data and ARAMIS strain were compared to finite element (FE) analysis predictions. Both the forward skirt and tooling were modeled. This allows the analysis to simulate the loading as close as possible to actual test configuration. FSTA-1 and FSTA-2 were instrumented with over 200 strain gages to ensure all possible failure modes could be captured. However, it turned out that three gages provided critical strain data. One was located in the post bore and two on the post radius. More gages were not specified due to space limitations and the desire to not interfere with the use of the ARAMIS system on the post radius. Measured strains were compared to analysis results for the load cycle to failure. Note that FSTA-1 gages were lost before failure was reached. FSTA-2 gages made it to the failure load but one of the radius gages was lost before testing began. This gage was not replaced because of the time and cost associated with disassembly of the test structure. Correlation to analysis was excellent for FSTA-1. FSTA-2 was not quite as good because there was more residual strain from previous load cycles. FSTA-2 was loaded and unloaded with 12 liftoff cases and two ascent cases before taking the skirt to failure. FSTA-1 only had two liftoff cases and one ascent case before taking the skirt to failure. The ARAMIS system was used to determine strain at the post radius by processing digital images of a speckled paint pattern. Digital cameras recorded images of the speckled paint pattern. ARAMIS strain results for FSTA-2 just prior to failure. Note a high strain location develops near the left side. This high strain compares well to analysis prediction for both FSTA-1 and FSTA-2. The strain at this location was also plotted versus limit load. Both FSTA-1 and FSTA-2 had excellent correlation between ARAMIS and analysis strains. Acoustic emission (AE) sensors were used to monitor for damage formation that may occur during testing (e.g., crack formation and growth or propagation). AE was very important because after disassembly of FSTA-1, a crack was observed in the ball fitting radius. The ball fitting did not crack on FSTA-2. AE data was used to reconstruct when the crack occurred. The AE energy versus time plot for FSTA. The energy increased considerably at 850 seconds (152% limit load), indicating a crack could have formed at this point. The only visual evidence found that could have corresponded to this was the crack that initiated in the ball fitting. The cracks in the forward skirt aluminum structures would likely have been lower energy due to a lower modulus and all that were found after failure correlated to occurring after the initial crack in the post radius. This was verified by high-speed cameras used to record the failure.

Lohrer, J. D.↗
DOE OSTI2016

Controlled release of cortisone drugs from block copolymers synthetized by ATRP

Diseases affecting posterior eye segment, like macular edema, infection and neovascularization, may cause visual impairment. Traditional treatments, such as steroidal-drugs intravitreal injections, involve chronic course of therapy usually over a period of years. Moreover, they can require frequent administrations of drug in order to have an adequately disease control. This dramatically reduce patient’s compliance. Efforts have been made to develop implantable devices that offer an alternative therapeutic approach to bypass many challenges of conventional type of therapy. Implantable drug delivery systems (DDS) have been developed to optimize therapeutic properties of drugs and ensure their slow release in the specific site. Polymeric materials can play an essential role in modulating drug delivery and their use in such field has become indispensable. During last decades, acrylic polymers have obtained growing interest. Biocompatibility and chemical properties make them extremely versatile, allowing their use in many field such as biomedical. In particular, block methacrylate copolymer with a balance of hydrophilic and hydrophobic properties can be suitable for prolonged DDS in biomedical devices. In this work, we focused on the realization of a system for controlled and long term release of betamethasone 17,21-dipropionate (BDP), a cortisone drug, from methacrylic block copolymers, to be tested in the treatment of the posterior eye’s diseases. Different series of methyl methacrylate/hydroxyethyl methacrylate (MMA/HEMA) block and random copolymers, with different monomer compositions (10–60% HEMA), were synthetized by Atom Transfer Radical Polymerization (ATRP) to find the best hydrophilic/hydrophobic ratio, able to ensure optimal kinetic release. Copolymer samples were characterized by NMR spectroscopy ({sup 1}H-NMR, {sup 13}C-NMR, CosY), SEC, TGA and DSC. Monitoring of drug release from films loaded with BDP was carried out by HPLC analysis. Evaluation of different kinetic models allowed to deduce that release of BDP is controlled over time from PMMA-b-PHEMA 53/47. In particular, PMMA-b-PHEMA 53/47 showed the best release profile to achieve the therapeutic reference dose of 3 µg/die, employed in treatment of posterior eye disease, up to four months. Accordingly, PMMA-b-PHEMA 53/47 has been tested to prepare ocular inserts. Ocular inserts with different shape and the same area of polymer films have been obtained using silicon moulds made by a 3D printer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Coastal Ohio Wind Project

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

17 WIND ENERGY↗