Search NASA⌕ Search

SEARCH · Search NASA

Results for “transient storage”

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 289 records · Page 16

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.↗

Time-Dependent Parabolic Finite Difference Formulation for Harmonic Sound Propagation in a Two-Dimensional Duct with Flow

An explicit finite difference real time iteration scheme is developed to study harmonic sound propagation in aircraft engine nacelles. To reduce storage requirements for future 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 for a harmonic monochromatic sound field, a parabolic (in time) approximation is introduced to reduce the order of the governing equation. The analysis begins with a harmonic sound source radiating into a quiescent duct. This fully explicit iteration method then calculates stepwise in time to obtain the 'steady state' harmonic solutions of the acoustic field. For stability, applications of conventional impedance boundary conditions requires coupling to explicit hyperbolic difference equations at the boundary. 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.

Kreider, Kevin L.↗

Finite Difference Time Marching in the Frequency Domain: A Parabolic Formulation for the Convective Wave Equation

An explicit finite difference iteration scheme is developed to study harmonic sound propagation in ducts. 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, K. J.↗

Reversion to the mean of the Solar Wind Input to the Magnetosphere

Earth’s magnetosphere responds to solar wind input over time scales ranging from minutes to days. On minute times scales, researchers seek to identify discontinuities in the solar wind that may trigger transient magnetospheric phenomena such as flux transfer events or substorm onsets. On hour time scales, researchers seek to determine the solar wind phenomena that drive the growth phase of geomagnetic substorms and the storage of magnetic flux within the magnetotail. On longer time scales, researchers seek to determine the solar wind features than propel geomagnetic storms through their initial, main, and recovery phases. Researchers must, perforce, rely upon spacecraft located near the L1 libration point for observations of the solar wind input into the magnetosphere over periods ranging from hours to days. Some of the hazards of doing this are well known. Features seen far upstream at locations away from the Sun-Earth line may never arrive at Earth due to structure in the solar wind transverse to the Sun-Earth line. Alternatively, they may arrive at unexpected times due to errors in calculating their orientation and velocity. Other hazards when predicting the magnetospheric response to varying solar wind conditions are less well understood. When observations of large solar wind drivers at spacecraft far upstream and off the Sun-Earth line are used to select events for study, the ‘reversion to the mean’ effect associated with spatial structure and errors in arrival timing for solar wind features at Earth demands that the magnetosphere will actually be subjected to and respond to solar wind drivers with much lower strengths. Consequently, the magnetospheric response must saturate for the highest solar wind drivers.

David Sibeck↗

Enhanced Dynamics and Charge Transport at the Eutectic Point: A New Paradigm for the Use of Deep Eutectic Solvent Systems

Deep eutectic solvents (DESs) are a class of versatile solvents with promise for a wide range of applications, from separation processes to electrochemical energy storage technologies. A fundamental understanding of the correlation among the structure, thermodynamics, and dynamics of these materials necessary for targeted rational design for specific applications is still nascent. Here, we employ differential scanning calorimetry (DSC), broadband dielectric spectroscopy (BDS), and femtosecond transient absorption spectroscopy (fs-TAS) to investigate the correlation among thermodynamics, dynamics, and charge transport in mixtures comprising a wide range of compositions of choline chloride (ChCl) and ethylene glycol (EG). Detailed analyses reveal that (i) the eutectic composition of this prototypical DES occurs in the 15–20 mol % ChCl in the EG range rather than the previously assumed 33 mol %, and (ii) both rotational dynamics and charge transport at the eutectic composition are enhanced in this composition range. These findings highlight the fundamental interplay between thermodynamics and dynamics in determining the properties of DESs that are relevant to many applications.

25 ENERGY STORAGE↗

Thermal Conductivity of a Novel Salt-Hydrogel Complex as a Phase Change Material for Building Thermal Management

