Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct numerical simulation”

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 1,585 records · Page 88

Numerical Analysis of a Rotating Detonation Engine in the Relative Reference Frame

A two-dimensional, computational fluid dynamic (CFD) simulation of a semi-idealized rotating detonation engine (RDE) is described. The simulation operates in the detonation frame of reference and utilizes a relatively coarse grid such that only the essential primary flow field structure is captured. This construction yields rapidly converging, steady solutions. Results from the simulation are compared to those from a more complex and refined code, and found to be in reasonable agreement. The performance impacts of several RDE design parameters are then examined. Finally, for a particular RDE configuration, it is found that direct performance comparison can be made with a straight-tube pulse detonation engine (PDE). Results show that they are essentially equivalent.

pulse detonation engines↗

Numerical Analysis of a Rotating Detonation Engine in the Relative Reference Frame

A two-dimensional, computational fluid dynamic (CFD) simulation of a semi-idealized rotating detonation engine (RDE) is described. The simulation operates in the detonation frame of reference and utilizes a relatively coarse grid such that only the essential primary flow field structure is captured. This construction yields rapidly converging, steady solutions. Results from the simulation are compared to those from a more complex and refined code, and found to be in reasonable agreement. The performance impacts of several RDE design parameters are then examined. Finally, for a particular RDE configuration, it is found that direct performance comparison can be made with a straight-tube pulse detonation engine (PDE). Results show that they are essentially equivalent.

propulsion↗

Step Bunching: Influence of Impurities and Solution Flow

Step bunching results in striations even at relatively early stages of its development and in inclusions of mother liquor at the later stages. Therefore, eliminating step bunching is crucial for high crystal perfection. At least 5 major effects causing and influencing step bunching are known: (1) Basic morphological instability of stepped interfaces. It is caused by concentration gradient in the solution normal to the face and by the redistribution of solute tangentially to the interface which redistribution enhances occasional perturbations in step density due to various types of noise; (2) Aggravation of the above basic instability by solution flowing tangentially to the face in the same directions as the steps or stabilization of equidistant step train if these flows are antiparallel; (3) Enhanced bunching at supersaturation where step velocity v increases with relative supersaturation s much faster than linear. This v(s) dependence is believed to be associated with impurities. The impurities of which adsorption time is comparable with the time needed to deposit one lattice layer may also be responsible for bunching; (4) Very intensive solution flow stabilizes growing interface even at parallel solution and step flows; (5) Macrosteps were observed to nucleate at crystal corners and edges. Numerical simulation, assuming step-step interactions via surface diffusion also show that step bunching may be induced by random step nucleation at the facet edge and by discontinuity in the step density (a ridge) somewhere in the middle of a face. The corresponding bunching patterns produce the ones observed in experiment. The nature of step bunching generated at the corners and edges and by dislocation step sources, as well as the also relative importance and interrelations between mechanisms 1-5 is not clear, both from experimental and theoretical standpoints. Furthermore, several laws controlling the evolution of existing step bunches have been suggested, though unambiguous conclusions are still missing. Addressing these issues is the major goal of the present project. The theory addressing the above problem, experimental methods, several figures which include: (1) the spatial wave numbers at which the system is neutrally stable as a function of growth velocity for linear kinetics and supersaturation for nonlinear kinetics; (2) a schematic of the experiment of lysozyme crystal growing under conditions of natural convection; (3) fluctuations in time, t, of the normal growth rate, R(t), vicinal slope, p(t) and Fourier Spectra of R(t), discussions and conclusions are presented.

Chernov, A. A.↗

A numerical model of magnetosphere-ionosphere coupling Preliminary results

A three-dimensional simulation model was developed to study magnetosphere-ionosphere coupling in the auroral region. One-fluid MHD equations are adopted to model the magnetosphere, and current density continuity equations were solved consistently to model the ionosphere. In the preliminary simulation runs described here, the electrodynamics of region 1 field-aligned currents were modeled. Initially, the electric field is taken to be in the magnetospheric equatorial plane. Alfven waves then propagate down to the ionosphere, accompanying a field-aligned current and exciting the electrostatic potential in the ionosphere by electrodynamic coupling. The results of the preliminary runs directly correspond with the fundamental characteristics of global magnetosphere-ionosphere coupling. These characteristics include an ionospheric electrostatic potential which varies in its development in time, depending upon the ratio of ionospheric resistance to magnetospheric impedance, the flowing of field-aligned currents into the ionosphere on the dawnside and out on the duskside, and the distribution along the geomagnetic field of the amplitude of the field-aligned current density in proportion to the field intensity.

Watanabe, K.↗

