SEARCH · Search NASA

Results for “POINT MATCHING 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

TPSAS-NF1676L-32020-DND

The Earth Polychromatic Imaging Camera (EPIC) instrument aboard the Deep Space Climate Observatory (DSCOVR) satellite has a constant unique view of the sunlit disk of the Earth from the Lagrange-1 (L1) point nearly a million miles away from the Earth. The DSCOVR EPIC team released the version 2 dataset last year which included many enhancements, such as improved navigation, stray-light correction, and flat-fielding. Since there is no on-board calibration for EPIC, its 10 spectral channels must be calibrated using other methods, such as inter-calibrating with sensors with on-board calibration hardware. VIIRS is a well-calibrated instrument with similar spectral channels as EPIC. Spectral band induced radiance differences between VIIRS and EPIC are mitigated using Spectral Band Adjustment Factor (SBAF). We will extend the previous calibrations from prior studies to include more recent data by ray-matching EPIC with VIIRS over all-sky tropical ocean (ATO) and deep convective cloud (DCC) conditions. In addition to ATO and DCC ray-matching, we will validate these results by utilizing DCC as invariant targets to monitor sensor degradation. DCC behave as solar diffusers in the tropics and when collectively analyzed can be used to monitor the degradation of the EPIC sensor. Improvement to our dynamic navigation adjustment algorithm, as well as exploring the view angle response of the EPIC CCD array, will be examined for their effectiveness and presented at the conference.

C O Haney↗

NASA NTRS1991

Unified aeroacoustics analysis for high speed turboprop aerodynamics and noise. Volume 5: Propagation of propeller tone noise through a fuselage boundary layer

An analysis of tone noise propagation through a boundary layer and fuselage scattering effects was derived. This analysis is a three dimensional and the complete wave field is solved by matching analytical expressions for the incident and scattered waves in the outer flow to a numerical solution in the boundary layer flow. The outer wave field is constructed analytically from an incident wave appropriate to the source and a scattered wave in the standard Hankel function form. For the incident wave, an existing function - domain propeller noise radiation theory is used. In the boundary layer region, the wave equation is solved by numerical methods. The theoretical analysis is embodied in a computer program which allows the calculation of correction factors for the fuselage scattering and boundary layer refraction effects. The effects are dependent on boundary layer profile, flight speed, and frequency. Corrections can be derived for any point on the fuselage, including those on the opposite side from the source. The theory was verified using limited cases and by comparing calculations with available measurements from JetStar tests of model prop-fans. For the JetStar model scale, the boundary layer refraction effects produce moderate fuselage pressure reinforcements aft of and near the plane of rotation and significant attenuation forward of the plane of rotation at high flight speeds. At lower flight speeds, the calculated boundary layer effects result in moderate amplification over the fuselage area of interest. Apparent amplification forward of the plane of rotation is a result of effective changes in the source directivity due to boundary layer refraction effects. Full scale effects are calculated to be moderate, providing fuselage pressure amplification of about 5 dB at the peak noise location. Evaluation using available noise measurements was made under high-speed, high-altitude flight conditions. Comparisons of calculations made of free field noise, using a current frequency-domain propeller noise prediction method, and fuselage effects using this new procedure show good agreement with fuselage measurements over a wide range of flight speeds and frequencies. Correction factors for the JetStar measurements made on the fuselage are provided in an Appendix.

Magliozzi, B.↗

NASA NTRS1990

PAN AIR: A computer program for predicting subsonic or supersonic linear potential flows about arbitrary configurations using a higher order panel method. Volume 2: User's manual (version 3.0)

A comprehensive description of user problem definition for the PAN AIR (Panel Aerodynamics) system is given. PAN AIR solves the 3-D linear integral equations of subsonic and supersonic flow. Influence coefficient methods are used which employ source and doublet panels as boundary surfaces. Both analysis and design boundary conditions can be used. This User's Manual describes the information needed to use the PAN AIR system. The structure and organization of PAN AIR are described, including the job control and module execution control languages for execution of the program system. The engineering input data are described, including the mathematical and physical modeling requirements. Version 3.0 strictly applies only to PAN AIR version 3.0. The major revisions include: (1) inputs and guidelines for the new FDP module (which calculates streamlines and offbody points); (2) nine new class 1 and class 2 boundary conditions to cover commonly used modeling practices, in particular the vorticity matching Kutta condition; (3) use of the CRAY solid state Storage Device (SSD); and (4) incorporation of errata and typo's together with additional explanation and guidelines.

