Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Normal form methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 343 records · Page 19

Use of T7 RNA polymerase to direct expression of outer Surface Protein A (OspA) from the Lyme disease Spirochete, Borrelia burgdorferi

The OspA gene from a North American strain of the Lyme disease Spirochete, Borrelia burgdorferi, was cloned under the control of transciption and translation signals from bacteriophage T7. Full-length OspA protein, a 273 amino acid (31kD) lipoprotein, is expressed poorly in Escherichia coli and is associated with the insoluble membrane fraction. In contrast, a truncated form of OspA lacking the amino-terminal signal sequence which normally would direct localization of the protein to the outer membrane is expressed at very high levels (less than or equal to 100 mg/liter) and is soluble. The truncated protein was purified to homogeneity and is being tested to see if it will be useful as an immunogen in a vaccine against Lyme disease. Circular dichroism and fluorescence spectroscopy was used to characterize the secondary structure and study conformational changes in the protein. Studies underway with other surface proteins from B burgdorferi and a related spirochete, B. hermsii, which causes relapsing fever, leads us to conclude that a strategy similar to that used to express the truncated OspA can provide a facile method for producing variations of Borrelia lipoproteins which are highly expressed in E. coli and soluble without exposure to detergents.

Dunn, John J.↗

90 GHz and 150 GHz Observations of the Orion M42 Region. A Submillimeter to Radio Analysis

We have used the new 90GHz MUSTANG camera on the Robert C. Green Bank Telescope (GBT)to map the bright Huygens region of the star-forming region M42 with a resolution of 9" and a sensitivity of 2.8 mJy/beam. Ninety GHz is an interesting transition frequency, as MUSTANG detects both the free-free emission characteristic of the H II region created by the Trapezium stars, normally seen at lower frequencies, and thermal dust emission from the background OMCI molecular cloud, normally mapped at higher frequencies. We also present similar data from the 150 GHz GISMO camera taken on the IRAM 30 m telescope. This map has 15" resolution. By combining the MUSTANG data with 1.4, 8. and 31 GHz radio data from the VLA and GBT, we derive a new estimate of the emission measure averaged electron temperature of T(sub e) = 11376+/-1050 K by an original method relating free-free emission intensities at optically thin and optically thick frequencies. Combining Infrared Space Observatory-long wavelength spectrometer (ISO-LWS) data with our data, we derive a new estimate of the dust temperature and spectral emissivity index within the 80" ISO-LWS beam toward Orion KL/BN, T(sub d) = 42+/-3 K and Beta(sub d) = 1.3+/-0.1. We show that both T(sub d) and Beta(sub d) decrease when going from the H II region and excited OMCI interface to the denser UV shielded part OMCI (Orion KL/BN, Orion S). With a model consisting of only free-free and thermal dust emission, we are able to fit data taken at frequencies from 1.5 GHz to 854 GHz (350 micrometers).

Dicker, S. R.↗

Thermal images collected during large-scale 3D printing with Ingersoll MasterPrint

This data set contains images produced to test the performance of anomaly and fault detection methods in the context of additive manufacturing. Two print jobs were executed using the Ingersoll MasterPrint with the Model 30 Strangpresse extruder. The material used was Techmer compounded polylactic acid (PLA) with wood flour as a filler (80/20 PLA/WF by weight). The first print job consists of the 3D printing of a hexagonal cylinder with a two-bead wall. This print was sliced at a gantry velocity of 3000 millimeter per minute, and an extruder screw speed of 68.14 rotations per minute. This are considered the normal operating conditions. A second hexagonal cylinder was printed using a reduced extruder screw speed, 15% lower than under normal operating conditions. The images were collected with a Teledyne FLIR Lepton 3.5 infra-red camera, a small form factor radiometric long-wave infrared camera with a spectral range of 8 µm to 14 µm. Sensor resolution was 160x120 pixels, with a pixel size of 12 µm, a temperature range of -10 - 450°C, and an accuracy of +/- 10°C in its low gain configuration. Images in this data set were collected with cameras oriented at the printer nozzle. The nozzle camera setup consisted of two cameras located 12.5cm from the nozzle center. These were mounted directly to the print head, so that the camera positions relative to the print direction would remain constant as it rotated around its C-axis to follow the print path. One camera was placed ahead of the nozzle to capture the previous layer immediately before being covered by the new layer of material, while the second camera was placed behind the nozzle and captured the freshly extruded bead. Images are collecting during each phase of the print: (a) idle (i.e., no material deposited), (b) extrusion (i.e., to prime the extruder), and (c) printing (deposition of material to manufacture the hexagonal cylinder).

