Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational methods”

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 361 records · Page 20

Computational Relativistic Astrophysics Using the Flowfield-Dependent Variation Theory

Theoretical models, observations and measurements have preoccupied astrophysicists for many centuries. Only in recent years, has the theory of relativity as applied to astrophysical flows met the challenges of how the governing equations can be solved numerically with accuracy and efficiency. Even without the effects of relativity, the physics of magnetohydrodynamic flow instability, turbulence, radiation, and enhanced transport in accretion disks has not been completely resolved. Relativistic effects become pronounced in such cases as jet formation from black hole magnetized accretion disks and also in the study of Gamma-Ray bursts (GRB). Thus, our concern in this paper is to reexamine existing numerical simulation tools as to the accuracy and efficiency of computations and introduce a new approach known as the flowfield-dependent variation (FDV) method. The main feature of the FDV method consists of accommodating discontinuities of shock waves and high gradients of flow variables such as occur in turbulence and unstable motions. In this paper, the physics involved in the solution of relativistic hydrodynamics and solution strategies of the FDV theory are elaborated. The general relativistic astrophysical flow and shock solver (GRAFSS) is introduced, and some simple example problems for Computational Relativistic Astrophysics (CRA) are demonstrated.

Richardson, G. A.↗

Parallelization of the Flow Field Dependent Variation Scheme for Solving the Triple Shock/Boundary Layer Interaction Problem

A parallelized version of the Flowfield Dependent Variation (FDV) Method is developed to analyze a problem of current research interest, the flowfield resulting from a triple shock/boundary layer interaction. Such flowfields are often encountered in the inlets of high speed air-breathing vehicles including the NASA Hyper-X research vehicle. In order to resolve the complex shock structure and to provide adequate resolution for boundary layer computations of the convective heat transfer from surfaces inside the inlet, models containing over 500,000 nodes are needed. Efficient parallelization of the computation is essential to achieving results in a timely manner. Results from a parallelization scheme, based upon multi-threading, as implemented on multiple processor supercomputers and workstations is presented.

Schunk, Richard Gregory↗

Atmospheric studies with the Multi-angle Imaging SpectroRadiometer (MISR)

Multi-angle Imaging SpectroRadiometer observations facilitate new atmospheric remote sensing methods. Angular variation of radiometric and geometric scene attributes, in conjunction with new algorithms, enables retrieval of aerosol properties, cloud heights, and cloud-traced winds.

MISR atmospheric studies↗

Verification and Validation of Requirements on the CEV Parachute Assembly System Using Design of Experiments

The Crew Exploration Vehicle Parachute Assembly System (CPAS) project conducts computer simulations to verify that flight performance requirements on parachute loads and terminal rate of descent are met. Design of Experiments (DoE) provides a systematic method for variation of simulation input parameters. When implemented and interpreted correctly, a DoE study of parachute simulation tools indicates values and combinations of parameters that may cause requirement limits to be violated. This paper describes one implementation of DoE that is currently being developed by CPAS, explains how DoE results can be interpreted, and presents the results of several preliminary studies. The potential uses of DoE to validate parachute simulation models and verify requirements are also explored.

Schulte, Peter Z.↗

Robotics

Earth's upper atmosphere is an extreme environment: dry, cold, and irradiated. It is unknown whether our aerobiosphere is limited to the transport of life, or there exist organisms that grow and reproduce while airborne (aerophiles); the microenvironments of suspended particles may harbor life at otherwise uninhabited altitudes[2]. The existence of aerophiles would significantly expand the range of planets considered candidates for life by, for example, including the cooler clouds of a hot Venus-like planet. The X project is an effort to engineer a robotic exploration and biosampling payload for a comprehensive survey of Earth's aerobiology. While many one-shot samples have been retrieved from above 15 km, their results are primarily qualitative; variations in method confound comparisons, leaving such major gaps in our knowledge of aerobiology as quantification of populations at different strata and relative species counts[1]. These challenges and X's preliminary solutions are explicated below. X's primary balloon payload is undergoing a series of calibrations before beginning flights in Spring 2012. A suborbital launch is currently planned for Summer 2012. A series of ground samples taken in Winter 2011 is being used to establish baseline counts and identify likely background contaminants.

Rothschild, Lynn J.↗

Climate Change and Sounder Radiometric Stability

