SEARCH · Engineering Papers
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.
Using the dissipation-fluctuation approach, the brightness temperature of a stratified medium with inhomogeneous permittivities and nonuniform temperature profiles is solved. The solution is expressed in closed form ready for machine computation. Numerical results are illustrated and compared with closed-form analytical solutions and results obtained from the WKB method for simple profiles.
A significant amount of research has been conducted to select valve timing and area profiles that create efficient and quiet hydraulic pumps and motors. Numerous active valve architectures have been modeled and optimized, but the rationale for the final solution is often unclear. The solution is usually highly dependent on the modeled valve geometry constraints and the duty cycle of the pump or motor for which the valve was optimized. This paper first presents an ideal valve area profile and what compromises are made to design for a working pressure and speed range. The ideal valve area profile is resolved using a piston-cylinder model with a pressure dependent bulk modulus and variable valve openings. The piston-cylinder is first modeled as a pump with passive (check) valves and a specified pressure drop across the valves to determine the ideal area profile. A generalization of the ideal timing is demonstrated for positive pressure differential and positive rotation direction, also known as the first quadrant. The effects on valve timing when varying pressure are shown for the first quadrant, but the trends can be extrapolated for all quadrants of operation. The ideal valve profile is implemented as fixed valve timing, meaning the valve area is only a function of the timing angle, in a pump-motor. The pump-motor is simulated in one rotation direction through a pressure range. Performance is high in pumping operation, but when the pressure differential is reversed, cylinder pressure spikes ensue. Two strategies to modify an ideal valve area profile are presented: timing grooves and a pressure shifted valve timing. Timing grooves reduce pressure spikes and cavitation in the cylinder but generally increase throttling losses. A pressure shifted valve timing has lower throttling energy losses, making it the favored solution.
We numerically investigate the bandwidth and collimation characteristics of ultrasound beams generated by a simple collimated ultrasound beam source that consists of a piezoelectric disk operated near its radial mode resonances. We simulate the ultrasound beam generated in a fluid medium as a function of the excitation frequency for two cases: (1) free piezoelectric disk that corresponds to zero-traction along the lateral edge and (2) fixed piezoelectric disk that corresponds to zero-displacement along the lateral edge. We present and discuss the physical mechanism underpinning the frequency-dependent collimation and bandwidth properties of the ultrasound beams. We observe that the collimated beam generated by the free disk repeatedly lengthens/shortens and also extends/retracts sidelobes with increasing frequency. Alternatively, fixing the piezoelectric disk results in a consistent beam profile shape across a broad range of frequencies. This facilitates generating broadband signals such as a Gaussian pulse or chirp, which are common in ultrasound imaging. Thus, the fixed piezoelectric disk finds application as a collimated ultrasound beam source in a wide range of applications including medical ultrasound imaging, scanning acoustic microscopy, sonar detection, and other nondestructive ultrasound inspection techniques.
Drift-wave turbulence produces anomalous transport via cross-correlations between fluctuations. This transport has profound implications for confinement, structure formation, and virtually all aspects of the non-linear turbulent dynamics. Herein, we use a data-driven method based on deep learning in order to study turbulent transport in the 2D Hasegawa–Wakatani system and infer a reduced mean-field model from numerical solution. In addition to the usual turbulent diffusion, we find an effect which couples the particle flux to the local gradient of vorticity, which tends to modulate the density profile. The direct coupling to the shear is relatively weak. In addition, the deep learning method finds a model for spontaneous zonal flow generation by negative viscosity, stabilized by non-linear and hyperviscous terms. We compare these results to analytic calculations using quasilinear theory and wave kinetics, finding qualitative agreement, though the calculations miss certain higher-order effects. A simplified, 1-D model for the evolution of the profile, flow, and intensity based on the deep learning results is solved numerically and compared to previous models for staircasing based on bistability. We see that the physics uncovered by the deep learning method provided simple explanations for the formation of zonal structures in the density, flow, and turbulence fields. We highlight the important role of symmetry in the deep learning method and speculate on the portability of the method to other applications.
The direct problem is to compute reflection coefficients for an assumed electron-density profile, using a full-wave solution. The inverse problem is treated numerically, using derivatives of the reflection matrix with respect to model parameters. A technique for the rapid computation of the derivatives is described. An algorithm for determining resolution in the profile is discussed.
The Global Positioning System enables positioning in 3 dimensions about our planet. It has been operational since 1994. Twenty-four satellites are used to aclile\,e this performance. The signals sent by these satellites are electromagnetic waves travelling through our atmosphere down to the small receivers used by the civilian community and the military. Because of varying meteorological conditions (namely, temperature and humidity changes along the ray path), the rays do not travel in a straight line. They bend towards the surface. As a consequence, the ray path between two points is longer than a straight line, and the time it takes for a signal to travel this distance is longer. In 1995, a small GPS receiver was launched on a satellite (GPS/MET). It become possible to perform radio occultations around the Earth: the source - one of the 24 GPS satellites - is seen by the receiver as it rises or sets around the other side of the Earth. When the source disappears, the receiver progressively loses the signals. By measuring accurately the time delay between the emission and the reception of the signal, it is possible to infer which part of the delay is due to the atmosphere. We use GPS/MET data to retrieve temperature and humidity profiles simultaneously. A specific method is implemented: it combines information from numerical forecasts and GPS observations in an optimal way. Comparing the result with an independent source of observations (weather balloons), we demonstrate that GPS data have the potential to improve weather analyses. We also show that improved temperature and humidity profiles can be obtained using information from a forecast model. This confirms results obtained in this study using simulated data.
This paper presents the analytical formulation and numerical results for heat transfer in a high heat flux condenser that relies on capillary flow along shaped fins (Gregorig surfaces) and a drainage network embedded in the condenser walls. Results are shown for a variety of fin profile shapes in order to show the geometric trade-offs involved in seeking a maximum effective heat transfer coefficient for the fin. Predictions of the model show excellent agreement with previously reported measurements for steam. Based on this work, a profile has been selected for a 2 kW ammonia condenser currently under development for use in space. In that design the fin half width is 0.5 mm and the model predicts a heat transfer coefficient referred to the base of the fin of 9 W/sq cm deg C for a heat flux of 10/W sq cm at the base.
Radio occultation experiments have been used to study various properties of planetary atmospheres, including pressure and temperature profiles, and the abundance profiles of absorbing constituents in those planetary atmospheres. However, the reduction of amplitude data from such experiments to determine abundance profiles requires the application of the inverse Abel transform (IAT) and numerical differentiation of experimental data. These two operations preferentially amplify measurement errors above the true signal underlying the data. A new technique for processing radio occultation data has been developed that greatly reduces the errors in the derived absorptivity and abundance profiles. This technique has been applied to datasets acquired from Pioneer Venus Orbiter radio occultation studies and more recently to experiments conducted with the Magellan spacecraft. While primarily designed for radar studies of the Venus surface, the high radiated power (EIRP) from the Magellan spacecraft makes it an ideal transmitter for measuring the refractivity and absorptivity of the Venus atmosphere by such experiments. The longevity of the Pioneer Venus Orbiter has made it possible to study long-term changes in the abundance and distribution of sulfuric acid vapor, H2SO4(g), in the Venus atmosphere between 1979 and 1992. The abundance of H2SO4(g) can be inferred from vertical profiles of 13-cm absorptivity profiles retrieved from radio occultation experiments. Data from 1979 and 1986-87 suggest that the abundance of H2SO4(g) at latitudes northward of 70 deg decreased over this time period. This change may be due to a period of active volcanism in the late 1970s followed by a relative quiescent period, or some other dynamic process in the Venus atmosphere. While the cause is not certain, such changes must be incorporated into dynamic models of the Venus atmosphere. Potentially, the Magellan spacecraft will extend the results of Pioneer Venus Orbiter and allow the continued monitoring of the abundance of distribution of H2SO4(g) in the Venus atmosphere, as well as other interesting atmospheric properties. Without such measurements it will be difficult to address other issues such as the short-term spatial variability of the abundance of H2SO4(g) at similar latitudes in Venus atmosphere, and the identities of particles responsible for large-scale variations observed in NIR images.
For free-space optical communication systems, numerical wave optics simulations provide a useful technique for modeling turbulence-induced beam degradation when the analytical theory is insufficient for characterizing the atmospheric channel. Motivated by such applications we use a split-step method modeling the turbulence as a series of random phase screens using the Hufnagel-Valley turbulence profile. We employ a space-to-ground case study to examine the irradiance and phase statistics for uniformly and non-uniformly located screens and find better agreement with theory using a non-uniform discretization minimizing the contribution of each screen to the total scintillation. In this poster, we summarize the method and the results of the case study including a comparison to layered models used in astronomical imaging applications.
Description of experiments designed to test the sensitivity of a new method for retrieving vertical temperature profiles from satellite radiance measurements to the initial guess of temperature distribution. Emphasis is placed on initial guesses of the kind that would be provided by a numerical model of the atmospheric circulation, with the objective of developing a procedure for direct use of satellite radiation measurements as input to a numerical forecasting scheme.
As temperatures change worldwide, the pattern and competency of disease vectors will change, altering the global distribution of both the burden of infectious disease and the risk of the emergence of those diseases into new regions. To evaluate the risk of potential summer dengue outbreaks triggered by infected travelers under various climate scenarios, we develop an SEIR-type model, run numerical simulations, and conduct sensitivity analyses under a range of temperature profiles. Our model extends existing theoretical frameworks for studying dengue dynamics by introducing temperature dependence of two key parameters: the mosquito extrinsic incubation period and the lifespan of mosquitoes, which empirical data suggests are both highly temperature dependent. We find that changing temperature significantly alters dengue risk in an inverted U-shape, with temperatures in the range 27-31°C producing the highest risk. As temperatures increase beyond 31°C, the determinants of dengue risk begin to shift from mosquito biting rate and carrying capacity to the duration of the human infectious period, suggesting that changing temperatures not only alter dengue risk but also the potential efficacy of control measures. To illustrate the role of spatial and temporal temperature heterogeneity, we select five US cities where the primary dengue vector, the mosquito Aedes aegypti , has been observed, and which have had dengue cases in the past: Los Angeles, Houston, Miami, Brownsville, and Phoenix. Our analysis suggests that an increase of 3°C leads to an approximate doubling of the risk of dengue in Los Angeles and Houston, but a reduction of risk in Miami, Brownsville, and Phoenix due to extreme heat.
Thermal stresses developed during the growth of silicon crystal ribbon have been shown to be negligible, thus eliminating residual stresses and dislocations, if the temperature profile satisfies a second-order partial differential equation inside the ribbon. This has been numerically verified through a finite element model, an outline of which is presented here. This model shows that, for homogeneous isotropic material with temperature independent thermal expansion coefficients, thermal stresses will vanish if the temperature profile satisfies the Laplacian. A comparison of stresses due to uniform and nonuniform temperature gradients in the plane of the ribbon is also presented. The strategies employed to control the round-off error and to validate the computer model are discussed.
Explore the source record for details and available documents.
Abstract The classical problem of neutral beam ions slowing down in a magnetic mirror geometry is revisited to provide predictive capability for the new Wisconsin HTS Axisymmetric Mirror under development at the University of Wisconsin. A Fokker–Planck model named fast beam ion solver (FBIS) is developed to include the spatial non-uniformity of a physical mirror geometry. The mathematical framework allows for efficient orbit averaging of the pitch-angle scattering operator, and permits a determination of the axial profile of the ambipolar potential confining the electrons. The numerical results from FBIS are consistent with earlier work, but further show how mirror-ratio and a near square magnetic well optimizes the fusion gain. The numerical results are also applied to inform the conceptual design of WHAM++, a low capital-cost breakeven-class magnetic mirror device.
Abstract Design for cooling effectiveness in turbine blades relies on accurate models for dynamic losses and heat transfer of internal cooling passages. Metal additive manufacturing (AM) has expanded the design space for these configurations, but can give rise to large-scale roughness features. The range of roughness length scales in these systems makes morphology resolved computational fluid dynamics (CFD) impractical. However, volumetric roughness models can be leveraged, as they have computational costs orders of magnitude lower. In this work, a discrete element roughness model (DERM), based on the double-averaged Navier–Stokes equations, is presented and applied to additively manufactured rough channels, representative of gas turbine blade cooling passages. Unique to this formulation of DERM is a generalized sheltering-based treatment of drag, a two-layer model for spatially averaged Reynolds stresses, and explicit treatment of dispersion. Six different AM rough surface channel configurations are studied, with roughness trough to peak sizes ranging from 15% to 60% nominal channel passage half-width, and the roughness Reynolds number ranges from Rek = 60 to 300. DERM predictions for spatially and temporally averaged mean flow quantities are compared to previously reported direct numerical simulation results. Good agreement in the mean velocity profiles, stress balances, and drag partitions are observed. While DERM models are typically calibrated to specific deterministic roughness morphologies at comparatively small roughness Reynolds numbers, the present more generalized DERM formulation has wider applicability. Here, it is demonstrated that the model can accommodate random roughness of large scale, typical of AM.
The DOE-GTO-funded project, award number 5.1.2.12, entitled "EXERGETIC - De-risking Next Generation Resilient Geothermal Hybrids via at-Scale Evaluation Using a Virtual Emulation Digital Twin Environment for Efficient Operation," advances the solution to these challenges by developing and validating a geothermal co-emulation environment implemented at the National Laboratory of the Rockies (NLR)'s Advanced Research on Integrated Energy Systems (ARIES) platform. This framework enables the de-risking of next-generation geothermal and geothermal hybrid systems through high-fidelity modeling, real-time digital emulation, advanced control strategies, and techno-economic assessment. The project focused on geothermal hybrid configurations that integrate geothermal power plants with concentrated solar power and underground thermal energy storage, enabling enhanced efficiency, flexibility, and grid support capabilities. The main goal of this project was the development of a geothermal digital co-emulation environment to demonstrate the technical and economic value of geothermal hybrid systems and their contribution to grid stability and flexibility. The EXERGETIC framework combined physics-based models, controls, and real assets at ARIES, including digital real-time simulators (DRTS), a 20-MW-scale controllable grid interface (CGI), and a 2-MW conventional generator. Detailed transient models were developed for the key subsystems of a hybrid geothermal plant, including parabolic trough solar collectors, reservoir thermal energy storage (RTES), and a binary Organic Rankine Cycle (ORC) power plant. The ORC model explicitly captured thermal inertia and off-design operation and integrated control strategies to dynamically respond to electric load profiles. The models were validated against published experimental and numerical studies, demonstrating strong agreement and confirming the accuracy and robustness of the modeling approach. The resulting digital twin represents geothermal-solar-storage systems at multiple scales (1 MW to 100 MW) and enables realistic emulation of grid-connected operation. The control architecture allows the geothermal resource to provide stable baseload generation, while solar and stored thermal energy supply flexible, dispatchable support during periods of high demand or variable grid conditions. A key contribution of the EXERGETIC project is the demonstration that geothermal hybrid systems can be designed to be active grid assets rather than passive baseload generators. Using the ARIES platform, the digital twin was evaluated under multiple grid scenarios, including load following, voltage support at the distribution level, and frequency response at the transmission level. Results show that hybrid geothermal systems can respond effectively to dynamic grid conditions, providing inertia-like behavior, primary frequency support, and voltage regulation through coordinated control. In addition to the performance and grid services capability analysis of geothermal and hybrid geothermal systems, the EXERGETIC project also focused on scalability and techno-economic analysis of geothermal hybrid plants. In particular, for the scalability analysis, machine-learning (ML)-based surrogate models were trained using data generated from the geothermal digital twin under different grid-connected scenarios and plant capacities. These ML models demonstrated strong interpolation and extrapolation capabilities across plant sizes, accurately reproducing both steady-state and transient responses with very low errors. Regarding the techno-economic analysis, plant performance results were integrated with cost models for hybrid geothermal systems, and the levelized cost of electricity (LCOE) was used as the main economic metric to evaluate system performance across a range of system capacities, solar shares, solar multiples, and storage durations. Results indicate that economies of scale significantly reduce geothermal LCOE as plant capacity increases, with large-scale systems (25-100 MW) achieving substantially lower costs than small plants. Hybridization with solar thermal energy and storage further improves economic performance by increasing capacity utilization and enabling flexible dispatch. In addition, thermal storage plays a critical role in reducing LCOE by maximizing geothermal, solar, and stored energy resources. In summary, the results from this project demonstrate that geothermal hybrid systems represent a promising alternative for increasing the energy conversion efficiency of geothermal technologies, contributing to the preservation of geothermal resources, and supporting the transition of geothermal plants from traditional baseload resources into flexible, resilient, and cost-competitive energy conversion technologies.
Numerical study of nonlinear inversion applied to atmospheric temperature profiles
Four W Ursae Majoris-type binaries were investigated using the simple model of blackbody solid spheres. Numerical computations were carried out for the rotationally broadened profiles of spectral lines Sr II (4077.7) and Ca I (4226.7) and the mean radial velocities. These mean velocities at different phases outside eclipses were compared with sinusoidal values of the radial velocity curve for point masses. Within the present accuracy of observations, the reflection effect is small and may be neglected in the analysis of radial velocity curves.