SEARCH · Engineering Papers

Results for “numerical modelling”

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,657 records · Page 92

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Theoretical overview and modeling of the sodium and potassium atmospheres of mercury

A general theoretical overview for the sources, sinks, gas-surface interactions, and transport dynamics of sodium and potassium in the exospheric atmsophere of Mercury is given. Information for these four factors, which control the spatial distribution of these two alkali-group gases about the planet, is incorporated in numerical models. The spatial nature and relative importance of the initial source atom atmosphere and the ambient (ballistic hopping) atom atmosphere are then examined and are shown to be controlled and coupled to a great extent by the extremely large and variable solar radiation acceleration experienced by sodium and potassium as they resonantly scatter solar photons. The lateral (antisunward) transport rate of thermally accommodated sodium and potassium ambient atoms is shown to be driven by the solar radiation acceleration and, over a significant portion of Mercury's orbit about the Sun, is sufficiently rapid to be competitive with the short photoionization lifetimes for these atoms when they are located on the summit surface near or within about 30 deg of the terminator. The lateral transport rate is characterized by a migration time determined by model calculations for an ensemble of atoms initially starting at a point source on the surface (i.e., a numerical spacetime dependent Green's function). Four animations for the spacetime evolution of the sodium (or potassium) atmosphere produced by a point source on the surface are presented on a videotape format. For extended surface sources for sodium and potassium, the local column density is determined by competition between the photoionization lifetimes and the lateral transport times of atoms originating from different surface source locations. Sodium surface source fluxes (referenced to Mercury at perihelion) that are required on the sunlit hemisphere to reproduce the typically observed several megarayleighs of D2 emission-line brightness and the inferred column densities of 1-2 x 10(exp 11) atoms per sq cm range from approximately 2-5 x 10(exp 7) atoms/sq cm/sec. The sodium model is applied to study observational data that document an anticorrelation in the average sodium column density and solar radiation acceleration. Lateral transport driven by the solar radiation acceleration is shown to produce this behavior for combinations of different sources and surface accomodation coefficients. The best fit model fits to the observational data require a significant degree of thermal accommodation of the ambient sodium atoms to the surface and a source rate that decreases as an inverse power of 1.5 to 2 in heliocentric distance.

Smyth, William H.↗

DOE OSTI2020

Clouds and Radiation Processes in Regional Climate Models Evaluated Using Observations Over the Ice‐free Arctic Ocean

Abstract The presence of clouds in the Arctic regulates the surface energy budget (SEB) over the sea‐ice surface and the ice‐free ocean. Following several previous field campaigns, the cloud‐radiation relationship, including cloud vertical structure and phase, has been elucidated; however, modeling of this relationship has matured slowly. In recognition of the recent decline in the Arctic sea‐ice extent, representation of the cloud system in numerical models should consider the effects of areas covered by sea ice and ice‐free areas. Using an in situ stationary meteorological observation data set obtained over the ice‐free Arctic Ocean by the Japanese Research Vessel Mirai (September 2014), coordinated evaluation of six regional climate models (RCMs) with nine model runs was performed by focusing on clouds and the SEB. The most remarkable findings were as follows: (1) reduced occurrence of unstable stratification with low‐level cloud water in all models in comparison to the observations, (2) significant differences in cloud water representations between single‐ and double‐moment cloud schemes, (3) extensive differences in partitioning of hydrometeors including solid/liquid precipitation, and (4) pronounced lower‐tropospheric air temperature biases. These issues are considered as the main sources of SEB uncertainty over ice‐free areas of the Arctic Ocean. The results from a coupled RCM imply that the SEB is constrained by both the atmosphere and the ocean (and sea ice) with considerable feedback. Coordinated improvement of both stand‐alone atmospheric and coupled RCMs would promote a more comprehensive and improved understanding of the Arctic air‐ice‐sea coupled system.

Inoue, Jun↗

DOE OSTI2022

Convective Updraft and Downdraft Characteristics of Continental Mesoscale Convective Systems in the Model Gray Zone

