Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “delta method”

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 649 records · Page 36

Estimating Mangrove Aboveground Biomass from Airborne LiDAR Data: a Case Study from the Zambezi River Delta

Mangroves are ecologically and economically important forested wetlands with the highest carbon (C) density of all terrestrial ecosystems. Because of their exceptionally large C stocks and importance as a coastal buffer, their protection and restoration has been proposed as an effective mitigation strategy for climate change. The inclusion of mangroves in mitigation strategies requires the quantification of C stocks (both above and belowground) and changes to accurately calculate emissions and sequestration. A growing number of countries are becoming interested in using mitigation initiatives, such as REDD+ (reducing emissions from deforestation and forest degradation), in these unique coastal forests. However, it is not yet clear how methods to measure C traditionally used for other ecosystems can be modified to estimate biomass in mangroves with the precision and accuracy needed for these initiatives. Airborne Lidar (ALS) data has often been proposed as the most accurate way for larger scale assessments but the application of ALS for coastal wetlands is scarce, primarily due to a lack of contemporaneous ALS and field measurements. Here, we evaluated the variability in field and Lidar-based estimates of aboveground biomass (AGB) through the combination of different local and regional allometric models and standardized height metrics that are comparable across spatial resolutions and sensor types, the end result being a simplified approach for accurately estimating mangrove AGB at large scales and determining the uncertainty by combining multiple allometric models.We then quantified wall-to-wall AGB stocks of a tall mangrove forest in the Zambezi Delta, Mozambique. Our results indicate that the Lidar H100 height metric correlates well with AGB estimates, with R(exp 2) between 0.80 and 0.88 and RMSE of 33% or less. When comparing Lidar H100 AGB derived from three allometric models, mean AGB values range from 192 Mg ha(exp −1) up to 252 Mg ha(exp −1).We suggest the best model to predict AGB was based on the East Africa specific allometry and a power-based regression that used Lidar H100 as the height input with an R(exp 2) of 0.85 and an RMSE of 122 Mg ha(exp −1) or 33%. The total AGB of the Lidar inventoried mangrove area (6654 ha) was 1 350 902 Mg with a mean AGB of 203 Mg ha(exp −1) +/-166 Mg ha(exp −1). Because the allometry suggested here was developed using standardized height metrics, it is recommended that the models can generate AGB estimates using other remote sensing instruments that are more readily accessible over other mangrove ecosystems on a large scale, and as part of future carbon monitoring efforts in mangroves.

Fatoyinbo, Temilola↗

The CFL condition for spectral approximations to hyperbolic initial-boundary value problems

The stability of spectral approximations to scalar hyperbolic initial-boundary value problems with variable coefficients are studied. Time is discretized by explicit multi-level or Runge-Kutta methods of order less than or equal to 3 (forward Euler time differencing is included), and spatial discretizations are studied by spectral and pseudospectral approximations associated with the general family of Jacobi polynomials. It is proved that these fully explicit spectral approximations are stable provided their time-step, delta t, is restricted by the CFL-like condition, delta t less than Const. N(exp-2), where N equals the spatial number of degrees of freedom. We give two independent proofs of this result, depending on two different choices of approximate L(exp 2)-weighted norms. In both approaches, the proofs hinge on a certain inverse inequality interesting for its own sake. The result confirms the commonly held belief that the above CFL stability restriction, which is extensively used in practical implementations, guarantees the stability (and hence the convergence) of fully-explicit spectral approximations in the nonperiodic case.

Gottlieb, David↗

The CFL condition for spectral approximations to hyperbolic initial-boundary value problems

The stability of spectral approximations to scalar hyperbolic initial-boundary value problems with variable coefficients are studied. Time is discretized by explicit multi-level or Runge-Kutta methods of order less than or equal to 3 (forward Euler time differencing is included), and spatial discretizations are studied by spectral and pseudospectral approximations associated with the general family of Jacobi polynomials. It is proved that these fully explicit spectral approximations are stable provided their time-step, delta t, is restricted by the CFL-like condition, delta t less than Const. N(exp-2), where N equals the spatial number of degrees of freedom. We give two independent proofs of this result, depending on two different choices of approximate L(exp 2)-weighted norms. In both approaches, the proofs hinge on a certain inverse inequality interesting for its own sake. The result confirms the commonly held belief that the above CFL stability restriction, which is extensively used in practical implementations, guarantees the stability (and hence the convergence) of fully-explicit spectral approximations in the nonperiodic case.