Hybrid-finite-element analysis of some nonlinear and 3-dimensional problems of engineering fracture mechanics

In this paper, efficient numerical methods for the analysis of crack-closure effects on fatigue-crack-growth-rates, in plane stress situations, and for the solution of stress-intensity factors for arbitrary shaped surface flaws in pressure vessels, are presented. For the former problem, an elastic-plastic finite element procedure valid for the case of finite deformation gradients is developed and crack growth is simulated by the translation of near-crack-tip elements with embedded plastic singularities. For the latter problem, an embedded-elastic-singularity hybrid finite element method, which leads to a direct evaluation of K-factors, is employed.

Atluri, S. N.↗

NASA Ames Research Center's parabolized Navier-Stokes code - A critical evaluation of heat-transfer predictions

The viscous supersonic/hypersonic flow over a biconic configuration was numerically simulated using the NASA Ames Research Center's Parabolized Navier-Stokes (PNS) code, and results obtained for the effects of various computational parameters on the heat transfer are compared with experimental results. The PNS code is found to provide accurate results for heat transfer parameters using low smoothing coefficient values, assuming that the radial spacing is at some small value and that the body is simple in cross section. For the cases considered, the nominal value of the spatial step size was 0.05, and the grid density was 45 points in the meridional direction and 30 points in the radial direction.

Chaussee, Denny S.↗

Improving Rain/No-Rain Detection Skill by Merging Precipitation Estimates from Different Sources

Rain/no-rain detection error is a key source of uncertainty in regional and global precipitation products that propagates into offline hydrological and land surface modeling simulations. Such detection error is difficult to evaluate and/or filter without access to high-quality reference precipitation datasets. For cases where such access is not available, this study proposes a novel approach for improved rain/no-rain detection. Based on categorical triple collocation (CTC) and a probabilistic framework, a weighted merging algorithm (CTC-M) is developed to combine noisy, but independent, precipitation products into an optimal binary rain/no-rain time series. Compared with commonly used approaches that directly apply the best parent product for rain/no-rain detection, the superiority of CTC-M is demonstrated analytically and numerically using spatially dense precipitation measurements over Europe. Our analysis also suggests that CTC-M is tolerant to a range of cross-correlated rain/no-rain detection errors and detection biases of the parent products. As a result, CTC-M will benefit global precipitation estimation by improving the representation of precipitation occurrence in gauge-based and multisource merged precipitation products.

Jianzhi Dong↗

Introduction to the JPSS-2 Advanced Technology Microwave Sounder (ATMS) Government Calibration Data Book (GCDB)

The third Advanced Technology Microwave Sounder (ATMS) is an instrument onboard the Joint Polar Satellite System (JPSS), JPSS-2 (renamed NOAA-21 in orbit) mission. This report is to introduce the JPSS-2 Government Calibration Data Book (J2 GCDB) for ATMS, SN 304. This J2 GCDB document contains key information generated during the calibration testing campaign that is driving parameters for radiometric performance. This document also contains supporting data that augments the calibration results. The values in this document are utilized by ATMS’s calibration packet which is, in turn, an integral component in the interpretation of science data. The calibration data in this report was collected from tests such as shelf-level testing, antenna testing, instrument thermal vacuum (TVAC) testing; satellite TVAC testing; and JPSS-2 post-launch tests. JPSS-2 was launched on November 10, 2022. In the subsequent years, the Government will release an ATMS GCDB for each JPSS mission. We expect that all public users can download these ATMS GCDBs from the NOAA operational Integrated Calibration and Validation System (ICVS) website, see more discussions below. The goal of this GCDB is to demonstrate how to characterize ATMS measurements using JPSS-2 ATMS on-orbit operational data and to provide relevant explanations. This document serves as a primary public domain reference for calibrating operational ATMS Raw Data Records (RDR) science data, as used in the current operational Interface Data Processing Segment (IDPS) system. This same RDR science data is distributed through direct broadcast (DB) to DB users for use in their ground processing systems. This J2 GCDB provides the results of the ATMS system radiometric calibration, the antenna flat reflector emissivity [1], the antenna pattern measurements, the antenna pattern corrected brightness temperature [2], the brightness temperature of the lunar disk [3], Lunar Intrusion (LI) correction algorithm [4], receiver spectral parameters, and mechanical alignment on-orbit pointing results, and the striping effect appeared significantly in S-NPP on-orbit radiance data when the data are compared to the Radiative Transfer Model (RTM) simulation in numerical weather prediction (NWP) system [5]. It also provides the parameters required for conversion of telemetry counts to engineering units, for radiometric calibration, and for antenna beam geo-location. Moreover, it provides JPSS-2 ATMS Spectral Response Functions data, some additional information related to ATMS on-orbit performance, on-orbit lunar intrusion correction parameters and Earth contamination bias, and on how to derive ATMS RDR, antenna Temperature Data Records (TDR), and Sensor Data Records (SDR). Furthermore, an introduction of NOAA operational Integrated Calibration and Validation System (ICVS) website and services is added in this J2 GCDB. This ICVS hosts a long-term monitoring system which allows to visualization and comparison of data from JPSS missions, NOAA legacy Polar Operational Environmental Satellites (POES), and Geostationary Operational Environmental Satellites (GOES). From NOAA Comprehensive Large Array-data Stewardship System (CLASS), the public users can download all JPSS ATMS data products for all JPSS missions.