Sidwell, Kenneth W.↗

NASA NTRS1995

Analysis of the two-point velocity correlations in turbulent boundary layer flows

The general objective of the present work is to explore the use of Rapid Distortion Theory (RDT) in analysis of the two-point statistics of the log-layer. RDT is applicable only to unsteady flows where the non-linear turbulence-turbulence interaction can be neglected in comparison to linear turbulence-mean interactions. Here we propose to use RDT to examine the structure of the large energy-containing scales and their interaction with the mean flow in the log-region. The contents of the work are twofold: First, two-point analysis methods will be used to derive the law-of-the-wall for the special case of zero mean pressure gradient. The basic assumptions needed are one-dimensionality in the mean flow and homogeneity of the fluctuations. It will be shown that a formal solution of the two-point correlation equation can be obtained as a power series in the von Karman constant, known to be on the order of 0.4. In the second part, a detailed analysis of the two-point correlation function in the log-layer will be given. The fundamental set of equations and a functional relation for the two-point correlation function will be derived. An asymptotic expansion procedure will be used in the log-layer to match Kolmogorov's universal range and the one-point correlations to the inviscid outer region valid for large correlation distances.

Oberlack, M.↗

NASA NTRS1992

Diffraction by a symmetric material junction simulated with generalized sheet transistion conditions

When diffraction by a material discontinuity in a thick dielectric/ferrite slab is treated for the dual integral equation approach (a variation of the Weiner-Hopf method), the resulting solution displays unknown constants, characteristic of the incompleteness associated with the generalized transition conditions employed in modeling a thick slab. The constants can be shown to depend on the geometry and properties of the discontinuity. Their explicit determination by enforcing field continuity across the slab junction is treated. The field internal to the slab is determined by first expressing the field as an expansion of discrete and continuous eigenmodes valid in the exterior and interior slab region. The expansion or mode coefficients are then identified by recasting the Weiner-Hopf solution in a suitable form, enabling the solution constants to be numerically computed by matching the field expansions to the left and right of the slab discontinuity at a discrete set of points.

Ricoy, Mark A.↗

NASA NTRS1997

Turbulent Premixed Flame Propagation in Microgravity

A facility in which turbulent Couette flow could be generated in a microgravity environment was designed and built. To fit into the NASA Lewis drop tower the device had to be very compact. This means that edge effects and flow re-circulation were expected to affect the flow. The flow was thoroughly investigated using LDV and was found to be largely two dimensional away from the edges with constant turbulence intensities in the core. Slight flow asymmetries are introduced by the non symmetric re-circulation of the fluid outside the test region. Belt flutter problems were remedied by adding a pair of guide plates to the belt. In general, the flow field was found to be quite similar to previously investigated Couette flows. However, turbulence levels and associated shear stresses were higher. This is probably due to the confined re-circulation zone reintroducing turbulence into the test section. An estimate of the length scales in the flow showed that the measurements were able to resolve nearly all the length scales of interest. Using a new LES method for subgrid combustion it has been demonstrated that the new procedure is computational feasible even on workstation type environment. It is found that this model is capable of capturing the propagation of the premixed names by resolving the flame in the LES grid within 2-3 grid points. In contrast, conventional LES results in numerical smearing of the flame and completely inaccurate estimate of the turbulent propagation speed. Preliminary study suggests that there is observable effect of buoyancy in the 1g environment suggesting the need for microgravity experiments of the upcoming experimental combustion studies. With the cold flow properties characterized, an identical hot flow facility is under construction. It is assumed that the turbulence properties ahead of the flame in this new device will closely match the results obtained here. This is required since the hot facility will not enable LDV measurements. The reacting flow facility is also being constructed with planned drop tower experiments in mind. Therefore, issues related to safety and structural integrity are being take into account. Further development of the numerical model will also be carried out to include finite-rate kinetics for representative premixed cases. More detail analysis of the flame structure and propagation nature will be investigated. Simulations will also be compared to the flame properties observed in the experiments.

Menon, S.↗

NASA NTRS2010

Efficient Web Services Policy Combination