Gottlieb, David↗

Path planning for assembly of strut-based structures

A path planning method with collision avoidance for a general single chain nonredundant or redundant robot is proposed. Joint range boundary overruns are also avoided. The result is a sequence of joint vectors which are passed to a trajectory planner. A potential field algorithm in joint space computes incremental joint vectors delta-q = delta-q(sub a) + delta-q(sub c) + delta-q(sub r). Adding delta-q to the robot's current joint vector leads to the next step in the path. Delta-q(sub a) is obtained by computing the minimum norm solution of the underdetermined linear system J delta-q(sub a) = x(sub a) where x(sub a) is a translational and rotational force vector that attracts the robot to its goal position and orientation. J is the manipulator Jacobian. Delta-q(sub c) is a collision avoidance term encompassing collisions between the robot (links and payload) and obstacles in the environment as well as collisions among links and payload of the robot themselves. It is obtained in joint space directly. Delta-q(sub r) is a function of the current joint vector and avoids joint range overruns. A higher level discrete search over candidate safe positions is used to provide alternatives in case the potential field algorithm encounters a local minimum and thus fails to reach the goal. The best first search algorithm A* is used for graph search. Symmetry properties of the payload and equivalent rotations are exploited to further enlarge the number of alternatives passed to the potential field algorithm.

Muenger, Rolf↗

The Effect of Infrastructure Sharing in Estimating Operations Cost of Future Space Transportation Systems

NASA and the aerospace industry are extremely serious about reducing the cost and improving the performance of launch vehicles both manned or unmanned. In the aerospace industry, sharing infrastructure for manufacturing more than one type spacecraft is becoming a trend to achieve economy of scale. An example is the Boeing Decatur facility where both Delta II and Delta IV launch vehicles are made. The author is not sure how Boeing estimates the costs of each spacecraft made in the same facility. Regardless of how a contractor estimates the cost, NASA in its popular cost estimating tool, NASA Air force Cost Modeling (NAFCOM) has to have a method built in to account for the effect of infrastructure sharing. Since there is no provision in the most recent version of NAFCOM2002 to take care of this, it has been found by the Engineering Cost Community at MSFC that the tool overestimates the manufacturing cost by as much as 30%. Therefore, the objective of this study is to develop a methodology to assess the impact of infrastructure sharing so that better operations cost estimates may be made.

Sundaram, Meenakshi↗

Small Launch Vehicle Trade Space Definition: Development of a Zero Level Mass Estimation Tool with Trajectory Validation

Recent high level interest in the capability of small launch vehicles has placed significant demand on determining the trade space these vehicles occupy. This has led to the development of a zero level analysis tool that can quickly determine the minimum expected vehicle gross liftoff weight (GLOW) in terms of vehicle stage specific impulse (Isp) and propellant mass fraction (pmf) for any given payload value. Utilizing an extensive background in Earth to orbit trajectory experience a total necessary delta v the vehicle must achieve can be estimated including relevant loss terms. This foresight into expected losses allows for more specific assumptions relating to the initial estimates of thrust to weight values for each stage. This tool was further validated against a trajectory model, in this case the Program to Optimize Simulated Trajectories (POST), to determine if the initial sizing delta v was adequate to meet payload expectations. Presented here is a description of how the tool is setup and the approach the analyst must take when using the tool. Also, expected outputs which are dependent on the type of small launch vehicle being sized will be displayed. The method of validation will be discussed as well as where the sizing tool fits into the vehicle design process.

Waters, Eric D.↗

Four-body trajectory optimization

A collection of typical three-body trajectories from the L1 libration point on the sun-earth line to the earth is presented. These trajectories in the sun-earth system are grouped into four distinct families which differ in transfer time and delta V requirements. Curves showing the variations of delta V with respect to transfer time, and typical two and three-impulse primer vector histories, are included. The development of a four-body trajectory optimization program to compute fuel optimal trajectories between the earth and a point in the sun-earth-moon system are also discussed. Methods for generating fuel optimal two-impulse trajectories which originate at the earth or a point in space, and fuel optimal three-impulse trajectories between two points in space, are presented. A brief qualitative comparison of these methods is given. An example of a four-body two-impulse transfer from the Li libration point to the earth is included.

Pu, C. L.↗

Utility of correlation techniques in gravity and magnetic interpretation

