Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “factorization 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 199 records · Page 11

Quantifying VOC-Reaction Tracers, Ozone Production, and Continuing Aerosol Production Rates in Urban and Far-Downwind Atmospheres

We have found a surprisingly informative decomposition of the complex question of smoggy ozone production (basically, [HO2] in a more locally determined field of [NO]) in the process of linked investigations of modestly smoggy Eastern North America (by NASA aircraft, July 2004) and rather polluted Flushing, NYC (Queens College, July, 2001). In both rural and very polluted situations, we find that a simple contour graph parameterization of the local principal ozone production rate can be estimated using only the variables [NO] and j(sub rads) [HCHO]: Po(O3) = c (j(sub rads) [HCHO])(sup a) [HCHO](sup b). Here j(sub rads) is the photolysis of HCHO to radicals, presumably capturing many harder-UV photolytic processes and the principle ozone production is that due to HO2; mechanisms suggest that ozone production due to RO2 is closely correlated, often suggesting a limited range of different proportionality factors. The method immediately suggests a local interpretation for concepts of VOC limitation and NOx limitation. We believe that the product j(sub rads) [HCHO] guages the oxidation rate of observed VOC mixtures in a way that also provides [HO2] useful for the principle ozone production rate k [HO2] [NO], and indeed, all ozone chemical production. The success of the method suggests that dominant urban primary-HCHO sources may transition to secondary plume-HCHO sources in a convenient way. Are there other, simple, near-terminal oxidized VOC's which help guage ozone production and aerosol particle formation? Regarding particles, we report on, to the extent NASA Research resources allow, on appealing relationships between far-downwind (Atlantic PBL) HCHO and very fine aerosol (including sulfate. Since j(sub rads) [HCHO] provides a time-scale, we may understand distant-plume particle production in a more quantitative manner. Additionally we report on a statistical search in the nearer field for relationships between glyoxals (important near-terminal aromatic and isoprene reaction products) and aerosol production, looking for VOC's that might be most implicated. All three variables j(sub rads), [HCHO], and [NO] are relatively easily measured in widespread air pollution monitoring networks, and all are deducible form space-borne observations, though estimation of [NO] from [NO2] (the species observable from space) may require care. We report also on airborne and surface observations of HCHO, suggesting that concentrated (urban) and more diffuse (forest) sources may be distinguishable from space. The use of the 3.58 micron microwindow for HCHO remote sensing should allow much sharper resolution of HCHO than the UV. UV sensing requires large and expensive instruments, but even these seem justified since formaldehyde is so informative.

Chatfield, Robert↗

Microphysics and Radiation Effect of Dust on Saharan Air Layer: An HS3 Case Study

A Saharan air layer (SAL) event associated with a nondeveloping African easterly wave (AEW) over the main development region of the eastern Atlantic was sampled by the NASA Global Hawk aircraft on 24-25 August 2013 during the NASA Hurricane and Severe Storm Sentinel (HS3) campaign and was simulated with the NASA Unified Weather Research and Forecasting (NU-WRF) Model. Airborne, ground-based, and spaceborne measurements were used to evaluate the model performance. The microphysical and radiative effects of dust and other aerosols on the SAL structure and environment were investigated with the factor-separation method. The results indicate that relative to a simulation without dust-radiative and microphysical impacts, Saharan dust and other aerosols heated the SAL air mainly through shortwave heating by the direct aerosol-radiation (AR) effect, resulting in a warmer (up to 0.6 K) and drier (up to 5% RH reduction) SAL and maintaining the strong temperature inversion at the base of the SAL in the presence of predominant longwave cooling. Radiative heating of the dust accentuated a vertical circulation within the dust layer, in which air rose (sank) in the northern (southern) portions of the dust layer. Furthermore, above and to the south of the dust layer, both the microphysical and radiative impacts of dust tended to counter the vertical motions associated with the Hadley circulation, causing a small weakening and southward shift of convection in the intertropical convergence zone (ITCZ) and reduced anvil cloud to the north. Changes in moisture and cloud/precipitation hydrometeors were largely driven by the dust induced changes in vertical motion. Dust strengthened the African easterly jet by up to ~1ms(exp -1) at the southern edge of the jet, primarily through the AR effect, and produced modest increases in vertical wind shear within and in the vicinity of the dust layer. These modulations of the SAL and AEW environment clearly contributed to the nondevelopment of this AEW.

Tao, Zhining↗

Implicit calculations of transonic flows using monotone methods

Implicit approximate-factorization algorithms have been developed that use monotone methods for the calculation of steady and unsteady transonic flows governed by the small-disturbance-potential equation. These algorithms use the new Engquist-Osher switch in the type-dependent differencing in place of the standard Murman-Cole switch. The resulting algorithms are more stable; hence, calculations can be done more efficiently. For steady flows, the convergence rate is about 35% faster, and for unsteady flows the allowable time step is about 10 times larger. These improvements are achieved with no increase in computer storage and with only minor modifications in codes that use the Murman-Cole switch. Also an implicit algorithm has been developed for the steady full-potential equation in one-dimension, which uses monotone methods.

