SEARCH · Engineering Papers

Results for “Boundary value problem”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,171 records · Page 65

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Airborne lidar and radiometric observations of PBL- and low clouds

Boundary layer- and low altitude clouds over open ocean and continent areas have been studied during several field campaigns since mid-1990 using the French airborne backscatter lidar LEANDRE in conjunction with on-board IR and visible radiometers. LEANDRE is an automatic system, and a modification of the instrumental parameters, when airborne, is computer controlled through an operator keyboard. The vertical range squared lidar signals and instrument status are displayed in real time on two dedicated monitors. The lidar is used either down- or up-looking while the aircraft is flying above or below clouds. A switching of the viewing configuration takes about a minute. The lidar measurements provide a high resolution description of cloud morphology and holes in cloud layers. The flights were conducted during various meteorological conditions on single or multilayer stratocumulus and cumulus decks. Analysis on a single shot basis of cloud top (or bottom) altitude and a plot of the corresponding histogram allows one to determine a probability density function (PDF). The preliminary results show the PDFs for cloud top are not Gaussian and symmetric about the mean value. The skewness varies with atmospheric conditions. An example of results recorded over the Atlantic ocean near Biarritz is displayed, showing: (1) the range squared lidar signals as a function of time (here 100 s corresponds to about 8 km, 60 shots are averaged on horizontal); the Planetary Boundary Layer (PBL) - up to 600 m - is observed at the beginning of the leg as well as on surface returns, giving an indication of the porosity; (2) the cloud top altitude variation between 2.4 to 2.8 km during the 150 to 320 s section; and (3) the corresponding PDF. Similar results are obtained on stratocumulus over land. Single shot measurements can be used also to determine an optical porosity at a small scale as well as a fractional cloudiness at a larger scale. A comparison of cloud top altitude retrieved from lidar and narrowbeam IR radiometer is conducted to study the scale integration problem. A good agreement within less than 100 m relies on spatial uniformity and an optically thick layer. In the presence of holes, a discrepancy is observed. This is illustrated in figure 2, displaying as a function of time (1) the lidar signals; (2) the target temperature (either clouds or sea surface) retreived from a narrowbeam IR radiometer, 17 C is the sea surface temperature on that day; and (3) the visible flux, linked to cloud albedo, measured by a pyranometer. In preparation of ASTEX, down- and up-looking measurements where conducted on stratocumulus clouds over the Atlantic Ocean near Quimper in Brittany. Depending on the flight pattern orientation with respect to the wind, the top and bottom cloud morphologies are different. Preliminary results are given on cloud morphology, cloud top PDFs, optical porosity, fractional cloudiness, and comparison of lidar and radiometric measurements.

Flamant, P. H.↗

NASA NTRS1997

Multidisciplinary Optimization of Tilt Rotor Blades Using Comprehensive Composite Modeling Technique

An optimization procedure is developed for addressing the design of composite tilt rotor blades. A comprehensive technique, based on a higher-order laminate theory, is developed for the analysis of the thick composite load-carrying sections, modeled as box beams, in the blade. The theory, which is based on a refined displacement field, is a three-dimensional model which approximates the elasticity solution so that the beam cross-sectional properties are not reduced to one-dimensional beam parameters. Both inplane and out-of-plane warping are included automatically in the formulation. The model can accurately capture the transverse shear stresses through the thickness of each wall while satisfying stress free boundary conditions on the inner and outer surfaces of the beam. The aerodynamic loads on the blade are calculated using the classical blade element momentum theory. Analytical expressions for the lift and drag are obtained based on the blade planform with corrections for the high lift capability of rotor blades. The aerodynamic analysis is coupled with the structural model to formulate the complete coupled equations of motion for aeroelastic analyses. Finally, a multidisciplinary optimization procedure is developed to improve the aerodynamic, structural and aeroelastic performance of the tilt rotor aircraft. The objective functions include the figure of merit in hover and the high speed cruise propulsive efficiency. Structural, aerodynamic and aeroelastic stability criteria are imposed as constraints on the problem. The Kreisselmeier-Steinhauser function is used to formulate the multiobjective function problem. The search direction is determined by the Broyden-Fletcher-Goldfarb-Shanno algorithm. The optimum results are compared with the baseline values and show significant improvements in the overall performance of the tilt rotor blade.

Chattopadhyay, Aditi↗

