Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model based definitions”

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 937 records · Page 52

A GPU-based Approach for Turbomachinery Application

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Gas Turbine↗

Cirrus Parcel Model Comparison Project: The Critical Components to Simulate Cirrus Initiation Explicitly - Phase 1

The Cirrus Parcel Model Comparison Project, a project of the GCSS (GEWEX Cloud System Studies) Working Group on Cirrus Cloud Systems, involves the systematic comparison of current models of ice crystal nucleation and growth for specified, typical, cirrus cloud environments. In Phase I of the project reported here, simulated cirrus cloud microphysical properties are compared for situations of "warm" (40 C) and "cold" (-60 C) cirrus, both subject to updrafts of 4, 20 and 100 centimeters per second. Five models participated. The various models employ explicit microphysical schemes wherein the size distribution of each class of particles (aerosols and ice crystals) is resolved into bins or treated separately. Simulations are made including both the homogeneous and heterogeneous ice nucleation mechanisms. A single initial aerosol population of sulfuric acid particles is prescribed for all simulations. To isolate the treatment of the homogeneous freezing (of haze droplets) nucleation process, the heterogeneous nucleation mechanism is disabled for a second parallel set of simulations. Qualitative agreement is found for the homogeneous-nucleation- only simulations, e.g., the number density of nucleated ice crystals increases with the strength of the prescribed updraft. However, significant quantitative differences are found. Detailed analysis reveals that the homogeneous nucleation rate, haze particle solution concentration, and water vapor uptake rate by ice crystal growth (particularly as controlled by the deposition coefficient) are critical components that lead to differences in predicted microphysics. Systematic bias exists between results based on a modified classical theory approach and models using an effective freezing temperature approach to the treatment of nucleation. Each approach is constrained by critical freezing data from laboratory studies, but each includes assumptions that can only be justified by further laboratory research. Consequently, it is not yet clear if the two approaches can be made consistent. Large haze particles may deviate considerably from equilibrium size in moderate to strong updrafts (20-100 centimeters per second) at -60 C when the commonly invoked equilibrium assumption is lifted. The resulting difference in particle-size- dependent solution concentration of haze particles may significantly affect the ice particle formation rate during the initial nucleation interval. The uptake rate for water vapor excess by ice crystals is another key component regulating the total number of nucleated ice crystals. This rate, the product of particle number concentration and ice crystal diffusional growth rate, which is particularly sensitive to the deposition coefficient when ice particles are small, modulates the peak particle formation rate achieved in an air parcel and the duration of the active nucleation time period. The effects of heterogeneous nucleation are most pronounced in weak updraft situations. Vapor competition by the heterogeneously nucleated ice crystals may limit the achieved ice supersaturation and thus suppresses the contribution of homogeneous nucleation. Correspondingly, ice crystal number density is markedly reduced. Definitive laboratory and atmospheric benchmark data are needed for the heterogeneous nucleation process. Inter-model differences are correspondingly greater than in the case of the homogeneous nucleation process acting alone.

Lin, Ruei-Fong↗

Geophysical and Environmental Monitoring Data, and Subsurface Flow Modelling Results for Chicken Bone Meadow, Mt. Snodgrass, Crested Butte, CO

This dataset includes geoelectrical monitoring data acquired between October 2021 and November 2022, soil moisture and temperature data, groundwater data obtained from borehole SNIB covering the period from June 2021 to September 2022, and hydrological modelling results. The data were acquired to investigate how variations in bedrock type and topography, and vegetation cover control subsurface flow dynamics. To provide insights into the subsurface flow dynamics and their controls, a monitoring transect was installed at the Chicken Bone Meadow, Mt. Snodgrass, Crested Butte, CO, measuring the spatio-temporal variations of soil moisture, soil and snow temperature, subsurface electrical resistivity variations, and groundwater dynamics. Field data are organized in a folder structure, with Electrical Resistivity Tomography (ERT) data being provided as one file per measurement, and data of the soil moisture and temperature sensors being provided as text files covering the entire monitoring period. The ‘Locations.csv’ file contains the location of all sensors, given in NAD83 – UTM Zone 13N. ERT monitoring data has been processed to filter data based on reciprocal errors (data with errors > 30% were removed), a linear error model was fitted to each survey, and to ensure a constant set of measurements for time-lapse inversion, filtered data were interpolated and assigned a 100% measurement error. Soil moisture and temperature data were acquired at 15 min intervals, and averaged to provide 1h data. Weather data and borehole data (groundwater depth, conductivity and temperature) were acquired at 30 min intervals, and are provided as daily measurements; all measurements are averaged, except of precipitation values, which are given as daily accumulation. The hydrological model was set up along the ERT monitoring transect, and net infiltration was used as surface boundary condition and derived from the weather data. Four different results are provided, (1) results for a parameterization using hydraulic permeability and porosity as derived from the ERT data through petrophysical relationships, and (2) three simplified model results, using 1 to 3 geological layers above the bedrock. Modelling was performed using PFLOTRAN, and for each model the PFLOTRAN input files are provided. The result files include weekly hydrological modelling results (e.g., saturation, velocities, pressures), as well as the model parameterization. The dataset additionally includes a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; and a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type.

