Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistent discretization”

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 127 records · Page 7

Electron and Proton Auroral Dynamics

Data from the Wide-band Imaging Camera (WIC) sensitive to far ultraviolet auroras and from the Spectrographic Imager (SI) channel SI12, sensitive to proton precipitation induced Lyman alpha were analyzed during a high altitude orbit segment of the IMAGE spacecraft. This segment began during the expansive phase of a substorm. The aurora changed into a double oval configuration, consisting of a set of discrete pole-ward forms and a separate diffuse auroral oval equatorwards, Although IMF Bz was strongly southward considerable activity could be seen poleward of the discrete auroras in the region that was considered to be the polar cap. The SI12 Doppler shifted Lyman alpha signature of precipitating protons show that the proton aurora is on the equatorward side of the diffuse aurora. In the following several hours the IMF Bz field changed signed. Although the general character of the proton and electron aurora did not change, the dayside aurora moved equatorward when the Bz was negative and more bright aurora was seen in the central polar cap during periods of positive Bz.

Mende, S. B.↗

Assessing Organic Preservation and the Implications for Potential Biosignatures in the Bastide Member of the Séítah Formation, Jezero Crater

Introduction: Olivine has the highest CO2 trapping potential of ultramafic minerals, due to its rapid rate of dissolution and high percentage of divalent cations/unpolymerized silicate tetrahedra [1]. It generates divalent carbonates from CO2 and sequesters CO2 into the mineral matrix of the target lithology. The co-occurrence of olivine and carbonate within abraded targets from the Bastide Member of the Séítah Formation (Fm) of the Jezero crater floor [2] suggests the carbonation of olivine occurred within the mineral matrix hosted in a subsurface system. Hydrothermal origins for the subsurface system are hypothesized from orbital data [3]. This is supported by the detection of hydration features within the rock as well as the carbonate features are solely detected within the abrasion patch but not the rock’s surface. Organic preservation potential of abrasion patches: We incorporated the SHERLOC/WATSON results acquired from the Dourbes, Garde, and Quartier abrasion patches in the Séítah Fm to investigate the organomineral associations, and determine the biosignature preservation potential of these rocks. Dourbes is dominated by olivine and has minor amounts of carbonate, hydrated Ca-sulfate, and amorphous or microcrystalline silicate. Fluorescence features (330-340 nm) are detected in discrete locales and could be consistent with double ring aromatic organic molecules; yet, these features do not appear to be associated with an identified mineral phase. Dark subhedral to euhedral olivine grains within the Garde abrasion patch often co-occur with carbonate-consistent spectral signatures in all analyzed scans. The availability of Fe2+ is a known influence on olivine dissolution rates [1] and SuperCam estimates of the olivine composition (Forsterite-60 average for Sols 202-234) may thereby provide a constraint on the carbonation extent. Carbonated olivine phenocrysts within the matrix may be due to aqueous alteration, as the carbonation of nodules is consistent with observations in other hydrothermal systems and within Martian meteorites (ALH84001) [4]. In comparison, the Quartier abrasion patch represents an extensively altered endmember within the Bastide Member of the Séítah Fm. It contains a fluorescence doublet at 305/325 nm that coexists with multiple species of Na-sulfate and Mg-sulfate, Mg-carbonates, olivine [5]. Aqueous alteration and implications for habitability: The identification of primary and secondary mineral phases observed in the Bastide Member suggests the release of cations from primary ultramafic minerals through aqueous alteration. Within Garde, the carbonate detected appears to be Mg-rich and likely formed from the in situ alteration of Mg-rich olivine within the region as carbonate has only been detected within the rock via in situ analysis. Carbonates derived from abiotic and biotic reactions preserve biosignatures (i.e. indices of habitability) on Earth. The detected carbonate phase found in association with fluorescent features within the rock matrix may also indicate potential organic compounds preserved in a putative hydrothermal system. Fluorescent features (~330 nm) are unique to Garde and Dourbes, though they are co-located to carbonate signatures solely within Garde, and between the light toned minerals. Identification of these fluorescent features may be consistent with 1-2 ring aromatic compounds. The limit of detection for Raman is multiple orders of magnitude greater than the limit of detection required for fluorescence [6]. Implications for provenance: Hydrothermal systems represent disequilibrium chemical conditions that are hypothesized to have supported the emergence of life and also preserve ancient carbon within precipitated carbonates [7]. Hydrothermal system associated carbonates are capable of preserving biosignatures up to an estimated ~3.77 – 4.28 Gya [8]. Thus, carbonated olivine found within the Séítah Fm may represent a high-potential biosignature preserving environment on Mars. The formation of carbonates by an aqueous alteration process, such as carbonation of olivine is also consistent with hypotheses for carbonate within the greater regional-olivine bearing unit [2,3], which contains Garde, Dourbes, and Quartier. Acknowledgments: This work was carried out at the Jet Propulsion Laboratory, The California Institute of Technology under a contract from NASA. References: [1] Wood et al., (2019) ES&T, 6, 10. [2] Stack, K. et al., (2020) Space Sci Rev, 216, 127. [3] Tarnas, J. et al., (2021) JGR: Planets, 126, 11. [4] Steele et al., (2007) Meteorit. Planet. Sci., 42, 9. [5] Murphy, A.E. et al., (2022) LPSC [6] Bhartia et al., (2021) Space Sci Rev, 217, 58. [7] Luther (2021) GRL, 48, e2021GL094869. [8] Dodd et al., (2017) Nature, 543, 60-64