Large-scale Web security systems usually involve cooperation between domains with non-identical policies. The network management and Web communication software used by the different organizations presents a stumbling block. Many of the tools used by the various divisions do not have the ability to communicate network management data with each other. At best, this means that manual human intervention into the communication protocols used at various network routers and endpoints is required. Developing practical, sound, and automated ways to compose policies to bridge these differences is a long-standing problem. One of the key subtleties is the need to deal with inconsistencies and defaults where one organization proposes a rule on a particular feature, and another has a different rule or expresses no rule. A general approach is to assign priorities to rules and observe the rules with the highest priorities when there are conflicts. The present methods have inherent inefficiency, which heavily restrict their practical applications. A new, efficient algorithm combines policies utilized for Web services. The method is based on an algorithm that allows an automatic and scalable composition of security policies between multiple organizations. It is based on defeasible policy composition, a promising approach for finding conflicts and resolving priorities between rules. In the general case, policy negotiation is an intractable problem. A promising method, suggested in the literature, is when policies are represented in defeasible logic, and composition is based on rules for non-monotonic inference. In this system, policy writers construct metapolicies describing both the policy that they wish to enforce and annotations describing their composition preferences. These annotations can indicate whether certain policy assertions are required by the policy writer or, if not, under what circumstances the policy writer is willing to compromise and allow other assertions to take precedence. Meta-policies are specified in defeasible logic, a computationally efficient non-monotonic logic developed to model human reasoning. One drawback of this method is that at one point the algorithm starts an exhaustive search of all subsets of the set of conclusions of a defeasible theory. Although the propositional defeasible logic has linear complexity, the set of conclusions here may be large, especially in real-life practical cases. This phenomenon leads to an inefficient exponential explosion of complexity. The current process of getting a Web security policy from combination of two meta-policies consists of two steps. The first is generating a new meta-policy that is a composition of the input meta-policies, and the second is mapping the meta-policy onto a security policy. The new algorithm avoids the exhaustive search in the current algorithm, and provides a security policy that matches all requirements of the involved metapolicies.

Vatan, Farrokh↗

NASA NTRSDate not supplied

Comparing Tropopause-Penetrating Convection Identifications Derived from NEXRAD and GOES over the Contiguous United States

Overshooting cloud tops (OTs) are a result of powerful updrafts that rapidly transport air from the troposphere to the lower stratosphere. Aside from the well-known relationships between strong updrafts and severe weather, the boundaries of updrafts, including the OT, are turbulent, which leads to irreversible mixing of air between the stratosphere and troposphere, altering the composition of both layers. In order to better understand the impact of OTs, knowing precisely when and where they occur is vital. Within the last decade remote-sensing observations from satellites and radars have been used to identify OTs, but the results thus far have not been entirely consistent due to challenges with accurately identifying OTs with these sensor datasets. This study compares OTs identified using NEXRAD reflectivity data from the GridRad dataset with an approach that uses GOES infrared (IR) data. The GridRad method estimates the altitude of echo tops directly from the radar pointing geometry, merging a combination of radars that enables products every 5 mins, with 2 km spatial and 1 km vertical grid spacing. The GOES satellite method uses spatial patterns in tropopause-relative IR brightness temperature to derive an OT likelihood. The study region covers a large part of the contiguous United States during selected active convection dates in 2017 and includes data from GOES-13/16. A large sample of OTs observed by NEXRAD and GOES-13/16 satellites are analyzed to better understand NEXRAD and IR observations of OTs, quantify agreement between the OT detection methods, and demonstrate how the 4x increase in spatial resolution from GOES-13 to GOES-16 impacts OT algorithm detection performance. The results show that for nearly time-matched scenes, with many severe storms, GOES-13 and -16 detection accuracy is quite comparable; however, some missed detections are attributed to poorer GOES-13 image resolution. It was found that GOES-13 OT regions are ~2 K warmer than GOES-16 on average, but such differences exceed 5 K for the most intense OT regions. It is found that false GOES IR-based detections are unavoidable, given that cold cloud pixels significantly colder than the tropopause from nearby or recently decayed OTs can persist within anvils for 15 minutes or more.

John W. Cooney↗

NASA NTRS1990

Joint US/USSR study: Comparison of effects of horizontal and head-down bed rest