36 MATERIALS SCIENCE↗

Quantitative form and fit of N95 filtering facepiece respirators are retained after dry and humid heat treatments for coronavirus deactivation

Re-use of filtering facepiece respirators (FFRs, commonly referred to as N95s) normally meant for single use only is becoming common in healthcare facilities due to shortages caused by the COVID19 pandemic. Here we report that mouse hepatitis virus (MHV) initially seeded on FFR filter material is inactivated (6 log reduction as measured by 50% tissue culture infective dose (TCID50)) after dry heating at 75 °C for 30 minutes. We also find that the quantitative fit of FFRs after heat treatment at this temperature, under dry conditions or at 90% relative humidity, is not affected by single or ten heating cycles. Previous studies have reported that the filtration efficiency of FFR filters is not negatively impacted by these heating conditions. These results suggest that thermal inactivation of coronaviruses is a potentially rapid and widely deployable method to re-use N95 FFRs in emergency situations where reusing FFRs is a necessity and broad-spectrum sterilization is unavailable. However, we also found that a heat source that emits radiation (e.g., an exposed heating element) results in rapid qualitative degradation of the FFR. Finally, we discuss differences in the results reported here and other recent studies investing heat as a means to recycle FFRs and suggest that overall wear time and donning/doffing cycles are important factors that need to be considered.

36 MATERIALS SCIENCE↗

Roll-to-Roll, Dual-Layer Slot Die Coating of Chitin and Cellulose Oxygen Barrier Films for Renewable Packaging

Cellulose and chitin are the two most abundant naturally produced biopolymers and are being extensively studied as candidates for renewable oxygen barrier films used in packaging. It has been shown that bilayers formed from cellulose nanocrystals (CNCs) and chitin nanofibers (ChNFs) exhibit oxygen barrier properties similar to polyethylene terephthalate (PET). However, this prior work explored only coating each layer individually in sequence through techniques such as spray coating. Here, we demonstrate the viability of dual-layer slot die coating of CNC/ChNF bilayers onto cellulose acetate (CA) substrates. The dual-layer slot die method enables significantly lower oxygen permeability versus spray coating while using a roll-to-roll system that applies the bilayer in a single pass. This work discusses suspension properties, wetting, and drying conditions required to achieve well-controlled ChNF/CNC bilayers. Spray-coated bilayer films were on average 25% thinner than the dual-layer bilayer film; however, the thickness-normalized oxygen permeability (OP) of the dual-layer-coated ChNF/CNC bilayer film on CA was 20 times better than that of the spray-coated bilayers. Additionally, it has been shown that ChNF contributes to the wetting and barrier properties. Values of OP for the slot die-coated bilayers under optimized drying conditions were as low as 1.2 cm 3 ∙μm∙m -2 ∙d -1 ∙kPa -1 , corresponding to a normalized oxygen transmission rate of 0.32 cm 3 ∙m -2 ∙d -1 at 23 °C and 50% relative humidity. It is also noted that the adhesive properties of the dual-layer coating are also improved when films are air-dried and that ChNF contributes to the wetting and barrier properties.

09 BIOMASS FUELS↗

Microhardness Distribution of Long Magnesium Block Processed through Powder Metallurgy

Powder metallurgy is a popular method of making raw powders into specific shaped samples. However, the pressure distribution and the microhardness difference within the sample are nonnegligible and unclear when the sample is long or exceeds a specific size. In this study, the long magnesium blocks, with a ratio of about 2.8 between the sample height and the sample side length, are successfully synthesized under three uniaxial and two biaxial conditions. Then, the sample hardness values on the outer surface and the center plane are tested to study the microhardness distribution. The modified analytical expression indicates that the normal pressure exponentially decreases along the compression direction, which is consistent with the hardness distribution trend. Because higher pressure leads to a more compact arrangement of the powders, more metal bonds are formed after sintering. During the first pressing, the sidewall pressure makes the surface hardness higher. The secondary reverse compression mainly improves the bottom and core hardness due to the re-orientation and re-location of the powders. The obtained relationship between the applied pressure and the hardness distribution is instructive in predicting and improving the sample quality.

