Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical forward model”

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 145 records · Page 8

Modeling and numerical simulation of microwave pulse propagation in air breakdown environment

It is shown that ionization occurs wherever the field intensity of the pulse exceeds the local breakdown threshold field of the background air. The produced plasma then attenuates the pulse and gives rise to a tail erosion phenomenon that plays the primary role in limiting the energy transfer of the pulse from source to destination. A theoretical model describing the propagation of an intense microwave pulse in an air breakdown environment is developed that includes the possible focusing effect introduced by either using phase array antennas or the other arrangements. The self-consistent description of the propagation process is provided by a set of two modal equations. These include a continuity equation (Poynting's equation) for the energy density of the pulse and a rate equation of the electron density. A forward wave approximation is used to simplify Poynting's equation, and a semiempirical formula is used for the ionization frequency, nu(sub i). This frequency provides the coupling between the two modal equations, and is used to express the electron rate equation explicitly. In terms of the relevant parameters of the atmosphere, these two equations are normalized for numerical analysis of pulse propagation in the atmosphere. The dependencies of the propagation characteristics of the pulse on intensity, frequency, width, and shape of the pulse are determined. The numerical simulations lead to a useful empirical relation p(exp 3)w = alpha = constant, where p and w are the incident power and width of the pulse and alpha depends on the percentage of pulse energy transferred from the source point to a destined position. The density distribution of the pulse's self-generated plasma is also evaluated. The results also show that for ionization caused by a single unfocused microwave pulse transmitted upwards from the ground, the maximum electron density produced at, for example, 50 km altitude is limited by the tail erosion effect to below 10(exp 6) cm(exp -3). Repetitive pulse and focused beam approaches are also examined. Both approaches can increase the maximum electron density by no more than an order of magnitude. A scheme using two obliquely propagating pulses intersecting at the destined height, e.g. 50 km, is considered. It is shown that the electron density generated at the lowest intersecting position can easily reach a value of 6.6 x 10(exp 8) cm(exp -3), which is considered to be high enough for artificial ionospheric mirror (AIM) application.

Kim, J.↗

Propagation of propeller tone noise through a fuselage boundary layer

In earlier experimental and analytical studies, it was found that the boundary layer on an aircraft could provide significant shielding from propeller noise at typical transport airplane cruise Mach numbers. In this paper a new three-dimensional theory is described that treats the combined effects of refraction and scattering by the fuselage and boundary layer. The complete wave field is solved by matching analytical expressions for the incident and scattered waves in the outer flow to a numerical solution in the boundary layer flow. The model for the incident waves is a near-field frequency-domain propeller source theory developed previously for free field studies. Calculations for an advanced turboprop (Prop-Fan) model flight test at 0.8 Mach number show a much smaller than expected pressure amplification at the noise directivity peak, strong boundary layer shielding in the forward quadrant, and shadowing around the fuselage. Results are presented showing the difference between fuselage surface and free-space noise predictions as a function of frequency and Mach number. Comparison of calculated and measured effects obtained in a Prop-Fan model flight test show good agreement, particularly near and aft of the plane of rotation at high cruise Mach number.

Hanson, D. B.↗

Integral Velocity Sampling for Unsteady Rotor Models on Cartesian Meshes

An unsteady body-force rotor model utilizing integral velocity sampling is developed for Cartesian meshes with a focus on accurate performance predictions. In this approach, the rotors are modeled as body-force source terms, and integral sampling is used to define the local freestream velocity vector. The method is verified in two dimensions for airfoils, providing accurate angle of attack and freestream Mach number estimates even on relatively coarse meshes, and demonstrating excellent mesh convergence across a variety of source region thicknesses. Extensions of the model to three dimensions for both fixed and rotating wings are considered, with numerical experiments demonstrating the need for a spanwise-continuous representation of the induced velocity or angle of attack. The model provides reasonable spanwise lift distributions without any explicit tip-loss model. Comparisons to high-fidelity simulations and experimental data show excellent agreement for predictions of an isolated DJI Phantom 3 propeller in hover. Finally, simulations of the full aircraft in edgewise forward flight also compare well with wind tunnel data across a variety of pitch angles and propeller speeds.

ARMD↗

Model Simulations of CO2 Transport Using Assimilated Meteorological Fields

The numerical simulation of CO2 transport (and other tracers such as CO, CH4, and biomass burning tracers) in the atmosphere is required to determine the fate of anthropogenic source gases. Estimation of the CO2 exchange between the ocean surface, the terrestrial biosphere, and the atmosphere is of first-order importance to understanding the global carbon cycle and the processes that are most crucial in determining the atmospheric CO2 concentration. Forward transport simulations have been conducted using two-dimensional, time-dependent grids of average surface fluxes (from TRANSCOM) and three-dimensional wind data from a prototype data assimilation system (FV-DAS) run by the Goddard Data Assimilation Office. The objective is to better understand the contribution of meteorological variability to changes in CO2 and other constituents, By accurately accounting for meteorological variability, through use of assimilated winds, we hope to better characterize the distribution of surface sources and sinks (and chemistry where applicable). With assimilated meteorology such chemistry/transport runs provide the basic framework to analyze existing (and proposed) measurement data on a point-by-point, real-time basis. We compare with measured CO2 concentration gradients on a daily, seasonal, regional, and interhemispheric basis to examine the consistency of sources, sinks, and transport formulation. We will also examine the inter-annual variability of atmospheric CO2 due to atmospheric circulation changes using longer runs with assimilated winds.

Kaea, S. R.↗

initMIP-Antarctica: an Ice Sheet Model Initialization Experiment of ISMIP6

Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMIP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMIP-Greenland, initMIP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMIP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue.

ice sheet model↗

An on-line expert system for diagnosing environmentally induced spacecraft anomalies using CLIPS

A new rule-based, expert system for diagnosing spacecraft anomalies is under development. The knowledge base consists of over two-hundred rules and provide links to historical and environmental databases. Environmental causes considered are bulk charging, single event upsets (SEU), surface charging, and total radiation dose. The system's driver translates forward chaining rules into a backward chaining sequence, prompting the user for information pertinent to the causes considered. The use of heuristics frees the user from searching through large amounts of irrelevant information (varying degrees of confidence in an answer) or 'unknown' to any question. The expert system not only provides scientists with needed risk analysis and confidence estimates not available in standard numerical models or databases, but it is also an effective learning tool. In addition, the architecture of the expert system allows easy additions to the knowledge base and the database. For example, new frames concerning orbital debris and ionospheric scintillation are being considered. The system currently runs on a MicroVAX and uses the C Language Integrated Production System (CLIPS).

Lauriente, Michael↗

Flame Spread Across Liquids: Experimental Results

The goal of our research on flame spread across a pool of liquid fuel is the quantitative identification of the mechanisms that control the rate and nature of flame spread when the initial temperature of the liquid pool is below the fuel's flash point temperature. Besides numerous experiments in drop towers and 1 g laboratories, we have flown five microgravity (mu-g) experiments on sounding rockets. As described in earlier papers, the first three flights examined the effect of forced opposed airflow over a 2.5 cm deep x 2 cm wide x 30 cm long pool of 1-butanol in mu-g. It was found that the flame spread is much slower and steadier than in 1 g where flame spread has a pulsating character. It was speculated that the flame spread in mu-g resembled the character of pseudo-uniform spread in 1 g; Ito et al later confirmed this conclusively in 1 g experiments. Much of the mu-g flame is also farther from the surface, dimmer, and with less soot, when compared to the 1 g flame. Three-dimensional liquid-phase flow patterns that control the liquid preheating were discovered in both 1 g and mu-g. Our numerical model, restricted to two dimensions, had predicted faster, pulsating flame spread in mu-g for opposed airflow. In examining the differences in the dimensionality of the model and experiment, it was noted that the experiment allowed gas expansion in the lateral direction (across the width of the pool), for which the model could not account. Such lateral expansion could reduce the expansion in the forward and upward directions. Because only these latter directions could be modeled, it was decided to artificially reduce the gas thermal expansion in the predictions. When this was done, satisfactory agreement could be obtained between the predicted and observed spread rates and the steadiness of the spread in microgravity. In 1 g, however, the predicted flame spread character also changed to pseudo-uniform, which disagreed with our 1 g experiments where the spread is pulsating. It was then speculated that gas-phase buoyant flow might oppose the lateral gas expansion, so that the 1 g experiments retained their pulsating flame spread character. If this speculation was valid, a difference in lateral gas expansion should be observable when comparing 1 g and mu-g experiments. Specifically, it was anticipated that greater flow divergence caused by lateral expansion would be measured in mu-g in the absence of a buoyant flow directed towards the flame.

Ross, H. D.↗

Three-dimensional Finite Element Formulation and Scalable Domain Decomposition for High Fidelity Rotor Dynamic Analysis

This paper has two objectives. The first objective is to formulate a 3-dimensional Finite Element Model for the dynamic analysis of helicopter rotor blades. The second objective is to implement and analyze a dual-primal iterative substructuring based Krylov solver, that is parallel and scalable, for the solution of the 3-D FEM analysis. The numerical and parallel scalability of the solver is studied using two prototype problems - one for ideal hover (symmetric) and one for a transient forward flight (non-symmetric) - both carried out on up to 48 processors. In both hover and forward flight conditions, a perfect linear speed-up is observed, for a given problem size, up to the point of substructure optimality. Substructure optimality and the linear parallel speed-up range are both shown to depend on the problem size as well as on the selection of the coarse problem. With a larger problem size, linear speed-up is restored up to the new substructure optimality. The solver also scales with problem size - even though this conclusion is premature given the small prototype grids considered in this study.

Datta, Anubhav↗

Toward Coupled Data Assimilation in NASA’s GEOS: Developments in the Ocean Context

The Global Modeling & Assimilation Office (GMAO) at NASA GSFC produces analyses and predictions of the Earth system using various configurations of the Goddard Earth Observing System (GEOS) model and assimilation system. The current sub-seasonal-to-seasonal prediction system (GEOS-S2S) is based on a coupled atmosphere-ocean-land-ice configuration of GEOS which includes the Modular Ocean Model version 5 (MOM5) run at approximately 50-km resolution and a de-coupled OI-based ocean analysis that uses an initialization of MOM5 forced by the MERRA-2 reanalysis. GMAO will soon implement an updated GEOS-S2S system that will run at 25-km resolution and adopt aspects of the hybrid four-dimensional ensemble-variational (H4DEnVar) system already running in the production-version atmospheric analysis system, including a Local Ensemble Transform Kalman Filter (LETKF) to provide initial conditions for the oceanic state. This presentation will focus on developments to sustain the GMAO's systems on longer time horizons, where more radical transformations will be required to adapt to advanced computing environments, higher resolution and more diverse model components, and new observations for the Earth system. Results will describe progress toward a version of the GEOS coupled system that will be based around the Joint Effort for Data assimilation Integration (JEDI) framework being developed within Joint Center for Satellite Data Assimilation (JCSDA) and include an updated ocean model, MOM6. Discussion will focus specifically on the use of a Unified Forward Operator (UFO) for simulating observations and the Object Oriented Prediction System (OOPS) for providing the state estimate. These features are being developed as a multi-agency effort under the auspices of the JCSDA and are being adopted in the GMAO for all its applications of coupled data assimilation including S2S, numerical weather prediction, and reanalysis.

Mahajan, Rahul↗

TPSAS-NF1676L-11045-DND

The International TOVS Working Group (ITWG) is convened as a sub-group of the International Radiation Commission (IRC) of the International Association of Meteorology and Atmospheric Physics (IAMAP). The ITWG continues to organize international TOVS study conferences which have met approximately every 18 months since 1983. The most recent conference, the seventeenth International TOVS Study Conference (ITSC-XVII) was held in Monterey, California from 12 to 20 April 2010. This report will cover the current ITWG status, objectives, activities, along with recommendations formulated during the recent ITSC conference to be forwarded to space agencies, operational NWP centers and the scientific community on issues ranging from data processing methods, derived products, and the impacts of radiances and inferred atmospheric temperature and moisture fields on numerical weather prediction, and weather and climate studies. Also to be reported on are activities of the technical sub-groups which meet informally to coordinate ATOVS processing software, radiative transfer models, sounding data for climate studies, use of sounding data in data assimilation/NWP, international issues, and future systems and frequency protection issues relevant to ATOVS. As the result of ITSC-XVII, a paper report and electronic proceedings are being formulated and will soon be published and distributed. The conference Working Group Report will summarize the recommendations and actions of these working sub-groups. Technical Proceedings of the scientific presentations and posters will also published. The ITWG web site (http://cimss.ssec.wisc.edu/itwg/) contains electronic versions of the conference papers, presentations and posters from earlier ITSCs. Together, these documents and web pages reflect the conduct of highly successful international collaborations.

Steve English↗

A differential game solution to the Coplanar tail-chase aerial combat problem

Numerical results obtained in a simplified version of the one on one aerial combat problem are presented. The primary aim of the data is to specify the roles of pursuer and evader as functions of the relative geometry and of the significant physical parameters of the problem. Numerical results are given in a case in which the slower aircraft is more maneuverable than the faster aircraft. A third order dynamic model of the relative motion is described, for which the state variables are relative range, bearing, and heading. The ranges at termination are arbitary in the present version of the problem, so the weapon systems of both aircraft can be visualized as forward firing high velocity weapons, which must be aimed at the tail pipe of the evader. It was found that, for the great majority of the ralative geometries, each aircraft can evade the weapon system of the other.

Merz, A. W.↗

An implicit method for the nonlinear modelling and simulation of piezoceramic actuators displaying hysteresis

Experiments have shown that piezoceramic materials display a nonlinear relationship between the applied electric field and the actuation strain. This relationship also displays a substantial hysteresis upon reversal of the applied field. In this paper, piezoceramic actuator models are incorporated into the structural equations of motion to arrive at a set of nonlinear actively controlled structural equations of motion. A new implicit algorithm for determining the time history of the actively controlled structure is presented. The algorithm employs the trapezoidal rule for stepping the equations forward in time. The algorithm is compared to an explicit algorithm and is shown to provide greater numerical accuracy. In addition, the numerical stability and convergence characteristics are presented via example. The new algorithm displays convergence and stability properties that are comparable to the standard trapezoidal rule used for dynamic systems without hysteresis.

Leigh, Timothy D.↗

Study of curved and planar frequency-selective surfaces with nonplanar illumination

A locally planar technique (LPT) is investigated for determining the forward-scattered field from a generally shaped inductive frequency-selective surface (FSS) with nonplanar illumination. The results of an experimental study are presented to assess the LPT accuracy. The effects of a nonplanar incident field are determined by comparing the LPT numerical results with a series of experiments with the feed source placed at varying distances from the planar FSS. The limitations of the LPT model due to surface curvature are investigated in an experimental study of the scattered fields from a set of hyperbolic cylinders of different curvatures. From these comparisons, guidelines for applying the locally planar technique are developed.

Caroglanian, Armen↗

Interaction of an interplanetary shock with the heliospheric plasma sheet

Interplanetary shocks often propagate along the heliospheric plasma sheet (HPS) where the interplanetary magnetic field (IMF) changes its polarity. This problem is investigated by the time-dependent 2.5-D MHD numerical model in the meridional plane. An example of computation is shown in the figure using density (log) contours and IMF vectors. Values of plasma parameters along the HPS fluctuate in time due to the Kelvin-Helmholtz instability. The HPS with its decreased intensity of the IMF as well as with its increased mass density causes a dimple in the shock structure (relatively weak for the forward shock, significant for the reverse shock, and very large for the contact discontinuity). Beyond the forward shock, the HPS is slightly compressed due to the post-shock increase of the azimuthal IMF component. Then follows expansion of the HPS surrounded by the highly-deformed contact discontinuity. A significant draping of IMF lines occurs around this structure that increases the meridional component of the IMF. This can cause a favorable condition for initiation of a geomagnetic storm.

Odstrcil, D.↗

A Comparative Study of Three Methodologies for Modeling Dynamic Stall

During the past two decades, there has been an increased reliance on the use of computational fluid dynamics methods for modeling rotors in high speed forward flight. Computational methods are being developed for modeling the shock induced loads on the advancing side, first-principles based modeling of the trailing wake evolution, and for retreating blade stall. The retreating blade dynamic stall problem has received particular attention, because the large variations in lift and pitching moments encountered in dynamic stall can lead to blade vibrations and pitch link fatigue. Restricting to aerodynamics, the numerical prediction of dynamic stall is still a complex and challenging CFD problem, that, even in two dimensions at low speed, gathers the major difficulties of aerodynamics, such as the grid resolution requirements for the viscous phenomena at leading-edge bubbles or in mixing-layers, the bias of the numerical viscosity, and the major difficulties of the physical modeling, such as the turbulence models, the transition models, whose both determinant influences, already present in static maximal-lift or stall computations, are emphasized by the dynamic aspect of the phenomena.

Sankar, L.↗

GC13I-0857: Designing a Frost Forecasting Service for Small Scale Tea Farmers in East Africa

Kenya is the third largest tea exporter in the world, producing 10% of the world's black tea. Sixty percent of this production occurs largely by small scale tea holders, with an average farm size of 1.04 acres, and an annual net income of $1,075. According to a recent evaluation, a typical frost event in the tea growing region causes about $200 dollars in losses which can be catastrophic for a small holder farm. A 72-hour frost forecast would provide these small-scale tea farmers with enough notice to reduce losses by approximately 80 USD annually. With this knowledge, SERVIR, a joint NASA-USAID initiative that brings Earth observations for improved decision making in developing countries, sought to design a frost monitoring and forecasting service that would provide farmers with enough lead time to react to and protect against a forecasted frost occurrence on their farm. SERVIR Eastern and Southern Africa, through its implementing partner, the Regional Centre for Mapping of Resources for Development (RCMRD), designed a service that included multiple stakeholder engagement events whereby stakeholders from the tea industry value chain were invited to share their experiences so that the exact needs and flow of information could be identified. This unique event allowed enabled the design of a service that fit the specifications of the stakeholders. The monitoring service component uses the MODIS Land Surface Temperature product to identify frost occurrences in near-real time. The prediction component, currently under testing, uses the 2-m air temperature, relative humidity, and 10-m wind speed from a series of high-resolution Weather Research and Forecasting (WRF) numerical weather prediction model runs over eastern Kenya as inputs into a frost prediction algorithm. Accuracy and sensitivity of the algorithm is being assessed with observations collected from the farmers using a smart phone app developed specifically to report frost occurrences, and from data shared through our partner network developed at the stakeholder engagement meeting. This presentation will illustrate the efficacy of our frost forecasting algorithm, and a way forward for incorporating these forecasts in a meaningful way to the key decision makers - the small-scale farmers of East Africa.

frost↗

Two space scatterer formalism calculation of bulk parameters of thunderclouds

In a previous study, we used a modified two-space scatterer formalism of Twersky to establish for a cloud modeled as a statistically homogeneous distribution of spherical water droplets, the dispersion relations that determine its bulk propagation numbers and bulk indexes of refraction in terms of the vector equivalent scattering amplitude and the dyadic scattering amplitude of the single water droplet in isolation. The results were specialized to the forward direction of scattering while demanding that the scatterers preserve the incident polarization. We apply this approach to obtain specific numerical values for the macroscopic parameters of the cloud. We work with a cloud of density rho = 100 cm(exp -3), a wavelength lambda = 0.7774 microns, and with spherical water droplets of common radius alpha = 10 microns. In addition, the scattering medium is divided into three parts, the medium outside the cloud, moist air (the medium inside the cloud but outside the droplets), and the medium inside the spherical water droplets. The results of this report are applicable to a cloud of any geometry since the boundary does not interfere with the calculations. Also, it is important to notice the plane wave nature of the incidence wave in the moist atmosphere.

Phanord, Dieudonne D.↗

Technical Challenges Associated with In-Air Wingtip Docking of Aircraft in Forward Flight

Autonomous in-air wingtip docking of aircraft offers significant opportunity for system level performance gains for numerous aircraft applications. Several of the technical challenges facing wingtip docking of fixed-wing aircraft are addressed in this paper, including: close proximity aerodynamic coupling; mechanisms and operations for robust docking; and relative state estimation methods. A simulation framework considering the aerodynamics, rigid-body dynamics, and vehicle controls is developed and used to perform docking sensitivity studies for a system of two 5.5% scale NASA Generic Transport Model aircraft. Additionally, proof of- concept testing of a candidate docking mechanism designed to move the primary wingtip vortex inboard suggests the viability of such an approach for achieving robust docking.

Quinlan, Jesse R.↗