An account is given of the results of the first joint U.S./U.S.S.R. bed rest study. The study was accomplished in two parts: A soviet part (May to June 1979) and an American part (July to August 1979). Both studies were conducted under identical conditions and provided a basis for comparison of physiologic reactions and standardizing procedures and methods. Each experiment consisted of three periods: 14 days of pre-bed rest control, 7 days of bed rest, and a 10 to 14 day recovery period. Ten males participated in each study, with five subjects experiencing horizontal bed rest and five subjects a -6 deg head-down body position. Biochemical and hormonal measurements were made of blood and urine, with particular attention to electrolyte metabolism and kidney function; cardio-pulmonary changes at rest and exercise; influence of Lower Body Negative Pressure (LBNP); and incremental exercise using a bicyle ergometer while supine and sitting. Expected moderate changes were noted to occur for various physiologic parameters. Clinical evidence pointed to the fact that head-down bed rest when compared to horizontal conditions more closely matched the conditions seen after manned spaceflight. For the most part, statistically significant differences between the two body positions were not observed.

Sandler, Harold↗

NASA NTRS2019

Time Series Analysis in the Search for Other Worlds Through Transit Photometry

The Kepler Mission launched in June 2009 to commence NASA's first mission to search for potentially habitable, Earth-size planets orbiting Sun-like stars. Kepler discovered explanets via the transit method: searching for minute (100 ppm) drops in brightness lasting 1 - 13 hours corresponding to occasions where the planet crosses the face of its host star from Kepler's point of view. The exquisite precision required to carry out the Kepler mission (20 ppm in 6.5 hours) pushed astronomical time series analysis to the limits, and motivated the development of novel algorithmic approaches. Transit signatures of rocky planets are often dwarfed by the intrinsic stellar variability, which is not white noise, and often is non-stationary, and by instrumental systematic effects, which can include transients and electronic artifacts. Surmounting this challenging regime of weak, temporally compact, periodic signals in observation noise with strong systematics and other sources of variability motivated the development of 1) an overcomplete, non-decimated, wavelet-based matched filter to jointly estimate the properties of the non-stationary, non-white observation noise process, and 2) a multi-scale, maximum a posteriori (msMAP) approach to identifying and removing instrumental systematic effects. After over nine years of observations, the Kepler spacecraft finally ran out of fuel in November 2018, ending its data collection activities. Over 2300 planets were discovered by Kepler in its primary mission, and over 355 have been discovered by K2, the repurposed mission that followed Kepler's primary mission after the loss of a second reaction wheel in May 2013. We have ported the Kepler science pipeline for the Transiting Exoplanet Survey Satellite (TESS) Mission, which began science observations in July 2019, and report initial results and performance of the modified science pipeline.The Kepler and TESS Missions are supported by NASA's Science Mission Directorate.

transit surveys↗

NASA NTRS1961

Estimates of Minimum Energy Requirements for Range-Controlled Return of a Nonlifting Satellite from a Circular Orbit

Existing expressions are used to obtain the minimum propellant fraction required for return from a circular orbit as a function of vacuum trajectory range. trajectory are matched to those of the atmospheric trajectory to obtain a complete return from orbit to earth. The results are restricted by the assumptions of (1) impulsive velocity change, (2) nearly circular transfer trajectory, ( 3) spherical earth, atmosphere, and gravitational field, (4) exponential atmospheric density variation with attitude, and (5) a nonrotating atmosphere. The solutions for the parameters of the vacuum Calculations are made t o determine the effects of longitudinal and lateral range on required propeUant fraction and reentry loading for a nonrotating earth and for several orbital altitudes. the single- and two-impulse method of return is made and the results indicate a "trade off" between propellant fraction required and landing- position accuracy. A comparison of An example of a return mission from a polar orbit is discussed where the initial deorbit point is the intersection of the North Pole horizon with the satellite orbit. Some effects of a rotating earth are also considered. It is found that, for each target-orbital-plane longitudinal difference, there exists a target latitude for which the required propellant fraction is a minimum.

Jackson, Charlie M., Jr.↗

NASA NTRS2009

Axial and Centrifugal Compressor Mean Line Flow Analysis Method