DOE OSTI2022

FROSch Preconditioners for Land Ice Simulations of Greenland and Antarctica

Numerical simulations of Greenland and Antarctic ice sheets involve the solution of large-scale highly nonlinear systems of equations on complex shallow geometries. This work is concerned with the construction of Schwarz preconditioners for the solution of the associated tangent problems, which are challenging for solvers mainly because of the strong anisotropy of the meshes and wildly changing boundary conditions that can lead to poorly constrained problems on large portions of the domain. In this study, two-level generalized Dryja--Smith--Widlund (GDSW)--type Schwarz preconditioners are applied to different land ice problems, i.e., a velocity problem, a temperature problem, as well as the coupling of the former two problems. We employ the message passing interface (MPI)--parallel implementation of multilevel Schwarz preconditioners provided by the package FROSch (fast and robust Schwarz) from the Trilinos library. The strength of the proposed preconditioner is that it yields out-of-the-box scalable and robust preconditioners for the single physics problems. To the best of our knowledge, this is the first time two-level Schwarz preconditioners have been applied to the ice sheet problem and a scalable preconditioner has been used for the coupled problem. The preconditioner for the coupled problem differs from previous monolithic GDSW preconditioners in the sense that decoupled extension operators are used to compute the values in the interior of the subdomains. Additionally, several approaches for improving the performance, such as reuse strategies and shared memory OpenMP parallelization, are explored. In our numerical study we target both uniform meshes of varying resolution for the Antarctic ice sheet as well as nonuniform meshes for the Greenland ice sheet. We present several weak and strong scaling studies confirming the robustness of the approach and the parallel scalability of the FROSch implementation. Among the highlights of the numerical results are a weak scaling study for up to 32,K processor cores (8,K MPI ranks and 4 OpenMP threads) and 566,M degrees of freedom for the velocity problem as well as a strong scaling study for up to 4,K processor cores (and MPI ranks) and 68,M degrees of freedom for the coupled problem.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1996

Heat and Momentum Transfer Studies in High Reynolds Number Wavy Films at Normal and Reduced Gravity Conditions

We examined the effect of the gas flow on the liquid film when the gas flows in the countercurrent direction in a vertical pipe at normal gravity conditions. The most dramatic effect of the simultaneous flow of gas and liquid in pipes is the greatly increased transport rates of heat, mass, and momentum. In practical situations this enhancement can be a benefit or it can result in serious operational problems. For example, gas-liquid flow always results in substantially higher pressure drop and this is usually undesirable. However, much higher heat transfer coefficients can be expected and this can obviously be of benefit for purposes of design. Unfortunately, designers know so little of the behavior of such two phase systems and as a result these advantages are not utilized. Due to the complexity of the second order boundary model as well as the fact that the pressure variation across the film is small compared to the imposed gas phase pressure, the countercurrent gas flow affect was studied for the standard boundary layer model. A different stream function that can compensate the shear stress affect was developed and this stream function also can predict periodic solutions. The discretized model equations were transformed to a traveling wave coordinate system. A stability analysis of these sets of equations showed the presence of a Hopf bifurcation for certain values of the traveling wave velocity and the shear stress. The Hopf celerity was increased due to the countercurrent shear. For low flow rate the increases of celerity are more than for the high flow rate, which was also observed in experiments. Numerical integration of a traveling wave simplification of the model also predicts the existence of chaotic large amplitude, nonperiodic waves as observed in the experiments. The film thickness was increased by the shear.

Balakotaiah, V.↗

DOE OSTI2023

Sensitivity analysis of a layered piezoelectric system using ZFEM

The complex variable finite element method (ZFEM) is a numerical technique which aims to find the partial derivatives of the independent variables with respect to variation in dependent parameters declared in the physics. This is done by combining the complex Taylor series expansion within the weak formulation of the governing equation in a coupled system of linear equations forming a complex valued block matrix given by the Cauchy–Riemann matrix representation. In this work, two-dimensional linear first-order elements have been implemented in ZFEM to predict the design derivatives of the mechanical displacement field and the voltage potential field for a layered piezoelectric system in a steady-state study with Dirichlet boundary condition applied at the top and bottom edges of the geometry. This approach allows the standard FEM solution to quantify the sensitivity of the mechanical displacement and voltage potential fields with respect to small variations in the material properties through the information obtained from the computation of the derivatives. The domain is formed by a layered body with PZT-4 and PZT-5 stacked together. For result verification, the numerical solution obtained with ZFEM was compared to results from a commercial FEM package and the solution from the imaginary part was compared to the exact solution of a well-known benchmark problem. In conclusion, comparison of the results showed good agreement for both the real and imaginary parts of the solution and the largest sensitivities were found in PZT-5 specifically in C 13 , C 33 , and ε 33 .

