Search NASA⌕ Search

SEARCH · Search NASA

Results for “tangent linear”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Understanding and Utilizing PBL Height Data from Multiple Observing Systems in the GEOS System

The accuracy of PBL height simulation is a key issue in many applications including forecasting near surface meteorology and air quality, however, it is a very challenging problem due to the lack of not only comprehensive, global Planetary Boundary Layer (PBL) observations but also a strategy and infrastructure to utilize PBL height data from a variety of sensors. Following the designation of PBL as an incubation class observable in the 2017 Decadal Survey, the PBL Incubation Study Team Report [14] made clear that “a future global PBL observing system requires modeling and data assimilation as essential components.” There is an urgent need for global modeling development in order to utilize Program of Record (POR) observations, assess their impacts, and identify gaps to be filled by future PBL missions. Our overall objective is to develop PBL data assimilation capabilities in the NASA Global Earth Observing System (GEOS), focusing on PBL height from multiple observing systems, to support the assessment and use of future PBL observations. The NASA GEOS system is composed of the GEOS global atmospheric general circulation model (AGCM) and the atmospheric data assimilation system (ADAS). The PBL parameterizations include the “Lock” K-profile scheme driven by surface and cloud-top buoyancy fluxes ([4]), and the “Louis” local scheme for stable conditions based on the Richardson number ([5]). Above the mixed layer defined by the Lock surface plume, shallow cumulus convection is represented by the mass flux scheme of [9]. Additional parameterizations are summarized in [1]. The ADAS employs the hybrid 4D Ensemble- Variational (EnVar) configuration ([15]), with the ensemble providing flow-dependent background error covariance information. The resultant analysis increments are fed back to the forecast model through the 4D incremental analysis update (IAU) approach ([11]). In this study, PBL height data are being or have been generated from radiosondes, GNSS RO, satellite (CATS, CALIPSO and ICESat-2) and ground-based (MPLNET) lidars, and wind profiler. Investigations have been conducted to specify quality marks for PBL height retrievals for the data assimilation purpose. These PBL height data have different strengths and weaknesses ([2], [3], [6], [7], [8], [10]), and the satellite PBL height data provide better global coverage and complement in-situ PBL height data. Radiosondes offer high accuracy and in situ measurement of temperature and humidity profiles, but with poor spatio-temporal sampling. The in-situ observing systems like MPLNET and wind profiler provide long history of PBL height records at each station. The GNSS RO based PBL height is retrieved based on the sharp gradients in refractivity profile that represent the fine vertical structure of temperature and moisture changes above the PBL. However, not all RO refractivity profiles reach the surface depending on location and regime, and RO refractivity retrievals can be negatively biased below 2km. The PBL height data from satellite lidars provide high resolution along track PBL height retrievals, but over land they are affected by previous day convective PBL aerosol and strongly associated with mixing layer and retrievals cannot be made below thick, attenuating clouds. A successful assimilation of PBL height data requires a thorough understanding of the observing method and the retrieval algorithm for each observing system in order to use the PBL height data from multiple observing systems properly. Due to the sensitivity of PBL height data to the observing method and choice of algorithm, it is important to use a model definition appropriate for each observation type to compute differences between PBL height data and model PBL height (OmFs). The GEOS model currently includes two PBL height definitions suitable for direct comparison with observed PBL height, and additional definitions are being added in this study. Evaluation of different model PBL height definitions is underway. Meanwhile, efforts have been made in the GEOS data assimilation system to develop PBL height data assimilation capability. PBL height data can be assimilated using two different approaches. The traditional approach is to construct an observation operator and its tangent linear and adjoint, which link control variables to PBL height data from each observing system. This observation operator can be very complicated, e.g., the lidar-based PBL height observation operator includes the backscatter lidar forward observation operator, the algorithm to derive PBL height from attenuated total backscatter, interpolation, and calculations handling the mismatch between observed and model scales. The other approach is to augment PBL height to the control variable vector, and it is adopted in this study. The latter approach was also used in previous studies, e.g., the assimilation of PBL height data from radiosonde and aircraft in the Real Time Mesoscale Analysis (RTMA) system for a dispersion modelling study ([13]); the PBL height assimilation study using lidar PBL height data at Greensburg, Kansas for a field campaign ([12]). The PBL height assimilation from multiple observing systems in this study allows us to take advantage of the diverse PBL height data that provide much better global coverage collectively under different meteorological conditions and with different temporal and spatial scales. As all the PBL heights are tightly coupled with the PBL thermodynamic variables, the strong correlations, which are provided by the 4D ensemble forecast, enable PBL height data from various sources to interact and combine coherently and provide additional information for PBL temperature and moisture fields. The results of comparisons among PBL height data from different sources and the evaluation of the model PBL height definitions with the PBL height data will be presented, and the PBL height data synergy strategies and preliminary results will also be discussed at the conference.

