Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time integrator”

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 217 records · Page 12

Local-in-Time Adjoint-Based Method for Optimal Control/Design Optimization of Unsteady Compressible Flows

.We study local-in-time adjoint-based methods for minimization of ow matching functionals subject to the 2-D unsteady compressible Euler equations. The key idea of the local-in-time method is to construct a very accurate approximation of the global-in-time adjoint equations and the corresponding sensitivity derivative by using only local information available on each time subinterval. In contrast to conventional time-dependent adjoint-based optimization methods which require backward-in-time integration of the adjoint equations over the entire time interval, the local-in-time method solves local adjoint equations sequentially over each time subinterval. Since each subinterval contains relatively few time steps, the storage cost of the local-in-time method is much lower than that of the global adjoint formulation, thus making the time-dependent optimization feasible for practical applications. The paper presents a detailed comparison of the local- and global-in-time adjoint-based methods for minimization of a tracking functional governed by the Euler equations describing the ow around a circular bump. Our numerical results show that the local-in-time method converges to the same optimal solution obtained with the global counterpart, while drastically reducing the memory cost as compared to the global-in-time adjoint formulation.

Yamaleev, N. K.↗

Visible-Infrared Hyperspectral Image Projector

The VisIR HIP generates spatially-spectrally complex scenes. The generated scenes simulate real-world targets viewed by various remote sensing instruments. The VisIR HIP consists of two subsystems: a spectral engine and a spatial engine. The spectral engine generates spectrally complex uniform illumination that spans the wavelength range between 380 nm and 1,600 nm. The spatial engine generates two-dimensional gray-scale scenes. When combined, the two engines are capable of producing two-dimensional scenes with a unique spectrum at each pixel. The VisIR HIP can be used to calibrate any spectrally sensitive remote-sensing instrument. Tests were conducted on the Wide-field Imaging Interferometer Testbed at NASA s Goddard Space Flight Center. The device is a variation of the calibrated hyperspectral image projector developed by the National Institute of Standards and Technology in Gaithersburg, MD. It uses Gooch & Housego Visible and Infrared OL490 Agile Light Sources to generate arbitrary spectra. The two light sources are coupled to a digital light processing (DLP(TradeMark)) digital mirror device (DMD) that serves as the spatial engine. Scenes are displayed on the DMD synchronously with desired spectrum. Scene/spectrum combinations are displayed in rapid succession, over time intervals that are short compared to the integration time of the system under test.

Bolcar, Matthew↗

High-Order Implicit-Explicit Multi-Block Time-stepping Method for Hyperbolic PDEs

This work seeks to explore and improve the current time-stepping schemes used in computational fluid dynamics (CFD) in order to reduce overall computational time. A high-order scheme has been developed using a combination of implicit and explicit (IMEX) time-stepping Runge-Kutta (RK) schemes which increases numerical stability with respect to the time step size, resulting in decreased computational time. The IMEX scheme alone does not yield the desired increase in numerical stability, but when used in conjunction with an overlapping partitioned (multi-block) domain significant increase in stability is observed. To show this, the Overlapping-Partition IMEX (OP IMEX) scheme is applied to both one-dimensional (1D) and two-dimensional (2D) problems, the nonlinear viscous Burger's equation and 2D advection equation, respectively. The method uses two different summation by parts (SBP) derivative approximations, second-order and fourth-order accurate. The Dirichlet boundary conditions are imposed using the Simultaneous Approximation Term (SAT) penalty method. The 6-stage additive Runge-Kutta IMEX time integration schemes are fourth-order accurate in time. An increase in numerical stability 65 times greater than the fully explicit scheme is demonstrated to be achievable with the OP IMEX method applied to 1D Burger's equation. Results from the 2D, purely convective, advection equation show stability increases on the order of 10 times the explicit scheme using the OP IMEX method. Also, the domain partitioning method in this work shows potential for breaking the computational domain into manageable sizes such that implicit solutions for full three-dimensional CFD simulations can be computed using direct solving methods rather than the standard iterative methods currently used.

Nielsen, Tanner B.↗

Solar System Studies with the Space Infrared Telescope Facility (SIRTF)