E L Cardarelli↗

Volatile Abundances, Extended Coma Sources, and Nucleus Ice Associations in Comet C/2014 Q2 (Lovejoy)

High-resolution infrared spectra of comet C/2014 Q2 Lovejoy were acquired with NIRSPEC at the W. M. Keck Observatory on two post-perihelion dates(UT 2015 February 2 and 3).H2O was measured simultaneously with CO,CH3OH, H2CO, CH4, C2H6,C2H4, C2H2, HCN, and NH3 on both dates, and rotational temperatures, production rates, relative abundances, H2O ortho-to-para ratios, and spatial distributions in the coma were determined. The first detection of C2H4in a comet from ground-based observations is reported. Abundances relative to H2O for all species were found to be in the typical range compared with values for other comets in the overall population to date. There is evidence of variability in rotational temperatures and production rates on timescales that are small compared with the rotational period of the comet. Spatial distributions of volatiles in the coma suggest complex outgassing behavior.CH3OH, HCN, C2H6, and CH4 spatial distributions in the coma are consistent with direct release from associated ices in the nucleus and are peaked in a more sunward direction compared with co-measured dust. H2O spatial profiles are clearly distinct from these other four species, likely due to a sizable coma contribution from icy grain sublimation.Spatial distributions for C2H2,H2CO, and NH3suggest substantial contributions from extended coma sources, providing further evidence for distinct origins and associations for these species in comets. CO shows a different spatial distribution compared with other volatiles, consistent with jet activity from discrete nucleus ice sources.

Neil Dello Russo↗

Goal-Oriented Discretization Error Control in Coupled Nearfield-Farfield Low-Boom Simulations

The method of adjoint weighted residuals is used to determine the level of discretization error in loudness predictions of sonic booms on the ground. We analyze the standard nearfied-farfield domain decomposition approach. In the nearfield domain, the three-dimensional Euler equations are solved to obtain a pressure signature generated by the aircraft. In the farfield, this waveform is propagated through the atmosphere to the ground by solving the augmented Burgers’ equation. Loudness is characterized using weighted sound-exposure-level metrics. We formulate discretization error estimates for the ground signature and loudness metrics for this one-way coupled system. Although the nearfield solution is independent of the farfield, the adjoint formulation for the coupled system provides feedback from the farfield to identify high-error regions in the nearfield. The results demonstrate that the discrete adjoint implementation is asymptotically consistent and provides reliable error estimates. Furthermore, we show how the error can be controlled through adaptive refinement of the nearfield mesh. The approach is evaluated on two- and three-dimensional problems, including the X-59 flight demonstration aircraft.

CST↗

A Multi-Disciplinary Analysis Framework for the Design of Small Launch Vehicles