54 ENVIRONMENTAL SCIENCES↗

Coupling Satellite and Ground-Based Instruments to Map Climate Forcing by Anthropogenic Aerosols

Climate forcing by anthropogenic aerosols is a significant but highly uncertain factor in global climate change. Only satellites can offer the global coverage essential to reducing this uncertainty; however, satellite measurements must be coupled with correlative, in situ measurements both to constrain the aerosol optical properties required in satellite retrieval algorithms and to provide chemical identification of aerosol sources. This grant funded the first two years of a three-year project which seeks to develop methodologies for combining spaceborne lidar with in-situ aerosol data sets to improve estimates of direct aerosol climate forcing. Progress under this two-year grant consisted in the development and deployment of a new in-situ capability for measuring aerosol 180' backscatter and the extinction-to-backscatter ratio. This new measurement capacity allows definitive lidar/in-situ comparisons and improves our ability to interpret lidar data in terms of climatically relevant quantities such as the extinction coefficient and optical depth. Measurements were made along the coast of Washington State, in Central Illinois, over the Indian Ocean, and in the Central Pacific. Thus, this research, combined with previous measurements by others, is rapidly building toward a global data set of extinction-to-backscatter ratio for key aerosol types. Such information will be critical to interpreting lidar data from the upcoming PICASSO-CENA, or P-C, satellite mission. Another aspect of this project is to investigate innovative ways to couple the lidar-satellite signal with targeted in-situ measurements toward a direct determination of aerosol forcing. This aspect is progressing in collaboration with NASA Langley's P-C lidar simulator and radiative transfer modeling by the University of Lille, France.

Charlson, Robert J.↗

Toward Better Understanding of Microphysical Processes and Resulting Precipitation Physics: A Merger Of Observations and Cloud Models (Final Report)

A large database of global disdrometer observations and a diverse set of simulations from the Regional Atmospheric Modeling System (RAMS) were compared using principal component analysis (PCA) in order to better understand warm and ice-based precipitation processes. The analysis demonstrated that six distinct precipitation groups (PGs) with common characteristics were revealed in both the observations and model simulations. These PGs were defined by similar co-variability of rain parameters. However, the model showed large concentrations of drops with a given size compared to the observations. Detailed investigations determined the parameterization of drop breakup was forcing the drops to an equilibrium size with a stronger constraint than was indicated by the observations. A series of sensitivity studies permutating the rain shape parameter in the RAMS 2-moment bulk scheme demonstrated how the assumed shape of the rain distribution influences the microphysics as well as precipitation characteristics. Results were compared to the same simulations performed with a bin microphysics scheme where the rain shape parameter is freely evolving. In the bin scheme, a wide range of shape parameters were found, with horizontal and vertical variability, in addition to being a function of storm lifetime. The results of this study highlight the limitations of a fixed assumed rain distribution in 2-moment microphysics schemes. Model simulations were used to probe the microphysical origins of the six PGs. Rain budgets from a supercell case showed that several of the groups hypothesized to be associated with strong ice microphysical processes had significant contributions from melting hail, but the complexity of ice and warm-rain processes in the supercell precluded definitive determination of microphysical origins for many of the PGs. A case study from the Mid-latitude Continental Convection and Clouds Experiment (MC3E) compared the model PGs to disdrometer observations, and demonstrated the spatial and temporal variability of the PGs which is not possible from disdrometer point measurements alone. The analyses performed during this project demonstrated how statistical analysis can provide a robust method for comparing model simulations and observations, which ultimately resulted identification of limitations in some current microphysics parameterizations as well as insights into precipitation variability which ultimately benefits both observations and modeling efforts.

54 ENVIRONMENTAL SCIENCES↗

High dopant activation in arsenic doped single-crystal CdTe thin films: Insights from MBE growth and rapid thermal processing