This paper describes a method to estimate key aerodynamic parameters of single and multistage axial and centrifugal compressors. This mean-line compressor code COMDES provides the capability of sizing single and multistage compressors quickly during the conceptual design process. Based on the compressible fluid flow equations and the Euler equation, the code can estimate rotor inlet and exit blade angles when run in the design mode. The design point rotor efficiency and stator losses are inputs to the code, and are modeled at off design. When run in the off-design analysis mode, it can be used to generate performance maps based on simple models for losses due to rotor incidence and inlet guide vane reset angle. The code can provide an improved understanding of basic aerodynamic parameters such as diffusion factor, loading levels and incidence, when matching multistage compressor blade rows at design and at part-speed operation. Rotor loading levels and relative velocity ratio are correlated to the onset of compressor surge. NASA Stage 37 and the three-stage NASA 74-A axial compressors were analyzed and the results compared to test data. The code has been used to generate the performance map for the NASA 76-B three-stage axial compressor featuring variable geometry. The compressor stages were aerodynamically matched at off-design speeds by adjusting the variable inlet guide vane and variable stator geometry angles to control the rotor diffusion factor and incidence angles.

Veres, Joseph P.↗

NASA NTRS2019

Scaling Evaluation of Ice-Crystal Icing on a Modern Turbofan Engine in PSL Using the COMDES-MELT Code

This paper presents preliminary ice-crystal icing (ICI) altitude scaling evaluation results of a Honeywell Uncertified Research Engine (HURE) that was tested in the NASA Glenn Research Center Propulsion Systems Laboratory (PSL) during January of 2018. This engine geometry features a hidden core design to keep the core less exposed. The engine was fitted with internal video cameras to observe various ice buildup processes at multiple selected locations within the engine core flow path covering the fan stator, the splitter-lip/shroud/strut, and the high pressure compressor (HPC) variable inlet guide vane (IGV) regions. The potential ice accretion risk was pre-determined to occur by using NASA’s in-house 1D Engine Icing Risk assessment code, COMDES-MELT. The code was successful in predicting the risk of ice accretion in adiabatic regions like the fan-stator of the HURE at specific engine operating points. However at several operating points during the test, liquid water was observed running along the shroud toward the variable IGV of the HPC regions with an air temperature well below freezing, thus no particle melting could have occurred due to heating from the air alone. It was reasoned that other sources of heat were present in that region. To account for these heat sources the inlet total temperature was adjusted to give a wet bulb temperature of 24 F below the standard minimum wet bulb temperature of 492 R to allow ice to accrete in the splitter-lip/shroud/strut region, which was determined from a reference case where hard ice was observed in that region. With that adjustment the COMDES-MELT code was successful in providing operating points where there was a risk of ice accretion during the test campaign. In addition to calculating possible conditions at different selected lower altitudes, simulations were run to determine potential inlet conditions that could lead to ice-crystal accretion along the prescribed stations where the cameras were available. From there, scaled test conditions were determined by best matching the following three icing related parameters of the reference condition: (1) the local air total wet bulb temperature, (2) the local ice crystal cloud melt ratio and (3) the engine fan face ice/water to air mass flux ratio of the ice crystal cloud. Instantaneous images taken from the time-lapsed movies of ice buildup were used along with the relevant thermodynamic data of air, water vapor and local icing condition to help evaluate how closely the proposed altitude scaling method could be used in ground based test facility to duplicate selected reference ICI features observed at specific location inside this engine at different scale altitudes. Discussions on observed limitation for engine icing scaling application from this test campaign and needed improvement are provided. A scaling test procedure to help identify potential ICI risk conditions and possible ice accretion locations of a new turbofan engine is evaluated in PSL.

Engine Testing in PSL↗

NASA NTRS1992

Simultaneous mapping of the unsteady flow fields by Particle Displacement Velocimetry (PDV)

Current experimental and computational techniques must be improved in order to advance the prediction capability of the longitudinal vortical flows shed by underwater vehicles. The generation, development, and breakdown mechanisms of the shed vortices at high Reynolds numbers are not fully understood. The ability to measure hull separated vortices associated with vehicle maneuvering does not exist at present. The existing point-by-point measurement techniques can only capture approximately the large 'mean' eddies but fail to meet the dynamics of small vortices during the initial stage of generation. A new technique, which offers a previously unavailable capability to measure the unsteady cross-flow distribution in the plane of the laser light sheet, is called Particle Displacement Velocimetry (PDV). PDV consists of illuminating a thin section of the flowfield with a pulsed laser. The water is seeded with microscopic, neutrally buoyant particles containing imbedded fluorescing dye which responds with intense spontaneous fluorescence with the illuminated section. The seeded particles in the vortical flow structure shed by the underwater vehicle are illuminated by the pulse laser and the corresponding particle traces are recorded in a single photographic frame. Two distinct approaches were utilized for determining the velocity distribution from the particle traces. The first method is based on matching the traces of the same particle and measuring the distance between them. The direction of the flow can be identified by keeping one of the pulses longer than the other. The second method is based on selecting a small window within the image and finding the mean shift of all the particles within that region. The computation of the auto-correlation of the intensity distribution within the selected sample window is used to determine the mean displacement of particles. The direction of the flow is identified by varying the intensity of the laser light between pulses. Considerable computational resources are required to compute the auto-correction of the intensity distribution. Parallel processing will be employed to speed up the data reduction. A few examples of measured unsteady vortical flow structures shed by the underwater vehicles will be presented.