Wang, Jiaying↗

Influence of Turbulence on the Restraint of Liquid Jets by Surface Tension in Microgravity Investigated

Microgravity poses many challenges to the designer of spacecraft tanks. Chief among these are the lack of phase separation and the need to supply vapor-free liquid or liquidfree vapor to the spacecraft processes that require fluid. One of the principal problems of phase separation is the creation of liquid jets. A jet can be created by liquid filling, settling of the fluid to one end of the tank, or even closing a valve to stop the liquid flow. Anyone who has seen a fountain knows that jets occur in normal gravity also. However, in normal gravity, the gravity controls and restricts the jet flow. In microgravity, with gravity largely absent, surface tension forces must contain jets. To model this phenomenon, a numerical method that tracks the fluid motion and the surface tension forces is required. Jacqmin has developed a phase model that converts the discrete surface tension force into a barrier function that peaks at the free surface and decays rapidly away. Previous attempts at this formulation were criticized for smearing the interface. This can be overcome by sharpening the phase function, double gridding the fluid function, and using a higher order solution for the fluid function. The solution of this equation can be rewritten as two coupled Poisson equations that also include the velocity. After the code was implemented in axisymmetric form and verified by several test cases at the NASA Glenn Research Center, the drop tower runs of Aydelott were modeled. Work last year with a laminar model was found to overpredict Aydelott's results, except at the lowest Reynolds number conditions of 400. This year, a simple turbulence model was implemented by adding a turbulent viscosity based on the mixing-length hypothesis and empirical measurements of previous works. Predictions made after this change was implemented have been much closer to experimentally observed flow patterns and geyser heights. Two model runs is shown. The first, without any turbulence correction, breaks through the free surface and strikes the far end of the tank. In the second, the turbulence spreads the jet momentum over more of the free surface, enabling the surface tension forces to turn the jet back into the bulk liquid. The model geyser height with the second model is 1.1 cm. This is quite close to the 1.5-cm geyser height measured by Aydelott.

Chato, David J.↗

CTF Validation and Verification: Version 4.3

Coolant-Boiling in Rod Arrays- Two Fluids (COBRA-TF) is a thermal/hydraulic (T/H) simulation code designed for light water reactor (LWR) analysis. It uses a two-fluid, three-field (i.e., fluid film, fluid drops, and vapor) modeling approach. Both subchannel and 3D Cartesian forms of nine conservation equations are available for LWR modeling. The code was originally developed by Pacific Northwest Laboratory in 1980 and has been used and modified by several institutions over the last several decades. COBRA-TF is also used at the Pennsylvania State University (PSU) by the Reactor Dynamics and Fuel Modeling Group (RDFMG) and has been improved, updated, and subsequently became the CTF. One part of the improvement process includes validating the methods in CTF. This document seeks to provide a certain level of certainty and confidence in the predictive capabilities of the code for the scenarios it was designed to model—rod bundle geometries with operating conditions that are representative of prototypical pressurized water reactor (PWR)s and boiling water reactor (BWR)s in both normal and accident conditions. This is done by modeling a variety of experiments that simulate these scenarios and then presenting a qualitative and quantitative analysis of the results that demonstrates the accuracy to which CTF is capable of capturing specific quantities of interest.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

CTF Validation and Verification: Version 4.4

Coolant-Boiling in Rod Arrays- Two Fluids (COBRA-TF) is a thermal/hydraulic (T/H) simulation code designed for light water reactor (LWR) analysis. It uses a two-fluid, three-field (i.e., fluid film, fluid drops, and vapor) modeling approach. Both subchannel and 3D Cartesian forms of nine conservation equations are available for LWR modeling. The code was originally developed by Pacific Northwest Laboratory in 1980 and has been used and modified by several institutions over the last several decades. COBRA-TF is also used at the Pennsylvania State University (PSU) by the Reactor Dynamics and Fuel Modeling Group (RDFMG) and has been improved, updated, and subsequently became the CTF. One part of the improvement process includes validating the methods in CTF. This document seeks to provide a certain level of certainty and confidence in the predictive capabilities of the code for the scenarios it was designed to model—rod bundle geometries with operating conditions that are representative of prototypical pressurized water reactor (PWR)s and boiling water reactor (BWR)s in both normal and accident conditions. This is done by modeling a variety of experiments that simulate these scenarios and then presenting a qualitative and quantitative analysis of the results that demonstrates the accuracy to which CTF is capable of capturing specific quantities of interest.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Disentangling AGN and Star Formation in Soft X-Rays