Single-crystal model systems are valuable tools to investigate fundamental material properties. In this work, we use molecular beam epitaxy to deposit in situ arsenic (As) doped single-crystal CdTe films on large area Si substrates to better understand As doping for photovoltaic applications. We found that As incorporation is highly temperature dependent: a substrate temperature difference of 50 °C can lead to several orders of magnitude difference in As concentration. Cd overpressure during in situ doping may limit out-diffusion of As but decrease As incorporation, especially at lower growth temperatures. Carrier concentrations greater than 10 16 cm −3 can be achieved with or without Cd overpressure when annealed at temperatures above 500 °C. However, unlike the low (∼1% to 5%) dopant activation commonly observed in polycrystalline CdTe, our films achieve significantly higher activation ratios—exceeding 50%, and in some cases approaching 80%. These values are consistent with or exceed prior reports in single-crystal CdTe systems. In addition to as-deposited arsenic concentrations, we also consider arsenic distribution after different rapid thermal processing temperatures. We propose a detailed definition and description of how arsenic incorporation is considered and calculated. Due to carrier concentration saturation, As incorporation also needs to be controlled to average levels of 10 17 cm −3 to achieve high activation. These findings suggest that higher annealing temperature regimes may be beneficial to polycrystalline CdTe based PV devices.

36 MATERIALS SCIENCE↗

Stochastic Continuous-time Flexibility Scheduling and Pricing in Wholesale Electricity Markets

Large-scale integration of intermittent renewable energy sources (RES) is calling for additional flexibility resources as well as more advanced modeling and optimization techniques to account for the increasing uncertainty and variability in power systems operation. As the RES integration gains momentum, the magnitude and frequency of their variations increase, which may trigger ramping scarcity events in real-time power systems operation. This necessitates revisiting the present definition of power systems flexibility and reserve services to reflect their robustness and adequacy towards sub-interval variations of the load and RES, as well as adjusting the operation models to accommodate the new reserve services. This project took a fundamental approach and aimed at developing continuous-time scheduling and pricing model that accurately models the continuous-time variations of load and RES and efficiently deploys the ramping capability of flexible resources to compensate the sources of variability and uncertainty in the market. In this regard, this project pursued the following goals: Developing stochastic multi-fidelity continuous-time optimization models for scheduling of energy storage (ES) systems and flexible loads in wholesale energy markets; Developing the theory and practices of continuous-time locational marginal pricing for valuating energy storage systems and flexible loads in wholesale energy markets; Developing function space solution approach to convert the proposed stochastic multi-fidelity continuous-time optimization models into tractable mixed-integer linear optimization models; and Defining flexibility reserve as a new type of reserve in markets that would enable ultimate participation of energy storage devices in provision of services to compensate the variability and uncertainty of RES in electricity markets. This project successfully completed all five major tasks defined in the SOPO, and produced 8 high-impact journal papers, 6 conference papers, 3 published U.S. patents, and one web-based software for continuous-time operation optimization of power systems. The application of the proposed flexibility reserve and the stochastic multi-fidelity continuous-time operation scheduling models would modify the forward commitment and schedule of generating units, ES devices and flexible loads, and would line up the resources in such a way that the composition of available resources is better prepared to respond to the sub-hourly variations of the load and renewable resources in real-time operation. Therefore, this project paves the way to sustainable, reliable, and economic integration of renewable energy resources in power system, supporting the progress towards reaching the national targets on energy independence. Even if the proposed models offers a radically different point of view as compared to existing models, it does not alter fundamentally the architecture of power systems operations, nor the complexity of the scheduling problem, so the integration of this project in power systems is extremely practical.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Formal methods in the design of Ada 1995