Satellite instrument radiometric stability is critical for climate studies. The Atmospheric Infrared Sounder (AIRS) radiances are of sufficient stability and accuracy to serve as a climate data record as evidenced by comparisons with the global network of buoys. In this paper we examine the sensitivity of derived geophysical products to potential instrument radiometric stability issues due to diurnal, orbital and seasonal variations. Our method is to perturb the AIRS radiances and examine the impact to retrieved parameters. Results show that instability in retrieved temperature products will be on the same order of the brightness temperature error in the radiances and follow the same time dependences. AIRS excellent stability makes it ideal for examining impacts of instabilities of future systems on geophysical parameter performance.

Sounding↗

Temporal Gain Correction for X-Ray Calorimeter Spectrometers

Calorimetric X-ray detectors are very sensitive to their environment. The boundary conditions can have a profound effect on the gain including heat sink temperature, the local radiation temperature, bias, and the temperature of the readout electronics. Any variation in the boundary conditions can cause temporal variations in the gain of the detector and compromise both the energy scale and the resolving power of the spectrometer. Most production X-ray calorimeter spectrometers, both on the ground and in space, have some means of tracking the gain as a function of time, often using a calibration spectral line. For small gain changes, a linear stretch correction is often sufficient. However, the detectors are intrinsically non-linear and often the event analysis, i.e., shaping, optimal filters etc., add additional non-linearity. Thus for large gain variations or when the best possible precision is required, a linear stretch correction is not sufficient. Here, we discuss a new correction technique based on non-linear interpolation of the energy-scale functions. Using Astro-HSXS calibration data, we demonstrate that the correction can recover the X-ray energy to better than 1 part in 104 over the entire spectral band to above 12 keV even for large-scale gain variations. This method will be used to correct any temporal drift of the on-orbit per-pixel gain using on-board calibration sources for the SXS instrument on the Astro-H observatory.

X-ray detector X-ray spectrometer X-ray calorimete↗

An Integrated Examination of AMSR2 Products over Ocean

Integrated processing and analysis are presented for oceanic retrievals of precipitation, sea ice concentration, columnar water vapor and cloud water, sea surface temperature, and near surface wind from the advanced microwave scanning radiometer-2 (AMSR2) sensor. By developing a common algorithmic framework and permitting iterative interaction between historically separate science algorithms, ambiguous and contradictory retrieval results are minimized. The integration also serves to improve each algorithm individually, decreasing reliance on ancillary datasets. Case studies are presented that exemplify both ongoing challenges for these retrievals and potential research uses of such an integrated satellite product. Two cases are presented each for precipitation at high latitudes and marginal sea ice detection, with analysis supplemented by satellite data from CloudSat profiles and Himawari-8 imagery. Detection and retrieval of snowfall remains a challenge, while light rainfall detection can be aided by a variational algorithm. Sea ice concentrations below 20% cause disagreements, but the algorithms otherwise agree well on the sea ice edge. Potential mitigation strategies for ambiguous areas of light rainfall and marginal sea ice are discussed. The analysis demonstrates potential avenues for future algorithm development, but also some physical limitations of remote sensing with the AMSR2 frequencies.

variational methods↗

UAS Integration in the NAS: Encounter Rate Simulation Study with UAS Missions

This study investigates the effects that different parameters of a given well clear metric has on the rate of well clear violations (WCV) in Class E airspace. The proposed metrics are taken from variations on methods used by TCAS IIand also considers different alerting schemes against the well clear metric. This presentation is divided into three analysis. The first analysis presented investigates the effect of well clear definitions on the rate of WCV occurrence. The second analysis characterises encounters at the WCV and the third analysis evaluates the alerting threshold definition and criteria. This study should inform refinement of an appropriate definition of well clear for UAS sense and avoid systems and motivate further research on exploring alerting definitions and criteria for a given well clear definition.

Johnson, Marcus↗

Comparison of Equilibrium Climate Sensitivity Estimates From Slab Ocean, 150‐Year, and Longer Simulations

We compare equilibrium climate sensitivity (ECS) estimates from pairs of long (≥ 800‐year) control and abruptly quadrupled CO2 simulations with shorter (150, 300 year) coupled atmosphere‐ocean simulations and Slab Ocean Models (SOM). Consistent with previous work, ECS estimates from shorter coupled simulations based on annual averages for years 1‐150 underestimate those from SOM (‐8% ± 13%) and long (‐14% ± 8%) simulations. Analysis of only years 21‐150 improved agreement with SOM (‐2% ± 14%) and long (‐8% ± 10%) estimates. Use of pentadal averages for years 51‐150 results in improved agreement with long simulations (‐4% ± 11%). While ECS estimates from current generation US models based on SOM and coupled annual averages of years 1‐150 range from 2.6°C to 5.3°C, estimates based longer simulations of the same models range from 3.2°C to 7.0°C. Such variations between methods argues for caution in comparison and interpretation of ECS estimates across models.