42 ENGINEERING↗

NASA NTRS2016

Avoiding Stair-Step Artifacts in Image Registration for GOES-R Navigation and Registration Assessment

In developing software for independent verification and validation (IVV) of the Image Navigation and Registration (INR) capability for the Geostationary Operational Environmental Satellite R Series (GOES-R) Advanced Baseline Imager (ABI), we have encountered an image registration artifact which limits the accuracy of image offset estimation at the subpixel scale using image correlation. Where the two images to be registered have the same pixel size, subpixel image registration preferentially selects registration values where the image pixel boundaries are close to lined up. Because of the shape of a curve plotting input displacement to estimated offset, we call this a stair-step artifact. When one image is at a higher resolution than the other, the stair-step artifact is minimized by correlating at the higher resolution. For validating ABI image navigation, GOES-R images are correlated with Landsat-based ground truth maps. To create the ground truth map, the Landsat image is first transformed to the perspective seen from the GOES-R satellite, and then is scaled to an appropriate pixel size. Minimizing processing time motivates choosing the map pixels to be the same size as the GOES-R pixels. At this pixel size image processing of the shift estimate is efficient, but the stair-step artifact is present. If the map pixel is very small, stair-step is not a problem, but image correlation is computation-intensive. This paper describes simulation-based selection of the scale for truth maps for registering GOES-R ABI images.

stair-step artifact↗

NASA NTRS2016

Avoiding Stair-Step Artifacts in Image Registration for GOES-R Navigation and Registration Assessment

image registration↗

DOE OSTI2021

Deep Learning Approaches to Surrogates for Solving the Diffusion Equation for Mechanistic Real-World Simulations

In many mechanistic medical, biological, physical, and engineered spatiotemporal dynamic models the numerical solution of partial differential equations (PDEs), especially for diffusion, fluid flow and mechanical relaxation, can make simulations impractically slow. Biological models of tissues and organs often require the simultaneous calculation of the spatial variation of concentration of dozens of diffusing chemical species. One clinical example where rapid calculation of a diffusing field is of use is the estimation of oxygen gradients in the retina, based on imaging of the retinal vasculature, to guide surgical interventions in diabetic retinopathy. Furthermore, the ability to predict blood perfusion and oxygenation may one day guide clinical interventions in diverse settings, i.e., from stent placement in treating heart disease to BOLD fMRI interpretation in evaluating cognitive function (Xie et al., 2019; Lee et al., 2020). Since the quasi-steady-state solutions required for fast-diffusing chemical species like oxygen are particularly computationally costly, we consider the use of a neural network to provide an approximate solution to the steady-state diffusion equation. Machine learning surrogates, neural networks trained to provide approximate solutions to such complicated numerical problems, can often provide speed-ups of several orders of magnitude compared to direct calculation. Surrogates of PDEs could enable use of larger and more detailed models than are possible with direct calculation and can make including such simulations in real-time or near-real time workflows practical. Creating a surrogate requires running the direct calculation tens of thousands of times to generate training data and then training the neural network, both of which are computationally expensive. Often the practical applications of such models require thousands to millions of replica simulations, for example for parameter identification and uncertainty quantification, each of which gains speed from surrogate use and rapidly recovers the up-front costs of surrogate generation. We use a Convolutional Neural Network to approximate the stationary solution to the diffusion equation in the case of two equal-diameter, circular, constant-value sources located at random positions in a two-dimensional square domain with absorbing boundary conditions. Such a configuration caricatures the chemical concentration field of a fast-diffusing species like oxygen in a tissue with two parallel blood vessels in a cross section perpendicular to the two blood vessels. To improve convergence during training, we apply a training approach that uses roll-back to reject stochastic changes to the network that increase the loss function. The trained neural network approximation is about 1000 times faster than the direct calculation for individual replicas. Because different applications will have different criteria for acceptable approximation accuracy, we discuss a variety of loss functions and accuracy estimators that can help select the best network for a particular application. We briefly discuss some of the issues we encountered with overfitting, mismapping of the field values and the geometrical conditions that lead to large absolute and relative errors in the approximate solution.