Use of electricity for space heating, ventilation and air conditioning in buildings may be better managed through temporary energy storage. Compared with electrochemical storage, the levelized cost of thermal energy storage systems can be lower depending on geographical location and thermophysical properties [1]. In any phase change material (PCM), the rate of absorption and release of thermal energy increases with higher thermal conductivity [2]. However, leading PCMs like paraffin waxes have thermal conductivities on the order of 0.1 W/mK. Glauber’s salt, a salt hydrate, possesses a relatively high thermal conductivity but suffers from supercooling and phase segregation. Research has shown that polymer-based composites with salt hydrates can prevent phase segregation and limit supercooling over hundreds of cycles [3] though the thermal conductivity of the composite may be reduced leading to low energy storage rates. Here, we report thermal conductivity measurements on polymer hydrogels which are stable in solutions of Glauber’s salt. To measure thermal conductivity of soft gels, we use a 3-omega method [4] where heat flows bidirectionally through a substrate below and through a gel placed on top. This suits measurements of both gels and liquids unlike transient plane heat source or transient hot wire methods. We report measurements of thermal conductivity of salt-hydrogel complexes as a function of temperature, cross-linking and salt concentration. In addition, we present modeling of thermal conductivity of sodium sulfate solutions with available theories. We finally employ polarized optical microscopy to observe inside the transparent hydrogel networks, showing that precipitating crystals may influence thermal conductivity measurements for high salt concentrations. Our thermal conductivity results combined with enthalpy of phase change, mass density and specific heat are essential to accurately design thermal storage systems for energy-efficient buildings.

composite, thermal conductivity↗

Deep learning model for fast, science-based forecasting of fluid migration along faults in geologic carbon storage scenarios

Effective long-term geologic storage depends on robust site selection and credible, science-based forecasting of subsurface behavior to ensure storage integrity. For this work, we develop a deep learning–based reduced-order model (ROM) to quantify potential carbon dioxide (CO₂) and brine migration through geological faults. The ROM combines a Transformer model for binary classification and a Stacked Ensemble for regression, trained on a comprehensive dataset generated from 1400 physics-based reservoir simulations. Key geologic and operational parameters—including fault geometry, reservoir structure, and injection conditions—were systematically varied to capture a wide range of fluid migration scenarios. The ROM accurately predicts the onset of migration, cumulative migration volumes of both CO₂ and brine, and associated migration rates, as compared to an independent set of validation simulations, while significantly reducing computational cost compared to traditional simulation methods. Model performance was evaluated across diverse fault configurations, revealing that shallow reservoir geometry and fault angle are among the most influential factors governing migration behavior. Sensitivity analysis using SHapley Additive exPlanations (SHAP) provided interpretability, revealing distinct patterns in how geological and operational features drive transient versus cumulative migration outcomes. The ROM’s ability to rapidly simulate fault migration scenarios enables efficient sensitivity analyses, scenario evaluations, and decision support for site selection and monitoring design. This approach enhances the safety, scalability, and long-term operational performance of geologic carbon storage (GCS) systems by providing a robust, interpretable tool for predicting subsurface fluid migration and assessing fault-related migration potential.

42 ENGINEERING↗

Overview of the distributed image processing infrastructure to produce the Legacy Survey of Space and Time

The Vera C. Rubin Observatory is preparing to execute the most ambitious astronomical survey ever attempted, the Legacy Survey of Space and Time (LSST). Currently the final phase of construction is under way in the Chilean Andes, with the Observatory’s ten-year science mission scheduled to begin in 2025. Rubin’s 8.4-meter telescope will nightly scan the southern hemisphere collecting imagery in the wavelength range 320–1050 nm covering the entire observable sky every 4 nights using a 3.2 gigapixel camera, the largest imaging device ever built for astronomy. Automated detection and classification of celestial objects will be performed by sophisticated algorithms on high-resolution images to progressively produce an astronomical catalog eventually composed of 20 billion galaxies and 17 billion stars and their associated physical properties. In this article we present an overview of the system currently being constructed to perform data distribution as well as the annual campaigns which reprocess the entire image dataset collected since the beginning of the survey. These processing campaigns will utilize computing and storage resources provided by three Rubin data facilities (one in the US and two in Europe). Each year a Data Release will be produced and disseminated to science collaborations for use in studies comprising four main science pillars: probing dark matter and dark energy, taking inventory of solar system objects, exploring the transient optical sky and mapping the Milky Way. Also presented is the method by which we leverage some of the common tools and best practices used for management of large-scale distributed data processing projects in the high energy physics and astronomy communities. We also demonstrate how these tools and practices are utilized within the Rubin project in order to overcome the specific challenges faced by the Observatory.

