A mathematical analysis of the suppression of high frequency combustion instabilities by a sound absorbing liner.
Sound absorbing liner effect on harmonic oscillations at cylindrical boundary of rocket combustion chamber
SEARCH · Engineering Papers
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.
Sound absorbing liner effect on harmonic oscillations at cylindrical boundary of rocket combustion chamber
Inorganic CsPbIBr 2 perovskites have recently attracted enormous attention as a viable alternative material for optoelectronic applications due to their higher efficiency, thermal stability, suitable bandgap, and proper optical absorption. However, the CsPbIBr 2 perovskite films fabricated using a one-step deposition technique is usually comprised of small grain size with a large number of grain boundaries and compositional defects. In this work, silver iodide (AgI) will be incorporated as an additive into the CsPbIBr 2 perovskite precursor solution to prepare the unique perovskite CsI(PbBr 2 ) 1-x (AgI) x . The AgI additive in the precursor solution works as a nucleation promoter which will help the perovskite to grow and merge into a continuous film with reduced defects. With detailed characterizations, we found that incorporating AgI additive resulted in a uniform perovskite film with fewer grain boundaries, increased grain size, crystallinity, optical absorption while decreasing carrier recombination and trap density. Using the AgI in an optimum amount, we fabricated CsPbIBr 2 perovskite solar cells (PSCs) with a simple structure and achieved a power conversion efficiency (PCE) of 7.2% with a reduced hysteresis index. This work offers an alternative approach towards preparing high-quality CsPbIBr 2 perovskite films for solar cells with higher stability and other optoelectronic applications.
The U.S. Department of Energy (DOE) has selected vitrification for stabilizing legacy tank waste at the Hanford site, where radioactive waste from plutonium production was historically stored in underground tanks. This waste will be separated into low-activity waste (LAW) and high-level waste (HLW) fractions and processed at the Waste Treatment and Immobilization Plant (WTP). At WTP, glass melters are used for the vitrification of radioactive tank waste, transforming it into a stable borosilicate glass form for safe long-term storage. The melter vessel is constructed from highly durable and heat-resistant materials, where the vitrification process occurs. The main regions that are modeled are the melt pool, plenum, cold cap, riser/discharge chamber, and surrounding structure with insulation layers. Forced convection induced by air bubblers at the base of the melter ensure uniform temperature distribution and provide heat to the cold cap layer. The cold cap is a region of reacting batch feed that floats on top of the molten glass and is where the batch-to-glass reactions occur. Joule heating provided by electrodes mounted along the vertical walls of the melter and immersed directly in the glass, generates the necessary heat for the net endothermic conversion processes that occur in the cold cap. The high temperatures, radioactivity, and opaque nature of the glass prevent direct observation inside the melters. Therefore, computational models are essential for providing insight into factors that affect melter throughput. Thermocouples in the plenum provide operators with plenum temperature measurements. Operational adjustments include bubbling rate, voltage supplied to the electrodes, feed adjustments, and glass removal rate. Different computational fluid dynamics (CFD) models have been developed, each serving a specific purpose. There are CFD models of different scale melters, as well as models that capture the two-phase flow interfaces of rising bubbles in the molten glass or models with a simplified molten glass region so that the surrounding structure and plenum can be feasibly incorporated. Pilot-scale melter models have been developed to serve as validation of the methods employed in the simulation of the full-scale WTP melters. Models incorporating resolved bubbling are used to develop momentum source terms to implement into a single phase, multi-region, steady-state flow model that is being validated by measured process parameters such as glass production rate, voltage, input power, plenum temperatures, etc. The resolved bubbling model uses the multiphase volume of fluid approach to model the system with a high-resolution interface capturing scheme to maintain sharp interfaces between the molten glass and the air phase. The suite of CFD models is continually being improved to incorporate more realistic physics and achieve faster turnaround time. For example, an incremental controller is implemented to automatically adjust electrode voltage within the simulation to a molten glass set point temperature of 1150°C. Newer models feature improved meshes to ensure conformal meshes between regions and eliminate unnecessary mesh refinement in areas that are not of interest (such as boundary layers in offgas ports). Instead of explicitly modeling the structural, refractory, and insulation layers of the melter, a thermal resistance approach is used with published correlations used for boundary conditions. The development of robust and efficient CFD models will be instrumental in enabling the WTP to successfully fulfill its mission of safely stabilizing legacy nuclear waste.
We propose a spectral method for the 1D-1V Vlasov–Poisson system where the discretization in velocity space is based on asymmetrically-weighted Hermite functions, dynamically adapted via a scaling α and shifting u of the velocity variable. Specifically, at each time instant an adaptivity criterion selects new values of α and u based on the numerical solution of the discrete Vlasov–Poisson system obtained at that time step. Once the new values of the Hermite parameters α and u are fixed, the Hermite expansion is updated and the discrete system is further evolved for the next time step. The procedure is applied iteratively over the desired temporal interval. The key aspects of the adaptive algorithm are: the map between approximation spaces associated with different values of the Hermite parameters that preserves total mass, momentum and energy; and the adaptivity criterion to update α and u based on physics considerations relating the Hermite parameters to the average velocity and temperature of each plasma species. For the discretization of the spatial coordinate, we rely on Fourier functions and use the implicit midpoint rule for time stepping. The resulting numerical method possesses intrinsically the property of fluid-kinetic coupling, where the low-order terms of the expansion are akin to the fluid moments of a macroscopic description of the plasma, while kinetic physics is retained by adding more spectral terms. Moreover, the scheme features conservation of total mass, momentum and energy associated in the discrete, for periodic boundary conditions. A set of numerical experiments confirms that the adaptive method outperforms the non-adaptive one in terms of accuracy and stability of the numerical solution.
The rare earth neodymium arguably exhibits the most complex magnetic ordering and series of magnetic phase transitions of the elements. Here we report the results of small-angle neutron scattering (SANS) measurements as a function of temperature and applied magnetic field to study magnetic correlations on nanometer length scales in Nd. We support these results with additional neutron diffraction, thermodynamic and synchrotron x-ray measurements. The SANS measurements reveal the presence of modulation vectors characterizing the ordered spin configuration which exhibit changes in magnitude and direction that are phase dependent. Between 5.9 and 7.6 K an additional modulation vector is observed with magnitude Q = 0.12 Å -1 . This scattering appears to originate from order of the Nd layers which contain a center of inversion. Tracking this modulation vector as a function of magnetic field indicates a phase boundary at ≈1 T. The modulation vectors observed in this study establish the presence of nanometer length scale spin textures which are likely stabilized by frustrated Ruderman–Kittel–Kasuya–Yosida (RKKY) interactions.
Compared with their three-dimensional (3D) counterparts, low-dimensional metal halide perovskites (2D and quasi-2D; B 2 A n-1 M n X 3n+1 , such as B = R-NH 3 + , A = HC(NH 2 ) 2 + , Cs + ; M = Pb 2+ , Sn 2+ ; X = Cl - , Br - , I - ) with periodic inorganic-organic structures have shown promising stability and hysteresis-free electrical performance. However, their unique multiple-quantum-well structure limits the device efficiencies because of the grain boundaries and randomly oriented quantum wells in polycrystals. In single crystals, the carrier transport through the thickness direction is hindered by the layered insulating organic spacers. Furthermore, the strong quantum confinement from the organic spacers limits the generation and transport of free carriers. Also, lead-free metal halide perovskites have been developed but their device performance is limited by their low crystallinity and structural instability 11 . Here, in this paper, we report a low-dimensional metal halide perovskite BA 2 MA n-1 Sn n I 3n+1 (BA, butylammonium; MA, methylammonium; n = 1, 3, 5) superlattice by chemical epitaxy. The inorganic slabs are aligned vertical to the substrate and interconnected in a criss-cross 2D network parallel to the substrate, leading to efficient carrier transport in three dimensions. A lattice-mismatched substrate compresses the organic spacers, which weakens the quantum confinement. The performance of a superlattice solar cell has been certified under the quasi-steady state, showing a stable 12.36% photoelectric conversion efficiency. Moreover, an intraband exciton relaxation process may have yielded an unusually high open-circuit voltage (V OC ).
One of the largest uncertainties in the terrestrial carbon cycle is the timing and magnitude of soil organic carbon (SOC) response to climate and vegetation change. This uncertainty prevents models from adequately capturing SOC dynamics and challenges the assessment of management and climate change effects on soils. Reducing these uncertainties requires simultaneous investigation of factors controlling the amount (SOC abundance) and duration (SOC persistence) of stored C. We present a global synthesis of SOC and radiocarbon profiles (n Profile = 597) to assess the timescales of SOC storage. We use a combination of statistical and depth–resolved compartment models to explore key factors controlling the relationships between SOC abundance and persistence across pedo–climatic regions and with soil depth. This allows us to better understand (i) how SOC abundance and persistence covary across pedo–climatic regions and (ii) how the depth dependence of SOC dynamics relates to climatic and mineralogical controls on SOC abundance and persistence. We show that SOC abundance and persistence are differently related; the controls on these relationships differ substantially between major pedo–climatic regions and soil depth. For example, large amounts of persistent SOC can reflect climatic constraints on soils (e.g., in tundra/polar regions) or mineral absorption, reflected in slower decomposition and vertical transport rates. In contrast, lower SOC abundance can be found with lower SOC persistence (e.g., in highly weathered tropical soils) or higher SOC persistence (e.g., in drier and less productive regions). We relate variable patterns of SOC abundance and persistence to differences in the processes constraining plant C input, microbial decomposition, vertical C transport and mineral SOC stabilization potential. This process–oriented grouping of SOC abundance and persistence provides a valuable benchmark for global C models, highlighting that pedo–climatic boundary conditions are crucial for predicting the effects of climate change and soil management on future C abundance and persistence.
It has previously been inferred that small isolated flux tubes appearing in supergranule boundaries are compressed to 1500 gauss or more. This paper considers whether some dynamic condition within a flux tube exists which provides both stability and a 'mechanical advantage' so that a small force over a small period of time can accomplish the enormous compression from the weak-field to the strong-field state. It is found that the equipartition solutions to the hydromagnetic equations apparently may have the desired property of permitting an infinitesimal external pressure to convert a gentle flow of gas along a weak field into a very intense field through a succession of equipartition states. An illustrative example is presented, and field compression by convective forces is analyzed.
Graduate student engineering research in aeronautics at Old Dominion University is surveyed. Student participation was facilitated through a NASA sponsored university program which enabled the students to complete degrees. Research summaries are provided and plans for the termination of the grant program are outlined. Project topics include: Failure modes for mechanically fastened joints in composite materials; The dynamic stability of an earth orbiting satellite deploying hinged appendages; The analysis of the Losipescu shear test for composite materials; and the effect of boundary layer structure on wing tip vortex formation and decay.
A new hybrid-stress finite element algorithm, suitable for analyses of large, quasistatic, inelastic deformations, is presented. The algorithm is based upon a generalization of de Veubeke's (1972) complementary energy principle. The principal variables in the formulation are the nominal stress rate and spin, and the resulting finite element equations are discrete versions of the equations of compatibility and angular momentum balance. The algorithm produces true rates, time derivatives, as opposed to 'increments'. There results a boundary value problem (for stress rate and velocity) and an initial value problem (for total stress and deformation). A discussion of the numerical treatment of the boundary value problem is followed by a detailed examination of the numerical treatment of the initial value problem, covering the topics of efficiency, stability, and objectivity. The paper is closed with a set of examples, finite homogeneous deformation problems, which serve to bring out important aspects of the algorithm.
A new hybrid-stress finite element algorithm, suitable for analyses of large, quasistatic, inelastic deformations, is presented. The algorithm is base upon a generalization of de Veubeke's complementary energy principle. The principal variables in the formulation are the nominal stress rate and spin, and thg resulting finite element equations are discrete versions of the equations of compatibility and angular momentum balance. The algorithm produces true rates, time derivatives, as opposed to 'increments'. There results a complete separation of the boundary value problem (for stress rate and velocity) and the initial value problem (for total stress and deformation); hence, their numerical treatments are essentially independent. After a fairly comprehensive discussion of the numerical treatment of the boundary value problem, we launch into a detailed examination of the numerical treatment of the initial value problem, covering the topics of efficiency, stability and objectivity. The paper is closed with a set of examples, finite homogeneous deformation problems, which serve to bring out important aspects of the algorithm.
The 3-D boundary layer on a swept wing can support different types of hydrodynamic instability. Attention is focused on the so-called spanwise contamination problem, which occurs when the attachment line boundary layer on the leading edge becomes unstable to Tollmien-Schlichting waves. In order to gain insight into the interactions important in that problem, a simplified basic state is considered. This simplified flow corresponds to the swept attachment line boundary layer on an infinite flat plate. The basic flow here is an exact solution of the Navier-Stokes equations and its stability to 2-D waves propagating along the attachment can be considered exactly at finite Reynolds number. This has been done in the linear and weakly nonlinear regimes. The corresponding problem is studied for oblique waves and their interaction with 2-D waves is investigated. In fact, oblique modes cannot be described exactly at finite Reynolds number so it is necessary to make a high Reynolds number approximation and use triple deck theory. It is shown that there are two types of oblique wave which, if excited, cause the destabilization of the 2-D mode and the breakdown of the disturbed flow at a finite distance from the leading edge. First, a low frequency mode related to the viscous stationary crossflow mode is a possible cause of breakdown. Second, a class of oblique wave with frequency comparable with that of the 2-D mode is another cause of breakdown. It is shown that the relative importance of the modes depends on the distance from the attachment line.
The three-dimensional boundary layer on a swept wing can support different types of hydrodynamic instability. Here attention is focused on the so-called 'spanwise instability' problem which occurs when the attachment-line boundary layer on the leading edge becomes unstable to Tollmien-Schlichting waves. In order to gain insight into the interactions that are important in that problem a simplified basic state is considered. This simplified flow corresponds to the swept attachemnt-line boundary layer on an infinite flat plate. The basic flow here is an exact solution of the Navier-Stokes equations and its stability to two-dimensional waves propagating along the attachment line can be considered exactly at finite Reynolds number. The corresponding problem is studied for oblique waves, and their interaction with two-dimensional waves is investigated. It is shown that there are two types of oblique wave which, if excited, cause the destabilization of the two-dimensional mode and the breakdown of the disturbed flow at a finite distance from the leading edge.
The results are presented of a study to identify potential real time remote computational applications to support monitoring HRV flight test experiments along with definitions of preliminary requirements. A major expansion of the support capability available at Ames-Dryden was considered. The focus is on the use of extensive computation and data bases together with real time flight data to generate and present high level information to those monitoring the flight. Six examples were considered: (1) boundary layer transition location; (2) shock wave position estimation; (3) performance estimation; (4) surface temperature estimation; (5) critical structural stress estimation; and (6) stability estimation.
The periodic streamwise modulation of the supersonic and hypersonic boundary layers by a two dimensional first mode or second mode wave makes the resulting base flow susceptible to a broadband spanwise-periodic three dimensional type of instability. The principal parametric resonance of this instability (subharmonic) was analyzed using Floquet theory. The effect of Mach number and the effectiveness of wall cooling or wall suction in controlling the onset, the growth rate, and the vortical nature of the subharmonic secondary instability are assessed for both a first mode and a second mode primary wave. Results indicate that the secondary subharmonic instability of the insulated wall boundary layer is weakened as Mach number increases. Cooling of the wall destabilizes the secondary subharmonic of a second mode primary wave, but stabilizes it when the primary wave is a first mode. Suction stabilizes the secondary subharmonic at all Mach numbers.
A flexible fuselage capability has been developed and implemented within version 1.2 of the CAP-TSD code. The capability required adding time dependent terms to the fuselage surface boundary conditions and the fuselage surface pressure coefficient. The new capability will allow modeling the effect of a flexible fuselage on the aeroelastic stability of complex configurations. To assess the flexible fuselage capability several steady and unsteady calculations have been performed for slender fuselages with circular cross-sections. Steady surface pressures are compared with experiment at transonic flight conditions. Unsteady cross-sectional lift is compared with other analytical results at a low subsonic speed and a transonic case has been computed. The comparisons demonstrate the accuracy of the flexible fuselage modifications.
It is known that the incompressible, infinite swept attachment line flow is unstable to streamwise disturbances that originate in the boundary layer when the cross-flow exceeds a critical magnitude. Furthermore, a small degree of suction at the surface has a significant stabilizing influence while a small degree of blowing has a considerable destabilizing influence. This paper investigates the stabilizing and destabilizing effects of, respectively, cooling or heating the plate and the competing or enhancing effects of suction or blowing. A nonorthogonal flow with respect to the attachment line is also considered by adding a component of shear to the mean flow.
A new procedure seeks to combine the thin-layer Navier-Stokes solver LAURA with the parabolized Navier-Stokes solver UPS for the aerothermodynamic solution of chemically-reacting air flowfields. The interface protocol is presented and the method is applied to two slender, blunted shapes. Both axisymmetric and three dimensional solutions are included with surface pressure and heat transfer comparisons between the present method and previously published results. The case of Mach 25 flow over an axisymmetric six degree sphere-cone with a noncatalytic wall is considered to 100 nose radii. A stability bound on the marching step size was observed with this case and is attributed to chemistry effects resulting from the noncatalytic wall boundary condition. A second case with Mach 28 flow over a sphere-cone-cylinder-flare configuration is computed at both two and five degree angles of attack with a fully-catalytic wall. Surface pressures are seen to be within five percent with the present method compared to the baseline LAURA solution and heat transfers are within 10 percent. The effect of grid resolution is investigated and the nonequilibrium results are compared with a perfect gas solution, showing that while the surface pressure is relatively unchanged by the inclusion of reacting chemistry the nonequilibrium heating is 25 percent higher. The procedure demonstrates significant, order of magnitude reductions in solution time and required memory for the three dimensional case over an all thin-layer Navier-Stokes solution.