60 APPLIED LIFE SCIENCES↗

NASA NTRS2008

Automated Generation and Assessment of Autonomous Systems Test Cases

This slide presentation reviews some of the issues concerning verification and validation testing of autonomous spacecraft routinely culminates in the exploration of anomalous or faulted mission-like scenarios using the work involved during the Dawn mission's tests as examples. Prioritizing which scenarios to develop usually comes down to focusing on the most vulnerable areas and ensuring the best return on investment of test time. Rules-of-thumb strategies often come into play, such as injecting applicable anomalies prior to, during, and after system state changes; or, creating cases that ensure good safety-net algorithm coverage. Although experience and judgment in test selection can lead to high levels of confidence about the majority of a system's autonomy, it's likely that important test cases are overlooked. One method to fill in potential test coverage gaps is to automatically generate and execute test cases using algorithms that ensure desirable properties about the coverage. For example, generate cases for all possible fault monitors, and across all state change boundaries. Of course, the scope of coverage is determined by the test environment capabilities, where a faster-than-real-time, high-fidelity, software-only simulation would allow the broadest coverage. Even real-time systems that can be replicated and run in parallel, and that have reliable set-up and operations features provide an excellent resource for automated testing. Making detailed predictions for the outcome of such tests can be difficult, and when algorithmic means are employed to produce hundreds or even thousands of cases, generating predicts individually is impractical, and generating predicts with tools requires executable models of the design and environment that themselves require a complete test program. Therefore, evaluating the results of large number of mission scenario tests poses special challenges. A good approach to address this problem is to automatically score the results based on a range of metrics. Although the specific means of scoring depends highly on the application, the use of formal scoring - metrics has high value in identifying and prioritizing anomalies, and in presenting an overall picture of the state of the test program. In this paper we present a case study based on automatic generation and assessment of faulted test runs for the Dawn mission, and discuss its role in optimizing the allocation of resources for completing the test program.

Testing challenges↗

NASA NTRS1994

Effect of nonzero surface admittance on receptivity and stability of compressible boundary layer

The effect of small-amplitude short-scale variations in surface admittance on the acoustic receptivity and stability of two-dimensional compressible boundary layers is examined. In the linearized limit, the two problems are shown to be related both physically and mathematically. This connection between the two problems is used, in conjunction with some previously reported receptivity results, to infer the modification of stability properties due to surface permeability. Numerical calculations are carried out for a self-similar flat-plate boundary layer at subsonic and low supersonic speeds. Variations in mean suction velocity at the perforated admittance surface can also induce receptivity to an acoustic wave. For a subsonic boundary layer, the dependence of admittance-induced receptivity on the acoustic-wave orientation is significantly different from that of the receptivity produced via mean suction variation. The admittance-induced receptivity is generally independent of the angle of acoustic incidence, except in a relatively narrow range of upstream-traveling waves for which the receptivity becomes weaker. However, this range of angles is precisely that for which the suction-induced receptivity tends to be large. At supersonic Mach numbers, the admittance-induced receptivity to slow acoustic models is relatively weaker than that in the case of the fast acoustic modes. We also find that purely real values for the surface admittance tend to have a destabilizing effect on the evolution of an instability wave over a slightly permeable surface. The limits on the validity of the linearized approximation are also assessed in one specific case.

Choudhari, Meelan↗

NASA NTRS2001

Hydrodynamic Stability of Liquid-Propellant Combustion: Landau's Problem Revisited

