Search NASA⌕ Search

SEARCH · Search NASA

Results for “splitting 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 1,351 records · Page 75

Cover crops and poultry litter impact on soil structural stability in dryland soybean production in southeastern United States

Abstract This study explored the efficacy of soil aggregate indices in quantifying soil structural development, utilizing 5‐year field experiment data from the Southeastern United States. The experiment utilized a split‐plot design with cover crops (native vegetation as control, cereal rye (Secale cerealeL.), winter wheat (Triticum aestivum), hairy vetch (Vicia villosa), and mustard (Brassica rapa) plus cereal rye as the main factor and fertilizer source (no fertilizer as control, inorganic fertilizer with phosphorus, potassium, and elemental sulfur, and poultry litter) as the secondary factor. Aggregate size fractions were determined using the wet‐sieving method, and aggregate stability index (ASI), mean weight diameter (MWD), geometric mean diameter (GMD), and fractal dimension (FD) were calculated to assess soil structural stability. Main effects results indicated that cereal rye (55.11%) and poultry litter (50.97%) exhibited the highest ASI values. The highest MWD, GMD, and FD were observed under mustard plus cereal rye (1.187 mm), cereal rye (0.462 mm), and hairy vetch (2.573), respectively. Principal component analysis revealed that cover crops significantly improved soil aggregate structure and stability, overcoming limitations of sole fertilization practices. Regression analysis suggested that ASI, MWD, and GWD positively correlated with soil organic carbon, whereas FD negatively correlated with MWD, GMD, and ASI. Principal component analysis exhibited that FD decreased with increasing soil organic carbon, ASI, MWD, and GMD, demonstrating that lower FD values indicate enhanced soil aggregation and structure. Assessed indices, FD included, effectively gauged soil structural stability. These metrics should be prioritized in managerial decisions to support soil productivity and health in agricultural systems.

Agriculture↗

Formation of monomeric Sn(ii) and Sn(iv) perfluoropinacolate complexes and their characterization by 119 Sn Mössbauer and 119 Sn NMR spectroscopies

The synthesis and characterization of a series of Sn(II) and Sn(IV) complexes supported by the highly electron-withdrawing dianionic perfluoropinacolate (pinF) ligand are reported herein. Three analogs of [Sn IV (pin F ) 3 ] 2- with NEt 3 H + (1), K + (2), and {K(18C6)} + (3) counter cations and two analogs of [SnII(pinF)2]2- with K+ (4) and {K(15C5) 2 } + (5) counter cations were prepared and characterized by standard analytical methods, single-crystal X-ray diffraction, and 119 Sn Mössbauer and NMR spectroscopies. The six-coordinate SnIV(pinF) complexes display 119 Sn NMR resonances and 119 Sn Mössbauer spectra similar to SnO 2 (cassiterite). In contrast, the four-coordinate SnII(pinF) complexes, featuring a stereochemically-active lone pair, possess low 119 Sn NMR chemical shifts and relatively high quadrupolar splitting. Furthermore, the Sn(II) complexes are unreactive towards both Lewis bases (pyridine, NEt 3 ) and acids (BX 3 , Et 3 NH + ). Calculations confirm that the Sn(II) lone pair is localized within the 5s orbital and reveal that the Sn 5p x LUMO is energetically inaccessible, which effectively abates reactivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Narrow band laser source for efficient 3 He polarization

Polarized 3 He is widely in polarized neutron scattering as neutron polarizers or analyzers. Increasing 3 He polarization can improve both neutron polarization and transmission which are important for any polarized neutron experiments. Here at the US, 3 Hee is usually polarized using spin-exchange optical pumping (SEOP), where 3 He gains polarization via spin-exchange collisions with laser pumped alkali metals, typically rubidium (Rb) and potassium (K). Over the last two decades, 3 He polarization has improved from around 50% to around 80% due to the application of spectral narrowed diode lasers which narrows the laser linewidth from 2 nm to between 0.2 nm and 0.3 nm. However, the highest 3 He polarization achieved so far is limited to about 85%. The mechanism behind this is still not well understood, but it is suspected that a laser with an ultra-narrow linewidth and lower broadband background may further improve the 3 He polarization. In addition, if we define the spin-exchange efficiency as the ratio of the rate at which the polarized 3 He is produced to the rate at which the alkali atom polarization is lost through various collisions, K has much higher spin-exchange efficiency than Rb does. In principle, direct pumping of K should improve both the pumping efficiency and the production rate of polarized 3 He. However, because of the small fine structure splitting in K (D1 770 nm, D2 766.7 nm), the laser used to pump K needs to have a ultra-narrow spectrum tuned at 770 nm and minimal broadband background as the pumping of the K D2 transition will reduce the pumping efficiency. Such a laser is not readily available. Therefore, the current practice of polarizing 3 He is through a hybrid SEOP method in which the laser tuned at Rb D1 transition (795 nm) polarizes Rb first and then Rb transfers its polarization to K. This hybrid method takes advantages of the more accessible 795 nm laser and gains some pumping efficiency. But to realize the full potential in SEOP, direct pumping of K would be necessary. The new Optigrate laser provides a tool to investigate the impact on the pumping efficiency and achievable 3 He polarization under the condition of ultra-narrow spectrum and very low background.

3He polarizing filters↗

IM3/HyperFACETS Thermodynamic Global Warming (TGW) Simulation Datasets

Publication For a thorough description of the methods, see the peer-reviewed paper: Jones, A.D., Rastogi, D., Vahmani, P. et al. Continental United States climate projections based on thermodynamic modification of historical weather. Sci Data 10, 664 (2023). https://doi.org/10.1038/s41597-023-02485-5 Overview The IM3 / HyperFACETS climate simulations provide 40-year historical (1980-2019) as well as four 80-year future simulations (2020-2099) over the U.S. The future simulations are split into near (2020-2059) and far future (2060-2099) segments. The future scenarios span a range of plausible changes in future climate (both Global Circulation Model (GCM) and Representative Concentration Pathways/Shared Socioeconomic Pathway (RCP/SSP) dimensions). The simulations provide climate variables with high spatiotemporal resolution (25 hourly variables and 207 3-hourly variables at 12 km2). The datasets are generated using dynamical downscaling with the WRF (Weather Research and Forecasting) model (version 4.2.1) and therefore preserve physical consistency across variables. WRF is a state-of-the-art, fully compressible, non-hydrostatic, mesoscale numerical weather prediction model. WRF is coupled with an urban canopy model (UCM), which resolves urban surfaces. The future scenarios were developed using a thermodynamic global warming approach where past events are replayed under a range of future warming conditions. These scenarios therefore provide a perspective on potential increases in extreme event intensity, geographic scope, and duration, with previously non-extreme conditions potentially crossing new thresholds to be considered extreme by today's standards. This approach is not intended to estimate future changes in extreme event frequency that might result from changes in large-scale atmospheric dynamics. This dataset has NOT been bias corrected. A bias corrected version of selected variables is under development and will be released here when available. Scenarios Files Data for each scenario is provided in weekly NetCDF files. 25 variables are available at hourly resolution, and 207 variables are available at three-hourly resolution. Spatial resolution is 12km and spans the conterminous United States (CONUS), including some areas of Canada and Mexico, resulting in a grid of 424 by 299 cells. The spatial projection is a Lambert Conformal Conic with the following proj-string: "+proj=lcc +lat_0=40.0000076293945 +lon_0=-97 +lat_1=30 +lat_2=45 +x_0=0 +y_0=0 +R=6370000 +units=m +no_defs". The available scenarios and simulation periods are listed below: historical | 1980 - 2019 rcp45cooler | 2020 - 2059 rcp45cooler | 2060 - 2099 rcp45hotter | 2020 - 2059 rcp45hotter | 2060 - 2099 rcp85cooler | 2020 - 2059 rcp85cooler | 2060 - 2099 rcp85hotter | 2020 - 2059 rcp85hotter | 2060 - 2099 * The first year (1979, 2019, and 2059) of data within each scenario represents a model warmup period and should not be used. These are located in the `spinup_files` directory. Historical year 2020 is considered an extra year of data beyond the simulation period and can be found in the `additional_files` directory. For information on specific variables and a more in-depth discussion of methodology, please refer to the data landing page at https://tgw-data.msdlive.org. Delta Warming Files The global and CONUS warming deltas for each scenario are provided in degrees Celsius annually and monthly. Restart Files Yearly restart files are provided for each scenario which can be used to restart the WRF model at a particular point in time. Spinup Files The first year of data within each simulation period represents a model warmup period and should not be used. The files are provided here for the sake of reproducibility. Additional Files Additional years of data are provided as an extension of the historic simulation.

Jones, Andrew D.↗

Interlaminar fracture characterization: A current review

Interlaminar fracture characterization has been investigated for several years. Only now is it well enough understood for standardization organizations to attempt to write standard test methods. A review of the current philosophies in characterizing interlaminar fracture is presented. First, the mode 1, double cantilever beam test (DCB) for measuring G(sub 1c) and the end notched flexure (ENF) test for measuring G(sub 2c) are discussed. These tests have undergone the most extensive research throughout the years and are furthest towards standardization. In addition, the mode 2, end loaded split (ELS) specimen is discussed. Mixed mode fracture is also discussed and the recently developed mixed mode bending (MMB) test is detailed. Then, tests for evaluating mode 3 fracture toughness, including the split cantilever beam (SCB), are reviewed. Finally, the work done on interlaminar fracture characterization in fatigue is reviewed.

Martin, Roderick H.↗

High-frequency ECG

The standard ECG is by convention limited to 0.05-150 Hz, but higher frequencies are also present in the ECG signal. With high-resolution technology, it is possible to record and analyze these higher frequencies. The highest amplitudes of the high-frequency components are found within the QRS complex. In past years, the term "high frequency", "high fidelity", and "wideband electrocardiography" have been used by several investigators to refer to the process of recording ECGs with an extended bandwidth of up to 1000 Hz. Several investigators have tried to analyze HF-QRS with the hope that additional features seen in the QRS complex would provide information enhancing the diagnostic value of the ECG. The development of computerized ECG-recording devices that made it possible to record ECG signals with high resolution in both time and amplitude, as well as better possibilities to store and process the signals digitally, offered new methods for analysis. Different techniques to extract the HF-QRS have been described. Several bandwidths and filter types have been applied for the extraction as well as different signal-averaging techniques for noise reduction. There is no standard method for acquiring and quantifying HF-QRS. The physiological mechanisms underlying HF-QRS are still not fully understood. One theory is that HF-QRS are related to the conduction velocity and the fragmentation of the depolarization wave in the myocardium. In a three-dimensional model of the ventricles with a fractal conduction system it was shown that high numbers of splitting branches are associated with HF-QRS. In this experiment, it was also shown that the changes seen in HF-QRS in patients with myocardial ischemia might be due to the slowing of the conduction velocity in the region of ischemia. This mechanism has been tested by Watanabe et al by infusing sodium channel blockers into the left anterior descending artery in dogs. In their study, 60 unipolar ECGs were recorded from the entire ventricular surface and were signal-averaged and filtered in the 30-250 Hz frequency range. The results showed that the decrease noted in the HF-QRS correlated linearly with the local conduction delay. The results suggest that HF-QRS is a potent indicator of disturbed local conduction. An alternative theory is that HF-QRS reflect the shape of the original electrocardiographic signal. Bennhagen et al showed that root mean square (RMS) voltage values of the depolarization signal correlate poorly with the signal amplitude but highly with the first and second derivatives, i.e. the velocity and the acceleration of the signal. It has also been suggested that the autonomic nervous system affects HF-QRS. For example, sitting up causes significant changes in HF-QRS in some leads compared to the supine position [Douglas et al., 2006]. Unpublished results indicate that familial dysautonomic patients (both vagal and sympathetic degeneration) have very little Reduced Amplitude Zones (RAZ) formation . Athletic individuals, especially elite athletes, who have vagally-mediated changes on the conventional ECG (i.e. early repolarization, bradycardia) have increased RAZ formation. Further electrophysiological studies are needed, however, to better understand the underlying mechanisms of HF-QRS. Several investigators have studied HF-QRS in different cardiac conditions, including acute myocardial ischemia and myocardial infarction (MI). However, in order for clinicians to confidently use HF-QRS as an adjunct to standard ECG, more knowledge about the characteristics of HF-QRS is needed.

Tragardh, Elin↗

Metal-facilitated, sustainable nitroarene hydrogenation under ambient conditions

Hydrogenation is a critical reaction in the chemical industry, yielding a range of important compounds such as fine chemicals, pharmachemicals and agrochemicals. However, conventional hydrogenation typically requires pressurized hydrogen, high temperatures and involves noble metal catalysts. Here, we proposed a two-step hydrogenation process, utilizing water as the hydrogen source for the industrially important reduction of nitroarenes to anilines. A metal or reduced metal oxide, which can be obtained from solar thermal or electrochemical reduction, acts as the active site for nitrobenzene adsorption, H 2 O dissociation and in-situ hydrogen generation. Among the 15 metal and reduced metal oxides investigated, Zn and Sn emerged as highly efficient catalysts for the reduction of a broad range of organic nitro compounds under mild conditions, with H 2 utilization efficiency 1-2 orders of magnitude above the state-of-the-art. The presented protocol provides extra dimensions for designing and optimizing conventional hydrogenation process with an alternative pathway. The reactive hydrogen atoms generated in-situ effectively overcome the barriers associated with hydrogen gas dissolution and its subsequent dissociation on the catalyst surface, thereby greatly enhancing the overall effectiveness for the hydrogenation reaction. This research potentially establishes a sustainable, generally applicable alternative to conventional hydrogenation methods, simultaneously presenting a viable solution for renewable energy storage.

25 ENERGY STORAGE↗

Atomic and electronic structure of the CdTe(111)B–(2√3 × 4) orthogonal surface

The atomic and electronic structure of four variants of Te-terminated CdTe(111)B–(2√3 × 4) orthogonal polar surface (ideal, relaxed, reconstructed, and reconstructed with subsequent relaxation) are calculated ab initio for the first time. The surface is modeled by a film composed of 12 atomic layers with a vacuum gap of ~16 Å in the layered superlattice approximation. To close Cd dangling bonds on the opposite side of the film, 24 fictitious hydrogen atoms with a charge of 1.5 electrons each are added. Ab initio calculations are performed using the Quantum Espresso program based on density functional theory. It is demonstrated that relaxation leads to splitting of the four upper layers. The band energy structures and total and layer-by-layer densities of electronic states for the four surface variants are calculated and analyzed.

36 MATERIALS SCIENCE↗

A new relaxed splitting preconditioner for the generalized saddle point problems from the incompressible Navier–Stokes equations

Based on the modified relaxed splitting (MRS) preconditioner proposed by Fan and Zhu (Appl Math Lett 55:18–26, 2016), a new relaxed splitting preconditioner is proposed for the generalized saddle point problems arising from the incompressible Navier–Stokes equations. The eigenvalue distribution and an upper bound of the degree of the minimal polynomial of the preconditioned matrix are analyzed. Numerical results show that the proposed preconditioner is superior to the MRS preconditioner when they are used as preconditioners to accelerate the convergence rate of the Krylov subspace methods, such as GMRES.

97 MATHEMATICS AND COMPUTING↗

Electrochemical pH regulation in droplet microfluidics

Here we report a method for electrochemical pH regulation in microdroplets generated in a microfluidic device. The key finding is that controlled quantities of reagents can be generated electrochemically in moving microdroplets confined within a microfluidic channel. Additionally, products generated at the anode and cathode can be isolated within descendant microdroplets. Specifically, ~5 nL water-in-oil microdroplets are produced at a T-junction and then later split into two descendant droplets. During splitting, floor-patterned microelectrodes drive water electrolysis within the aqueous microdroplets to produce H + and OH - . This results in a change in the pHs of the descendant droplets. The droplet pH can be regulated over a range of 5.9 to 7.7 by injecting controlled amounts of charge into the droplets. When the injected charge is between -6.3 and 54.5 nC nL -1 , the measured pH of the resulting droplets is within ±0.1 pH units of that predicted based on the magnitude of the injected charge. This technique can likely be adapted to electrogeneration of other reagents within microdroplets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Space shuttle data handling and communications considerations.

Operational and development flight instrumentation, data handling subsystems and communication requirements of the space shuttle orbiter are discussed. Emphasis is made on data gathering methods, crew display data, computer processing, recording, and telemetry by means of a digital data bus. Also considered are overall communication conceptual system aspects and design features allowing a proper specification of telemetry encoders and instrumentation recorders. An adaptive bit rate concept is proposed to handle the telemetry bit rates which vary with the amount of operational and experimental data to be transmitted. A split-phase encoding technique is proposed for telemetry to cope with the excessive bit jitter and low bit transition density which may affect television performance.

Stoker, C. J.↗

An efficient nonlinear relaxation technique for the three-dimensional, Reynolds-averaged Navier-Stokes equations

An efficient implicit method for the computation of steady, three-dimensional, compressible Navier-Stokes flowfields is presented. A nonlinear iteration strategy based on planar Gauss-Seidel sweeps is used to drive the solution toward a steady state, with approximate factorization errors within a crossflow plane reduced by the application of a quasi-Newton technique. A hybrid discretization approach is employed, with flux-vector splitting utilized in the streamwise direction and central differences with artificial dissipation used for the transverse fluxes. Convergence histories and comparisons with experimental data are presented for several 3-D shock-boundary layer interactions. Both laminar and turbulent cases are considered, with turbulent closure provided by a modification of the Baldwin-Barth one-equation model. For the problems considered (175,000-325,000 mesh points), the algorithm provides steady-state convergence in 900-2000 CPU seconds on a single processor of a Cray Y-MP.

Edwards, Jack R.↗

Satellite Optical Remote Sensing of Clouds and Aerosols: From Particle Single-Scattering and Gaseous Absorption Through Radiative Transfer to Retrieval Products

Clouds and aerosols are fundamental regulators of Earth’s radiation budget and climate system, influencing both solar and terrestrial radiation through scattering, absorption, and emission processes. Accurate characterization of their physical and radiative properties from space requires a rigorous understanding of particle single-scattering, gaseous absorption, and radiative transfer in the atmosphere, as well as reliable inversion methods. This review synthesizes the physical foundations and algorithmic implementations of satellite-based passive optical remote sensing of clouds and aerosols, spanning the ultraviolet to thermal infrared spectral range. Beginning with electromagnetic scattering theory and state-of-the-art methods for computing single-scattering by nonspherical particles and computationally efficient methods for accounting for atmospheric absorption, we discuss the radiative transfer framework underpinning cloud and aerosol retrievals. The connection between single-scattering and multiple-scattering is rigorously formulated. We then summarize operational and research-grade retrieval techniques, including cloud masking and thermodynamic phase determination, CO₂ slicing for cloud-top pressure, the Nakajima-King shortwave bi-spectral, and infrared split-window approaches for cloud optical thickness and effective particle size, inversion algorithms for determining aerosol properties from multi-spectral and/or multi-angle radiometric and polarimetric measurements, and active-passive sensing synergy. Examples of the global cloud and aerosol climatologies are illustrated using observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Multi-angle Imaging SpectroRadiometer (MISR). Furthermore, the unique strengths of active remote sensing techniques based on spaceborne lidar observations are briefly elaborated in the context of studying ice clouds composed of randomly and horizontally oriented ice crystals, which is a significant challenge for conventional passive remote sensing techniques. By connecting physical theory to practical retrievals, this review highlights both the maturity of current methodologies and the remaining challenges in reducing uncertainties in particle morphology, vertical structure, absorption, and aerosol-cloud interactions. Furthermore, the impact of artificial intelligence (AI) on atmospheric remote sensing is briefly addressed.

Aerosols↗

Universal Reversible Hydrogen Potential for Electrocatalytic Ammonia Splitting Reactions in Nonaqueous Solvents from Unified pH Measurements

In this work, we introduce a new approach of using differential potentiometric measurements in four nonaqueous solvents─MeCN, THF, DMF, and PC─to determine the universal pH abs H 2 O values aligned to the aqueous pH scale for dilute NH 4 + /NH 3 solutions. Knowledge of the pH abs H 2 O values allows simple determination of the reversible hydrogen potential in any given solvent relative to the aqueous standard hydrogen electrode (SHE) and, most importantly, ensures comparability across different solvents. As an independent method, Open Circuit Potenial measurements were carried out in the same solvents titrated with NH 4 + /NH 3 to obtain alternative values for the reversible hydrogen potential in these solvents. The close agreement of these two methods, as well as calculated potentials from literature values when available, substantiates the new, simpler, and more robust approach to determine the reversible hydrogen potential introduced here. We further use the reversible hydrogen potential values established here to report, for the first time, the overpotential for ammonia oxidation as a function of solvent, with a recently discovered ruthenium catalyst.

ammonia↗

Improved Time-Stepping Methods in Global to Regional Ocean Modeling (Annual Status Report)

Time stepping algorithms are an important part of ocean models, and strongly influence both the accuracy of solution and performance. There have been a number of projects investigating various improvements for ocean time-stepping schemes in the Model for Prediction Across Scales-Ocean (MPAS-Ocean), a component of the DOE Energy Exascale Earth System Model. Ocean dynamics include fast surface gravity waves, which are two-dimensional, and slower internal waves, which are three-dimensional, so ocean models use a split time-stepping scheme that separates these barotropic and baroclinic modes for efficiency. MPAS-Ocean runs on variable-resolution horizontal meshes, and must scale to tens of thousands of cores and millions of horizontal gridcells. Ocean models require time stepping algorithms that are customized to these needs, and which are tuned for performance on various resolutions and architectures.

58 GEOSCIENCES↗

Cosmological constraints from the redshift dependence of the Alcock-Paczynski effect: application to the SDSS-III BOSS DR12 galaxies

We apply the methodology developed in Li et al. to BOSS DR12 galaxies and derive cosmological constraints from the redshift dependence of the Alcock–Paczynski (AP) effect. The apparent anisotropy in the distribution of observed galaxies arise from two main sources, the redshift-space distortion (RSD) effect due to the galaxy peculiar velocities, and the geometric distortion when incorrect cosmological models are assumed for transforming redshift to comoving distance, known as the AP effect. Anisotropies produced by the RSD effect are, although large, maintaining a nearly uniform magnitude over a large range of redshift, while the degree of anisotropies from the AP effect varies with redshift by a much larger magnitude. We split the DR12 galaxies into six redshift bins, measure the two-point correlation function in each bin, and assess the redshift evolution of anisotropies. We obtain constraints of Ω{sub m}=0.290±0.053,w=−1.07±0.15, which are comparable with the current constraints from other cosmological probes such as SNe Ia, cosmic microwave background, and baryon acoustic oscillation (BAO). Combining these cosmological probes with our method yield tight constraints of Ω{sub m}=0.301±0.006,w=−1.054±0.025. Our method is complementary to the other large-scale structure (LSS) probes like BAO and topology. We expect this technique will play an important role in deriving cosmological constraints from LSS surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Numerical modeling for dilute and dense sprays

We have successfully implemented a numerical model for spray-combustion calculations. In this model, the governing gas-phase equations in Eulerian coordinate are solved by a time-marching multiple pressure correction procedure based on the operator-splitting technique. The droplet-phase equations in Lagrangian coordinate are solved by a stochastic discrete particle technique. In order to simplify the calculation procedure for the circulating droplets, the effective conductivity model is utilized. The k-epsilon models are utilized to characterize the time and length scales of the gas phase in conjunction with turbulent modulation by droplets and droplet dispersion by turbulence. This method entails random sampling of instantaneous gas flow properties and the stochastic process requires a large number of computational parcels to produce the satisfactory dispersion distributions even for rather dilute sprays. Two major improvements in spray combustion modelings were made. Firstly, we have developed a probability density function approach in multidimensional space to represent a specific computational particle. Secondly, we incorporate the Taylor Analogy Breakup (TAB) model for handling the dense spray effects. This breakup model is based on the reasonable assumption that atomization and drop breakup are indistinguishable processes within a dense spray near the nozzle exit. Accordingly, atomization is prescribed by injecting drops which have a characteristic size equal to the nozzle exit diameter. Example problems include the nearly homogeneous and inhomogeneous turbulent particle dispersion, and the non-evaporating, evaporating, and burning dense sprays. Comparison with experimental data will be discussed in detail.

Chen, C. P.↗

Modelling Sonic Boom Propagation through Planetary Boundary Layer Turbulence near the Lateral Extent of the Carpet

Nonlinear acoustic propagation of sonic booms in the atmospheric boundary layer is considered in the context of a one-way solution to a third order wave equation. A split-step approach is utilized to compute different physical effects efficiently. Unlike previous approaches, the diffraction effects are computed exactly in the forward direction with no restriction on the angle of propagation. This results in more accurate modelling of sonic boom near the carpet edge. Heterogenous flow effects are incorporated with a wide-angle parabolic approximation. Nonlinear propagation is computed with a Burgers-Hayes method. Turbulence in the medium is constructed by the method of random Fourier modes. The turbulence spectra are constructed using an altitude dependent 3-dimensional von Karman spectrum. The turbulence field is considered frozen, as the eddy turnover time in the atmospheric boundary layer is generally much larger than the propagation time of the acoustic wave. Comparisons with benchmark predictions from the PCBoom software are conducted. The benchmark cases provide insight into the accuracy of the current prediction code for non-turbulent atmospheres. Future work is discussed regarding the prediction of shaped booms, boundary conditions, and predictions near the lateral extent of the boom carpet where previous approaches fail to capture the diffraction effects.

sonic boom↗