Formal, mathematical methods are most useful when applied early in the design and implementation of a software system--that, at least, is the familiar refrain. I will report on a modest effort to apply formal methods at the earliest possible stage, namely, in the design of the Ada 95 programming language itself. This talk is an 'experience report' that provides brief case studies illustrating the kinds of problems we worked on, how we approached them, and the extent (if any) to which the results proved useful. It also derives some lessons and suggestions for those undertaking future projects of this kind. Ada 95 is the first revision of the standard for the Ada programming language. The revision began in 1988, when the Ada Joint Programming Office first asked the Ada Board to recommend a plan for revising the Ada standard. The first step in the revision was to solicit criticisms of Ada 83. A set of requirements for the new language standard, based on those criticisms, was published in 1990. A small design team, the Mapping Revision Team (MRT), became exclusively responsible for revising the language standard to satisfy those requirements. The MRT, from Intermetrics, is led by S. Tucker Taft. The work of the MRT was regularly subject to independent review and criticism by a committee of distinguished Reviewers and by several advisory teams--for example, the two User/Implementor teams, each consisting of an industrial user (attempting to make significant use of the new language on a realistic application) and a compiler vendor (undertaking, experimentally, to modify its current implementation in order to provide the necessary new features). One novel decision established the Language Precision Team (LPT), which investigated language proposals from a mathematical point of view. The LPT applied formal mathematical analysis to help improve the design of Ada 95 (e.g., by clarifying the language proposals) and to help promote its acceptance (e.g., by identifying a verifiable subset that would meet the needs of safety-critical applications). The first LPT project, which ran from the fall of 1990 unti the end of 1992, produced studies of several language issues: optimization, sharing and storage, tasking and protected records, overload resolution, the floating point model, distribution, program erros, and object-oriented programming. The second LPT project, in 1994, formally modeled the dynamic semantics of a large part of the (almost) final language definition, looking especially for interactions between language features.

Guaspari, David↗

From LDEF to a national Space Environment and Effects (SEE) program: A natural progression

As the LDEF program draws to a close, it leaves in place the fundamental building blocks for a Space Environment and Effects (SEE) program. Results from LDEF data analyses and investigations now form a substantial core of knowledge on the long term effects of the space environment on materials, system and structures. In addition, these investigations form the basic structure of a critically-needed SEE archive and database system. An agency-wide effort is required to capture all elements of a SEE program to provide a more comprehensive and focused approach to understanding the space environment, determining the best techniques for both flight and ground-based experimentation, updating the models which predict both the environments and those effects on subsystems and spacecraft, and, finally, ensuring that this multitudinous information is properly maintained, and inserted into spacecraft design programs. Many parts and pieces of a SEE program already exist at various locations to fulfill specific needs. The primary purpose of this program, under the direction of the Office of Advanced Concepts and Technology (OACT) in NASA Headquarters, is to take advantage of these parts; apply synergisms where possible; identify and when possible fill-in gaps; coordinate and advocate a comprehensive SEE program. The SEE program must coordinate and support the efforts of well-established technical communities wherein the bulk of the work will continue to be done. The SEE program will consist of a NASA-led SEE Steering Committee, consisting of government and industry users, with the responsibility for coordination between technology developers and NASA customers; and Technical Working Groups with primary responsibility for program technical content in response to user needs. The Technical Working Groups are as follows: Materials and Processes; Plasma and Fields; Ionizing Radiation; Meteoroids and Orbital Debris; Neutral External Contamination; Thermosphere, Thermal, and Solar Conditions; Electromagnetic Effects; Integrated Assessments and Databases. Specific technology development tasks will be solicited through a NASA Research Announcement to be released in May of 1994. The areas in which tasks are solicited include: (1) engineering environment definitions, (2) environments and effects design guidelines, (3) environments and effects assessment models and databases, and (4) flight/ground simulation/technology assessment data.

Bowles, David E.↗

Should Pruning be a Pre-Processor of any Linear System?