Here, the "gray zone" of convective modeling is defined as the range of horizontal grid spacings (Δx) at which turbulent transport processes are only partially resolved by the dynamics of the numerical model. This zone typically covers Δx from a few kilometers to several hundred meters, wherein the realistic representation of convective cloud processes can be challenging. This study characterizes the convective draft behaviors at multiple Δx across the gray zone and determines the appropriate Δx that can reliably capture these salient convective properties. We perform an ensemble of idealized simulations of mesoscale convective systems (MCS) using the Weather Research and Forecasting model at various Δx from 4 km to 250 m over the central U.S. An evaluation of key MCS kinematic properties is constrained using unique, long-term vertical velocity estimates obtained by radar wind profilers deployed by the Department of Energy Atmospheric Radiation Measurement user facility. MCS simulations for all Δx tested overestimate (underestimate) the probabilities of convective updrafts (downdrafts) compared to the observations. In terms of the convective draft intensity, finer-Δx models overestimate the updraft intensity, while the opposite is found for downdrafts. Moving from Δx = 4 km to 250 m, downdrafts become stronger and more frequent especially at middle and upper levels, attributed to additional drag from an increasing graupel frequency and compensation of enhanced updraft velocity. Simulated draft characteristics, including core size, intensity, and probability of occurrence, exhibit pronounced changes at Δx greater than 500 m, but suggest less sensitivity when Δx is reduced below 500 m.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

A Methodology for Estimating the Energy and Moisture Budget of the Convective Boundary Layer Using Continuous Ground-Based Infrared Spectrometer Observations

Abstract Land–atmosphere interactions play a critical role in both the atmospheric water and energy cycles. Changes in soil moisture and vegetation alter the partitioning of surface water and energy fluxes, influencing diurnal evolution of the planetary boundary layer (PBL). The mixing-diagram framework has proven useful in understanding the evolution of the heat and moisture budget within the convective boundary layer (CBL). We demonstrate that observations from the Department of Energy Atmospheric Radiation Measurement (ARM) Southern Great Plains site provide all of the needed inputs needed for the mixing-diagram framework, allowing us to quantify the impact from the surface fluxes, advection, radiative heating, encroachment, and entrainment on the evolution of the CBL. Profiles of temperature and humidity retrieved from the ground-based infrared spectrometer [Atmospheric Emitted Radiance Interferometer (AERI)] are a critical component in this analysis. Large-eddy simulation results demonstrate that mean mixed-layer values derived are shown to be critical to close the energy and moisture budgets. A novel approach demonstrated here is the use of network of AERIs and Doppler lidars to quantify the advective fluxes of heat and moisture. The framework enables the estimation of the entrainment fluxes as a residual, providing a way to observe the entrainment fluxes without using multiple lidar systems. The high temporal resolution of the AERI observations enables the morning, midday, and afternoon evolution of the CBL to be quantified. This work provides a new way to use observations in this framework to evaluate weather and climate models. Significance Statement The energy and moisture budget of the planetary boundary layer (PBL) is influenced by multiple sources, and accurately representing this evolution in numerical models is critical for weather forecasts and climate predictions. The mixing-diagram approach, driven by profiling observations as illustrated here, provides a powerful way to quantify the contributions from each of these sources. In particular, the energy and moisture mixed into the PBL from above the PBL can be determined accurately from ground-based remote sensors using this approach.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Vadose Zone Flow and Transport Parameters Data Package for the Hanford Site Composite Analysis and Cumulative Impact Evaluation

This report provides a description of the basis for the development and implementation of a conceptual model for vadose zone flow and transport for the composite analysis (CA) groundwater pathway analysis and the cumulative impact evaluation (CIE). The parameterization for a numerical model is intimately linked to the conceptual model framework. The report describes the basis for the selection of hydraulic and transport parameters for the hydrostratigraphic units (HSUs) identified in the 200 East and 200 West Areas. Whenever data are sparse or unavailable, surrogate hydraulic properties are chosen based on samples collected within the 200 Areas and nearby locations that are representative of sediments characteristic of the HSUs identified elsewhere.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS1975

