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.
The development and characteristics of heat pipes are examined. The subjects discussed are: (1) principles of operation, (2) heat pipe theory, (3) pressure gradient effects, (4) variable conductance, (5) design procedure, and (6) performance limit evaluation.
The atmospheric general circulation experiment (AGCE) uses a rotating fluid flow cell assembly. The key technical areas affecting the feasibility of the design and operation of the AGCE are investigated. The areas investigated include materials for the flow cell assembly, thermal design, high voltage power supply design, effective retrieval and handling of experiment data and apparatus configuration. Several materials, DMSO and m-tolunitrile, were selected as candidate fluids for the flow cell principally for their high dielectric constant which permits the high voltage power supply design to be held to 15 kV and still simulate terrestrial gravity. Achievement of a low dissipation factor in the fluid to minimize internal heating from the applied electrical field depends strongly on purification and handling procedures. The use of sapphire as the outer hemisphere for the flow cell provides excellent viewing conditions without a significant impact on attaining the desired thermal gradients. Birefringent effects from sapphire can be held to acceptably low limits. Visualization of flow fluid is achieved through the motion of a dot matrix formed by photochromic dyes. Two dyes found compatible with the candidate fluids are spiropyran and triarylmethane. The observation of the dot motion is accomplished using a flying spot scanner.
The equations for the single-particle orbits in a nonneutral high density plasma in the presence of inhomogeneous crossed fields are obtained. Using these orbits, the linearized Vlasov equation is solved as an expansion in the orbital radii in the presence of inhomogeneities and density gradients. A model distribution function is introduced whose cold-fluid limit is exactly the same as that used in many previous studies of the cold-fluid equations. This model function is used to reduce the linearized Vlasov-Poisson equations to a second-order ordinary differential equation for the linearized electrostatic potential whose eigenvalue is the perturbation frequency.
Wing leading edge heat pipes were conceptually designed for three types of vehicle: an entry research vehicle, aero-space plane, and advanced shuttle. A full scale, internally instrumented sodium/Hastelloy X heat pipe was successfully designed and fabricated for the advanced shuttle application. The 69.4 inch long heat pipe reduces peak leading edge temperatures from 3500 F to 1800 F. It is internally instrumented with thermocouples and pressure transducers to measure sodium vapor qualities. Large thermal gradients and consequently large thermal stresses, which have the potential of limiting heat pipe life, were predicted to occur during startup. A test stand and test plan were developed for subsequent testing of this heat pipe. Heat pipe manufacturing technology was advanced during this program, including the development of an innovative technique for wick installation.
The problem of day-to-day variability in the occurrence of equatorial spread F (ESF) is addressed using multidiagnostic observations and semiempirical modeling. The observational results are derived from a two-night case study of ESF onset conditions observed at Kwajalein Atoll (Marshall Islands) using the ALTAIR incoherent scatter radar and all-sky optical imaging techniques. The major difference between nights when ESF instabilities did not occur (August 14, 1988) and did occur (August 15, 1988) in the Kwajalein sector was that the northern meridional gradient of 6300-A airglow was reduced on the night of limited ESF activity. Modeling results suggest that this unusual airglow pattern is due to equatorward neutral winds. Previous researchers have shown that transequatorial thermospheric winds can exert a control over ESF seasonal and longitudinal occurrence patterns by inhibiting Rayleigh-Taylor instability growth rates. Evidence is presented to suggest that this picture can be extended to far shorter time scales, namely, that 'surges' in transequatorial winds acting over characteristic times of a few hours to a day can result in a stabilizing influence upon irregularity growth rates. The seemingly capricious nature of ESF onset may thus be controlled, in part, by the inherent variability of low-latitude thermospheric winds.
Passive depth and displacement map determinations have become an important part of computer vision processing. Applications that make use of this type of information include autonomous navigation, robotic assembly, image sequence compression, structure identification, and 3-D motion estimation. With the reliance of such systems on visual image characteristics, a need to overcome image degradations, such as random image-capture noise, motion, and quantization effects, is clearly necessary. Many depth and displacement estimation algorithms also introduce additional distortions due to the gradient operations performed on the noisy intensity images. These degradations can limit the accuracy and reliability of the displacement or depth information extracted from such sequences. Recognizing the previously stated conditions, a new method to model and estimate a restored depth or displacement field is presented. Once a model has been established, the field can be filtered using currently established multidimensional algorithms. In particular, the reduced order model Kalman filter (ROMKF), which has been shown to be an effective tool in the reduction of image intensity distortions, was applied to the computed displacement fields. Results of the application of this model show significant improvements on the restored field. Previous attempts at restoring the depth or displacement fields assumed homogeneous characteristics which resulted in the smoothing of discontinuities. In these situations, edges were lost. An adaptive model parameter selection method is provided that maintains sharp edge boundaries in the restored field. This has been successfully applied to images representative of robotic scenarios. In order to accommodate image sequences, the standard 2-D ROMKF model is extended into 3-D by the incorporation of a deterministic component based on previously restored fields. The inclusion of past depth and displacement fields allows a means of incorporating the temporal information into the restoration process. A summary on the conditions that indicate which type of filtering should be applied to a field is provided.
Observations and numerical simulations show that winds near Gibraltar Strait cause an Atlantic Ocean to Mediterranean Sea sea level difference of 20 cm peak to peak with a 3-cm standard deviation for periods of days to years. Theoretical arguments and numerical experiments establish that this wind-driven sea level difference is caused in part by storm surges due to alongshore winds near the North African coastline on the Atlantic side of Gibraltar. The fraction of the Moroccan coastal current offshore of the 284-m isobath is deflected across Gibraltar Strait, west of Camarinal Sill, resulting in a geostrophic surface pressure gradient that contributes to a sea level difference at the stationary limit. The sea level difference is also caused in part by the along-strait wind setup, with a contribution proportional to the along-strait wind stress and to the length of Gibraltar Strait and adjoining regions and inversely proportional to its depth.
Using Multi-angle Imaging SpectroRadiometer (MISR) and sub-orbital measurements from the 2006 INTEX-B/MILAGRO field campaign, in this study we explore MISR's ability to map different aerosol air mass types over the Mexico City metropolitan area. The aerosol air mass distinctions are based on shape, size and single scattering albedo retrievals from the MISR Research Aerosol Retrieval algorithm. In this region, the research algorithm identifies dust-dominated aerosol mixtures based on non-spherical particle shape, whereas spherical biomass burning and urban pollution particles are distinguished by particle size. Two distinct aerosol air mass types based on retrieved particle microphysical properties, and four spatially distributed aerosol air masses, are identified in the MISR data on 6 March 2006. The aerosol air mass type identification results are supported by coincident, airborne high-spectral-resolution lidar (HSRL) measurements. Aerosol optical depth (AOD) gradients are also consistent between the MISR and sub-orbital measurements, but particles having single-scattering albedo of approx. 0.7 at 558 nm must be included in the retrieval algorithm to produce good absolute AOD comparisons over pollution-dominated aerosol air masses. The MISR standard V22 AOD product, at 17.6 km resolution, captures the observed AOD gradients qualitatively, but retrievals at this coarse spatial scale and with limited spherical absorbing particle options underestimate AOD and do not retrieve particle properties adequately over this complex urban region. However, we demonstrate how AOD and aerosol type mapping can be accomplished with MISR data over complex urban regions, provided the retrieval is performed at sufficiently high spatial resolution, and with a rich enough set of aerosol components and mixtures.
This paper will describe the results that have been obtained to date concerning MMS formation flying. The MMS spacecraft spin at a rate of 3.1 RPM, with spin axis roughly aligned with Ecliptic North. Several booms are used to deploy instruments: two 5 m magnetometer booms in the spin plane, two rigid booms of length 12.5 m along the positive and negative spin axes, and four flexible wire booms of length 60 m in the spin plane. Minimizing flexible motion of the wire booms requires that reorientation of the spacecraft spin axis be kept to a minimum: this is limited to attitude maneuvers to counteract the effects of gravity-gradient and apparent solar motion. Orbital maneuvers must therefore be carried out in essentially the nominal science attitude. These burns make use of a set of monopropellant hydrazine thrusters: two (of thrust 4.5 N) along the spin axis in each direction, and eight (of thrust 18 N) in the spin plane; the latter are pulsed at the spin rate to produce a net delta-v. An on-board accelerometer-based controller is used to accurately generate a commanded delta-v. Navigation makes use of a weak-signal GPS-based system: this allows signals to be received even when MMS is flying above the GPS orbits, producing a highly accurate determination of the four MMS orbits. This data is downlinked to the MMS Mission Operations Center (MOC) and used by the MOC Flight Dynamics Operations Area (FDOA) for maneuver design. These commands are then uplinked to the spacecraft and executed autonomously using the controller, with the ground monitoring the burns in real time.
This paper presents an analysis and simulation of evaporation and condensation at a motionless liquid/vapor interface. A 1-D model equation, emphasizing heat and mass transfer at the interface, is solved in two ways, and incorporated into a subgrid interface model within a CFD simulation. Simulation predictions are compared with experimental data from the CPST Engineering Design Unit tank, a cryogenic fluid management test tank in 1-g. The numerical challenge here is the physics of the liquid/vapor interface; pressurizing the ullage heats it by several degrees, and sets up an interfacial temperature gradient that transfers heat to the liquid phase-the rate limiting step of condensation is heat conducted through the liquid and vapor. This physics occurs in thin thermal layers O(1 mm) on either side of the interface which is resolved by the subgrid interface model. An accommodation coefficient of 1.0 is used in the simulations which is consistent with theory and measurements. This model is predictive of evaporation/condensation rates, that is, there is no parameter tuning.
Recently, Tao and Mo developed a semilocal exchange-correlation density functional. The exchange part of this functional is derived from a density-matrix expansion corrected to reproduce the fourth-order gradient expansion of the exchange energy in the slowly-varying-density limit, while the correlation part is based on the Tao-Perdew-Staroverov-Scuseria (TPSS) correlation functional, with a modification for the low-density limit. In this work, the Tao-Mo (TM) functional is assessed by computing various properties of solids and jellium surfaces. This includes 22 lattice constants and bulk moduli, 30 band gaps, seven cohesive energies, and jellium surface exchange and correlation energies for the density parameter r s in the range from 2 to 3 bohr. Our calculations show that the TM approximation can yield consistently high accuracy for most properties considered here, with mean absolute errors (MAEs) of 0.025 Å for lattice constants, 7.0 GPa for bulk moduli, 0.08 eV/atom for cohesive energies, and 35 erg / cm 2 for surface exchange-correlation energies. The MAE in band gaps is larger than that of TPSS, but slightly smaller than the errors of the local spin-density approximation, Perdew-Burke-Ernzerhof generalized gradient approximation, and revised TPSS. However, band gaps are still underestimated, particularly for large-gap semiconductors, compared to the Heyd-Scuseria-Ernzerhof nonlocal screened hybrid functional.
Numerical weather prediction (NWP) and climate models equipped with wind-farm parameterizations (WFPs) can simulate cluster wake effects affecting downstream wind farms in both onshore and offshore environments. This study evaluates wake recovery behind a wind farm represented by the NWP-WFP approach in the Weather Research and Forecasting (WRF) model using either the Fitch et al. (2012) or Ma et al. (2022a, b) WFPs. Results are benchmarked against large-eddy simulations (LES) of an idealized offshore wind farm with aligned and staggered layouts under neutral atmospheric stability. Near-farm wake recovery is underestimated in NWP-WFP simulations due to its representation on a coarse mesoscale grid. This limitation leads to slow wake recovery through two interconnected mechanisms: (i) spatial gradients in the wind velocity field are weaker compared to LES and (ii) turbulence kinetic energy (TKE) remains low not because of excessive dissipation but due to insufficient shear production caused by these weakened gradients. For the scenario considered here, a wind-speed bias develops in the near-farm wake and persists into the far wake. Differences between the NWP-WFP simulations and LES emerge within a short distance downstream of the farm exit, where the mesoscale simulations recover too slowly. This reduced recovery contributes approximately 0.15-0.50 m s-1 to the near-farm wind-speed bias. The bias established in this region is not subsequently compensated for downstream but instead propagates into the far wake, where wind-speed differences of approximately 0.4-0.6 m s-1 remain up to 50 km downstream. Higher-resolution mesoscale simulations partially reduce this bias. Increasing turbine-added TKE or including subgrid wake effects provides additional improvement, but neither fully addresses the underlying cause. The slow wake recovery is not caused by limitations of the WFPs themselves, as it also occurs outside their region of influence, and adding subgrid wake effects does not significantly impact recovery. Rather, the slow wake recovery is a consequence of mesoscale flow representation. This behavior is not limited to regions downstream of the wind farm but is less visible within the farm, where wake recovery occurs simultaneously with turbine-induced momentum extraction. These results highlight the need for improved representations of wake recovery both within and downstream of wind farms. While enhanced subgrid modeling, shear-driven TKE production, and refined WFP formulations may improve intra-farm dynamics, accurately capturing near-farm wake recovery downstream remains challenging, as WFPs do not act in this region.
Deep learning has been shown to be an effective tool in solving partial differential equations (PDEs) through physics-informed neural networks (PINNs). PINNs embed the PDE residual into the loss function of the neural network, and have been successfully employed to solve diverse forward and inverse PDE problems. However, one disadvantage of the first generation of PINNs is that they usually have limited accuracy even with many training points. Here, we propose a new method, gradient-enhanced physics-informed neural networks (gPINNs), for improving the accuracy of PINNs. gPINNs leverage gradient information of the PDE residual and embed the gradient into the loss function. We tested gPINNs extensively and demonstrated the effectiveness of gPINNs in both forward and inverse PDE problems. Our numerical results show that gPINN performs better than PINN with fewer training points. Additionally, we combined gPINN with the method of residual-based adaptive refinement (RAR), a method for improving the distribution of training points adaptively during training, to further improve the performance of gPINN, especially in PDEs with solutions that have steep gradients.
The minimization of cost and energy consumption of future particle accelerators, both large and small, depends upon the development of new materials for the active surfaces of superconducting RF (SRF) accelerating structures. SRF properties are inherently a surface phenomenon as the RF only penetrates the London penetration depth ?, typically between 20 and 400 nm depending on the material. The greatest potential for dramatic new performance capabilities lies with methods and materials which deliberately produce the sub-micron-thick critical surface layer in a controlled way. Thus, two opportunities arise for the use of SRF thin ?lms as single layer superconductor or multilayer Superconductor-Insulator-Superconductor structures: SRF properties are inherently a surface phenomenon, and when other technological processes are optimized, the fundamental limit to the maximum supportable RF ?eld amplitude is understood to be the ?eld at which the magnetic ?ux ?rst penetrates into the surface, Hc1. Niobium, the material most exploited for SRF accelerator applications, has Hc1 ? 170 mT, which yields a maximum acceler?ating gradient of less than 50 MV/m. Thus, new materials must be investigated to push this limit to the desired range. Although many compounds have shown higher critical superconducting transition temperature (Tc), their critical ?eld (Hc1) is lower than niobium. To circumvent this issue, a Super?conductor/Insulator/Superconductor (S-I-S) multilayer ?lms approach has been proposed to achieve performance in excess of that of bulk Nb. In this system, the critical ?eld is enhanced and the ?uxon penetration delayed by using, on top of the Nb surface, a superposition of superconductor layers, each less than a penetration depth in thickness, separated by a dielectric layer a few nanometers thick. For bulk type-II superconducting materials, vortex entry occurs at surface magnetic ?eld Hc1. The key aspect of Gurevich's analysis is the recognition that vortex entry is inhibited if the surface layer thick?ness (d) of the superconductor is less than the penetration depth, ?. The overlayer provides magnetic screening of the underlying Nb with the attenuation factor exp(d/ ?0). The Nb surface can remain in the Meissner state at ?elds much higher than in bulk due to the increase of the parallel Hc1 in a thin ?lm. This could potentially lead to further improvement in RF cavity performance using the bene?t of the higher-Tc superconductors without being limited by their lower Hc1. Niobium on copper (Nb/Cu) technology for superconducting cavities has proven over the years to be a viable alternative to bulk niobium. However the deposition techniques used for cavities, mainly magnetron sputtering, have not yielded, so far, SRF surfaces suitable for high ?eld performance. High quality ?lms can be grown using methods of energetic condensation, such as Electron Cy?clotron Resonance (ECR) Nb ion source in UHV and High Power Impulse Magnetron Sputtering (HiP?IMS) with energetic ions only and neutral atoms signi?cantly blocked. The main advantages of these techniques are the production of higher ?ux of ions with controllable incident angle and kinetic energy and the absence of macroparticle production. The project in general, is focusing onto understanding the physics underlying the RF superconduc?tivity in layered structures and developing SRF material ?lms suitable as substrate and S-I-S multilayer structures. In ?rst hand, the primary focus is to understand the physics governing the relationship be?tween growth conditions, ?lm microstructure and RF performance for SRF materials produced by energetic condensation. The deposition by energetic condensation in UHV via ECR and HiPIMS are studied for Nb and B1 superconducting compounds on substrates with various lattice match and sur?face conditions (in-situ etching, high energy assisted atomic annealing, deposition temperature). The general energy-temperature zone model developed by André Anders is explored to de?ne the optimized coating conditions for ?lms suitable for single layer and multilayer SRF structures. The contribution of both the ion incident energy and the substrate lattice energy on the nucleation of the SRF materials considered, and some insulating materials (AlN and MgO) and the in?uence of the coating conditions on the resulting materials, superconducting and/or RF properties will be investigated. Attention is be given to minimize the residual stress (intrinsic and thermal) with the choice of lattice match between the ?lm and substrate materials, the coating temperature, and to evaluate the coating methods used in terms of conformality for 3D structures. In second hand, as a proof of concept, the coating condi?tions developed will be used to prepare multi-layer superconductor/insulator/superconductor (S-I-S) ?lms having high quality factor and high gradient breakthrough potential. A thick Nb ?lm (> 1µm) could be used as a substrate. Particular attention is given to the interface between superconducting and insulator ?lms. For this work, various material analyses techniques have been used to character?ize each type of thin ?lm material: electron backscatter di?raction (EBSD), X-ray di?raction (XRD), high-resolution scanning electron microscopy (SEM), transmission electron microscopy (TEM), stylus and AFM pro?lometry, secondary ion mass spectroscopy (SIMS), scanning auger microscopy (SAM). Superconducting and RF measurements are also performed (RRR, Tc, and low-temperature SRF sur- face impedance, SQUID-magnetometry). The insight gained from this project on the physics involved in RF superconductivity of layered structures open the way to enable major reductions in both capital and operating costs associated with future particle accelerators across the spectrum from low foot?print compact machines to energy frontier facilities. Realization of quality thin ?lm SRF cavities may dramatically change the cost structure of acceleration systems and signi?cantly improve performance consistency and reliability.
Quiescent H-mode (QH-mode) is a naturally ELM-stable high performance operation mode which has been obtained in DIII-D, ASDEX Upgrade, JET, and JT-60U [1 and references within]. The additional transport to maintain constant density and radiation is provided by the benign coherent edge harmonic oscillations mode (EHO) [2] which is thought to be a kink/peeling mode [3]. A new QH-mode regime with enhanced pedestal has been discovered in DIII-D at reactor-relevant low torque and collisionality [4]. The regime was originally found in conventional QH-mode when the counter-Ip neutral beam torque drops to ~2Nm in double null plasma shapes. It is referred to as ‘wide-pedestal QH-mode’ because the pedestal width exceeds the EPED kinetic ballooning mode (KBM) limit [5]. Across the transition from QH to wide-pedestal QH, the pedestal electron pressure generally increases by 60% and widens by 50% and the plasma confinement (ITER Hy98y2) rises by 40%. The onset of broadband edge MHD modes and micro-turbulence accompanied with a lower ExB shear in this region is a common feature of the wide-pedestal QH, instead of EHO. It is conjectured that the increased transport provided by these edge modes reduces the pedestal gradients resulting in a higher pedestal while still remaining below the ELM-limit.
We show that Streptomyces biogeography in soils across North America is influenced by the regional diversification of microorganisms due to dispersal limitation and genetic drift. Streptomyces spp. form desiccation-resistant spores, which can be dispersed on the wind, allowing for a strong test of whether dispersal limitation governs patterns of terrestrial microbial diversity. We employed an approach that has high sensitivity for determining the effects of genetic drift. Specifically, we examined the genetic diversity and phylogeography of physiologically similar Streptomyces strains isolated from geographically distributed yet ecologically similar habitats. We found that Streptomyces beta diversity scales with geographic distance and both beta diversity and phylogenetic diversity manifest in a latitudinal diversity gradient. This pattern of Streptomyces biogeography resembles patterns seen for diverse species of plants and animals, and we therefore evaluated these data in the context of ecological and evolutionary hypotheses proposed to explain latitudinal diversity gradients. The data are consistent with the hypothesis that niche conservatism limits dispersal, and historical patterns of glaciation have limited the time for speciation in higher-latitude sites. Most notably, higher-latitude sites have lower phylogenetic diversity, higher phylogenetic clustering, and evidence of range expansion from lower latitudes. In addition, patterns of beta diversity partition with respect to the glacial history of sites. Furthermore, the data support the hypothesis that extant patterns of Streptomyces biogeography have been driven by historical patterns of glaciation and are the result of demographic range expansion, dispersal limitation, and regional diversification due to drift.
Ground-based Doppler-lidar instrumentation provides atmospheric wind data at dramatically improved accuracies and spatial/temporal resolutions. These capabilities have provided new insights into atmospheric flow phenomena, but they also should have a strong role in NWP model improvement. Insight into the nature of model errors can be gained by studying recurrent atmospheric flows, here a regional summertime diurnal sea breeze and subsequent marine-air intrusion into the arid interior of Oregon-Washington, where these winds are an important windenergy resource. These marine intrusions were sampled by three scanning Doppler lidars in the Columbia-River Basin as part of the Second Wind Forecast Improvement Project (WFIP2), using data from summer 2016. Lidar time-height cross sections of wind speed identified eight days when the diurnal-flow cycle, (peak wind speeds at midnight, afternoon minima) was obvious and strong. Eight-day composite time-height cross sections of lidar wind speeds are used to validate those generated by the operational NCEP-HRRR model. HRRR simulated the diurnal wind cycle, but produced errors in the timing of onset and significant errors due to a premature nighttime demise of the intrusion flow, producing low-bias errors of 6 m s-1 . Day-to-day and in the composite, whenever a marine intrusion occurred, HRRR made these same errors. The errors occurred under a range of gradient wind conditions indicating that they resulted from the misrepresentation of physical processes within a limited region around the measurement locations. Because of their generation within a limited geographical area, field-measurement programs can be designed to find and address the sources of these NWP errors.
Hamiltonian Dielectric Solvent (HADES) is a recent method [S. Bauer et al., J. Chem. Phys. 140, 104103 (2014)] which enables atomistic Hamiltonian molecular dynamics (MD) simulations of peptides and proteins in dielectric solvent continua. Such simulations become rapidly impractical for large proteins, because the computational effort of HADES scales quadratically with the number N of atoms. If one tries to achieve linear scaling by applying a fast multipole method (FMM) to the computation of the HADES electrostatics, the Hamiltonian character (conservation of total energy, linear, and angular momenta) may get lost. Here, we show that the Hamiltonian character of HADES can be almost completely preserved, if the structure-adapted fast multipole method (SAMM) as recently redesigned by Lorenzen et al. [J. Chem. Theory Comput. 10, 3244-3259 (2014)] is suitably extended and is chosen as the FMM module. By this extension, the HADES/SAMM forces become exact gradients of the HADES/SAMM energy. Their translational and rotational invariance then guarantees (within the limits of numerical accuracy) the exact conservation of the linear and angular momenta. Also, the total energy is essentially conserved—up to residual algorithmic noise, which is caused by the periodically repeated SAMM interaction list updates. These updates entail very small temporal discontinuities of the force description, because the employed SAMM approximations represent deliberately balanced compromises between accuracy and efficiency. The energy-gradient corrected version of SAMM can also be applied, of course, to MD simulations of all-atom solvent-solute systems enclosed by periodic boundary conditions. However, as we demonstrate in passing, this choice does not offer any serious advantages.