79 ASTRONOMY AND ASTROPHYSICS↗

Transient Convection Due to Imposed Heat Flux: Application to Liquid-Acquisition Devices

A model problem is considered that addresses the effect of heat load from an ambient laboratory environment on the temperature rise of liquid nitrogen inside an enclosure. This model has applications to liquid acquisition devices inside the cryogenic storage tanks used to transport vapor-free propellant to the main engine. We show that heat loads from Q = 0.001 to 10 W, with corresponding Rayleigh numbers from Ra = 109 to 1013, yield a range of unsteady convective states and temperature rise in the liquid. The results show that Q = 1 to 10 W (Ra = 1012 to 1013) yield temperature distributions along the enclosure height that are similar in trend to experimental measurements. Unsteady convection, which shows selfsimilarity in its planforms, is predicted for the range of heat-load conditions. The onset of convection occurs from a free-convection-dominated base flow that becomes unstable against convective instability generated at the bottom of the enclosure while the top of the enclosure is convectively stable. A number of modes are generated with small-scale thermals at the bottom of the enclosure in which the flow selforganizes into two symmetric modes prior to the onset of the propagation of the instability. These symmetric vertical modes transition to asymmetric modes that propagate as a traveling-wave-type motion of convective modes and are representative of the asymptotic convective state of the flow field. Intense vorticity production is created in the core of the flow field due to the fact that there is shear instability between the vertical and horizontal modes. For the higher Rayleigh numbers, 1012 to 1013, there is a transition from a stationary to a nonstationary response time signal of the flow and temperature fields with a mean value that increases with time over various time bands and regions of the enclosure.

Duval, Walter M. B.↗

Analysis and Design of Crew Sleep Station for ISS

This paper details the analysis and design of the Temporary Sleep Station (TeSS) environmental control system for International Space Station (ISS). The TeSS will provide crewmembers with a private and personal space, to accommodate sleeping, donning and doffing of clothing, personal communication and performance of recreational activities. The need for privacy to accommodate these activities requires adequate ventilation inside the TeSS. This study considers whether temperature, carbon dioxide, and humidity within the TeSS remain within crew comfort and safety levels for various expected operating scenarios. Evaluation of these scenarios required the use and integration of various simulation codes. An approach was adapted for this study, whereby results from a particular code were integrated with other codes when necessary. Computational Fluid Dynamics (CFD) methods were used to evaluate the flow field inside the TeSS, from which local gradients for temperature, velocity, and species concentration such as CO (sub 2) could be determined. A model of the TeSS, containing a human, as well as equipment such as a laptop computer, was developed in FLUENT, a finite-volume code. Other factors, such as detailed analysis of the heat transfer through the structure, radiation, and air circulation from the TeSS to the US Laboratory Aisle, where the TeSS is housed, were considered in the model. A complementary model was developed in G189A, a code which has been used by NASA/JSC for environmental control systems analyses since the Apollo program. Boundary conditions were exchanged between the FLUENT and G189A TeSS models. G189A provides human respiration rates to the FLUENT model, while the FLUENT model provides local convective heat transfer coefficients to G189A model. An additional benefit from using an approach with both a systems simulation and CFD model, is the capability to verify the results of each model by comparison to the results of the other model. The G189A and FLUENT models were used to evaluate various ventilation designs for the TeSS over a range of operating conditions with varying crew metabolic load, equipment operating modes, ventilation flow rates, and with the TeSS doors open and closed. Results from the study were instrumental in the optimization of a design for the TeSS ventilation hardware. A special case was considered where failure of the TeSS ventilation system occurred. In this case, a study was conducted in order to determine the time required for the CO (sub 2) concentration inside the TeSS to increase to ISS limit values under transient conditions. A lumped-capacitance code, SINDA-FLUINT was used in this case to provide accurate predictions of the human reaction to the TeSS cabin conditions including core and skin temperatures and body heat storage. A simple two-dimensional CFD model of a crewmember inside the TeSS was developed in FLUENT in order to determine the volume envelope of the respired air from the human, which maintained a minimum velocity profile. This volume was then used in the SINDA-FLUINT model to facilitate the calculations of CO (sub 2) concentrations, dry bulb temperatures and humidity levels inside the TeSS.