Microwave Sounder↗

A new numerical method for the simulation of three dimensional flow in a pipe

A new numerical technique for simulating three dimensional, unsteady, incompressible pipe flows is presented and its utility and accuracy is shown. Each vector function in the expansion of the velocity field is divergence free and satisfies the boundary conditions for viscous flow. Some of the benefits of the expansion technique are that pressure is eliminated from the dynamics, only two unknowns per mesh point are required, implicit treatment of the viscous terms is provided at no extra computational cost, and no fractional time steps are required. The method uses spectral expansions: Fourier series in the azimuthal and streamwise directions, and Jacobi polynominals in the radial direction.

Leonard, A.↗

A new numerical method for the simulation of three-dimensional flow in a pipe

A new numerical technique for simulating three dimensional, unsteady, incompressible pipe flows is presented and its utility and accuracy is shown. Each vector function in the expansion of the velocity field is divergence free and satisfies the boundary conditions for viscous flow. Some of the benefits of the expansion technique are that pressure is eliminated from the dynamics, only two unknowns per mesh point are required, implicit treatment of the viscous terms is provided at no extra computational cost, and no fractional time steps are required. The method uses spectral expansions: Fourier series in the azimuthal and streamwise directions, and Jacobi polynominals in the radial direction. Previously announced in STAR as N82-31644

Leonard, A.↗

The mesoscale forcing of a midlatitude upper-tropospheric jet streak by a simulated convective system. 2: Kinetic energy and resolution analysis

A kinetic energy (KE) analysis of the forcing of a mesoscale upper-tropospheric jet streak by organized diabatic processes within the simulated convective system (SCS) that was discussed in Part 1 is presented in this study. The relative contributions of the ageostrophic components of motion to the generation of KE of the convectively generated jet streak are compared, along with the KE generation by the rotational (nondivergent) and irrotational (divergent) mass transport. The sensitivity of the numerical simulations of SCS development to resolution is also briefly examined. Analysis within isentropic coordinates provides for an explicit determination of the influence of the diabatic processes on the generation of KE. The upper-level production of specific KE is due predominatly to the inertial advective ageostrophic component (IAD), and as such represents the primary process through which the KE of the convectively generated jet streak is realized. A secondary contribution by the inertial diabatic (IDI) term is observed. Partitioning the KE generation into its rotational and irrotational components reveals that the latter, which is directly linked to the diabatic heating within the SCS through isentropic continuity requirements, is the ultimate source of KE generation as the global area integral of generation by the rotational component vanishes. Comparison with an identical dry simulation reveals that the net generation of KE must be attributed to latent heating. Both the IAD and IDI ageostrophic components play important roles in this regard. Examination of results from simulations conducted at several resolutions supports the previous findings in that the effects of diabatic processes and ageostrophic motion on KE generation remain consistent. Resolution does impact the location and timing of SCS development, a result that has important implications in forecasting the onset of convection that develops from evolution of the large-scale flow and moisture transport. Marked differences are observed in the momentum field aloft subsequent to the life cycle of the SCS in the 1 deg, 30-level base case (MP130) simulation discussed in Part 1 versus its 2 deg counterparts in that the MP130 simulation with higher spatial resolution contains from 14% to 30% more total KE.

Wolf, Bart J.↗

Center for Advanced Computational Technology

The Center for Advanced Computational Technology (ACT) was established to serve as a focal point for diverse research activities pertaining to application of advanced computational technology to future aerospace systems. These activities include the use of numerical simulations, artificial intelligence methods, multimedia and synthetic environments, and computational intelligence, in the modeling, analysis, sensitivity studies, optimization, design and operation of future aerospace systems. The Center is located at NASA Langley and is an integral part of the School of Engineering and Applied Science of the University of Virginia. The Center has four specific objectives: 1) conduct innovative research on applications of advanced computational technology to aerospace systems; 2) act as pathfinder by demonstrating to the research community what can be done (high-potential, high-risk research); 3) help in identifying future directions of research in support of the aeronautical and space missions of the twenty-first century; and 4) help in the rapid transfer of research results to industry and in broadening awareness among researchers and engineers of the state-of-the-art in applications of advanced computational technology to the analysis, design prototyping and operations of aerospace and other high-performance engineering systems. In addition to research, Center activities include helping in the planning and coordination of the activities of a multi-center team of NASA and JPL researchers who are developing an intelligent synthesis environment for future aerospace systems; organizing workshops and national symposia; as well as writing state-of-the-art monographs and NASA special publications on timely topics.