There are many real-world problems whose mathematical models turn out to be linear systems Ax = b, where A is an m x n matrix. Each equation of the linear system is an information. An information, in a physical problem, such as 4 mangoes, 6 bananas, and 5 oranges cost $10, is mathematically modeled as an equation 4x(sub 1) + 6x(sub 2) + 5x(sub 3) = 10 , where x(sub 1), x(sub 2), x(sub 3) are each cost of one mango, that of one banana, and that of one orange, respectively. All the information put together in a specified context, constitutes the physical problem and need not be all distinct. Some of these could be redundant, which cannot be readily identified by inspection. The resulting mathematical model will thus have equations corresponding to this redundant information and hence are linearly dependent and thus superfluous. Consequently, these equations once identified should be better pruned in the process of solving the system. The benefits are (i) less computation and hence less error and consequently a better quality of solution and (ii) reduced storage requirements. In literature, the pruning concept is not in vogue so far although it is most desirable. It is assumed that at least one information, i.e. one equation is known to be correct and which will be our first equation. In a numerical linear system, the system could be slightly inconsistent or inconsistent of varying degree. If the system is too inconsistent, then we should fall back on to the physical problem (PP), check the correctness of the PP derived from the material universe, modify it, if necessary, and then check the corresponding mathematical model (MM) and correct it. In nature/material universe, inconsistency is completely nonexistent. If the MM becomes inconsistent, it could be due to error introduced by the concerned measuring device and/or due to assumptions made on the PP to obtain an MM which is relatively easily solvable or simply due to human error. No measuring device can usually measure a quantity with an accuracy greater that 0.005% or, equivalently with a relative error less than 0.005%. Hence measurement error is unavoidable in a numerical linear system when the quantities are continuous (or even discrete with extremely large number). Assumptions, though not desirable, are usually made when we find the problem sufficiently difficult to be solved within the available means/tools/resources and hence distort the PP and the corresponding MM. The . error thus introduced in the system could (not always necessarily though) make the system somewhat inconsistent. If the inconsistency (contradiction) is too much then one should definitely not proceed to solve the system in terms of getting a least-squares solution or the minimum-norm least-squares solution. All these solutions will be invariably of no real-world use. If, on the other hand, inconsistency is reasonably low, i.e. the system is near-consistent or, equivalently, has near-linearly-dependent rows, then the foregoing solutions are useful. Pruning in such a near-consistent system should be performed based on the desired accuracy and on the definition of near-linear dependence. In this article, we discuss pruning over various kinds of linear systems and strongly suggest its use as a pre-processor or as a part of an algorithm. Ideally pruning should (i) be a part of the solution process (algorithm) of the system, (ii) reduce both computational error and complexity of the process, and (iii) take into account the numerical zero defined in the context. These are precisely what we achieve through our proposed O(mn2) algorithm presented in Matlab, that uses a subprogram of solving a single linear equation and that has embedded in it the pruning.

Sen, Syamal K.↗

Lagrangian and Eulerian drag models that are consistent between Euler-Lagrange and Euler-Euler (two-fluid) approaches for homogeneous systems

The undisturbed flow of a particle is of fundamental importance since it controls both the undisturbed flow force and the perturbation force (which includes quasisteady, added-mass, and history forces). Here we use the pairwise interaction extended point particle framework to evaluate the undisturbed flow of each particle through superposition of the perturbation flow induced by all its neighbors. This approach allows calculation of various statistics related to undisturbed fluid velocity under conditions of both stationary and nonstationary particles. In a random distribution of stationary particles, while the macroscale undisturbed flow is slowly varying, the microscale undisturbed flow that arises due to the perturbation flow of neighbors varies substantially from one particle to another and this in turn leads to large variation in the hydrodynamic force exerted on the particles. The effect of particle motion is generally to increase the particle-to-particle variation in the undisturbed fluid velocity of the particles. We observe that this increase is greater for the transverse component than for the streamwise component. As a result, with increasing random particle motion, the distribution of undisturbed fluid-velocity fluctuation becomes isotropic. Three different normalized forces are defined: $Φ^L$ is the Lagrangian normalized force on an individual particle suitable for application in a microscale-informed Euler-Lagrange simulation, $Φ^E$ is the Eulerian normalized average force suitable for application in an Euler-Euler simulation, and $Φ^{LE}$ is the Lagrangian normalized force on an individual particle suitable for application in the standard Euler-Lagrange simulation. In this work, we establish precise relations between these different definitions. The drag laws developed based on particle-resolved direct numerical simulation results and experiments are appropriate for application only as the Eulerian normalized average force. We introduce the force consistency relation and use it to obtain an expression for $Φ^L$, which when applied to each particle and averaged over all the particles equals $Φ^E$. The results are first obtained in the limit of stationary particles and then extended to the general case of nonstationary particles.

42 ENGINEERING↗

Low Stretch Diffusion Flames Over a Solid Fuel