Between the years of 1995 and 2014 the number of small satellites (1-500 kg) went up from 20 to 180. [1] Out of the 180 launched in 2014 66% were Nano satellites (1-10 kg). [1] With this trend of smaller satellites, one would expect a rise in number of small launch vehicles (SLVs are defined by capability to carry 1-100 kg to orbit), but this has not happened: [3] only 8% of all small satellites are launched on SLVs and the others become secondary payloads on regular launch vehicles. [1] This results in small satellites being placed into either suboptimal orbits or waiting for launch schedules to align with bigger launches, resulting in long waiting times. Dedicated SLVs could improve the responsiveness of small satellite launches, but SLV design is complicated by the large architecture space needed to be explored in order to find efficient and affordable designs. The SLV architecture trade space contains many discrete options, e.g. solid vs. liquid fuels, air launch vs. ground launch, number of stages, etc. [2][3] Performing detailed analysis for all the options at once would be prohibitive. Thus, a sizing environment/framework that is capable of providing necessary information for conceptual-level trade studies and can rapidly explore the vast SLV architecture space is necessary. The framework, illustrated in Figure 1, consists of four disciplines central to the sizing and synthesis of launch vehicles: propulsion, aerodynamics, structures, and trajectory. For the aerodynamics and trajectory disciplines, Missile DATCOM and POST2 are used, respectively. The propulsion and structures disciplines are represented in the framework with tools developed at ASDL Georgia Tech. For propulsion, the Solid Motor Analysis Code (SMAC) is a physics-based conceptual design tool for solid rocket motors. SMAC is capable of geometric burn simulation, ballistic analysis, and prediction of thrust performance. [4] For structures, Launch Vehicle Structural Analysis (LVSA) tool is a physics-based tool that focuses on structural dynamic analysis with sizing capability. These tools are integrated into the framework illustrated in Figure 1 with the corresponding connections described in Table 1. The process flow is as follows: first, SMAC sizes insulation and calculates maximum operating pressures for each vehicle stage. The MEOP and insulation thickness values from SMAC are fed to LVSA which then utilizes this information to size the motor casing. This creates a feedback loop between LVSA and SMAC that converges on the radius available for fuel, casing, and insulation thickness. Once the stage sizing is converged on SMAC creates an engine deck that is passed to POST2. Next, the data from SMAC and LVSA goes into Missile DATCOM which generates an aerodynamics database for POST2. Finally, POST2 performs a targeting optimization while maximizing the payload mass to orbit. Within the framework, POST2 is automated in order to be robust to a wide variety of possible designs by performing a Monte Carlo simulation over the initialization vector of the POST2 optimization variables. To demonstrate the capability of this framework, a sample problem of exploring the design space of an SLV capable of placing satellites into a low Earth orbit (inclination=47 deg, 196.5 by 211.3 nm) is used. This sample problem is a ground-launched SLV consisting of four in-line SRM stages. The multidisciplinary design analysis (MDA) environment is utilized to explore a design space consisting of 22 continuous and 12 discrete variables, shown by Table 1 by blocks 1,2,3. Running a full factorial design of experiments (DOE) would have been prohibitively expensive even with this reduced design space, thus a space filling design with 3,502 and then additional expansion of 2,602 cases was used. The first DOE consisted of 3,502 cases, and all of the variables were varied. These input variables are listed in blocks 1, 2, and 3 in Table 1. Most of the variables are propulsion related with stage lengths determining the delta-V split of the SLV. The expansion consisted of 2,602 cases, and the continuous variables were set to be equal to the most promising designs from the sizing of first DOE. For each set of continuous variables, the discrete variables (grain type, star points, and propellant) were to varied. The results of the DOE can be seen in Figure 2; each of the points in this plot represents a closed launch vehicle that reaches the targeted orbit. For each of the cases, there is information on flown trajectory, structural, and propulsion properties of the SLV. For example, Figure 3 shows changes in altitude and velocity with time for a particular case. The right side of Figure 3 clearly shows the coasting (slow decrease) and burning phases (sharp increase) of the SLV mission. LV mass is positively correlated with the optimized payload mass to orbit because heavier LVs carry more fuel and thus have more stored chemical energy. Furthermore, for any given payload mass to orbit, the most efficient design would result in the smallest LV. The results as visualized in Figure 2 shows this tradeoff, and the Pareto frontier of the efficient designs can be seen along the dotted line. Figure 2 can be divided into regions with the lowest mass vehicles corresponding to the minimum bound on radius, and the highest mass vehicles corresponding to the maximum bound on radius. Within a mass region, the discrete variables, such as propellant type and propellant grain arrangement, have the most effect on payload mass. This paper presents an MDA framework that can perform an automated physics-based sizing of SLV designs and a corresponding methodology to utilize the MDA to explore the design space of SLVs. The proposed methodology was applied to a perform a design space exploration for a four stage SLV. The outputs show the expected pareto frontier forming and provide detailed information about the SLV performance and staging. Using the produced data, it will be possible to select a set of pareto optimal designs that can then be further explored in subsequent design cycles. This demonstration shows that automated design space exploration should be used in the early phase design of future SLV concepts.