Goorjian, P. M.↗

A History of Collapse Factor Modeling and Empirical Data for Cryogenic Propellant Tanks

One of the major technical problems associated with cryogenic liquid propellant systems used to supply rocket engines and their subassemblies and components is the phenomenon of propellant tank pressurant and ullage gas collapse. This collapse is mainly caused by heat transfer from ullage gas to tank walls and interfacing propellant, which are both at temperatures well below those of this gas. Mass transfer between ullage gas and cryogenic propellant can also occur and have minor to significant secondary effects that can increase or decrease ullage gas collapse. Pressurant gas is supplied into cryogenic propellant tanks in order to initially pressurize these tanks and then maintain required pressures as propellant is expelled from these tanks. The net effect of pressurant and ullage gas collapse is increased total mass and mass flow rate requirements of pressurant gases. For flight vehicles this leads to significant and undesirable weight penalties. For rocket engine component and subassembly ground test facilities this results in significantly increased facility hardware, construction, and operational costs. "Collapse Factor" is a parameter used to quantify the pressurant and ullage gas collapse. Accurate prediction of collapse factors, through analytical methods and modeling tools, and collection and evaluation of collapse factor data has evolved over the years since the start of space exploration programs in the 1950 s. Through the years, numerous documents have been published to preserve results of studies associated with the collapse factor phenomenon. This paper presents a summary and selected details of prior literature that document the aforementioned studies. Additionally other literature that present studies and results of heat and mass transfer processes, related to or providing important insights or analytical methods for the studies of collapse factor, are presented.

deQuay, Laurence↗

Loss Factor Estimation Using the Impulse Response Decay Method on a Stiffened Structure

High-frequency vibroacoustic modeling is typically performed using energy-based techniques such as Statistical Energy Analysis (SEA). Energy models require an estimate of the internal damping loss factor. Unfortunately, the loss factor is difficult to estimate analytically, and experimental methods such as the power injection method can require extensive measurements over the structure of interest. This paper discusses the implications of estimating damping loss factors using the impulse response decay method (IRDM) from a limited set of response measurements. An automated procedure for implementing IRDM is described and then evaluated using data from a finite element model of a stiffened, curved panel. Estimated loss factors are compared with loss factors computed using a power injection method and a manual curve fit. The paper discusses the sensitivity of the IRDM loss factor estimates to damping of connected subsystems and the number and location of points in the measurement ensemble.

Cabell, Randolph↗

Method and apparatus for lead-unity-lag electric power generation system

A method employing a lead-unity-lag adjustment on a power generation system is disclosed. The method may include calculating a unity power factor point and adjusting system parameters to shift a power factor angle to substantially match an operating power angle creating a new unity power factor point. The method may then define operation parameters for a high reactance permanent magnet machine based on the adjusted power level.

Ganev, Evgeni↗

Estimating electric field enhancement factors on an aircraft utilizing a small scale model: A method evaluation

A method for obtaining field enhancement factors at specific points on an aircraft utilizing a small scale model was evaluated by measuring several canonical shapes. Comparison of the form factors obtained by analytical means with measurements indicate that the experimental method has serious flaws. Errors of 200 to 300 percent were found between analytical values and measured values. As a result of the study, the analytical method is not recommended for calibration of field meters located on aircraft, and should not be relied upon in any application where the local spatial derivatives of the electric field on the model are large over the dimensions of the sensing probe.

Easterbrook, Calvin C.↗

Calculating parts factors for redundant systems

Method that is easily programmed simplifies calculation of parts factor. Individual module unreliabilities are computed as function of number of service intervals and service interval length. At each service interval, unreliability is sum of unreliabilities of replaced and original modules. It must be calculated for each module to obtain parts factor.

Derocher, W. L., Jr.↗

Stress-intensity factor calculations using the boundary force method

The Boundary Force Method (BFM) was formulated for the three fundamental problems of elasticity: the stress boundary value problem, the displacement boundary value problem, and the mixed boundary value problem. Because the BFM is a form of an indirect boundary element method, only the boundaries of the region of interest are modeled. The elasticity solution for the stress distribution due to concentrated forces and a moment applied at an arbitrary point in a cracked infinite plate is used as the fundamental solution. Thus, unlike other boundary element methods, here the crack face need not be modeled as part of the boundary. The formulation of the BFM is described and the accuracy of the method is established by analyzing a center-cracked specimen subjected to mixed boundary conditions and a three-hole cracked configuration subjected to traction boundary conditions. The results obtained are in good agreement with accepted numerical solutions. The method is then used to generate stress-intensity solutions for two common cracked configurations: an edge crack emanating from a semi-elliptical notch, and an edge crack emanating from a V-notch. The BFM is a versatile technique that can be used to obtain very accurate stress intensity factors for complex crack configurations subjected to stress, displacement, or mixed boundary conditions. The method requires a minimal amount of modeling effort.