A unique new way to study low gravity flames in normal gravity has been developed. To study flame structure and extinction characteristics in low stretch environments, a normal gravity low-stretch diffusion flame is generated using a cylindrical PMMA sample of varying large radii. Foutch and T'ien used the radiative loss as well as a densimetric Froude number to characterize the blowoff (small Da) and quenching extinction (large Da) boundaries in stagnation-point diffusion flames under various convective conditions. An important conclusion of this study was that the shape and location of the extinction boundary, as well as a number of important flame characteristics, were almost identical for the buoyant, forced, and mixed convective environments they modeled. This theory indicates it should be possible to understand a material's burning characteristics in the low stretch environment of spacecraft (induced by fans and crew movements) by understanding its burning characteristics in an equivalent Earth-based stretch environment (induced by normal gravity buoyancy). Similarly, the material's burning characteristics in Lunar or Martian stretch environments (induced by partial gravity buoyancy) can be assessed. Equivalent stretch rates can be determined as a function of gravity, imposed flow, and geometry. A generalized expression for stretch rate which captures mixed convection includes both buoyant and forced stretch is defined as a = a(sub f) ((1 + (a(sub b))exp 2/(a(sub b))exp 2))exp 1/2. For purely buoyant flow, the equivalent stretch rate is a(sub b) = [(rho(exp e)-rho(exp *)/rho(sub e)][g/R](exp 1/2). For purely forced flow, the equivalent stretch rate is characterized by either a(sub f)= 2U(sub infinity)/R for a cylinder, or a(sub f)=U(sub jet)/d(sub jet) for a jet impinging on a planar surface. In these experiments, the buoyant stretch is varied through R, the radius of curvature, but the buoyant stretch could also be varied through g, the gravity level. In this way the effect of partial gravity, such as those found on the Moon (1/6 g) or Mars (1/3 g) can be captured in the definition of flame stretch.

Olson, S. L.↗

LES-Based Analysis of Film Cooling in a High-Pressure Turbine Rotor Including Effects of Purge and Tip Flow

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Secondary Flows↗

Should Pruning be a Pre-Processor of any Linear System?

There are many real-world problems whose mathematical models turn out to be linear systems Ax = b , where A is an m by x n matrix. Each equation of the linear system is an information. An information, in a physical problem, such as 4 mangoes, 6 bananas, and 5 oranges cost $10, is mathematically modeled as 4x(sub 1) + 6x(sub 2) + 5x (sub 3) = 10, where x(sub 1), x(sub 2), x(sub 3) are each cost of one mango, that of one banana, and that of one orange, respectively. All the information put together in a specified context, constitutes the physical problem and need not be all distinct. Some of these could be redundant, which cannot be readily identified by inspection. The resulting mathematical model will thus have equations corresponding to this redundant information and hence are linearly dependent and thus superfluous. Consequently, these equations once identified should be better pruned in the process of solving the system. The benefits are (i) less computation and hence less error and consequently a better quality of solution and (ii) reduced storage requirements. In literature, the pruning concept is not in vogue so far although it is most desirable. In a numerical linear system, the system could be slightly inconsistent or inconsistent of varying degree. If the system is too inconsistent, then we should fall back on to the physical problem (PP), check the correctness of the PP derived from the material universe, modify it, if necessary, and then check the corresponding mathematical model (MM) and correct it. In nature/material universe, inconsistency is completely nonexistent. If the MM becomes inconsistent, it could be due to error introduced by the concerned measuring device and/or due to assumptions made on the PP to obtain an MM which is relatively easily solvable or simply due to human error. No measuring device can usually measure a quantity with an accuracy greater that 0.005% or, equivalently with a relative error less than 0.005%. Hence measurement error is unavoidable in a numerical linear system when the quantities are continuous (or even discrete with extremely large number). Assumptions, though not desirable, are usually made when we find the problem sufficiently difficult to be solved within the available means/tools/resources and hence distort the PP and the corresponding MM. The error thus introduced in the system could (not always necessarily though) make the system somewhat inconsistent. If the inconsistency (contradiction) is too much then one should definitely not proceed to solve the system in terms of getting a least-squares solution or a minimum norm solution or the minimum-norm least-squares solution. All these solutions will be invariably of no real-world use. If, on the other hand, inconsistency is reasonably low, i.e. the system is near-consistent or, equivalently, has near-linearly-dependent rows, then the foregoing solutions are useful. Pruning in such a near-consistent system should be performed based on the desired accuracy and on the definition of near-linear dependence. In this article, we discuss pruning over various kinds of linear systems and strongly suggest its use as a pre-processor or as a part of an algorithm. Ideally pruning should (i) be a part of the solution process (algorithm) of the system, (ii) reduce both computational error and complexity of the process, and (iii) take into account the numerical zero defined in the context. These are precisely what we achieve through our proposed O(mn2) algorithm presented in Matlab, that uses a subprogram of solving a single linear equation and that has embedded in it the pruning.

Sen, Syamal K.↗

Comparison of Artemis 2 and Artemis 5 Model Outcomes Using the Impact Probabilistic Risk Assessment Tool