Nikita S Birbasov↗

X-33 Hypersonic Boundary Layer Transition

Boundary layer and aeroheating characteristics of several X-33 configurations have been experimentally examined in the Langley 20-Inch Mach 6 Air Tunnel. Global surface heat transfer distributions, surface streamline patterns, and shock shapes were measured on 0.013-scale models at Mach 6 in air. Parametric variations include angles-of-attack of 20-deg, 30-deg, and 40-deg; Reynolds numbers based on model length of 0.9 to 6.6 million; and body-flap deflections of 0, 10 and 20-deg. The effects of discrete and distributed roughness elements on boundary layer transition, which included trip height, size, location, and distribution, both on and off the windward centerline, were investigated. The discrete roughness results on centerline were used to provide a transition correlation for the X-33 flight vehicle that was applicable across the range of reentry angles of attack. The attachment line discrete roughness results were shown to be consistent with the centerline results, as no increased sensitivity to roughness along the attachment line was identified. The effect of bowed panels was qualitatively shown to be less effective than the discrete trips; however, the distributed nature of the bowed panels affected a larger percent of the aft-body windward surface than a single discrete trip.

Berry, Scott A.↗

The Foggy EUV Corona and Coronal Heating by MHD Waves from Explosive Reconnection Events

In 0.5 arcsec/pixel TRACE coronal EUV images, the corona rooted in active regions that are at the limb and are not flaring is seen to consist of (1) a complex array of discrete loops and plumes embedded in (2) a diffuse ambient component that shows no fine structure and gradually fades with height. For each of two not-flaring active regions, found that the diffuse component is (1) approximately isothermal and hydrostatic and (2) emits well over half of the total EUV luminosity of the active-region corona. Here, from a TRACE Fe XII coronal image of another not-flaring active region, the large sunspot active region AR 10652 when it was at the west limb on 30 July 2004, we separate the diffuse component from the discrete loop component by spatial filtering, and find that the diffuse component has about 60% of the total luminosity. If under much higher spatial resolution than that of TRACE (e. g., the 0.1 arcsec/pixel resolution of the Hi-C sounding-rocket experiment proposed by J. W. Cirtain et al), most of the diffuse component remains diffuse rather being resolved into very narrow loops and plumes, this will raise the possibility that the EUV corona in active regions consists of two basically different but comparably luminous components: one being the set of discrete bright loops and plumes and the other being a truly diffuse component filling the space between the discrete loops and plumes. This dichotomy would imply that there are two different but comparably powerful coronal heating mechanisms operating in active regions, one for the distinct loops and plumes and another for the diffuse component. We present a scenario in which (1) each discrete bright loop or plume is a flux tube that was recently reconnected in a burst of reconnection, and (2) the diffuse component is heated by MHD waves that are generated by these reconnection events and by other fine-scale explosive reconnection events, most of which occur in and below the base of the corona where they are seen as UV explosive events, EUV blinkers, and type II spicules. These MHD waves propagate across field lines and dissipate, heating the plasma in the field between the bright loops and plumes.

Moore, Ron L.↗

Viewing hybrid systems as products of control systems and automata

The purpose of this note is to show how hybrid systems may be modeled as products of nonlinear control systems and finite state automata. By a hybrid system, we mean a network of consisting of continuous, nonlinear control system connected to discrete, finite state automata. Our point of view is that the automata switches between the control systems, and that this switching is a function of the discrete input symbols or letters that it receives. We show how a nonlinear control system may be viewed as a pair consisting of a bialgebra of operators coding the dynamics, and an algebra of observations coding the state space. We also show that a finite automata has a similar representation. A hybrid system is then modeled by taking suitable products of the bialgebras coding the dynamics and the observation algebras coding the state spaces.

Grossman, R. L.↗

Integrated Multilayer Insulation