Hydrodynamic, or Landau, instability in combustion is typically associated with the onset of wrinkling of a flame surface, corresponding to the formation of steady cellular structures as the stability threshold is crossed. As its name suggests, it stems from hydrodynamic effects connected with thermal expansion across the reaction region. In the context of liquid-propellant combustion, the classical models that originally predicted this phenomenon have been extended to include the important effects that arise from a dynamic dependence of the burning rate on the local pressure and temperature fields. Thus, the onset of Landau instability has now been shown to occur for sufficiently small negative values of the pressure sensitivity of the burning rate, significantly generalizing previous classical results for this problem that assumed a constant normal burning rate. It has also been shown that the onset of instability occurs for decreasing values of the disturbance wave number as the gravitational-acceleration parameter decreases. Consequently, in an appropriate weak-gravity limit, Landau instability becomes a long-wave phenomena associated with the formation of large cells on the liquid-propellant surface. Additionally, a pulsating form of hydrodynamic instability has been shown to occur as well, corresponding to the onset of temporal oscillations in the location of the liquid/gas interface. This instability occurs for sufficiently large negative values of the pressure sensitivity, and is enhanced by increasing values of the burning-rate temperature sensitivity. It is further shown that for sufficiently small values of this parameter, there exists a stable range of pressure sensitivities for steady, planar burning such that the classical cellular form of hydrodynamic instability and the more recent pulsating form of hydrodynamic instability can each occur as the corresponding stability threshold is crossed. For larger thermal sensitivities, however, the pulsating stability boundary evolves into a C-shaped curve in the (disturbance-wave number, pressure-sensitivity) plane, indicating loss of stability to pulsating perturbations for all sufficiently large disturbance wavelengths. It is thus concluded, based on characteristic parameter values that an equally likely form of hydrodynamic instability in liquid-propellant combustion is of a nonsteady, long-wave nature, distinct from the steady, cellular form originally predicted by Landau.

Margolis, S. B.↗

NASA NTRSDate not supplied

An Approach for Practical Grid-Resolved Roughness Aerothermodynamic Simulations for Woven Heatshield Surfaces

An approach for simulating the aerothermal environment over a patterned roughness surface within an otherwise axisymmetric flow is developed. The patterned roughness is simplified into a sinusoidal shape that captures the basic features of a charred 3MDCP ablator. Two orientations of this pattern, which are equally flight relevant due to the forming process, enable periodic boundary conditions along the edges of a grid spanning a single roughness element. This simplified three-dimensional grid represents the minimum problem size for a patterned roughness simulation, which enables computationally efficient grid-resolved roughness simulations. This approach is validated by simulating the Langley Mach 6 measurements made on a similar sinusoidal surface, which results in convective heating within 5% of the experimental data. Applying this approach to Mars Sample Return (MSR) Earth Entry System (EES) flight cases, considering simulations with both 11 species air and 30 species air with ablation products, results in heating augmentation that is slightly lower than the widely-used Dahm correlation approach, with the difference being dependent on pattern orientation. This result provides the only insight into the behavior of the roughness augmentation at the high temperature reacting flow conditions present for EES flight cases, which are not captured by available ground test measurements. For shear augmentation, the simulations revealed that the surface-parallel pressure component is dominant. For arc-jet tests targeted to flight values for shear, the higher arc-jet pressures required to match flight values for shear result in up to a 400% increase in the rough wall shear value, which is not captured by the Dahm correlation. This indicates that, to avoid significant over-testing, grid-resolved roughness simulations are required to determine arc-jet conditions that produce flight shear levels in rough-wall scenarios.

Christopher O Johnston↗

NASA NTRS1997

SYNTHESIS OF NOVEL ALL-DIELECTRIC GRATING FILTERS USING GENETIC ALGORITHMS

We are concerned with the design of inhomogeneous, all dielectric (lossless) periodic structures which act as filters. Dielectric filters made as stacks of inhomogeneous gratings and layers of materials are being used in optical technology, but are not common at microwave frequencies. The problem is then finding the periodic cell's geometric configuration and permittivity values which correspond to a specified reflectivity/transmittivity response as a function of frequency/illumination angle. This type of design can be thought of as an inverse-source problem, since it entails finding a distribution of sources which produce fields (or quantities derived from them) of given characteristics. Electromagnetic sources (electric and magnetic current densities) in a volume are related to the outside fields by a well known linear integral equation. Additionally, the sources are related to the fields inside the volume by a constitutive equation, involving the material properties. Then, the relationship linking the fields outside the source region to those inside is non-linear, in terms of material properties such as permittivity, permeability and conductivity. The solution of the non-linear inverse problem is cast here as a combination of two linear steps, by explicitly introducing the electromagnetic sources in the computational volume as a set of unknowns in addition to the material unknowns. This allows to solve for material parameters and related electric fields in the source volume which are consistent with Maxwell's equations. Solutions are obtained iteratively by decoupling the two steps. First, we invert for the permittivity only in the minimization of a cost function and second, given the materials, we find the corresponding electric fields through direct solution of the integral equation in the source volume. The sources thus computed are used to generate the far fields and the synthesized triter response. The cost function is obtained by calculating the deviation between the synthesized value of reflectivity/transmittivity and the desired one. Solution geometries for the periodic cell are sought as gratings (ensembles of columns of different heights and widths), or combinations of homogeneous layers of different dielectric materials and gratings. Hence the explicit unknowns of the inversion step are the material permittivities and the relative boundaries separating homogeneous parcels of the periodic cell.