We have explored the interplay of star formation and active galactic nucleus (AGN) activity in soft X-rays (0.5-2 keV) in two samples of Seyfert 2 galaxies (Sy2s). Using a combination of low-resolution CCD spectra from Chandra and XMM-Newton, we modeled the soft emission of 34 Sy2s using power-law and thermal models. For the 11 sources with high signal-to-noise Chandra imaging of the diffuse host galaxy emission, we estimate the luminosity due to star formation by removing the AGN, fitting the residual emission. The AGN and star formation contributions to the soft X-ray luminosity (i.e., L(sub x,AGN) and L(sub x,SF)) for the remaining 24 Sy2s were estimated from the power-law and thermal luminosities derived from spectral fitting. These luminosities were scaled based on a template derived from XSINGS analysis of normal star-forming galaxies. To account for errors in the luminosities derived from spectral fitting and the spread in the scaling factor, we estimated L(sub x,AGN) and L(sub x,SF))from Monte Carlo simulations. These simulated luminosities agree with L(sub x,AGN) and L(sub x,SF) derived from Chandra imaging analysis within a 3sigma confidence level. Using the infrared [Ne ii]12.8 micron and [O iv]26 micron lines as a proxy of star formation and AGN activity, respectively, we independently disentangle the contributions of these two processes to the total soft X-ray emission. This decomposition generally agrees with L(sub x,SF) and L(sub x,AGN) at the 3 sigma level. In the absence of resolvable nuclear emission, our decomposition method provides a reasonable estimate of emission due to star formation in galaxies hosting type 2 AGNs.

Disentangling↗

Method and apparatus for stable silicon dioxide layers on silicon grown in silicon nitride ambient

A method and apparatus for thermally growing stable silicon dioxide layers on silicon is disclosed. A previously etched and baked silicon nitride tube placed in a furnace is used to grow the silicon dioxide. First, pure oxygen is allowed to flow through the tube to initially coat the inside surface of the tube with a thin layer of silicon dioxide. After the tube is coated with the thin layer of silicon dioxide, the silicon is oxidized thermally in a normal fashion. If the tube becomes contaminated, the silicon dioxide is etched off thereby exposing clean silicon nitride and then the inside of the tube is recoated with silicon dioxide. As is disclosed, the silicon nitride tube can also be used as the ambient for the pyrolytic decomposition of silane and ammonia to form thin layers of clean silicon nitride.

Cohen, R. A.↗

VIPER project

The VIPER project has so far produced a formal specification of a 32 bit RISC microprocessor, an implementation of that chip in radiation-hard SOS technology, a partial proof of correctness of the implementation which is still being extended, and a large body of supporting software. The time has now come to consider what has been achieved and what directions should be pursued in the future. The most obvious lesson from the VIPER project was the time and effort needed to use formal methods properly. Most of the problems arose in the interfaces between different formalisms, e.g., between the (informal) English description and the HOL spec, between the block-level spec in HOL and the equivalent in ELLA needed by the low-level CAD tools. These interfaces need to be made rigorous or (better) eliminated. VIPER 1A (the latest chip) is designed to operate in pairs, to give protection against breakdowns in service as well as design faults. We have come to regard redundancy and formal design methods as complementary, the one to guard against normal component failures and the other to provide insurance against the risk of the common-cause failures which bedevil reliability predictions. Any future VIPER chips will certainly need improved performance to keep up with increasingly demanding applications. We have a prototype design (not yet specified formally) which includes 32 and 64 bit multiply, instruction pre-fetch, more efficient interface timing, and a new instruction to allow a quick response to peripheral requests. Work is under way to specify this device in MIRANDA, and then to refine the spec into a block-level design by top-down transformations. When the refinement is complete, a relatively simple proof checker should be able to demonstrate its correctness. This paper is presented in viewgraph form.