Noor, Ahmed K.↗

Development, Capabilities, and Impact on Wind Analyses of the Hurricane Imaging Radiometer (HIRAD)

The Hurricane Imaging Radiometer (HIRAD) is a new airborne microwave remote sensor for hurricane observations that is currently under development by NASA Marshall Space Flight Center in partnership with the NOAA Atlantic Oceanographic and Meteorological Laboratory/Hurricane Research Division, the University of Central Florida, the University of Michigan, and the University of Alabama in Huntsville. The instrument is being test flown in January and is expected to participate in the tropical cyclone experiment GRIP (Genesis and Rapid Intensification Processes) in the 2010 season. HIRAD is being designed to study the wind field in some detail within strong hurricanes and to enhance the real-time airborne ocean surface winds observation capabilities of NOAA and USAF Weather Squadron hurricane hunter aircraft currently using the operational Stepped Frequency Microwave Radiometer (SFMR). Unlike SFMR, which measures wind speed and rain rate along the ground track at a single point directly beneath the aircraft, HIRAD will provide images of the surface wind and rain field over a wide swath (approximately 3 x the aircraft altitude) with approximately 2 km resolution. This paper describes the HIRAD instrument and the physical basis for its operations, including chamber test data from the instrument. The potential value of future HIRAD observations will be illustrated with a summary of Observing System Simulation Experiments (OSSEs) in which measurements from the new instrument as well as those from existing instruments (air, surface, and space-based) are simulated from the output of a detailed numerical model, and those results are used to construct simulated H*Wind analyses. Evaluations will be presented on the impact on H*Wind analyses of using the HIRAD instrument observations to replace those of the SFMR instrument, and also on the impact of a future satellite-based HIRAD in comparison to instruments with more limited capabilities for observing strong winds through heavy rain. Potential impact on numerical prediction of hurricane intensity will also be discussed.

Miller, T.↗

Development, Capabilities, and Impact on Wind Analyses of the Hurricane Imaging Radiometer (HIRAD)

The Hurricane Imaging Radiometer (HIRAD) is a new airborne microwave remote sensor for hurricane observations that is currently under development by NASA Marshall Space Flight Center in partnership with the NOAA Atlantic Oceanographic and Meteorological Laboratory/Hurricane Research Division, the University of Central Florida, the University of Michigan, and the University of Alabama in Huntsville. The instrument is being test flown in January and is expected to participate in the tropical cyclone experiment GRIP (Genesis and Rapid Intensification Processes) in the 2010 season. HIRAD is being designed to study the wind field in some detail within strong hurricanes and to enhance the real-time airborne ocean surface winds observation capabilities of NOAA and USAF Weather Squadron hurricane hunter aircraft currently using the operational Stepped Frequency Microwave Radiometer (SFMR). Unlike SFMR, which measures wind speed and rain rate along the ground track at a single point directly beneath the aircraft, HIRAD will provide images of the surface wind and rain field over a wide swath (approximately 3 x the aircraft altitude) with approximately 2 km resolution. This paper describes the HIRAD instrument and the physical basis for its operations, including chamber test data from the instrument. The potential value of future HIRAD observations will be illustrated with a summary of Observing System Simulation Experiments (OSSEs) in which measurements from the new instrument as well as those from existing instruments (air, surface, and space-based) are simulated from the output of a detailed numerical model, and those results are used to construct simulated H*Wind analyses. Evaluations will be presented on the impact on H*Wind analyses of using the HIRAD instrument observations to replace those of the SFMR instrument, and also on the impact of a future satellite-based HIRAD in comparison to instruments with more limited capabilities for observing strong winds through heavy rain. Potential impact on numerical prediction of hurricane intensity will also be discussed.

Miller, T.↗

Computation of diffuse sky irradiance from multidirectional radiance measurements