BACKGROUND The Artemis campaign is a Moon exploration program with a series of six planned missions, five of which will be crewed. These five crewed missions will contain a single mission segment (space flight), or multiple mission segments involving space flight (Orion), lunar landing (LTV) and/or space habitat (Gateway). Each crewed segment faces the risk of unique medical conditions, necessitating medical sets/kits tailored to those specificities. To support and enable a data-driven and evidence-based decision-making process through out a mission’s life cycle, a software tool called IMPACT was developed. Using probabilistic risk assessment (PRA) methodologies, IMPACT (Informing Mission Planning via Analysis of Complex Tradespaces) is a novel tool built for analyzing the possibility of encountering complex medical risks during space flight, and for identifying the medical resources and capabilities needed to treat those potential at-risk medical conditions. This presentation will seek to compare IMPACT’s computational results upon potential complex space medical conditions (e.g., sprain/strain back or sleep disturbance) using IMPACT’s risk metrics and the associated optimized medical sets/kits between two Artemis missions: single segment Artemis 2 and multi-segmented Artemis 5. OVERVIEW By identifying potential medical conditions in space using input criteria such as crew quantity and composition, certain crew physical characteristics, mission duration and mission activities, IMPACT can produce analyses on the type of medical resources and capabilities needed to produce an optimized medical set/kit to address those medical conditions. IMPACT achieves this by performing hundreds of thousands of Monte Carlo simulations of missions to build aggregate pictures of medical risk. IMPACT’s risk metrics include loss of crew life (LOCL) – a measure of crew mortality due to medical conditions in space, return to definitive care (RTDC) – the need to perform crew evacuation, and task time lost (TTL) – a measure of the inability to perform activities due to crew disability. These risk metrics are applied to every medical condition identified by IMPACT’s computation analyses for every segment of the mission. Medical sets/kits are optimized to address these medical conditions but must fit within the stated Artemis Design Reference Mission (DRM) request for mass and volume physical size constraints. ANTICIPATED ANALYSIS AND CONCLUSION Using two Artemis missions, Artemis 2 and Artemis 5, IMPACT will provide the analyses for comparison of medical set/kit contents based upon mass and/or volume requirements and identify the at-risk medical conditions within both missions. This paper serves as an initial exploration of probabilistic risk assessment (PRA) medical risk calculations between two crewed Artemis missions and is not intended to be deemed the official medical response for the Artemis campaign.

probabilistic risk assessment↗

A Combined Observational and Modeling Approach to Study Modern Dust Transport from the Patagonia Desert to East Antarctica

The understanding of present atmospheric transport processes from Southern Hemisphere (SH) landmasses to Antarctica can improve the interpretation of stratigraphic data in Antarctic ice cores. In addition, long range transport can deliver key nutrients normally not available to marine ecosystems in the Southern Ocean and may trigger or enhance primary productivity. However, there is a dearth of observational based studies of dust transport in the SH. This work aims to improve current understanding of dust transport in the SH by showing a characterization of two dust events originating in the Patagonia desert (south end of South America). The approach is based on a combined and complementary use of satellite retrievals (detectors MISR, MODIS, GLAS ,POLDER, OMI,), transport model simulation (HYSPLIT) and surface observations near the sources and aerosol measurements in Antarctica (Neumayer and Concordia sites). Satellite imagery and visibility observations confirm dust emission in a stretch of dry lakes along the coast of the Tierra del Fuego (TdF) island (approx.54deg S) and from the shores of the Colihue Huapi lake in Central Patagonia (approx.46deg S) in February 2005. Model simulations initialized by these observations reproduce the timing of an observed increase in dust concentration at the Concordia Station and some of the observed increases in atmospheric aerosol absorption (here used as a dust proxy) in the Neumayer station. The TdF sources were the largest contributors of dust at both sites. The transit times from TdF to the Neumayer and Concordia sites are 6-7 and 9-10 days respectively. Lidar observations and model outputs coincide in placing most of the dust cloud in the boundary layer and suggest significant de- position over the ocean immediately downwind. Boundary layer dust was detected as far as 1800 km from the source and approx.800 km north of the South Georgia Island over the central sub-Antarctic Atlantic Ocean. Although the analysis suggests the presence of dust at approx.1500 km SW of South Africa five days after, the limited capabilities of existing satellite platforms to differentiate between aerosol types do not permit a definitive conclusion. In addition, the model simulations show dust lifting to the free troposphere as it travels south but it could not be confirmed by the satellite observations due to cloudiness. This work demonstrates that complementary information from existing transport models, satellite and surface data can yield a consistent picture of the dust transport from the Patagonia desert to Antarctica. It also illustrates the limitation of using any of these approaches individually to characterize the transport of dust in a heavily cloudy area.