Keener, John F.↗

Across‐model spread and shrinking in predicting peatland carbon dynamics under global change

Abstract Large across‐model spread in simulating land carbon (C) dynamics has been ubiquitously demonstrated in model intercomparison projects (MIPs), and became a major impediment in advancing climate change prediction. Thus, it is imperative to identify underlying sources of the spread. Here, we used a novel matrix approach to analytically pin down the sources of across‐model spread in transient peatland C dynamics in response to a factorial combination of two atmospheric CO 2 levels and five temperature levels. We developed a matrix‐based MIP by converting the C cycle module of eight land models (i.e., TEM, CENTURY4, DALEC2, TECO, FBDC, CASA, CLM4.5 and ORCHIDEE) into eight matrix models. While the model average of ecosystem C storage was comparable to the measurement, the simulation differed largely among models, mainly due to inter‐model difference in baseline C residence time. Models generally overestimated net ecosystem production (NEP), with a large spread that was mainly attributed to inter‐model difference in environmental scalar. Based on the sources of spreads identified, we sequentially standardized model parameters to shrink simulated ecosystem C storage and NEP to almost none. Models generally captured the observed negative response of NEP to warming, but differed largely in the magnitude of response, due to differences in baseline C residence time and temperature sensitivity of decomposition. While there was a lack of response of NEP to elevated CO 2 (eCO 2 ) concentrations in the measurements, simulated NEP responded positively to eCO 2 concentrations in most models, due to the positive responses of simulated net primary production. Our study used one case study in Minnesota peatland to demonstrate that the sources of across‐model spreads in simulating transient C dynamics can be precisely traced to model structures and parameters, regardless of their complexity, given the protocol that all the matrix models were driven by the same gross primary production and environmental variables.

Hou, Enqing↗

Nuclide inventory characterization of EBR-II MOX fuel test pins

Characterization of 3 mixed oxide (MOX) fuel pins from the Experimental Breeder Reactor II (EBR-II) SPA-2/-2B experiment program was performed in support of the Advanced Reactor Experiments for Sodium Fast Reactor Fuels (ARES) joint project between Idaho National Laboratory (INL) and the Japan Atomic Energy Agency (JAEA). The characterization efforts were performed to enable transient testing in the Transient Reactor Test (TREAT) facility at INL and will include post irradiation examination (PIE) and other analyses. MOX fuel pins UW02009, 02011, and 06024 were historically irradiated from 1989 to 1994 for a total of 822.42 effective full power days (EFPDs) and peak burnup of approximately 14.3 at.% (atomic percent of heavy metal atoms, ∼134.4 GWd/t). They had since been placed in storage until a future use was identified. Characterization of the fission products and activation of each fuel pin was necessary to enable removal from storage and subsequent irradiation testing or PIE activities. The initial fuel pin geometry and composition were established using historical documentation. Irradiation history details were utilized to compute irradiated and decayed fuel pin masses, atomic densities, and activities. A pellet-wise axial neutron flux was also computed to assist in evaluation of the axial distribution of {sup 235}U and {sup 239}Pu, per pellet, within each of the test fuel pins. The results computed herein supported experiment design for advanced high-burnup MOX fuel designs under slow transient overpower (MOXTOP) conditions, which was Phase I of the ARES project. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Solvation dynamics of choline fluoride in ethylene glycol – Water mixtures