A Method for Calculating the Induced Pressure Distribution Associated with a Jet in a Crossflow

A model is presented which can be used to study the loss of lift during hovering and horizontal flight of the VTOL aircraft. The model numerically predicts the pressure distribution induced by a round, turbulent, unheated, subsonic jet exhausting normally through a flat plate into a subsonic crossflow. The complete model assumes that the predominant features of the flow are jet entrainment and a pair of contrarotating vortices which form downstream of the jet. Experimentally determined vortex properties and a reasonable assumption concerning jet entrainment were used. Potential flow considerations were used except in the wake region, where a simple method for approximating the pressure distribution was suggested. The calculated pressure distribution, lift, and pitching moments on the flat plate are presented for a jet to crossflow velocity ratio of 8 and were compared with experimental results. A computer program is given which was used to calculate the pressure distribution across the flat plate.

Dietz, W. E., Jr.↗

NASA NTRS1976

Surface properties of ocean fronts

Background information on oceanic fronts is presented and the results of several models which were developed to study the dynamics of oceanic fronts and their effects on various surface properties are described. The details of the four numerical models used in these studies are given in separate appendices which contain all of the physical equations, program documentation and running instructions for the models.

Wolff, P. M.↗

NASA NTRS1979

Large eddy simulation - A predictive approach to turbulent flow computation

The filtering approach (Leonard, 1974) to large eddy simulation is described. Large eddy simulation is then surveyed with reference to subgrid scale models, numerical methods, channel flow, and free shear flows. It is concluded that the premise of LES (i.e., that the small eddies are much easier to model than the large ones) seems to be borne out by computational experience to date. The greatest obstacle in LES at present is the lack of computer power which effectively denies the ability to simulate all the significant eddies in flows of interest.

Ferziger, J. H.↗

NASA NTRS1988

Master oscillator power amplifier performance of Ti:Al2O3

An improved Ti:Al2O3 amplifier model is reported. Amplifier performance under various pump and probe fluence conditions is also reported. It is claimed to be in good agreement with the numerical model predictions.

Barnes, James C.↗

NASA NTRS1991

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program is to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative, manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, our immediate emphasis is on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we have developed unique global models that allow us to study the coupling between the different regions. These results are highlighted. Another important aspect of our NASA Theory Program concerns the effect that localized structure has on the macroscopic flow in the ionosphere, plasmasphere, thermosphere and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkeland current patterns) or time variations in these inputs due to storms and substorms. Also, some of the plasma flows that we predict with our macroscopic models may be unstable. Another one of our goals is to examine the stability of our predicted flows. Because time-dependent three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulations). Therefore, another long-range goal of our NASA Theory Program is to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This may involve a detailed comparison of kinetic, semikinetic, and hydrodynamic predictions for a given polar wind scenario or it may involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations will provide insight into when the various models can be used with confidence.

Schunk, Robert W.↗

NASA NTRS1992

Dilatation-dissipation corrections for advanced turbulence models

This paper analyzes dilatation-dissipation based compressibility corrections for advanced turbulence models. Numerical computations verify that the dilatation-dissipation corrections devised by Sarkar and Zeman greatly improve both the k-omega and k-epsilon model predicted effect of Mach number on spreading rate. However, computations with the k-gamma model also show that the Sarkar/Zeman terms cause an undesired reduction in skin friction for the compressible flat-plate boundary layer. A perturbation solution for the compressible wall layer shows that the Sarkar and Zeman terms reduce the effective von Karman constant in the law of the wall. This is the source of the inaccurate k-gamma model skin-friction predictions for the flat-plate boundary layer. The perturbation solution also shows that the k-epsilon model has an inherent flaw for compressible boundary layers that is not compensated for by the dilatation-dissipation corrections. A compressibility modification for k-gamma and k-epsilon models is proposed that is similar to those of Sarkar and Zeman. The new compressibility term permits accurate predictions for the compressible mixing layer, flat-plate boundary layer, and a shock separated flow with the same values for all closure coefficients.

Wilcox, David C.↗