Tan, P. W.↗

Time-fixed rendezvous by impulse factoring with an intermediate timing constraint

A method is presented for factoring a two-impulse orbital transfer into a three- or four-impulse transfer which solves the rendezvous problem and satisfies an intermediate timing constraint. Both the time of rendezvous and the intermediate time of a alinement are formulated as any element of a finite sequence of times. These times are integer multiples of a constant plus an additive constant. The rendezvous condition is an equality constraint, whereas the intermediate alinement is an inequality constraint. The two timing constraints are satisfied by factoring the impulses into collinear parts that vectorially sum to the original impulse and by varying the resultant period differences and the number of revolutions in each orbit. Five different types of solutions arise by considering factoring either or both of the two impulses into two or three parts with a limit for four total impulses. The impulse-factoring technique may be applied to any two-impulse transfer which has distinct orbital periods.

Green, R. N.↗

Modification of an impulse-factoring orbital transfer technique to account for orbit determination and maneuver execution errors

A method has previously been developed to satisfy terminal rendezvous and intermediate timing constraints for planetary missions involving orbital operations. The method uses impulse factoring in which a two-impulse transfer is divided into three or four impulses which add one or two intermediate orbits. The periods of the intermediate orbits and the number of revolutions in each orbit are varied to satisfy timing constraints. Techniques are developed to retarget the orbital transfer in the presence of orbit-determination and maneuver-execution errors. Sample results indicate that the nominal transfer can be retargeted with little change in either the magnitude (Delta V) or location of the individual impulses. Additonally, the total Delta V required for the retargeted transfer is little different from that required for the nominal transfer. A digital computer program developed to implement the techniques is described.

Kibler, J. F.↗

Comparison of Implicit Schemes for the Incompressible Navier-Stokes Equations

For a computational flow simulation tool to be useful in a design environment, it must be very robust and efficient. To develop such a tool for incompressible flow applications, a number of different implicit schemes are compared for several two-dimensional flow problems in the current study. The schemes include Point-Jacobi relaxation, Gauss-Seidel line relaxation, incomplete lower-upper decomposition, and the generalized minimum residual method preconditioned with each of the three other schemes. The efficiency of the schemes is measured in terms of the computing time required to obtain a steady-state solution for the laminar flow over a backward-facing step, the flow over a NACA 4412 airfoil, and the flow over a three-element airfoil using overset grids. The flow solver used in the study is the INS2D code that solves the incompressible Navier-Stokes equations using the method of artificial compressibility and upwind differencing of the convective terms. The results show that the generalized minimum residual method preconditioned with the incomplete lower-upper factorization outperforms all other methods by at least a factor of 2.

Rogers, Stuart E.↗

Estimation of peak winds from hourly observations

Two closely related methods to obtain estimates of the hourly peak wind at Cape Kennedy were compared by statistical tests. The methods evaluated the Monin-Obukhov stability length and the standard deviation of the hourly observed wind speed, so as to augment the latter quantity by F standard deviations. F is an optimized factor. A third method utilizing an optimized gust factor was also applied to the hourly wind. The latter procedure estimated 2952 peak winds with an rms error of 2.81 knots, an accuracy which was not surpassed by the other methods. Peak ground wind speed data were developed for use in space shuttle design operation analyses.

Graves, M. E.↗

Alternating direction implicit methods for parabolic equations with a mixed derivative

Alternating direction implicit (ADI) schemes for two-dimensional parabolic equations with a mixed derivative are constructed by using the class of all A sub 0-stable linear two-step methods in conjunction with the method of approximation factorization. The mixed derivative is treated with an explicit two-step method which is compatible with an implicit A sub 0-stable method. The parameter space for which the resulting ADI schemes are second order accurate and unconditionally stable is determined. Some numerical examples are given.

Beam, R. M.↗

Comparison of methods for prediction of transition by stability analysis

The paper compares several methods of transition prediction of linear stability analysis. The spectral stability analysis code SALLY is used to analyze flows over laminar flow control wings. It is shown that transition prediction by the envelope method and a new modified wave packet method are comparable in reliability, but that the envelope method is more efficient computationally. This is based on the results which show that the wave packet method provides N factors which are at best as consistent as those of the envelope method; since the wave packet method is at least 3 times as expensive as the envelope method, the latter is recommended for engineering design calculations.

Malik, M. R.↗

Alternating direction implicit methods for parabolic equations with a mixed derivative

Alternating direction implicit (ADI) schemes for two-dimensional parabolic equations with a mixed derivative are constructed by using the class of all A(0)-stable linear two-step methods in conjunction with the method of approximate factorization. The mixed derivative is treated with an explicit two-step method which is compatible with an implicit A(0)-stable method. The parameter space for which the resulting ADI schemes are second-order accurate and unconditionally stable is determined. Some numerical examples are given.

Beam, R. M.↗