Zuffada, Cinzia↗

NASA NTRS2015

A Module for Assimilating Hyperspectral Infrared Retrieved Profiles into the Gridpoint Statistical Interpolation System for Unique Forecasting Applications

Hyperspectral infrared sounder radiance data are assimilated into operational modeling systems however the process is computationally expensive and only approximately 1% of available data are assimilated due to data thinning as well as the fact that radiances are restricted to cloud-free fields of view. In contrast, the number of hyperspectral infrared profiles assimilated is much higher since the retrieved profiles can be assimilated in some partly cloudy scenes due to profile coupling other data, such as microwave or neural networks, as first guesses to the retrieval process. As the operational data assimilation community attempts to assimilate cloud-affected radiances, it is possible that the use of retrieved profiles might offer an alternative methodology that is less complex and more computationally efficient to solve this problem. The NASA Short-term Prediction Research and Transition (SPoRT) Center has assimilated hyperspectral infrared retrieved profiles into Weather Research and Forecasting Model (WRF) simulations using the Gridpoint Statistical Interpolation (GSI) System. Early research at SPoRT demonstrated improved initial conditions when assimilating Atmospheric Infrared Sounder (AIRS) thermodynamic profiles into WRF (using WRF-Var and assigning more appropriate error weighting to the profiles) to improve regional analysis and heavy precipitation forecasts. Successful early work has led to more recent research utilizing WRF and GSI for applications including the assimilation of AIRS profiles to improve WRF forecasts of atmospheric rivers and assimilation of AIRS, Cross-track Infrared and Microwave Sounding Suite (CrIMSS), and Infrared Atmospheric Sounding Interferometer (IASI) profiles to improve model representation of tropopause folds and associated non-convective wind events. Although more hyperspectral infrared retrieved profiles can be assimilated into model forecasts, one disadvantage is the retrieved profiles have traditionally been assigned the same error values as the rawinsonde observations when assimilated with GSI. Typically, satellitederived profile errors are larger and more difficult to quantify than traditional rawinsonde observations (especially in the boundary layer), so it is important to appropriately assign observation errors within GSI to eliminate potential spurious innovations and analysis increments that can sometimes arise when using retrieved profiles. The goal of this study is to describe modifications to the GSI source code to more appropriately assimilate hyperspectral infrared retrieved profiles and outline preliminary results that show the differences between a model simulation that assimilated the profiles as rawinsonde observations and one that assimilated the profiles in a module with the appropriate error values.

Berndt, Emily↗

NASA NTRS2016

Comparison of Turbulent Thermal Diffusivity and Scalar Variance Models

In this study, several variable turbulent Prandtl number formulations are examined for boundary layers, pipe flow, and axisymmetric jets. The model formulations include simple algebraic relations between the thermal diffusivity and turbulent viscosity as well as more complex models that solve transport equations for the thermal variance and its dissipation rate. Results are compared with available data for wall heat transfer and profile measurements of mean temperature, the root-mean-square (RMS) fluctuating temperature, turbulent heat flux and turbulent Prandtl number. For wall-bounded problems, the algebraic models are found to best predict the rise in turbulent Prandtl number near the wall as well as the log-layer temperature profile, while the thermal variance models provide a good representation of the RMS temperature fluctuations. In jet flows, the algebraic models provide no benefit over a constant turbulent Prandtl number approach. Application of the thermal variance models finds that some significantly overpredict the temperature variance in the plume and most underpredict the thermal growth rate of the jet. The models yield very similar fluctuating temperature intensities in jets from straight pipes and smooth contraction nozzles, in contrast to data that indicate the latter should have noticeably higher values. For the particular low subsonic heated jet cases examined, changes in the turbulent Prandtl number had no effect on the centerline velocity decay.

computational fluid dynamics↗

DOE OSTI2020