NASA NTRS1992

Collisional simulations of particles in a narrow planetary ring

A model is used to investigate how particle-particle collisions might modify some consequences of satellite perturbations relevant to short-term features of narrow planetary rings. A Monte Carlo-like simulation ring-model particle velocity alteration allows the simulation of collisions while individual particles are tracked. Periodic features visible in Voyager images are reproduced in a 2D numerical model which incorporates the collisional simulation as well as the perturbations of adjacent satellites. Collisions are noted to wash out periodic features within one collisional relaxation time.

Kolvoord, Robert A.↗

NASA NTRS1997

Global Transport of Aerosol and CO: Initial 3-D Simulations of MAPS, TOMS, and AVHRR

Carbon monoxide concentrations and aerosol properties provide the tracers of global tropospheric perturbation by humankind that are most easily observed from satellite platforms. Aircraft field observations are additionally needed for us to simulate and understand the patterns observed. We report on our accumulating experience in making detailed, situation-specific, 3-D simulations of these tropospheric constituents as observed from the MAPS CO sensor, the Advanced Very High Resolution Radiometer (AVHRR) aerosol scattering data, and the Total Ozone Mapping Spectrometer (TOMS) absorbing-aerosol information. We have found a strong complimentarity of satellite and aircraft data in this effort, and recommend modelers use: (1) satellite data to pose initial questions; (2) aircraft data for modelers' first detailed simulations; and (3) a return to satellite data for global generalization. This work provides one example. Our main tools are the MM5 numerical model for meteorological assimilation and GRACES, our NASA Ames tracer-chemistry model, which incorporates emissions estimates. We report on several six-week simulations of global biomass burning effects as observed in the MAPS October 1994 dataset, and show the usefulness of two aircraft datasets, the TRACE-A (1992) and PEM-Tropics missions of NASA's Global Tropospheric Experiment. Additional information is contained in the original extended abstract.

Chatfield, Robert B.↗

NASA NTRS2003

The Effect of Simulated Microgravity Environment of RWV Bioreactors on Surface Reactions and Adsorption of Serum Proteins on Bone-bioactive Microcarriers

Biomimetically modified bioactive materials with bone-like surface properties are attractive candidates for use as microcarriers for 3-D bone-like tissue engineering under simulated microgravity conditions of NASA designed rotating wall vessel (RWV) bioreactors. The simulated microgravity environment is attainable under suitable parametric conditions of the RWV bioreactors. Ca-P containing bioactive glass (BG), whose stimulatory effect on bone cell function had been previously demonstrated, was used in the present study. BG surface modification via reactions in solution, resulting formation of bone-like minerals at the surface and adsorption of serum proteins is critical for obtaining the stimulatory effect. In this paper, we report on the major effects of simulated microgravity conditions of the RWV on the BG reactions surface reactions and protein adsorption in physiological solutions. Control tests at normal gravity were conducted at static and dynamic conditions. The study revealed that simulated microgravity remarkably enhanced reactions involved in the BG surface modification, including BG dissolution, formation of bone-like minerals at the surface and adsorption of serum proteins. Simultaneously, numerical models were developed to simulate the mass transport of chemical species to and from the BG surface under normal gravity and simulated microgravity conditions. The numerical results showed an excellent agreement with the experimental data at both testing conditions.

Radin, Shula↗

NASA NTRS2011

Shunted Piezoelectric Vibration Damping Analysis Including Centrifugal Loading Effects