Kershaw, John↗

On a new method for calculating the potential flow past a body of revolution

A new method is presented for obtaining the velocity potential of the flow about a body of revolution moving uniformly in the direction of its axis of symmetry in a fluid otherwise at rest. This method is based essentially on the fact that the form of the differential equation for the velocity potential is invariant with regard to conformal transformation of the meridian plane. By means of the conformal transformation of the meridian profile into a circle a system of orthogonal curvilinear coordinates is obtained, the main feature of which is that one of the coordinate lines is the meridian profile itself. The use of this type of coordinate system yields a simple expression of the boundary condition at the surface of the solid and leads to a rational process of iteration for the solution of the differential equation for the velocity potential. It is shown that the velocity potential for an arbitrary body of revolution may be expressed in terms of universal functions which, although not normal, are obtainable by means of simple quadratures. The general results are applied to a body of revolution obtained by revolving a symmetrical Joukowski profile about its axis of symmetry. A numerical example further serves to illustrate the theory.

Kaplan, Carl↗

Thermophysical properties of simple liquid metals: A brief review of theory

In this paper, we review the current theory of the thermophysical properties of simple liquid metals. The emphasis is on thermodynamic properties, but we also briefly discuss the nonequilibrium properties of liquid metals. We begin by defining a 'simple liquid metal' as one in which the valence electrons interact only weakly with the ionic cores, so that the interaction can be treated by perturbation theory. We then write down the equilibrium Hamiltonian of a liquid metal as a sum of five terms: the bare ion-ion interaction, the electron-electron interaction, the bare electron-ion interaction, and the kinetic energies of electrons and ions. Since the electron-ion interaction can be treated by perturbation, the electronic part contributes in two ways to the Helmholtz free energy: it gives a density-dependent term which is independent of the arrangement of ions, and it acts to screen the ion-ion interaction, giving rise to effective ion-ion pair potentials which are density-dependent, in general. After sketching the form of a typical pair potential, we briefly enumerate some methods for calculating the ionic distribution function and hence the Helmholtz free energy of the liquid: monte Carlo simulations, molecular dynamics simulations, and thermodynamic perturbation theory. The final result is a general expression for the Helmholtz free energy of the liquid metal. It can be used to calculate a wide range of thermodynamic properties of simple metal liquids, which we enumerate. They include not only a range of thermodynamic coefficients of both metals and alloys, but also many aspects of the phase diagram, including freezing curves of pure elements and phase diagrams of liquid alloys (including liquidus and solidus curves). We briefly mention some key discoveries resulting from previous applications of this method, and point out that the same methods work for other materials not normally considered to be liquid metals (such as colloidal suspensions, in which the suspended microspheres behave like ions screened by the salt solution in which they are suspended). We conclude with a brief discussion of some non-equilibrium (i.e., transport) properties which can be treated by an extension of these methods. These include electrical resistivity, thermal conductivity, viscosity, atomic self-diffusion coefficients, concentration diffusion coefficients in alloys, surface tension and thermal emissivity. Finally, we briefly mention two methods by which the theory might be extended to non-simple liquid metals: these are empirical techniques (i.e., empirical two- and three-body potentials), and numerical many-body approaches. Both may be potentially applicable to extremely complex systems, such as nonstoichiometric liquid semiconductor alloys.

Stroud, David↗

Helium droplet infrared spectroscopy of the butyl radicals

Butyl radicals (n-, s-, i-, and tert-butyl) are formed from the pyrolysis of stable precursors (1-pentyl nitrite, 2-methyl-1-butyl nitrite, isopentyl nitrite, and azo-tert-butane, respectively). The radicals are doped into a beam of liquid helium droplets and probed with infrared action spectroscopy from 2700 to 3125 cm -1 , allowing for a low temperature measurement of the CH stretching region. The presence of anharmonic resonance polyads in the 2800–3000 cm -1 region complicates its interpretation. To facilitate spectral assignment, the anharmonic resonances are modeled with two model Hamiltonian approaches that explicitly couple CH stretch fundamentals to HCH bend overtones and combinations: a VPT2+K normal mode model based on coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] quartic force fields and a semi-empirical local mode model. Both of these computational methods provide generally good agreement with the experimental spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Using Fuzzy Clustering for Real-time Space Flight Safety