Internal correspondence uses Poisson's Theorem in a moving-window linear regression analysis between the anomalous first vertical derivative of gravity and total magnetic field reduced to the pole. The regression parameters provide critical information on source characteristics. The correlation coefficient indicates the strength of the relation between magnetics and gravity. Slope value gives delta j/delta sigma estimates of the anomalous source. The intercept furnishes information on anomaly interference. Cluster analysis consists of the classification of subsets of data into groups of similarity based on correlation of selected characteristics of the anomalies. Model studies are used to illustrate implementation and interpretation procedures of these methods, particularly internal correspondence. Analysis of the results of applying these methods to data from the midcontinent and a transcontinental profile shows they can be useful in identifying crustal provinces, providing information on horizontal and vertical variations of physical properties over province size zones, validating long wavelength anomalies, and isolating geomagnetic field removal problems.

Chandler, V. W.↗

Optimum two-impulse transfers for preliminary interplanetary trajectory design

A method for numerically determining the optimum two-impulse transfer between two positions in two different heliocentric orbits is described. The formulas are based on earlier work by McCue and Bender (1965) for optimum two-impulse transfers between two orbits about the earth. The contours of minimum total Delta-V(Delta-V to launch plus Delta-V to achieve a given final orbit at the target body) are plotted on axes of true anomaly of the launch body at launch and true anomaly of the target body at arrival. The optimum transfers obtained are 'time-open' and the Delta-V from the initial orbit to the final orbit is found using patched conics. These 'prime rib' plots are useful for selecting initial conditions in the early phases of mission design.

Hulkower, N. D.↗

An integrated CFD/experimental analysis of aerodynamic forces and moments

Aerodynamic analysis using computational fluid dynamics (CFD) is most fruitful when it is combined with a thorough program of wind tunnel testing. The understanding of aerodynamic phenomena is enhanced by the synergistic use of both analysis methods. A technique is described for an integrated approach to determining the forces and moments acting on a wind tunnel model by using a combination of experimentally measured pressures and CFD predictions. The CFD code used was FLO57 (an Euler solver) and the wind tunnel model was a heavily instrumented delta wing with 62.5 deg of leading-edge sweep. A thorough comparison of the CFD results and the experimental data is presented for surface pressure distributions and longitudinal forces and moments. The experimental pressures were also integrated over the surface of the model and the resulting forces and moments are compared to the CFD and wind tunnel results. The accurate determination of various drag increments via the combined use of the CFD and experimental pressures is presented in detail.

Melton, John E.↗

Constraints on the formation history of the elliptical galaxy NGC 3923 from the colors of its globular clusters

We present a study of the colors of globular clusters associated with the elliptical galaxy NGC 3923. Our final sample consists of Wasington system C and T(sub 1) photometry for 143 globular cluster candidates with an expected contamination of no more than 10%. We find that the color distribution of the NGC 3923 globular cluster system (GCS) is broad and appears to have at least two peaks. A mixture modeling analysis of the color distribution indicates that a two-component model is favored over a single-component one at a high level of confidence (greater than 99%). This evidence for more than one population in the GCS of NGC 3923 is similar to that previously noted for the four other elliptical galaxies for which similar data have been published. Furthermore, we find that the NGC 3923 GCS is redder than the GCSs of previously studed elliptical galaxies of similar luminosity. The median metallicity inferred from our (C-(T(sub 1)))(sub 0) colors is (Fe/H)(sub med) = -0.56, with an uncertainty of 0.14 dex arising from all sources of uncertainty in the mean color. This is more metal rich than the median metallicity found for the GCS of M87 using the same method, (Fe/H)(sub med) = -0.94. Since M87 is more luminous than NGC 3923, this result points to significant scatter about any trend of higher GCS metallicity with increasing galaxy luminosity. We also show that there is a color gradient in the NGC 3923 GCS corresponding to about -0.5 dex in Delta(Fe/H)/Delta(log r). We conclude that the shape of the color distribution of individual GCSs and the variation in mean color among the GCSs of ellipticals are difficult to understand if elliptical galaxies are formed in a single protogalactic collapse. Models in which ellipticals and their globular clusters are formed in more than one event, such as a merger scenario, are more successful in accounting for these observations.

Zepf, Stephen E.↗

The Practical Calculation of the Aerodynamic Characteristics of Slender Finned Vehicles