Huang, Thomas T.↗

NASA NTRS2001

Modeling of Fillets in Thin-Walled Structures for Dynamic Analysis

Relatively new developments in manufacturing methods have made it possible to produce machined parts with a wall thickness of less than 0.010". While the parts are being machined. Fillets of relatively large radius to thickness ratios are created by the milling tool. This paper discusses an accurate new technique for finite-element modeling of fillets for dynamic and stiffness analysis using efficient plate elements rather than more computationally intensive, high density solid meshes, A simple filleted cantilever beam of 0.040" thickness, 1.6" length, 1" depth, and 0.250" radius in the corner was modeled using both solid elements and several different plate element geometries. The finite element results were then compared with static and modal testing of a machined sample. The highest-scoring plate element technique uses a bridge of elements that is constructed through the tangent point of the Fillet radius such that the thickness of the element matches tile volume of the fillet. This model produces errors of less than 4.6 percent for static loading and less than 7.4 percent for modal analysis. This simple plate element technique will prove critical for efficient, timely, and accurate dynamic analysis of complex thin-walled structures.

Seugling, Richard M.↗

NASA NTRSDate not supplied

An Updated Secondary Lunar Meteoroid Ejecta Model for Engineering Design

Introduction: The surface of the Moon is constantly being bombarded by a flux of meteoroids of various sizes. Impacts due to these meteoroids produce secondary ejecta material at much lower speeds but with a total mass larger than the original impactor. Details about the secondary ejecta are important for planning missions on the lunar surface. In this work, an updated ejecta model is presented called the Meteoroid Model of Secondary Ejecta (MeMoSeE), to re-place the Apollo-era ejecta model, NASA SP-8013 [1], in the SLS-SPEC 159 Design Specification for Natural Environments (DSNE) [2]. The model pro-duces secondary ejecta flux environments for a user-specified location on the lunar surface, and sorts the incoming secondary flux by angular direction and speed. Methods: MeMoSeE is separated into three parts: the inputs, the conversion step, and the integration of fluxes. Inputs to the model include the primary meteoroid fluxes and the primary near-Earth object (NEO) fluxes. Meteoroid fluxes, both asteroidal and cometary, are calculated using the Meteoroid Engineering Model (MEM3) [3] for different locations on the Moon. For each surface location, an ephemeris is generated using the JPL HORIZONS System [4] that feeds into MEM3. The NEO fluxes are approximated by the high-density population of MEM3 (i.e., only the directionality), where the speed distribution of the NEO fluxes is renormalized to match observations [5]. The regolith properties are used as defined in the DSNE [2]. The conversion step utilizes scaling laws given by Housen & Holsapple 2011 [6] to convert the primary impactor flux to the total mass of secondary ejecta. The ejecta distribution, at the point-of-impact (POI), is separated into a zenith angle and azimuthal angle distribution. The zenith angle distribution follows a beta distribution where the peak depends on the impact altitude angle and the impact azimuth [7]. We employ an ejecta azimuth distribution that is based on Rival & Mandeville 1999 [8] which focuses ejecta in the downstream direction for more oblique impacts. Finally, during the integration step, we sum secondary ejecta number fluxes at a particular region-of-interest (ROI) that originated from many POI locations over the entire surface of the Moon. We keep track of both altitude and azimuth angle bins as well as a range of speed bins, following the igloo gridding as done in MEM3 [9]. The ejecta particle size distribution and density is assumed to be the same as the lunar regolith [10, 2]. Results: The primary fluxes are computed for one Metonic cycle (19 years) for various locations over the lunar surface with a fixed orientation. Both the angular and speed distributions of the primary fluxes are dependent on the latitude and longitude. The over-all primary fluxes show a roughly 13% increase from the eastern limb to the western limb. In general, the speed distribution of the secondary fluxes span from a user-defined minimum speed to the escape speed of the Moon (2.38 km/s), roughly following a power-law relation [6]. Different parts of the speed distribution come from different primary impact locations on the Moon. The secondary ejecta is dominated by the slowest speeds, where these particles originate nearby the ROI. For speeds around 71% of the escape speed, the secondary ejecta originates from locations near the antipodal point. On the other hand, for secondary ejecta speeds that exceed roughly 90% the escape speed, the ejecta particles come from all over the lunar surface to the ROI. Comparing the secondary ejecta fluxes from MeMoSeE with NASA SP-8013 [1], there is a reduction by about 2-3 orders of magnitude for secondary ejecta particles greater than 1 µg. These estimates agree with recent findings from Bjorkman & Christiansen 2019 [11]. The secondary ejecta fluxes are also compared with the primary fluxes, where the ejecta fluxes are roughly an order of magnitude greater than the primary fluxes. References: [1] Cour-Palais, B. G., (1969) NASA SP-8013. [2] NASA SLS-SPEC-159 Rev. H (2020). [3] Moorhead, A. V., et al. (2019) JS&R, 1-17. [4] Giorgini, J. D., (2015) IAUGA, 29, 2256293. [5] Moorhead, A. V., (2020) Memo OSMA/MEO/Lunar-001. [6] Housen, K. R., and Holsapple, K. A., (2011) Icarus, 211(1), 856-875. [7] Gault, D. E., and Wedekind, J. A., (1978) L&PSCP, 9, 3843-3875. [8] Rival, M., and Mandeville, J., (1999) Space Debris, 1(1), 45-57. [9] Moorhead, A. V., (2019) MEM3 User Guide. [10] Carrier, W. D., (1973) The Moon, 6(3-4), 250-263. [11] Bjorkman, M. D. & Christiansen, E. L., (2019) ODC