Y. Zhu↗
NASA NTRSDate not supplied

Preliminary Results Cycling GEOS-JEDI with GSI-based Background Errors

The first phase of transitioning the NASA GMAO GEOS atmospheric data assimilation capabilities to JEDI involves the replacement of the Grid-point Statistical Interpolation (GSI) with a corresponding JEDI analysis. This includes taking JEDI's Unified Observation Operator (UFO), its underlying dependencies, and the JEDI solver that enables a hybrid 4DEnVar strategy similar to what is used in the current GEOS-GSI system. Variational analysis involves at least two main components associated with the observation and background cost function terms. The first is directly related to the UFO, which is being carefully validated in a joint collaboration between GMAO and NCEP to demonstrate consistency with corresponding observations usage in GSI. The second component is the background term, which in a hybrid system involves the ability to set up both a climatologically-based term and an ensemble-based term. JEDI provides the means to implement both terms through its BUMP component. Use of BUMP would require a complete re-tune of both climatological and ensemble, which is a non-trivial exercise we would prefer to avoid. As an alternative, the work here studies the results of interfacing the GSI-background error capability (GSIBEC) into JEDI through SABER. With this, the exact same background error covariance formulation used in GSI can be employed in JEDI without need for re-tuning. This brief summary covers the work done to interface GSIBEC into JEDI and shows preliminary results where the background error covariances of the control (GEOS-GSI) and experiment (GEOS-JEDI) are identical in corresponding cycling experiments. The cycling exercise is obviously preliminary and so much can be expected from GEOS-JEDI when compared to GEOS-GSI. There is still a number of features that need closer attention and although in some cases in principle ready to cycle have been intentionally either turned off or not fully exercised (e.g., VarBC is applied but not cycled). Other features are still pending implementation, one such example is the implementation of the Tangent Linear Normal Mode Constraint. Still, results are quite encouraging as hopefully the discussion here illustrates.

Ricardo Todling↗
NASA NTRS2023

Developing a Radar Signal Simulator for the Community Radiative Transfer Model

Active radar instruments provide vertically resolved clouds and precipitation measurements that cannot be provided by the passive instruments. These active measurements are not conventionally assimilated into the data assimilation systems because of the lack of fast forward radiative transfer models and also difficulties in the error modelling of the measurements. This paper describes the development, evaluation, and sensitivity analysis for a forward radar model implemented in the Community Radiative Transfer Model (CRTM). The scattering properties required by the forward model are provided by the hydrometeor lookup tables that were generated using the discrete dipole approximation. The model is able to calculate both the reflectivity and the attenuated reflectivity for any given radar instrument at any given zenith angles as long as CRTM instrument specific coefficients are available. The evaluation using CloudSat measurements shows a very good agreement between the simulations and measurements as long as the input profiles of hydrometeors are consistent with the measured reflectivity profiles. Major sources contributing to the differences between the measured and simulated reflectivities are input hydrometeor profiles, scattering lookup tables, lack of melting layer in the forward model, CRTM scattering solvers, and attenuation calculations. In addition to the forward model, both Tangent Linear and Adjoint of the model are also implemented and tested within CRTM. These components may be required by some data assimilation systems for the assimilation of radar measurements.

radar↗
NASA NTRSDate not supplied

A Comprehensive Forward Model for Spaceborne Radar Instruments

We present the development and validation of a comprehensive forward model designed to enhance remote sensing capabilities of spaceborne radar instruments. To overcome limitations in existing models, we integrated a Discrete Dipole Approximation (DDA) cloud scattering database into our Radiative Transfer Model (RTM), focusing on microwave frequencies. By simulating the optical properties of non-spherical frozen hydrometeors, the DDA technique effectively reduced discrepancies between simulated and observed values, surpassing traditional Mie tables. The evaluation of DDA lookup tables involved comparisons with a collocated dataset comprising short-term forecasts and satellite microwave data, providing evidence of their superiority. Additionally, we address the challenges of assimilating active radar measurements, which offer vertically resolved insights into clouds and precipitation. We explored the assimilation of spaceborne radar measurements in Numerical Weather Prediction (NWP) models by integrating a forward radar model, along with its adjoint and tangent linear, into the data assimilation system. Evaluation using CloudSat measurements demonstrated promising agreement between simulations and observations, particularly when the input hydrometeor profiles aligned with the measured reflectivity profiles, showcasing the potential of the developed forward radar model. Moreover, we discuss other challenges in radar measurement assimilation within NWP models, including potential observation errors and biases.

Isaac Moradi↗
NASA NTRS2000