Integrated multilayer insulation (IMLI) is being developed as an improved alternative to conventional multilayer insulation (MLI), which is more than 50 years old. A typical conventional MLI blanket comprises between 10 and 120 metallized polymer films separated by polyester nets. MLI is the best thermal- insulation material for use in a vacuum, and is the insulation material of choice for spacecraft and cryogenic systems. However, conventional MLI has several disadvantages: It is difficult or impossible to maintain the desired value of gap distance between the film layers (and consequently, it is difficult or impossible to ensure consistent performance), and fabrication and installation are labor-intensive and difficult. The development of IMLI is intended to overcome these disadvantages to some extent and to offer some additional advantages over conventional MLI. The main difference between IMLI and conventional MLI lies in the method of maintaining the gaps between the film layers. In IMLI, the film layers are separated by what its developers call a micro-molded discrete matrix, which can be loosely characterized as consisting of arrays of highly engineered, small, lightweight, polymer (typically, thermoplastic) frames attached to, and placed between, the film layers. The term "micro-molded" refers to both the smallness of the frames and the fact that they are fabricated in a process that forms precise small features, described below, that are essential to attainment of the desired properties. The term "discrete" refers to the nature of the matrix as consisting of separate frames, in contradistinction to a unitary frame spanning entire volume of an insulation blanket.

Dye, Scott↗

HyDE Framework for Stochastic and Hybrid Model-Based Diagnosis

Hybrid Diagnosis Engine (HyDE) is a general framework for stochastic and hybrid model-based diagnosis that offers flexibility to the diagnosis application designer. The HyDE architecture supports the use of multiple modeling paradigms at the component and system level. Several alternative algorithms are available for the various steps in diagnostic reasoning. This approach is extensible, with support for the addition of new modeling paradigms as well as diagnostic reasoning algorithms for existing or new modeling paradigms. HyDE is a general framework for stochastic hybrid model-based diagnosis of discrete faults; that is, spontaneous changes in operating modes of components. HyDE combines ideas from consistency-based and stochastic approaches to model- based diagnosis using discrete and continuous models to create a flexible and extensible architecture for stochastic and hybrid diagnosis. HyDE supports the use of multiple paradigms and is extensible to support new paradigms. HyDE generates candidate diagnoses and checks them for consistency with the observations. It uses hybrid models built by the users and sensor data from the system to deduce the state of the system over time, including changes in state indicative of faults. At each time step when observations are available, HyDE checks each existing candidate for continued consistency with the new observations. If the candidate is consistent, it continues to remain in the candidate set. If it is not consistent, then the information about the inconsistency is used to generate successor candidates while discarding the candidate that was inconsistent. The models used by HyDE are similar to simulation models. They describe the expected behavior of the system under nominal and fault conditions. The model can be constructed in modular and hierarchical fashion by building component/subsystem models (which may themselves contain component/ subsystem models) and linking them through shared variables/parameters. The component model is expressed as operating modes of the component and conditions for transitions between these various modes. Faults are modeled as transitions whose conditions for transitions are unknown (and have to be inferred through the reasoning process). Finally, the behavior of the components is expressed as a set of variables/ parameters and relations governing the interaction between the variables. The hybrid nature of the systems being modeled is captured by a combination of the above transitional model and behavioral model. Stochasticity is captured as probabilities associated with transitions (indicating the likelihood of that transition being taken), as well as noise on the sensed variables.

Narasimhan, Sriram↗

Nonparametric probability density estimation by optimization theoretic techniques

Two nonparametric probability density estimators are considered. The first is the kernel estimator. The problem of choosing the kernel scaling factor based solely on a random sample is addressed. An interactive mode is discussed and an algorithm proposed to choose the scaling factor automatically. The second nonparametric probability estimate uses penalty function techniques with the maximum likelihood criterion. A discrete maximum penalized likelihood estimator is proposed and is shown to be consistent in the mean square error. A numerical implementation technique for the discrete solution is discussed and examples displayed. An extensive simulation study compares the integrated mean square error of the discrete and kernel estimators. The robustness of the discrete estimator is demonstrated graphically.

Scott, D. W.↗

On reliable control system designs

A mathematical model for use in the design of reliable multivariable control systems is discussed with special emphasis on actuator failures and necessary actuator redundancy levels. The model consists of a linear time invariant discrete time dynamical system. Configuration changes in the system dynamics are governed by a Markov chain that includes transition probabilities from one configuration state to another. The performance index is a standard quadratic cost functional, over an infinite time interval. The actual system configuration can be deduced with a one step delay. The calculation of the optimal control law requires the solution of a set of highly coupled Riccati-like matrix difference equations. Results can be used for off-line studies relating the open loop dynamics, required performance, actuator mean time to failure, and functional or identical actuator redundancy, with and without feedback gain reconfiguration strategies.

Birdwell, J. D.↗

Electrodynamics of the westward traveling surge