impacts↗

NASA NTRS2005

An Optimal Orthogonal Decomposition Method for Kalman Filter-Based Turbofan Engine Thrust Estimation

A new linear point design technique is presented for the determination of tuning parameters that enable the optimal estimation of unmeasured engine outputs such as thrust. The engine s performance is affected by its level of degradation, generally described in terms of unmeasurable health parameters related to each major engine component. Accurate thrust reconstruction depends upon knowledge of these health parameters, but there are usually too few sensors to be able to estimate their values. In this new technique, a set of tuning parameters is determined which accounts for degradation by representing the overall effect of the larger set of health parameters as closely as possible in a least squares sense. The technique takes advantage of the properties of the singular value decomposition of a matrix to generate a tuning parameter vector of low enough dimension that it can be estimated by a Kalman filter. A concise design procedure to generate a tuning vector that specifically takes into account the variables of interest is presented. An example demonstrates the tuning parameters ability to facilitate matching of both measured and unmeasured engine outputs, as well as state variables. Additional properties of the formulation are shown to lend themselves well to diagnostics.

Litt, Jonathan S.↗

NASA NTRS2007

An Optimal Orthogonal Decomposition Method for Kalman Filter-Based Turbofan Engine Thrust Estimation

A new linear point design technique is presented for the determination of tuning parameters that enable the optimal estimation of unmeasured engine outputs, such as thrust. The engine's performance is affected by its level of degradation, generally described in terms of unmeasurable health parameters related to each major engine component. Accurate thrust reconstruction depends on knowledge of these health parameters, but there are usually too few sensors to be able to estimate their values. In this new technique, a set of tuning parameters is determined that accounts for degradation by representing the overall effect of the larger set of health parameters as closely as possible in a least squares sense. The technique takes advantage of the properties of the singular value decomposition of a matrix to generate a tuning parameter vector of low enough dimension that it can be estimated by a Kalman filter. A concise design procedure to generate a tuning vector that specifically takes into account the variables of interest is presented. An example demonstrates the tuning parameters ability to facilitate matching of both measured and unmeasured engine outputs, as well as state variables. Additional properties of the formulation are shown to lend themselves well to diagnostics.

Litt, Jonathan S.↗