Climate sensitivity↗

Probability Bounds Analysis Applied to Multi-Purpose Crew Vehicle Nonlinearity

The Multi-Purpose Crew Vehicle (MPCV) Program Orion vehicle finite element model (FEM) was updated based on a modal test performed by Lockheed Martin. Due to nonlinearity observed in the test results, linear low force level (LL) and high force level (HL) FEMs were developed for use during various Space Launch System (SLS) flight regimes depending on expected forcing levels. Uncertainty models were derived for the combined MPCV and MPCV Stage Adaptor LL and HL Hurty/Craig-Bampton (HCB) components based on the MPCV structural test article Configuration 4 modal test-analysis correlation results. Subsequently, system-level uncertainty quantification analyses were performed using both models for various SLS flight configurations to determine the impact of the nonlinearity on important system metrics. The system metrics included both transfer functions associated with attitude control and dynamic loads associated with aerodynamic buffeting during ascent. In each case, an independent Monte Carlo (MC) analysis was performed, and no attempt was made to combine the results. The Hybrid Parametric Variation (HPV) method was used to develop the LL and HL MPCV HCB uncertainty models. The HPV method provides both parametric and non-parametric components of uncertainty. The non-parametric uncertainty accounts for the difference in model-form between the linearized analytical model and the corresponding linearized component test results in the form of mode shapes and frequencies at that force level. This linear model-form uncertainty is implemented in the HPV method using random matrix theory. However, the HPV uncertainty models developed for the linear LL and HL MPCV components do not account for the nonlinearity in the MPCV. With respect to the linearized models, this nonlinearity is also an uncertainty in model form, but in this case, it must be treated independently as an epistemic uncertainty. It represents a lack of knowledge, in contrast to an aleatory uncertainty due to the randomness of a variable. In the case of an epistemic variable, the true value is unknown, only the interval within which it lies is known. Epistemic uncertainty can be reduced with increased knowledge, while in general, aleatory uncertainty cannot. This work combines the epistemic uncertainty due to the MPCV nonlinearity with the parametric and non-parametric uncertainty within the HPV method using a second order propagation approach. The LL and HL test data is augmented with surrogate test data derived from a nonlinear MPCV representation. The impact of the MPCV nonlinearity on system response statistics is determined using a series of cumulative distribution functions in the form of a horsetail plot, or p-box. This results in an interval of probabilities for a specific response value, or an interval of response values at a specific probability.

Daniel C Kammer↗

Evaluation of a multivariate variational assimilation of conventional and satellite data for the diagnosis of cyclone systems

A variational data assimilation method for the study of cyclone-scale weather systems is described. The variational data assimilation method is to incorporate primitive equations for a moist, convectively unstable atmosphere and the radiative transfer equation. The variables to be adjusted include the three-dimensional vector wind, height, temperature, and moisture from rawinsonde data, and cloud-wind vectors, moisture, and radiance from satellite data. The development of variational model 1 which contains two nonlinear horizontal momentum equations, an integrated continuity equation, and a hydrostatic equation is examined. Examples applying the assimilation model to rawinsonde and satellite data are presented.

Achtemeier, Gary L.↗

Interplanetary Trajectories, Encke Method (ITEM)

Modified program has been developed using improved variation of Encke method which avoids accumulation of round-off errors and avoids numerical ambiguities arising from near-circular orbits of low inclination. Variety of interplanetary trajectory problems can be computed with maximum accuracy and efficiency.

Whitlock, F. H.↗

A unified approach for the application of general perturbation theories to the artificial satellite problem

An Encke-type method is developed as well as a variation of parameters method, both of which use a non-Keplerian reference orbit. The regularized time is used in the numerical integration and the optimum value of n is found for each type of orbit investigated. Accuracy and computation time comparisons are made with a classical Cowell method. It should be noted that the element formulations developed and tested were found to be exceptionally numerically stable in the sense that it was possible to achieve numerical consistency order of very long integration periods - a property not available with the Cowell formulation - and therefore these formulations may be helpful for high precision calculations.