Arc-Tangent Circuit for Continuous Linear Output

A device suitable for determining arc-tangent of an angle θ is provided. Circuitry generates a first square wave at a frequency ωt and a second square wave at the frequency ωt but shifted by a phase difference equal to the angle θ. A pulse width modulation signal generator processes the first and second square waves to generate a pulse width modulation signal having a frequency of rot and having a pulse width that is a function of the phase difference θ. The pulse width modulation signal is converted to a DC voltage that is a linear representation of the phase difference θ.

Dean C Alhorn↗
NASA NTRS2014

Recent Updates to the GEOS-5 Linear Model

Global Modeling and Assimilation Office (GMAO) is close to having a working 4DVAR system and has developed a linearized version of GEOS-5.This talk outlines a series of improvements made to the linearized dynamics, physics and trajectory.Of particular interest is the development of linearized cloud microphysics, which provides the framework for 'all-sky' data assimilation.

Adjoint↗
NASA NTRS1983

Linear approximations of nonlinear systems

The development of a method for designing an automatic flight controller for short and vertical take off aircraft is discussed. This technique involves transformations of nonlinear systems to controllable linear systems and takes into account the nonlinearities of the aircraft. In general, the transformations cannot always be given in closed form. Using partial differential equations, an approximate linear system called the modified tangent model was introduced. A linear transformation of this tangent model to Brunovsky canonical form can be constructed, and from this the linear part (about a state space point x sub 0) of an exact transformation for the nonlinear system can be found. It is shown that a canonical expansion in Lie brackets about the point x sub 0 yields the same modified tangent model.

Hunt, L. R.↗
NASA NTRS1983

Linear approximations of nonlinear systems

A method for designing an automatic flight controller for short and vertical takeoff aircraft is presently being developed at NASA Ames Research Center. This technique involves transformations of nonlinear systems to controllable linear systems and takes into account the nonlinearities of the aircraft. In general, the transformations cannot always be given in closed form. Using partial differential equations, an approximate linear system, called the modified tangent model, was recently introduced. A linear transformation of this tangent model to Brunovsky canonical form can be constructed, and from this an approximation of an exact transformation for the nonlinear system can be found. It is shown that a canonical expansion in Lie brackets about the point x(0) yields the same modified tangent model.

Hunt, L. R.↗
NASA NTRS1980

Computation of the shock-wave boundary layer interaction with flow separation

The boundary layer concept is used to describe the flow near the wall. The external flow is approximated by a pressure displacement relationship (tangent wedge in linearized supersonic flow). The boundary layer equations are solved in finite difference form and the question of the presence and unicity of the solution is considered for the direct problem (assumed pressure) or converse problem (assumed displacement thickness, friction ratio). The coupling algorithm presented implicitly processes the downstream boundary condition necessary to correctly define the interacting boundary layer problem. The algorithm uses a Newton linearization technique to provide a fast convergence.

Ardonceau, P.↗
NASA NTRS1990

Particle orbits in two-dimensional equilibrium models for the magnetotail

Assuming that there exist an equilibrium state for the magnetotail, particle orbits are investigated in two-dimensional kinetic equilibrium models for the magnetotail. Particle orbits in the equilibrium field are compared with those calculated earlier with one-dimensional models, where the main component of the magnetic field (Bx) was approximated as either a hyperbolic tangent or a linear function of z with the normal field (Bz) assumed to be a constant. It was found that the particle orbits calculated with the two types of models are significantly different, mainly due to the neglect of the variation of Bx with x in the one-dimensional fields.

Karimabadi, H.↗
NASA NTRS1982

An assessment of preliminary aerodynamic analysis methods for supersonic speeds including high angle-of-attack flow

The problems associated with predicting the aerodynamic loads for moderate supersonic Mach number, high angle-of-attack flows and for high supersonic Mach number flows are addressed. For these large disturbance-type flows, commonly used prediction methods based on linearized theory or on the tangent-wedge method are shown to be generally inadequate. For the preliminary analysis of aerodynamic shapes subjected to these type flows, two methods have recently been developed by combining features of the linearized theory and tangent-wedge method. These two new methods are discussed, and through extensive comparisons with experimental data, are shown to offer an improved prediction capability.

Pittman, J. L.↗
NASA NTRS1977

Assessment of existing analytic methods for prediction of high angle-of-attack loads on delta wings at supersonic speeds

An assessment of the applicability of four loading prediction methods to high angle-of-attack conditions for simplified wing-body configurations is provided. The methods are: The tangent wedge approximation, the linear theory methods of Middleton and Woodward, and a shock-fitting finite-difference technique. Estimates obtained by these methods were compared with experimental pressure data on delta wings to examine the effects of Mach number, camber, sweep angle, and angle of attack. Results indicate that all of the methods provided reasonable estimates at moderate angles of attack. At these moderate angles of attack, the methods of Middleton and Woodward provided good estimates at Mach numbers higher than those usually associated with linear theory. Only the finite-difference method provided reasonable load estimates at high angles of attack.