Due to their favorable physical properties, deep eutectic solvents (DESs) have been of interest in many applications in the energy and industry sectors. In this study, we investigate the new deep eutectic solvent ethalineF (a 1:2 M ratio of choline fluoride:ethylene glycol) and the effect of water addition on the solvation dynamics compared to the better known ethaline (a 1:2 M ratio of choline chloride:ethylene glycol). Femtosecond transient absorption spectroscopy and NMR diffusometry were used to study the solvation dynamics. 19 F NMR spectra were obtained for ethalineF solutions with different water additions. The 19 F NMR spectra show two fluorine signals; one can be assigned to fluoride anions solvated by EG and the other signal to fluoride anions associated with the much heavier and less mobile choline cation. Density, viscosity, conductivity, and E T (30) polarity measurements are also provided for ethalineF. The nuclear spin nature of fluorine atom permits to measure the self-diffusion coefficient of the negative ion in addition of the cation diffusion in a choline halide based deep eutectic solvent system. This allows to calculate the NMR-predicted conductivities σ NMR . The obtained ion dissociation coefficients α D at different water concentrations suggest a partial ion dissociation which was commensurate with the relatively low ion conductivity. Furthermore, the observed modest conductivity of ethalineF is explained in terms of ion association as well as the gradual increase in conductivity caused by the addition of water.

19F NMR↗

Failure Analysis for Molten Salt Thermal Energy Storage Tanks for In-Service CSP Plants

Thermal Energy Storage (TES) is a fundamental component in concentrating solar power (CSP) plants to increase the plant's dispatchability, capacity factor, while reducing the levelized cost of electricity. In central receivers CSP plants, nitrate molten salts have been used for several years for operation temperatures of up to 565 degrees C. Despite many efforts to advance nitrate salt to higher operation temperatures (even considering a replacement with molten chloride salts) to achieve higher energy conversion efficiencies, the 565 degrees C temperature is currently considered the state-of-the art. Although molten salt tanks have been broadly deployed in commercial CSP plants worldwide, several failures have been reported in these tanks after a few months or years of operation, causing significant economic loss and mistrust in CSP technologies. Most of these failures are associated with the infancy of the technology and multiple issues related to tank design, fabrication, commissioning, and aggressive operation. A technical standard dedicated to the design and fabrication of molten nitrate TES tanks does not exist today. Current in-service molten salt tanks have been generally designed based on the American Petroleum Institute's (API) 650 and ASME Section II standards. The API 650 code provides guidelines for dimensions and fabrication for oil storage tanks up to 260 degrees C. The ASME standard provides allowable stress values for various materials at a range of temperatures and conditions. Both standards seem to be limited for molten salt TES tanks where high temperatures, thermal cycling, and transient conditions are expected. In 2020, NREL released the Concentrating Solar Power Best Practices Study (NREL/TP-5500-75763) that summarized multiple issues in CSP plants, along with potential alternatives and recommendations to address those issues based on information collected from participants representing about 80% of operating CSP plants in the world. One of the recommendations from this study was the development of accurate and validated models to evaluate the plant's transient operation, capable of capturing the effect of short-term clouds and operator response, while being flexible in being adapted to various spatial and temporal resource data. The "Failure Analysis for Molten Salt Thermal Energy Tanks for In-Service CSP Plants" project was inspired on this recommendation and was focused on (1) the development and validation of a physics-based model for a representative, commercial-scale molten salt tank, (2) performing simulations to evaluate the behavior of the tank as a function of typical plant operation conditions, (3) understanding tank failures mechanisms, (4) determining the residual stress and distortion in the tank floor after welding fabrication and evaluating their impact in the stresses developed in the tank during operation, (5) assessing the impact of key operation parameters on the temperature and stress distribution, (6) conduct a preliminary evaluation of design features to reduce stress and improve tank's reliability, and (7) estimate the tank's service life based on the stress developed under diverse operation scenarios. From the analysis conducted in the project and presented in this report, it was found that maximum stresses surpassing the yield strength point of the stainless steel (SS) 347H are developed on the tank floor near the perimeter. These large stresses are strongly influenced by the initial residual stresses and distortion of the tank floor after welding fabrication. During operation, large stresses are developed in the tank floor at high operation temperatures with large salt inventory levels during transient operation. High stresses are also related to elevated temperature gradients in the tank floor that could be attributed to insufficient mixing within the salt inflow and the salt inventory. Based on the analysis, creep is the predominant failure mechanism. However, the large stress levels could favor the plastic deformation into buckles, and crack formation due to stress relaxation cracking during cycle operation. A lifetime below 3 years was estimated for the typical plant operation conditions and a specific initial residual stress and deformation distribution of the tank floor. The estimated life agrees with the service time to failure reported in several commercial molten salt tanks. Desing and operation guidelines can be extracted from the analysis presented in this report, which could be adopted by tank manufacturers and CSP operators to advance toward an ultimate solution for tank failures by reducing residual and operational stresses to achieve a tank service life of more than 30 years. Addressing failures in molten salt TES tanks is fundamental for the CSP industry's survivability, but it is also important for other industrial and power generation applications using this technology, including nuclear and concentrating solar thermal.