To ensure space flight safety, it is necessary to monitor myriad sensor readings on the ground and in flight. Since a space shuttle has many sensors, monitoring data and drawing conclusions from information contained within the data in real time is challenging. The nature of the information can be critical to the success of the mission and safety of the crew and therefore, must be processed with minimal data-processing time. Data analysis algorithms could be used to synthesize sensor readings and compare data associated with normal operation with the data obtained that contain fault patterns to draw conclusions. Detecting abnormal operation during early stages in the transition from safe to unsafe operation requires a large amount of historical data that can be categorized into different classes (non-risk, risk). Even though the 40 years of shuttle flight program has accumulated volumes of historical data, these data don t comprehensively represent all possible fault patterns since fault patterns are usually unknown before the fault occurs. This paper presents a method that uses a similarity measure between fuzzy clusters to detect possible faults in real time. A clustering technique based on a fuzzy equivalence relation is used to characterize temporal data. Data collected during an initial time period are separated into clusters. These clusters are characterized by their centroids. Clusters formed during subsequent time periods are either merged with an existing cluster or added to the cluster list. The resulting list of cluster centroids, called a cluster group, characterizes the behavior of a particular set of temporal data. The degree to which new clusters formed in a subsequent time period are similar to the cluster group is characterized by a similarity measure, q. This method is applied to downlink data from Columbia flights. The results show that this technique can detect an unexpected fault that has not been present in the training data set.

Lee, Charles↗

The Impact of Europium as a Fission Product Surrogate on Chromium Corrosion in Molten Chloride Salt

The rising interest in molten salts, both as thermal energy storage media and as high-temperature solvents for fissile materials in advanced nuclear reactors, has spurred significant growth in research on structural and functional salt-facing materials. Impurities, such as water, play a significant role in the degradation of materials in molten salt environments. Other impurities, such as fission products, are expected to accumulate in molten salts during normal reactor operation and can affect the properties of the salt and salt-facing materials. During normal reactor operation, rare-earth element species are expected to form as byproducts of fission reactions, potentially degrading or altering the performance of both structural and functional materials exposed to salts. Here, in this work, fission product surrogates, europium metal and EuCl 3 , were added to the KCl-MgCl 2 salt while exposing pure Cr specimens in the melt at 600 and 700 °C for 1000 h. Low reactivity was observed with the baseline salt, while both the Eu and EuCl 3 additions to the salt resulted in increased mass loss of Cr specimens. A capillary electrophoresis technique was developed to enable the detection of low parts-per-million (ppm) concentrations of europium (Eu) in chloride salt both before and after exposure of Cr metal. This technique provides a sensitive and reliable method for monitoring trace levels of Eu, which is critical for understanding the interactions and potential contamination effects in salt-facing materials during reactor operation. This quantification technique can be applied to other cationic and anionic contaminants in halide salts.

chloride salt↗

Finite Difference Time Marching in the Frequency Domain: A Parabolic Formulation for Aircraft Acoustic Nacelle Design

An explicit finite difference iteration scheme is developed to study harmonic sound propagation in aircraft engine nacelles. To reduce storage requirements for large 3D problems, the time dependent potential form of the acoustic wave equation is used. To insure that the finite difference scheme is both explicit and stable, time is introduced into the Fourier transformed (steady-state) acoustic potential field as a parameter. Under a suitable transformation, the time dependent governing equation in frequency space is simplified to yield a parabolic partial differential equation, which is then marched through time to attain the steady-state solution. The input to the system is the amplitude of an incident harmonic sound source entering a quiescent duct at the input boundary, with standard impedance boundary conditions on the duct walls and duct exit. The introduction of the time parameter eliminates the large matrix storage requirements normally associated with frequency domain solutions, and time marching attains the steady-state quickly enough to make the method favorable when compared to frequency domain methods. For validation, this transient-frequency domain method is applied to sound propagation in a 2D hard wall duct with plug flow.

Baumeister, Kenneth J.↗