Accurate determination of the diffuse solar spectral irradiance directly above the land surface is important in characterizing the reflectance properties of these surfaces, especially vegetation canopies. This determination is also needed to infer the net radiation budget of the earth-atmosphere system above these surfaces. An algorithm is developed here for the computation of hemispheric diffuse irradiance using the measurements from an instrument called PARABOLA, which rapidly measures upwelling and downwelling radiances in three selected wavelength bands. The validity of the algorithm is established from simulations. The standard reference data set of diffuse radiances of Dave (1978), obtained by solving the radiative transfer equation numerically for realistic atmospheric models, is used to simulate PARABOLA radiances. Hemispheric diffuse irradiance is estimated from a subset of simulated radiances by using the algorithm described. The algorithm is validated by comparing the estimated diffuse irradiance with the true diffuse irradiance of the standard data set. The validations include sensitivity studies for two wavelength bands (visible, 0.65-0.67 micron; near infrared, 0.81-0.84 micron), different atmospheric conditions, solar elevations, and surface reflectances. In most cases the hemispheric diffuse irradiance computed from simulated PARABOLA radiances and the true irradiance obtained from radiative transfer calculations agree within 1-2 percent. This technique can be applied to other sampling instruments designed to estimate hemispheric diffuse sky irradiance.

Ahmad, Suraiya P.↗

Three-dimensional in Situ Observations of Polycrystalline Microstructure Evolution During Directional Solidification of Transparent Alloys Aboard the ISS and Quantitative Comparison with Numerical Modeling

Cellular/dendritic microstructures formed during solidification have a crucial influence on the mechanical properties of a wide range of structural alloys. By minimizing the amount of gravity-induced convection in the liquid, directional solidification experiments using transparent organic alloys conducted in the DECLIC-DSI onboard the International Space Station have provided unique 3D in situ observations of the spatiotemporal evolution of the solid-liquid interface during the formation of cellular and dendritic microstructures under purely diffusive growth conditions in polycrystalline samples containing several grains with a small misorientation with respect to the temperature gradient. Those observations have made it possible to perform benchmark quantitative comparisons with the predictions of state-of-the-art phase-field simulations of microstructure formation in 3D on experimentally relevant length and time scales. This talk will report quantitative comparisons between microgravity experiments and phase-field simulations in succinonitrile-camphor alloys of two different compositions that shed new light on the role of subgrain boundaries in the spatiotemporal evolution of the primary cellular spacing and the selection of dendritic array structures.

Kaihua Ji↗

Assessment of Retrieved GMI Emissivity Over Land, Snow, and Sea Ice in the GEOS System

Directly assimilating microwave radiances over land, snow and sea ice remains a significant challenge for data assimilation systems. These data assimilation systems are critical to the success of global numerical weather prediction systems including the Global Earth Observing System-Atmospheric Data Assimilation System (GEOS-ADAS). Extending more surface sensitive microwave channels over land, snow and ice could provide a needed source of data for Numerical Weather Prediction particularly in the Planetary Boundary Layer (PBL). Unfortunately, the accuracy of emissivity models currently available within the GEOS-ADAS along with other data assimilation systems are insufficient to simulate and assimilate radiances. Recently, Munchak et al. (2020) published a 5-year climatological database for retrieved microwave emissivity from the GPM Microwave Imager (GMI) aboard the Global Precipitation Mission (GPM). In this work the database is utilized by modifying the GEOS-ADAS to use this emissivity database in place of the default emissivity value available in the Community Radiative Transfer Model (CRTM), which is the fast radiative transfer model used by the GEOS-ADAS. As a first step, the GEOS-ADAS is run in a so-called “stand-alone” mode to simulate radiances from GMI using the default CRTM emissivity, and replacing the default CRTM emissivity models with values from Munchak et al, 2020. The simulated observations using Munchak et al., 2020 agree more closely with observations from GMI. These results are presented along with a discussion of the implication for GMI observations within the GEOS-ADAS.

GPM↗

Numerical simulation of crossing/turbulent boundary layer interaction at Mach 8.3 comparison of zero and two-equation turbulence models

A 3D hypersonic crossing shock wave/turbulent boundary layer interaction is examined numerically. The test geometry consists of a pair of opposing sharp fins of angle alpha = 15 deg mounted on a flat plate. The freestream Mach number is 8.28. Two theoretical models are evaluated. The full 3D Reynolds-averaged Navier-Stokes equations are solved using the Baldwin-Lomax algebraic turbulent eddy viscosity model and the Rodi turbulence model. Computed results for both cases show good agreement with experiment for flat plate surface pressure and for pitot pressure and yaw angle profiles in the flowfield. General agreement is obtained for surface flow direction. Fair to poor agreement is obtained for surface heat transfer, indicating a need for more accurate turbulence models. The overall flowfield structure is similar to that observed in previous crossing shock interaction studies.

Narayanswami, N.↗