Excessive vibration of turbomachinery blades causes high cycle fatigue problems which require damping treatments to mitigate vibration levels. One method is the use of piezoelectric materials as passive or active dampers. Based on the technical challenges and requirements learned from previous turbomachinery rotor blades research, an effort has been made to investigate the effectiveness of a shunted piezoelectric for the turbomachinery rotor blades vibration control, specifically for a condition with centrifugal rotation. While ample research has been performed on the use of a piezoelectric material with electric circuits to attempt to control the structural vibration damping, very little study has been done regarding rotational effects. The present study attempts to fill this void. Specifically, the objectives of this study are: (a) to create and analyze finite element models for harmonic forced response vibration analysis coupled with shunted piezoelectric circuits for engine blade operational conditions, (b) to validate the experimental test approaches with numerical results and vice versa, and (c) to establish a numerical modeling capability for vibration control using shunted piezoelectric circuits under rotation. Study has focused on a resonant damping control using shunted piezoelectric patches on plate specimens. Tests and analyses were performed for both non-spinning and spinning conditions. The finite element (FE) shunted piezoelectric circuit damping simulations were performed using the ANSYS Multiphysics code for the resistive and inductive circuit piezoelectric simulations of both conditions. The FE results showed a good correlation with experimental test results. Tests and analyses of shunted piezoelectric damping control, demonstrating with plate specimens, show a great potential to reduce blade vibrations under centrifugal loading.

Min, James B.↗

DOE OSTI2023

Wind Turbine Noise Code Benchmark: A Comparison and Verification Exercise: Preprint

In a number of institutions and companies, researchers and engineers are developing modeling numerical frameworks that are used to predict the aerodynamic noise emission from a wind turbine rotor. The simulation codes vary from empirically-tuned engineering models to high-fidelity computational ones. Their common feature is the fact that they all specifically model the main aerodynamic noise mechanisms occurring at the rotating blades, namely trailing edge noise and turbulent inflow noise. Nevertheless, it can be expected that these different modeling frameworks will produce different results for a same given rotor design, and identical operating conditions. Trailing edge noise is put at the forefront of the study as it is recognized as the main audible source of noise from wind turbines. The present benchmark aims at comparing the results from the different modeling approaches, and drawing some conclusions from these comparisons. This effort, denoted as Wind Turbine Noise Code benchmark, was initiated in 2019 as a joint-activity between the IEA Wind Task 39 (Quiet Wind Turbine Technology) and Task 29 (Detailed Aerodynamics of Wind Turbines, now Task 47). In addition to the investigation of the noise emissions themselves, the rotor aerodynamic characteristics are investigated as they have a significant impact on the noise generation mechanisms. A number of test cases are defined and the aerodynamic and aeroacoustic predictions from the various models are compared. There exist some discrepancies between the different methods, but it is difficult to conclude if one methodology is better than another in order to design a wind turbine with noise as a constraint.

aeroacoustics↗

DOE OSTI2021

Modeling Sea Ice Effects for Wave Energy Resource Assessments

Wave-generated power has potential as a valuable coastal resource, but the wave climate needs to be mapped for feasibility before wave energy converters are installed. Numerical models are used for wave resource assessments to quantify the amount of available power and its seasonality. Alaska is the U.S. state with the longest coastline and has extensive wave resources, but it is affected by seasonal sea ice that dampens the wave energy and the full extent of this dampening is unknown. To accurately characterize the wave resource in regions that experience seasonal sea ice, coastal wave models must account for these effects. The aim of this study is to determine how the dampening effects of sea ice change wave energy resource assessments in the nearshore. Here, we show that by combining high-resolution sea ice imagery with a sea ice/wave dampening parameterization in an unstructured grid, the Simulating Waves Nearshore (SWAN) model improves wave height predictions and demonstrates the extent to which wave power decreases when sea ice is present. The sea ice parametrization decreases the bias and root mean square errors of wave height comparisons with two wave buoys and predicts a decrease in the wave power of up to 100 kW/m in areas around Prince William Sound, Alaska. The magnitude of the improvement of the model/buoy comparison depends on the coefficients used to parameterize the wave–ice interaction.

16 TIDAL AND WAVE POWER↗

NASA NTRS1989

Model atmospheres with periodic shocks

The pulsation of a long-period variable star generates shock waves which dramatically affect the structure of the star's atmosphere and produce conditions that lead to rapid mass loss. Numerical modeling of atmospheres with periodic shocks is being pursued to study the processes involved and the evolutionary consequences for the stars. It is characteristic of these complex dynamical systems that most effects result from the interaction of various time-dependent processes.

Bowen, G. H.↗