It is shown that the global convection pattern, the ionospheric current, and the field-aligned current associated with the westward traveling surge in the asymptotic state can be modeled quantitatively as consequences of a blockage of the Hall current from closure in the magnetosphere via field-aligned currents. The conductivity is allowed to increase self-consistently with increasing upward field-aligned current in the model. This inclusion of the self-consistent enhanced ionospheric conductivity due to discrete auroral precipitations is found to generate a localized intense westward electrojet on the poleward side of the Harang discontinuity. The westward electrojet is also found to rotate counterclockwise, merging into the eastward electrojet around the leading edge of the surge. Thus the major features of the westward traveling surge can be reproduced reasonably well in the model.

Kan, J. R.↗

Shock modelling in transonic and supersonic flow

Numerical simulation studies of problems in transonic flow are conducted in connection with the design of aircraft flying near the speed of sound, while recently there has also been a great deal of interest in supersonic and hypersonic flow calculations. Three important models for inviscid compressible gasdynamics are utilized in the study of these flow problems, taking into account the transonic small disturbance equation (TSD), the transonic full potential equation (FP), and the full Euler equations for compressible gasdynamics (EU). The present paper has the objective to provide a description of shock capturing methods recently developed for each model. The design principles behind the algorithms are also examined, giving attention to consistency, conservative form, monotone and sharp discrete shocks, and a bound on the variation, which prevents spurious oscillations.

Osher, S.↗

Optimal rejection of bounded persistent disturbances in periodic systems

The problem of optimal rejection of bounded persistent disturbances is solved in the case of linear discrete-time period systems. The solution consists of solving an equivalent time-invariant standard l1 optimization problem subject to an additional constraint. This constraint ensures the causality of the resulting periodic controller. By the duality theory, the problem is shown to be equivalent to a linear programming problem, which is no harder than the standard l1 problem.

Dahleh, Munther A.↗

Optimal and robust controllers for periodic and multirate systems

The problem of optimal rejection of bounded persistent disturbances is solved in the case of linear discrete-time periodic systems. The solution consists of solving an equivalent time-invariant standard l1 optimization problem subject to an additional constraint. This constraint assures the causality of the resulting periodic controller. By the duality theory, the problem is shown to be equivalent to a linear programming problem, which is no harder than the standard l1 problem. Also, it is shown that the method of solution presented applies exactly to the problem of disturbance rejection in the case of multirate sampled data systems. Finally, the results are applied to the problem of robust stabilization of periodic and multirate systems.

Dahleh, Munther A.↗

Preliminary design of an intermittent smoke flow visualization system

A prototype intermittent flow visualization system that was designed to study vortex flow field dynamics has been constructed and tested through its ground test phase. It produces discrete pulses of dense white smoke consisting of particles of terephthalic acid by the pulsing action of a fast-acting three-way valve. The trajectories of the smoke pulses can be tracked by a video imaging system without intruding in the flow around in flight. Two methods of pulsing the smoke were examined. The simplest and safest approach is to simply divert the smoke between the two outlet ports on the valve; this approach should be particularly effective if it were desired to inject smoke at two locations during the same test event. The second approach involves closing off one of the outlet ports to momentarily block the flow. The second approach requires careful control of valve dwell times to avoid excessive pressure buildup within the cartridge container. This method also increases the velocity of the smoke injected into the flow. The flow of the smoke has been blocked for periods ranging from 30 to 80 milliseconds, depending on the system volume and the length of time the valve is allowed to remain open between valve closings.

Ward, Donald T.↗

TBIEM3D: A Computer Program for Predicting Ducted Fan Engine Noise

This document describes the usage of the ducted fan noise prediction program TBIEM3D (Thin duct - Boundary Integral Equation Method - 3 Dimensional). A scattering approach is adopted in which the acoustic pressure field is split into known incident and unknown scattered parts. The scattering of fan-generated noise by a finite length circular cylinder in a uniform flow field is considered. The fan noise is modeled by a collection of spinning point thrust dipoles. The program, based on a Boundary Integral Equation Method (BIEM), calculates circumferential modal coefficients of the acoustic pressure at user-specified field locations. The duct interior can be of the hard wall type or lined. The duct liner is axisymmetric, locally reactive, and can be uniform or axially segmented. TBIEM3D is written in the FORTRAN programming language. Input to TBIEM3D is minimal and consists of geometric and kinematic parameters. Discretization and numerical parameters are determined automatically by the code. Several examples are presented to demonstrate TBIEM3D capabilities.

Dunn, M. H.↗