Risk to Ecological Resources on Burial Waste Sites on the Hanford Site, Washington - 20248

The US and other nations have a large nuclear and chemical waste remediation task remaining from World War II, the Cold War, and industrial activities. The US Department of Energy (DOE) has the largest task in the U.S., with Hanford (Washington State) being the most contaminated site, with the longest and most costly cleanup program projected to extend to the end of the century. One of the greatest remediation problems at Hanford are the burial waste facilities, some of which are inside the 200 Areas on the Central Plateau, while others are outside the 200 Areas. These areas are associated with waste sites, pipes, and contaminated trenches, posing technical cleanup challenges to meet DoE's remediation goals, including protection of human health and the environment commensurate with Tri- Party Agreement and current regulatory requirements. In this paper we select one of the legacy waste sites (the BC Cribs and Trenches) as a case study to explore how risk to ecological resources varies during and as a consequence of remediation on sites primarily considered legacy waste sites. We used the risk methodology developed by the Consortium for Risk Evaluation with Stakeholder Participation (CRESP) during the Hanford Site-wide Risk Review. The CRESP methodology included evaluating the risk to receptors, including human health, the Columbia River, groundwater, and ecological and cultural resources. The risk evaluation process for ecological receptors included defining the boundaries of each evaluation unit (EU) to be remediated, definition of a suitable buffer area that might be affected, and consideration of the known resource value based on field work. This information, coupled with known remediation options or strategies, allowed us to develop a risk rating for ecological resources. Site risks were evaluated on a scale of 'not discernible' to 'very high' risk. Consideration of the risks at the buried waste sites was complicated by uncertainties in remediation methods and the degree of restoration that would occur as part of the remediation process. We use the BC Cribs and Trenches as a case study to explore why the risks are low currently but increase during remediation. The current 'not discernible to low' risk rating is due to the lack of high quality ecological resources on the BC Cribs and Trenches EU, except for edge areas (with some state-listed species). During remediation the risks increase to 'medium', mainly due to high quality resources on adjacent buffer lands, which could be disrupted during active remediation. Medium and high quality resources on the edge of the site are continuous with higher quality native shrub-steppe habitat on the buffer. Ecological resources on some of the areas surrounding the 200 East and West Areas are some of the highest quality shrub-steppe areas on the Hanford Site, and risk to these resources could be prevented by planning before and during remediation to minimize disruption. For example, lay down areas for the BC Cribs and Trenches can be minimized by careful planning. Heavy truck and other vehicular traffic can be minimized. Spread of the noxious weeds now on some of the site can be prevented, and the spread of invasive species into the site (particularly near undisturbed shrub-steppe habitat) can be reduced. Further, the areas around and between the 200 East and West Areas have been largely undisturbed and can serve as a source of seeds and propagules for any revegetation that occurs as part of the remediation process. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS2024

Usage of ChatGPT for Engineering Design and Analysis Tool Development

ChatGPT, a generative AI large language model, has recently captured significant attention in both the computer science community and the broader public domain. It has demonstrated a wide range of capabilities, from answering simple questions to writing fully functional computer code. This study spotlights both the capabilities and limitations of ChatGPT when addressing engineering problems. The model's capacity to generate practical engineering tools is highlighted through an example of a prompt that leads to an interactive plotting tool, enabling the examination of the fluid boundary layer around a fan blade. Subsequently, the paper also uncovers potential pitfalls in ChatGPT’s application, shown through an unsuccessful attempt to use ChatGPT to automate a process in Ansys Workbench through scripting. The research further investigates ChatGPT's proficiency in addressing inquiries and providing explanations about the functionalities of OpenMDAO, an open-source, multidisciplinary design, analysis, and optimization tool developed at NASA Glenn Research Center. Finally, an optimization methodology, developed with ChatGPT’s help, is applied to the structural optimization of a fan blade. The developed optimization method utilizes T-Blade3 for geometry generation, Ansys Mechanical for meshing and finite element analysis, and sci-kit learn’s MLPRegressor method to generate a trained neural network model of the design space. OpenMDAO is then used to find the optimal point within the design space. The outcome is a significant reduction in stress in the optimized model—less than one-fifth of the stress value in the baseline model.

Design↗

NASA NTRSDate not supplied

Usage of ChatGPT for Engineering Design and Analysis Tool Development

Design↗