Landrum, E. J.↗
NASA NTRS1990

Numerical simulations of compressible mixing layers

Direct numerical simulations of two-dimensional temporally growing compressible mixing layers are presented. The Kelvin-Helmholtz instability is initially excited by a white-noise perturbation superimposed onto a hyperbolic tangent meanflow profile. The linear regime is studied at low resolution in the case of two flows of equal temperatures, for convective Mach numbers from 0.1 to 1 and for different values of the Reynolds number. At higher resolution, the complete evolution of a two-eddy mixing layer between two flows of different temperatures is simulated at moderate Reynolds number. Similarities and differences between flows of equal convective Mach numbers are discussed.

Normand, Xavier↗
NASA NTRS2006

Interferometer for Low-Uncertainty Vector Metrology

A simplified schematic diagram of a tilt-sensing unequal-path interferometer set up to measure the orientation of the normal vector of one surface of a cube mounted on a structure under test is described herein. This interferometer has been named a "theoferometer" to express both its interferometric nature and the intention to use it instead of an autocollimating theodolite. The theoferometer optics are mounted on a plate, which is in turn mounted on orthogonal air bearings for near-360 rotation in azimuth and elevation. Rough alignment of the theoferometer to the test cube is done by hand, with fine position adjustment provided by a tangent arm drive using linear inchwormlike motors.

Toland, Ronald W.↗
NASA NTRS1977

Lie theoretic aspects of the Riccati equation

Various features of the application of Lie theory to matrix Riccati equations, of basic importance in control and system theories, are discussed. Particular consideration is given to centralizer foliation, the Cartan decomposition, matrix Riccati equations as Lie systems on Grassmanians, local analysis near a zero point of a vector field, linearization in homogeneous space, the tangent bundle in terms of partitioned matrices, and stability properties of fixed points of Riccati vector fields.

Hermann, R.↗
NASA NTRS2011

Current Sheets in the Heliosheath: Voyager 1, 2009

We identified all of the current sheets for which we have relatively complete and accurate magnetic field (B) data from Voyager 1 (V1) from days of year (DOYs) 1 to 331, 2009, which were obtained deep in the heliosheath between 108.5 and 111.8 AU. Three types of current sheets were found: (1) 15 proton boundary layers (PBLs), (2) 10 and 3 magnetic holes and magnetic humps, respectively, and (3) 3 sector boundaries. The magnetic field strength changes across PBL, and the profile B(t) is linearly related to the hyperbolic tangent function, but the direction of B does not change. For each of the three sector boundaries, B rotated in a plane normal to the minimum variance direction, and the component of B along the minimum variance direction was zero within the uncertainties, indicating that the sector boundaries were tangential discontinuities. The structure of the sector boundaries was not as simple as that for PBLs. The average thickness of magnetic holes and humps (approx.30 RL) was twice that of the PBLs (approx.15 RL). The average thickness of the current sheets associated with sector boundaries was close to the thickness of the PBLs. Our observations are consistent with the hypothesis that magnetic holes and humps are solitons, which are initiated by the mirror mode instability, and evolve by nonlinear kinetic plasma processes to pressure balanced structures maintained by magnetization currents and proton drift currents in the gradients of B.

Burlaga, L. F.↗
NASA NTRS1992

Dielectric constant of liquid alkanes and hydrocarbon mixtures

The complex dielectric constants of n-alkanes with two to seven carbon atoms have been measured. The measurements were conducted using a slotted-line technique at 1.2 GHz and at atmospheric pressure. The temperature was varied from the melting point to the boiling point of the respective alkanes. The real part of the dielectric constant was found to decrease with increasing temperature and correlate with the change in the molar volume. An upper limit to all the loss tangents was established at 0.001. The complex dielectric constants of a few mixtures of liquid alkanes were also measured at room temperature. For a pentane-octane mixture the real part of the dielectric constant could be explained by the Clausius-Mosotti theory. For the mixtures of n-hexane-ethylacetate and n-hexane-acetone the real part of the dielectric constants could be explained by the Onsager theory extended to mixtures. The dielectric constant of the n-hexane-acetone mixture displayed deviations from the Onsager theory at the highest fractions of acetone. The dipole moments of ethylacetate and acetone were determined for dilute mixtures using the Onsager theory and were found to be in agreement with their accepted gas-phase values. The loss tangents of the mixtures exhibited a linear relationship with the volume fraction for low concentrations of the polar liquids.

NASA Program Exobiology↗