Alfriend, K. T.↗

Grid spacing control with variation diminishing splines

Methods used to specify and control two and three dimensional grids on which numerical solutions of partial differential equations may be obtained were studied. While initially focusing on grid generation, the research will evolve into a consideration of the interaction of grid generation with the solution of a partial differential equation. The multisurface method of grid generation was used to continuously patch a grid onto an existing grid. In the resulting grid the elements of the Jacobian matrix must be continuous across the boundary between the original grid and the patched grid. Programs were written which accept as input the coordinates of the original grid and the desired new boundary, and then use the three-surface or four-surface version of the multisurface method to extend the original grid out to a new boundary.

Smith, P. W.↗

Additional Trials & Tribulations in Synthesizing a Sulfide Standard

Well-characterized sulfide reference materials that can serve as matrix-matched calibrants for in-situ trace element analyses via laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) remain elusive. Here, we describe the creation an in-house sulfide standard at NASA JSC that will be used to measure siderophile and chalcophile trace elements in low pressure experimental products, specifically pentlandite ([FeNi]9S8) and pyrrhotite ([FeNi]1-xS), as well as natural sulfides in terrestrial and meteoritic mineral assemblages. Recent methods in creating a homogenous sulfide standard include pressed pellets, synthesized chips, or fused glass [e.g. 1-3]. Other studies, however, have had success in forming homogenous FeS standards via synthesization methods [e.g. 4-6]. We present a variation of the method described in [4], in which we create a pyrrhotite standard doped with a variety of trace elements (Zn, As, Se, Mo, Ru, Rh, Pd, Sn, Sb, Te, W, Os, Ir, Pt, Au) at ~10-40 ppm and Cu at ~200 ppm. Sulfide compositions were constructed using Fe and Ni metal and elemental sulfur powders. Trace elements were added to elemental sulfur from Atomic Absorption Spectroscopy (AAS) elemental standards as either nitrate or chloride solutions to prevent oxidation of the metal powders. The dried sulfur mixture was mixed with Fe and Ni powders and mechanically mixed before two aliquots were placed in separate SiO2 tubes. Each tube was held under vacuum for ~30 minutes, sealed under vacuum, and then heated at 800C for 48 hours. Like [4], the synthesis products were composed of pourous sulfide crystals. Major element analyses of both experimental aliquots, collected using an Electron Probe Microanalyser (EPMA), yield indistinguishable major element compositions (uncertainties in 2SE & 2RSE[%]), with an average of 57.90  0.09 (0.16 %), 4.95  0.03 (0.55 %), and 39.70  0.14 (0.34 %), for Fe, Ni, and S, respectively. Trace element data were measured using a Photon Machines 193nm laser ablation system coupled to a Thermo-Scientific Element-XR ICP-MS. Spot sizes were limited to 50 μm due to the porous nature of the sulfide target material. Trace element abundances, normalized to Fe as an internal standard, were also found to be homogenous between the two aliquots, with weighted mean 2RSE (%) values of <3.0 for all trace elements. Synthesized products were re-powdered and absolute concentrations measured via solution ICP-MS. Although the sulfide appears to be homogenous, sintering experiments will be performed to more closely match the standard density to natural sulfides and minimize differences in ablation behavior. Improved density also allows for higher sensitivity (i.e. more compact target material) and larger spot sizes or traverses, as void space is eliminated. Additional major and trace element analyses on the products of the sintering experiments will be undertaken.

Jacob B Setera↗

D-Wave Electron-H, -He+, and -Li2+ Elastic Scattering and Photoabsorption in P States of Two-Electron Systems

In previous papers [A. K. Bhatia, Phys. Rev. A 85, 052708 (2012); 86, 032709 (2012); 87, 042705 (2013)] electron-H, -He+, and -Li2+ P-wave scattering phase shifts were calculated using the variational polarized orbital theory. This method is now extended to the singlet and triplet D-wave scattering in the elastic region. The long-range correlations are included in the Schrodinger equation by using the method of polarized orbitals variationally. Phase shifts are compared to those obtained by other methods. The present calculation provides results which are rigorous lower bonds to the exact phase shifts. Using the presently calculated D-wave and previously calculated S-wave continuum functions, photoionization of singlet and triplet P states of He and Li+ are also calculated, along with the radiative recombination rate coefficients at various electron temperatures.

D-wave↗