The basic objective of this thesis is to provide a practical method of computing the aerodynamic characteristics of slender finned vehicles such as sounding rockets, high speed bombs, and guided missiles. The aerodynamic characteristics considered are the normal force coefficient derivative, c(sub N(sub alpha)); center of pressure, bar-X; roll forcing moment coefficient derivative, c(sub l(sub delta)); roll damping moment coefficient derivative, c(sub l(sub p)); pitch damping moment coefficient derivative, c(sub mq); and the drag coefficient, c (sub D). Equations are determined for both subsonic and supersonic flow. No attempts is made to analyze the transonic region. The general configuration to which the relations are applicable is a slender axisymmetric body having three or four fins.

Barrowman, James S.↗

A comparison of high resolution upwind solvers on 3-D inviscid hypersonic flows

A detailed comparison of numerical results obtained by solving the Euler equations for the inviscid flow over delta wings in reentry configurations is presented. The investigation involves a side-by-side comparison between independently developed upwind Euler solvers at VKI (M3D) and NASA Johnson Space Center (E3D) using identical grids. In both solvers the governing equations are integrated by means of time marching finite volume shock capturing methods, based on a cell centered upwind evaluation of the cell face fluxes and nonlinear limiters. High resolution schemes are obtained via MUSCL characteristic variable extrapolation to ensure total variation diminishing (TVD) properties and therefore monotonic discontinuity capturing. Further comparisons are made with several published results by other authors based on both upwind and central discretizations. The present results contribute in asserting the high resolution upwind TVD schemes as the most reliable numerical technique to handle the strong discontinuities typical of high speed flows. The performances of the two upwind solvers are satisfactory and the numerical results in good agreement. However, the important issue of reaching a grid converged solution for the present complex three-dimensional inviscid flow problems was not achieved.

Manna, M.↗

The Turbulent Boundary Layer on a Rough Curvilinear Surface