SIRTF (Space Infrared Telescope Facility) is the final element in NASA's 'Great Observatories' program. It consists of an 85-cm cryogenically-cooled observatory for infrared astronomy from space. SIRTF is scheduled for launch in late 2001 or early 2002 on a Delta rocket into a heliocentric orbit trailing the Earth. Data from SIRTF will be processed and disseminated to the community through the SIRTF Science Center (SSC) located at the Infrared Processing and Analysis Center (IPAC) at Caltech. Some 80/% of the total observing time (estimated at a minimum of 7500 hours of integration time per year for the mission lifetime of about 4 years) will be available to the scientific community at large through a system of refereed proposals. Three basic instruments are located in the SIRTF focal plane. The Multiband Imaging Photometer (MIPS), the Infrared Array Camera (IRAC), and the Infrared Spectrometer (IRS), taken together, provide imaging and spectroscopy from 3.5 to 160 microns. Among the solar system studies suited to SIRTF are the following: 1) spectroscopy and radiometry of small bodies from the asteroid main belt, through the Trojan clouds, to the Kuiper Disk; 2) dust distribution in the zodiacal cloud and the Earth's heliocentric dust ring; 3) spectroscopy and radiometry of comets; and 4) spectroscopy and radiometry of planets and their satellites. Searches for, and studies of dust disks around other stars, brown dwarfs, and superplanets will also be conducted with SIRTF. The SORTIE web site (http://ssc.ipac.caltech.edu/sirtf) contains important details and documentation on the project, the spacecraft, the telescope, instruments, and observing procedures. A community-wide workshop for solar system studies with SIRTF is in the planning stages by the author and Martha S. Hanner for the summer of 1999.

Cruikshank, Dale P.↗

Acoustic-optic spectrometer. 1: Noise contributions and system consideration

An acousto-optic spectrometer (AOS) used as an IF spectrometer to a heterodyne receiver is modeled as a total power multi-channel integrating receiver. Systematic noise contributions common to all total power, time integrating receivers, as well as noise terms unique to the use of optical elements and photo-detectors in an AOS are identified and discussed. In addition, degradation of signal-to-noise ratio of an unbalanced Dicke receiver compared to a balanced Dicke receiver is found to be due to gain calibration processing and is not an instrumental effect.

Chin, G.↗

Development of a Complete Landsat Evapotranspiration and Energy Balance Archive to Support Agricultural Consumptive Water Use Reporting and Prediction in the Central Valley, CA

Mapping evapotranspiration (ET) from agricultural areas in Californias Central Valley is critical for understanding historical consumptive use of surface and groundwater. In addition, long histories of ET maps provide valuable training information for predictive studies of surface and groundwater demands. During times of drought, groundwater is commonly pumped to supplement reduced surface water supplies in the Central Valley. Due to the lack of extensive groundwater pumping records, mapping consumptive use using satellite imagery is an efficient and robust way for estimating agricultural consumptive use and assessing drought impacts. To this end, we have developed and implemented an algorithm for automated calibration of the METRIC remotely sensed surface energy balance model on NASAs Earth Exchange (NEX) to estimate ET at the field scale. Using automated calibration techniques on the NEX has allowed for the creation of spatially explicit historical ET estimates for the Landsat archive dating from 1984 to the near present. Further, our use of spatial NLDAS and CIMIS weather data, and spatial soil water balance simulations within the NEX METRIC workflow, has helped overcome challenges of time integration between satellite image dates. This historical and near present time archive of agricultural water consumption for the Central Valley will be an extremely useful dataset for water use and drought impact reporting, and predictive analyses of groundwater demands.

valleys↗

Public Health Data Applications Using the CDC Tracking Network: Augmenting Environmental Hazard Information with Lower-latency NASA Data

Exposure to environmental hazards is an important determinant of health, and the frequency and severity of exposures is expected to be impacted by climate change. Through a partnership with the U.S. National Aeronautics and Space Administration, the U.S. Centers for Disease Control and Prevention’s National Environmental Public Health Tracking Network is integrating timely observations and model data of priority environmental hazards into its publicly accessible Data Explorer (https://ephtracking.cdc.gov/DataExplorer/). Newly integrated datasets over the contiguous U.S. (CONUS) include: daily 5-day forecasts of air quality based on the Goddard Earth Observing System Composition Forecast (GEOS-CF), daily historical (1980-present) concentrations of speciated PM2.5 based on the Modern Era Retrospective analysis for Research and Applications, version 2 (MERRA-2), and Moderate Resolution Imaging Spectroradiometer (MODIS) daily near real-time maps of flooding (MCDWD). Data integrated into the CDC Tracking Network are broadly intended to improve community health through action by informing both research and early warning activities, including (1) describing temporal and spatial trends in disease and potential environmental exposures, (2) identifying populations most affected, (3) generating hypotheses about associations between health and environmental exposures, and (4) developing, guiding, and assessing environmental public health policies and interventions aimed at reducing or eliminating health outcomes associated with environmental factors.

air quality↗

Stabilization by modification of the Lagrangian

In order to reduce the error growth during a numerical integration, a method for stabilizing the differential equations of Keplerian motion is offered. It is characterized by the use of the eccentric anomaly as independent variable in such a way that the time transformation is given by a generalized Lagrange formalism. The control terms in the equations of motion obtained by this modified Lagrangian give immediately a completely Lyapunov-stable set of differential equations. The equation of time integration is modified by a control term which leads to an integral which defined the time element for the perturbed Keplerian motion.

Baumgarte, J. W.↗

A numerical transport scheme which avoids negative mixing ratios

A new scheme for numerically integrating the transport equation has been developed. This square root method avoids the problem of negative mixing ratios by using the square root of the concentration instead of the concentration itself as an advective time step. Conservation of total mass at every time step is achieved by extending the concept of quadratic conservation to the numerical time integration method. A few time schemes that fulfill the requirement of step-by-step quadratic conservation are described, the simplest of which is the modified Lax-Wendroff method.

Schneider, H.-R.↗

Stabilization by modification of the Lagrangian

In order to reduce the error growth during a numerical integration, a method of stabilization of the differential equations of the Keplerian motion is offered. It is characterized by the use of the eccentric anomaly as an independent variable in such a way that the time transformation is given by a generalized Lagrange formalism. The control terms in the equations of motion obtained by this modified Lagrangian give immediately a completely Liapunov-stable set of differential equations. In contrast to other publications, here the equation of time integration is modified by a control term which leads to an integral which defined the time element for the perturbed Keplerian motion.

Baumgarte, J. W.↗

Full potential integral solution for transonic flows with and without embedded Euler domains

Two methods are presented to solve for the transonic airfoil flow problems. The first method is based on the integral equation solution of the full-potential equation in terms of the velocity field, and a Shock Capturing-Shock Fitting (SCSF) scheme has been developed. The SCSF-scheme consists of a shock-capturing part and a shock-fitting part in which shock panels are introduced at the shock location. The sock panels are fitted and crossed by using the Rankine-Hugoniot relations. The second method is based on coupling the integral equation of the full-potential equation with the pseudo time integration of Euler equations in a small embedded region around the shock. The integral solution provides the initial and boundary conditions for the Euler domain. This scheme is named as the Integral Equation-Embedded Euler (IEEE) scheme. The two methods are applied to NACA 0012 and NACA 64A010A over a wide range of Mach numbers, and the reults are in good agreement with the experimental data and other computational results. The schemes converge within a number of iterations which is one-order of magnitude less than the finite-difference schemes.

Kandil, Osama A.↗

A modified integrated NDVI for improving estimates of terrestrial net primary production

Logic is presented for a time-integrated NDVI that is modified by an AVHRR derived surface evaporation resistance factor sigma, and truncated by temperatures that cause plant dormancy, to improve environmental sensitivity. With this approach, NDVI observed during subfreezing temperatures is not integrated. Water stress-related impairment in plant activity is incorporated by reducing the effective NDVI at each integration with sigma, which is derived from the slope of the surface temperature to NDVI ratio for climatically similar zones of the scene. A comparison of surface resistance before and after an extended drought period for a 1200 sq km region of coniferous forest in Montana is presented.

Running, Steven W.↗

Estimating the Concentration of Large Raindrops from Polarimetric Radar and Disdrometer Observations

Estimation of rainfall integral parameters, including radar observables, and empirical relations between them are sensitive to the truncation of the drop size distribution (DSD), particularly at the large drop end. The sensitivity of rainfall integral parameters to the maximum drop diameter (D(sub max)) is exacerbated at C‐band since resonance effects are pronounced for large drops in excess of 5 mm diameter (D). Due to sampling limitations, it is often difficult to reliably estimate D(sub max) with disdrometers. The resulting uncertainties in D(sub max0 potentially increase errors in radar retrieval methods, particularly at C‐band, that rely on disdrometer observations for DSD input to radar models. In fact, D(sub max) is typically an assumed DSD parameter in the development of radar retrieval methods. Because of these very uncertainties, it is difficult to independently confirm disdrometer estimates of D(sub max) with polarimetric radar observations. A couple of approaches can be taken to reduce uncertainty in large drop measurement. Longer integration times can be used for the collection of larger disdrometer samples. However, integration periods must be consistent with a radar resolution volume (RRV) and the temporal and spatial scales of the physical processes affecting the DSD therein. Multiple co‐located disdrometers can be combined into a network to increase the sample size within a RRV. However, over a reasonable integration period, a single disdrometer sample volume is many orders of magnitudes less than a RRV so it is not practical to devise a network of disdrometers that has an equivalent volume to a typical RRV. Since knowledge of DSD heterogeneity and large drop occurrence in time and space is lacking, the specific accuracy or even general representativeness of disdrometer based D(sub max) and large drop concentration estimates within a RRV are currently unknown. To address this complex issue, we begin with a simpler question. Is the frequency of occurrence of large rain drops (D > 5 mm) in disdrometer observations, either stand alone or networked, generally representative and consistent with polarimetric radar observations? We first show from simulations that the concentration of large (D > 5 mm) rain drops (N(sub T5)) can be estimated from polarimetric observations of specific differential phase (K(sub dp)) and differential reflectivity (Z(sub dr)), N(sub T5)=F(K(sub dp),Z(sub dr)), or horizontal reflectivity (Z(sub h)) and Z(sub dr), N(sub T5)=(Z(sub h),Z(sub dr)). We assess the error associated with polarimetric retrieval of N(sub T5), including sensitivity to D(sub max) parameterization assumptions and measurement error in the radar simulations. Polarimetric measurements at S‐band and C‐band will then be used to retrieve estimates of N(sub T5) and compared to disdrometer estimates of N(sub T5). After careful consideration of retrieval error, we will check consistency between disdrometer and polarimetric radar estimates of N(sub T5) and the frequency of occurrence of large rain drops in a variety of precipitating regimes using data from NASA's Global Precipitation Measurement (GPM) Ground Validation (GV) program, including field campaigns such as MC3E (Oklahoma) and IFloodS (Iowa) and extended measurements over Huntsville, Alabama and NASA Wallops Flight Facility in coastal Virginia.

Carey, Lawrence D.↗

A New Formulation of Time Domain Boundary Integral Equation for Acoustic Wave Scattering in the Presence of a Uniform Mean Flow

It has been well-known that under the assumption of a constant uniform mean flow, the acoustic wave propagation equation can be formulated as a boundary integral equation, in both the time domain and the frequency domain. Compared with solving partial differential equations, numerical methods based on the boundary integral equation have the advantage of a reduced spatial dimension and, hence, requiring only a surface mesh. However, the constant uniform mean flow assumption, while convenient for formulating the integral equation, does not satisfy the solid wall boundary condition wherever the body surface is not aligned with the uniform mean flow. In this paper, we argue that the proper boundary condition for the acoustic wave should not have its normal velocity be zero everywhere on the solid surfaces, as has been applied in the literature. A careful study of the acoustic energy conservation equation is presented that shows such a boundary condition in fact leads to erroneous source or sink points on solid surfaces not aligned with the mean flow. A new solid wall boundary condition is proposed that conserves the acoustic energy and a new time domain boundary integral equation is derived. In addition to conserving the acoustic energy, another significant advantage of the new equation is that it is considerably simpler than previous formulations. In particular, tangential derivatives of the solution on the solid surfaces are no longer needed in the new formulation, which greatly simplifies numerical implementation. Furthermore, stabilization of the new integral equation by Burton-Miller type reformulation is presented. The stability of the new formulation is studied theoretically as well as numerically by an eigenvalue analysis. Numerical solutions are also presented that demonstrate the stability of the new formulation.

Hu, Fang↗

Finite element computation of a viscous compressible free shear flow governed by the time dependent Navier-Stokes equations

A finite element algorithm for solution of fluid flow problems characterized by the two-dimensional compressible Navier-Stokes equations was developed. The program is intended for viscous compressible high speed flow; hence, primitive variables are utilized. The physical solution was approximated by trial functions which at a fixed time are piecewise cubic on triangular elements. The Galerkin technique was employed to determine the finite-element model equations. A leapfrog time integration is used for marching asymptotically from initial to steady state, with iterated integrals evaluated by numerical quadratures. The nonsymmetric linear systems of equations governing time transition from step-to-step are solved using a rather economical block iterative triangular decomposition scheme. The concept was applied to the numerical computation of a free shear flow. Numerical results of the finite-element method are in excellent agreement with those obtained from a finite difference solution of the same problem.

Cooke, C. H.↗

Three-dimensional time-marching aeroelastic analyses using an unstructured-grid Euler method

The modifications to a 3D implicit, upwind, unstructured-grid Euler code for aeroelastic analysis of complete aircraft configurations described involve the addition of structural equations of motion for their simultaneous time integration with the governing flow equations. A detailed description of the time-marching aeroelastic procedure and comparisons with experimental data are presented, to provide an assessment of capabilities. Flutter results are shown for an isolated 45-deg sweptback wing and an SST with a fuselage, clipped delta wing, and two identical rearward-mounted nacelles. Comparisons between computed and experimental flutter characteristics show good agreement.

Rausch, Russ D.↗

Derivation of a New Smoke Emissions Inventory using Remote Sensing, and Its Implications for Near Real-Time Air Quality Applications

A new emissions inventory of particulate matter (PM) is being derived mainly from remote sensing data using fire radiative power (FRP) and aerosol optical depth (AOD) retrievals from the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument, as well as wind data from the Modern Era Retrospective-Analysis for Research and Applications (MERRA) reanalysis dataset, which spans the satellite era. This product is generated using a coefficient of emission, C(sub e), that has been produced on a 1x1 degree global grid such that, when it is multiplied with satellite measurements of FRP or its time-integrated equivalent fire radiative energy (FRE) retrieved over a given area and time period, the corresponding PM emissions are estimated. This methodology of using C(sub e) to derive PM emissions is relatively new and advantageous for near real-time air quality applications compared to current methods based on post-fire burned area that may not provide emissions in a timely manner. Furthermore, by using FRP to characterize a fire s output, it will represent better accuracy than the use of raw fire pixel counts, since fires in individual pixels can differ in size and strength by orders of magnitude, resulting in similar differences in emission rates. Here we will show examples of this effect and how this new emission inventory can properly account for the differing emission rates from fires of varying strengths. We also describe the characteristics of the new emissions inventory, and propose the process chain of incorporating it into models for air quality applications.

Ellison, Luke↗

Multigrid Approach to Incompressible Viscous Cavity Flows

Two-dimensional incompressible viscous driven-cavity flows are computed for Reynolds numbers on the range 100-20,000 using a loosely coupled, implicit, second-order centrally-different scheme. Mesh sequencing and three-level V-cycle multigrid error smoothing are incorporated into the symmetric Gauss-Seidel time-integration algorithm. Parametrics on the numerical parameters are performed, achieving reductions in solution times by more than 60 percent with the full multigrid approach. Details of the circulation patterns are investigated in cavities of 2-to-1, 1-to-1, and 1-to-2 depth to width ratios.

Wood, William A.↗