14 SOLAR ENERGY↗

Tundra SGS – Critical Importance of Site-Specific Injectivity Testing for Carbon Dioxide Injection Well Simulation Calibration

Paper presented at 16th International Conference on Greenhouse Gas Control Technologies (GHGT-16), October 23–27, 2022, Lyon, France. Brine injectivity tests were conducted on both the Broom Creek (BC) and Deadwood (DW) storage formations in one of the cored stratigraphic test wells prior to finalized simulations. The injectivity tests entailed injection of brine into the formations, followed by shutting in the wells. The shut-in duration represents the fall-off period which was used to carry out pressure transient analyses of the tests. The analysis of injectivity tests indicated that core-indicated and modeled permeability distributions should be conservatively multiplied by 2.5× and 5× in the BC and DW geologic models, respectively, prior to numerical simulations of CO 2 injection using Computer Modelling Group’s (CMG’s) GEM equation-of-state reservoir simulator. The scaler permeability multipliers allowed the fluid flow simulations to fully represent well test-measured injection well rates in the respective formations.

20 FOSSIL-FUELED POWER PLANTS↗

X-ray spectrometer experiment aboard the ISEE-C /Heliocentric/ spacecraft

This experiment is designed to provide continuous coverage of solar flare X-ray bursts and transient cosmic gamma-ray bursts. A proportional counter and a scintillation detector together cover the energy range from 5 to 228 keV with good sensitivity, large dynamic range, and high temporal resolution. This experiment provides data storage capability and good absolute timing so that in conjunction with similar experiments on other spacecraft, accurate source locations can be obtained for cosmic gamma-ray bursts.

Anderson, K. A.↗

Time dependent difference theory for sound propagation in axisymmetric ducts with plug flow

The time dependent governing acoustic-difference equations and boundary conditions are developed and solved for sound propagation in an axisymmetric (cylindrical) hard wall duct with a plug mean flow and spinning acoustic modes. The analysis begins with a harmonic sound source radiating into a quiescent duct. This explicit iteration method then calculates stepwise in real time to obtain the transient as well as the 'steady' state solutions of the acoustic field. The time dependent finite difference analysis has two advantages over the steady state finite difference and finite element techniques: (1)the elimination of large matrix storage requirements, and (2)shorter solution times under most conditions.

Baumeister, K. J.↗

Time dependent difference theory for sound propagation in axisymmetric ducts with plug flow

The time dependent governing/acoustic-difference equations and boundary conditions are developed and solved for sound propagation in an axisymmetric (cylindrical) hard wall duct with a plug mean flow and spinning acoustic modes. The analysis begins with a harmonic sound source radiating into a quiescent duct. This explicit iteration method then calculates stepwise in real time to obtain the transient as well as the 'steady' state solutions of the acoustic field. The time dependent finite difference analysis has two advantages over the steady state finite difference and finite element techniques: (1) the elimination of large matrix storage requirements, and (2) shorter solution times under most conditions.

Baumeister, K. J.↗