A number of semiempirical approximate methods exist for determining the characteristics of the turbulent boundary layer on a curvilinear surface. At present, among these methods, the one proposed by L. G. Loitsianskii is given frequent practical application. This method is sufficiently effective and permits, in the case of wing profiles with technically smooth surfaces, calculating the basic characteristics of the boundary layer and the values of the overall drag with an accuracy which suffices for practical purposes. The idea of making use of the basic integral momentum equation ((d delta(sup xx))/dx) + ((V' delta(sup xx))/V) (2 + H) = (tau(sub 0))/(rho V(exp 2)) proves to be fruitful also for the solution of the problems in the determination of the characteristics of the turbulent boundary layer on a rough surface.

Droblenkov, V. F.↗

Pressure-Sensitive Paint Investigation of Double-Delta Wing Vortex Flow Manipulation

A pressure-sensitive paint (PSP) technique was applied in a wind tunnel experiment in the NASA Langley Research Center 8-Foot Transonic Pressure Tunnel to quantify the effect of wing fillets on the global vortex-induced surface static pressure field about a sharp leading-edge 76 deg/40 deg double delta wing, or strake-wing, model at subsonic and transonic speeds. Global calibrations of the PSP were obtained at M = 0.50, 0.70, 0.85, 0.95, and 1.20, a Reynolds number per unit length of 2.0 million, and angles of attack from 10 degrees to 30 degrees using an in-situ method featuring the simultaneous acquisition of electronically-scanned pressures (ESP) at discrete locations on the model. The mean error in the PSP measurements relative to the ESP data was approximately 2 percent or less at M = 0.50 to 0.85 but increased to several percent at M = 0.95 and 1.20. The PSP pressure distributions and pseudo-colored planform view pressure maps clearly revealed the vortex-induced pressure signatures at all Mach numbers and angles of attack. Small fillets having a parabolic or diamond planform situated at the strake-wing intersection were designed to manipulate the vortical flows by, respectively, removing the leading-edge discontinuity or introducing additional discontinuities. The fillets caused global changes in the vortex-dominated surface pressure field that were effectively captured in the PSP measurements. The vortex surface pressure signatures were compared to available off-surface vortex cross-flow structures obtained using a laser vapor screen (LVS) flow visualization technique. The fillet effects on the PSP pressure distributions and the observed leading-edge vortex flow characteristics were consistent with the trends in the measured lift, drag, and pitching moment coefficients.

Erickson, Gary E.↗

Pressure-Sensitive Paint Investigation of Double-Delta Wing Vortex Flow Manipulation

A pressure-sensitive paint (PSP) technique was applied in a wind tunnel experiment in the NASA Langley Research Center 8-Foot Transonic Pressure Tunnel to quantify the effect of wing fillets on the global vortex-induced surface static pressure field about a sharp leading-edge 76o/40o double delta wing, or strake-wing, model at subsonic and transonic speeds. Global calibrations of the PSP were obtained at M = 0.50, 0.70, 0.85, 0.95, and 1.20, a Reynolds number per unit length of 2.0 million, and angles of attack from 10 degrees to 20 degrees using an in-situ method featuring the simultaneous acquisition of electronically-scanned pressures (ESP) at discrete locations on the model. The mean error in the PSP measurements relative to the ESP data was approximately 2 percent or less at M = 0.50 to 0.85 but increased to several percent at M =0.95 and 1.20. The PSP pressure distributions and pseudo-colored planform view pressure maps clearly revealed the vortex-induced pressure signatures at all Mach numbers and angles of attack. Small fillets having a parabolic or diamond planform situated at the strake-wing intersection were designed to manipulate the vortical flows by, respectively, removing the leading-edge discontinuity or introducing additional discontinuities. The fillets caused global changes in the vortex-dominated surface pressure field that were effectively captured in the PSP measurements. The vortex surface pressure signatures were compared to available off-surface vortex cross-flow structures obtained using a laser vapor screen (LVS) flow visualization technique. The fillet effects on the PSP pressure distributions and the observed leading-edge vortex flow characteristics were consistent with the trends in the measured lift, drag, and pitching moment coefficients.

Erickson, Gary E.↗

Wind Tunnel Application of a Pressure-Sensitive Paint Technique to a Double Delta Wing Model at Subsonic and Transonic Speeds

A pressure-sensitive paint (PSP) technique was applied in a wind tunnel experiment in the NASA Langley Research Center 8-Foot Transonic Pressure Tunnel to study the effect of wing fillets on the global vortex induced surface static pressure field about a sharp leading-edge 76 deg./40 deg. double delta wing, or strake-wing, model at subsonic and transonic speeds. Global calibrations of the PSP were obtained at M(sub infinity) = 0.50, 0.70, 0.85, 0.95, and 1.20, a Reynolds number per unit length of 2.0 million, and angles of attack from 10 degrees to 20 degrees using an insitu method featuring the simultaneous acquisition of electronically scanned pressures (ESP) at discrete locations on the model. The mean error in the PSP measurements relative to the ESP data was approximately 2 percent or less at M(sub infinity) = 0.50 to 0.85 but increased to several percent at M(sub infinity) =0.95 and 1.20. The PSP pressure distributions and pseudo-colored, planform-view pressure maps clearly revealed the vortex-induced pressure signatures at all Mach numbers and angles of attack. Small fillets having parabolic or diamond planforms situated at the strake-wing intersection were respectively designed to manipulate the vortical flows by removing the leading-edge discontinuity or introducing additional discontinuities. The fillets caused global changes in the vortex-dominated surface pressure field that were effectively captured in the PSP measurements. The vortex surface pressure signatures were compared to available off-surface vortex cross-flow structures obtained using a laser vapor screen (LVS) flow visualization technique. The fillet effects on the PSP pressure distributions and the observed leading-edge vortex flow characteristics were consistent with the trends in the measured lift, drag, and pitching moment coefficients.

Erickson, Gary E.↗

The Local Topological Free Energy of the SARS-CoV-2 Spike Protein

The novel coronavirus SARS-CoV-2 infects human cells using a mechanism that involves binding and structural rearrangement of its Spike protein. Understanding protein rearrangement and identifying specific amino acids where mutations affect protein rearrangement has attracted much attention for drug development. In this manuscript, we use a mathematical method to characterize the local topology/geometry of the SARS-CoV-2 Spike protein backbone. Our results show that local conformational changes in the FP, HR1, and CH domains are associated with global conformational changes in the RBD domain. The SARS-CoV-2 variants analyzed in this manuscript (alpha, beta, gamma, delta Mink, G614, N501) show differences in the local conformations of the FP, HR1, and CH domains as well. Finally, most mutations of concern are either in or in the vicinity of high local topological free energy conformations, suggesting that high local topological free energy conformations could be targets for mutations with significant impact of protein function. Namely, the residues 484, 570, 614, 796, and 969, which are present in variants of concern and are targeted as important in protein function, are predicted as such from our model.

60 APPLIED LIFE SCIENCES↗