Gasso, S.↗

2013-2014 Greater Fairbanks, Alaska, Transportation Survey

# 2013–2014 Greater Fairbanks, Alaska, Transportation Survey The 2013–2014 Greater Fairbanks Transportation Survey obtained behavior data for regional travel demand modeling. The planning region in Alaska comprised the North Star Borough—known as the PM2.5 nonattainment region—and included the cities of Fairbanks and North Pole. ## Data Collection Agency The Alaska Department of Transportation and Public Facilities conducted the survey. ## Methodology Data collection occurred in two phases: the first in fall 2013 and the second in winter 2014. The first phase employed address-based sampling to recruit more than 1,700 households for a one-day personal travel survey, and a sub-sample participated with global position system (GPS) and on-board diagnostic (OBD) loggers installed in their vehicles (282 vehicles) for one week. The purpose of phase one was to better understand the impact of vehicle emissions on air quality in the PM2.5 nonattainment region. Many of the households participating in the vehicle GPS/OBD portion of phase one were asked to participate in phase two. ## Drive Cycle Processing and Filtering NREL has developed a GPS data filtration routine to filter erroneous data points in individual drive cycles sourced from GPS devices mounted in vehicles. Second-by-second drive cycle data collected from GPS-instrumented vehicles during this survey have passed through NREL's drive cycle processing and filtering routines. ## Survey Records Survey records include 135 households. ## More Information For more information about the survey, see the [Greater Fairbanks Transportation Survey Final Report](https://www.nrel.gov/media/docs/libraries/tsdc/greater-fairbanks-transportation-survey-final-report.pdf?sfvrsn=16ecd45b_1). ## Transportation Data For details on available travel survey data and variable definitions, see the [data dictionary](https://www.nrel.gov/media/docs/libraries/tsdc/akdot_data_dictionary.pdf?sfvrsn=16778e11_1). NREL-generated drive cycle data are also available for this survey. For details on available data and variable definitions, see the [drive cycle data dictionary](https://www.nrel.gov/media/docs/libraries/tsdc/drive_cycles_data_dictionary.pdf?sfvrsn=7de7e888_1). Transportation data are available as zipped files. [Download Winzip](http://www.winzip.com/downwz.htm).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

2013-2014 Greater Fairbanks, Alaska, Transportation Survey

# 2013–2014 Greater Fairbanks, Alaska, Transportation Survey The 2013–2014 Greater Fairbanks Transportation Survey obtained behavior data for regional travel demand modeling. The planning region in Alaska comprised the North Star Borough—known as the PM2.5 nonattainment region—and included the cities of Fairbanks and North Pole. ## Data Collection Agency The Alaska Department of Transportation and Public Facilities conducted the survey. ## Methodology Data collection occurred in two phases: the first in fall 2013 and the second in winter 2014. The first phase employed address-based sampling to recruit more than 1,700 households for a one-day personal travel survey, and a sub-sample participated with global position system (GPS) and on-board diagnostic (OBD) loggers installed in their vehicles (282 vehicles) for one week. The purpose of phase one was to better understand the impact of vehicle emissions on air quality in the PM2.5 nonattainment region. Many of the households participating in the vehicle GPS/OBD portion of phase one were asked to participate in phase two. ## Drive Cycle Processing and Filtering NREL has developed a GPS data filtration routine to filter erroneous data points in individual drive cycles sourced from GPS devices mounted in vehicles. Second-by-second drive cycle data collected from GPS-instrumented vehicles during this survey have passed through NREL's drive cycle processing and filtering routines. ## Survey Records Survey records include 135 households. ## More Information For more information about the survey, see the [Greater Fairbanks Transportation Survey Final Report](https://www.nrel.gov/media/docs/libraries/tsdc/greater-fairbanks-transportation-survey-final-report.pdf?sfvrsn=16ecd45b_1). ## Transportation Data For details on available travel survey data and variable definitions, see the [data dictionary](https://www.nrel.gov/media/docs/libraries/tsdc/akdot_data_dictionary.pdf?sfvrsn=16778e11_1). NREL-generated drive cycle data are also available for this survey. For details on available data and variable definitions, see the [drive cycle data dictionary](https://www.nrel.gov/media/docs/libraries/tsdc/drive_cycles_data_dictionary.pdf?sfvrsn=7de7e888_1). Transportation data are available as zipped files. [Download Winzip](http://www.winzip.com/downwz.htm).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