Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Second moment method”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Impact of point spread function higher moments error on weak gravitational lensing

ABSTRACT Weak gravitational lensing is one of the most powerful tools for cosmology, while subject to challenges in quantifying subtle systematic biases. The point spread function (PSF) can cause biases in weak lensing shear inference when the PSF model does not match the true PSF that is convolved with the galaxy light profile. Although the effect of PSF size and shape errors – i.e. errors in second moments – is well studied, weak lensing systematics associated with errors in higher moments of the PSF model require further investigation. The goal of our study is to estimate their potential impact for LSST weak lensing analysis. We go beyond second moments of the PSF by using image simulations to relate multiplicative bias in shear to errors in the higher moments of the PSF model. We find that the current level of errors in higher moments of the PSF model in data from the Hyper Suprime-Cam survey can induce a ∼0.05 per cent shear bias, making this effect unimportant for ongoing surveys but relevant at the precision of upcoming surveys such as LSST.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2006

Exploratory Studies in Generalized Predictive Control for Active Gust Load Alleviation

The results of numerical simulations aimed at assessing the efficacy of Generalized Predictive Control (GPC) for active gust load alleviation using trailing- and leading-edge control surfaces are presented. The equations underlying the method are presented and discussed, including system identification, calculation of control law matrices, and calculation of commands applied to the control effectors. Both embedded and explicit feedforward paths for inclusion of disturbance effects are addressed. Results from two types of simulations are shown. The first used a 3-DOF math model of a mass-spring-dashpot system subject to user-defined external disturbances. The second used open-loop data from a wind-tunnel test in which a wing model was excited by sinusoidal vertical gusts; closed-loop behavior was simulated in post-test calculations. Results obtained from these simulations have been decidedly positive. In particular, results of closed-loop simulations for the wing model showed reductions in root moments by factors as high as 1000, depending on whether the excitation is from a constant- or variable-frequency gust and on the direction of the response.

Kvaternik, Raymond G.↗
NASA NTRSDate not supplied

Distribution of SLS Integrated Load Uncertainty to Surface Pressures and Sectional Loads

Aerodynamic loads that are important to launch vehicle programs such as NASA’s Space Launch System (SLS) include both integrated loads such as the force & moment on the entire vehicle and distributed loads. In this work two examples of distributed loads are considered: the pressure field on the surface of the vehicle and sectional loads, which are one-dimensional distributions along the axis of the launch vehicle. In some modern flight programs, the integrated loads, such as lift and drag, used to design the guidance and control laws for the vehicle come from wind tunnel testing, while distributed loads are produced using Computational Fluid Dynamics(CFD). The first task that is addressed in this paper, then, is to provide a formal method to adjust the distributed loads so that integrating them matches the prescribed integrated load. In addition, the integrated loads in a launch vehicle typically include an uncertainty estimate. The second task is to distribute this prescribed integrated uncertainty to each point in a distributed load. Both tasks are addressed using the same technique, which is to create distributed load profiles that isolate adjustments to one integrated load while leaving the others unaffected. These adjustments are informed by Proper Orthogonal Decomposition (POD) of the entire CFD-based distributed load database. Once applied, the adjusted distributed loads can be used to evaluate any scalar quantity of interest that might be needed by downstream users such as structural analysis or trajectory modelers.

SLS↗
NASA NTRS2017

Nearfield Summary and Statistical Analysis of the Second AIAA Sonic Boom Prediction Workshop

A summary is provided for the Second AIAA Sonic Boom Workshop held 8-9 January 2017 in conjunction with AIAA SciTech 2017. The workshop used three required models of increasing complexity: an axisymmetric body, a wing body, and a complete configuration with flow-through nacelle. An optional complete configuration with propulsion boundary conditions is also provided. These models are designed with similar nearfield signatures to isolate geometry and shock/expansion interaction effects. Eleven international participant groups submitted nearfield signatures with forces, pitching moment, and iterative convergence norms. Statistics and grid convergence of these nearfield signatures are presented. These submissions are propagated to the ground, and noise levels are computed. This allows the grid convergence and the statistical distribution of a noise level to be computed. While progress is documented since the first workshop, improvement to the analysis methods for a possible subsequent workshop are provided. The complete configuration with flow-through nacelle showed the most dramatic improvement between the two workshops. The current workshop cases are more relevant to vehicles with lower loudness and have the potential for lower annoyance than the first workshop cases. The models for this workshop with quieter ground noise levels than the first workshop exposed weaknesses in analysis, particularly in convective discretization.

Park, Michael A.↗
NASA NTRS1993

Magnetopause shape as a bivariate function of interplanetary magnetic field B(sub z) and solar wind dynamic pressure

We present a new method for determining the shape of the magnetopause as a bivariate function of the hourly averaged solar wind dynamic pressure (p) and the north-south component of the interplanetary magnetic field (IMF) B(sub z). We represent the magnetopause (for X(sub GSE) greater than -40 R(sub E)) as an ellipsoid of revolution in solar-wind-aberrated coordinates and express the (p, B(sub z)) dependence of each of the three ellipsoid parameters as a second-order (6-term) bivariate expansion in Inp and B(sub z). We define 12 overlapping bins in a normalized dimensionless (p, B(sub z)) `control space' and fit an ellipsoid to those magnetopause crossings having (p, B(sub z)) values within each bin. We also calculate the bivariate (Inp, B(sub z)) moments to second order over each bin in control space. We can then calculate the six control-space expansion coefficients for each of the three ellipsoid parameters in configuration space. From these coefficients we can derive useful diagnosis of the magnetopause shape as joint functions of p and B(sub z): the aspect ratio of the ellipsoid's minor-to-major axes; the flank distance, radius of curvature, and flaring angle (at X(sub GSE) = 0); and the subsolar distance and radius of curvature. We confirm and quantify previous results that during periods of southward B(sub z) the subsolar magnetopause moves inward, while at X(sub GSE) = 0 the flank magnetopause moves outward and the flaring angle increases.

Roelof, Edmond C.↗
NASA NTRS2018

Nearfield Summary and Statistical Analysis of the Second AIAA Sonic Boom PredictionWorkshop

A summary is provided for the Second AIAA Sonic Boom Workshop held 8–9 January 2017 in conjunction with AIAA SciTech 2017. The workshop used three required models of increasing complexity: an axisymmetric body, a wing body, and a complete configuration with flow-through nacelle. An optional complete configuration with propulsion boundary conditions is also provided. These models are designed with similar nearfield signatures to isolate geometry and shock/expansion interaction e ects. Eleven international participant groups submitted nearfield signatures with forces, pitching moment, and iterative convergence norms. Statistics and grid convergence of these nearfield signatures are presented. These submissions are propagated to the ground, and noise levels are computed. This allows the grid convergence and the statistical distribution of a noise level to be computed. While progress is documented since the first workshop, improvement to the analysis methods for a possible subsequent workshop are provided. The complete configuration with flow-through nacelle showed the most dramatic improvement between the two workshops. The current workshop cases are more relevant to vehicles with lower loudness and have the potential for lower annoyance than the first workshop cases. The models for this workshop with quieter ground noise levels than the first workshop exposed weaknesses in analysis, particularly in convective discretization.

Park, Michael A.↗
DOE OSTI2023

Structure and magnetism of AlCoCrCuFeNi high-entropy alloy

We report on the investigation of magnetic and structural properties of AlCoCrCuFeNi, which is known to crystallize in a dual phase solid solution: the face-centred cubic (FCC) or the body-centred cubic (BCC). The results of neutron (NPD) and synchrotron powder diffraction (SXRD) allow to partially resolve magnetic information coming from BCC and FCC phases, which is impossible in the bulk magnetic measurements. Electron diffraction (PED) revealed that AlCoCrCuFeNi forms dendritic microstructure with the Cu-rich FCC phase and the Ni-rich BCC phase. Lattice parameters obtained from PED method are in good agreement with parameters obtained after refinement on the basis of powder X-ray diffraction measurements. The local crystal and electronic structure around Co was studied using Co K X-ray Absorption Spectroscopy (XAS). The magnetic measurements show that AlCoCrCuFeNi reveal a ferromagnetic transition at about 330K and displays magnetic hysteresis loop at the room temperature. Results from NPD suggest that the magnetic moment is mostly located in the BCC subsystem. The alloy shows soft magnetic properties. Saturated magnetizations (M s ), remanence ratio (M r /M s ) and coercivity (H c ) of the cast are estimated to be 45.10emu/g, 5.1% and 56 Oe at 300K, respectively. Finally, the BCC-FCC phase transformation up to 673K (400°C) was investigated using temperature dependent NPD, where a possible second BCC phase was identified.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

A11I-2104 Evaluation of the NASA GEOS Chemistry-Climate Model Coupled Atmosphere-Ocean Configuration for its Suitability to Simulations of Atmospheric Composition in a Changing Climate

The NASA Goddard Earth Observing System Chemistry-Climate Model (GEOSCCM) is a full tropospheric-stratospheric chemistry enabled configuration of the NASA GEOS Earth system model. Among the objectives for research performed with GEOSCCM is to understand the impacts of climate change on recovery of the stratospheric ozone layer over the twenty-first century. Simulations performed in the past with GEOSCCM did not have an interactive ocean and were driven with projected sea surface temperature and sea ice boundary conditions provided by external models. This method limits the application of GEOSCCM for climate study because it doesn’t allow the feedback of composition changes to ocean. For a full representation of the chemistry-climate feedbacks in the Earth System, we have developed a configuration of GEOSCCM that makes use of the MOM5 ocean general circulation model coupled to the GEOS atmospheric general circulation model. This configuration leverages development of the GEOS Sub-seasonal-to-Seasonal (S2S) prediction system, and includes interactive, radiatively coupled aerosols and two-moment, aerosol-aware cloud microphysics. W e have assessed the baseline model climate sensitivity with an Monday, 11 December 2023 08:30 - 12:50 Poster Hall A-C - South (Exhibition Level, South, MC) experiment performed under pre-industrial (e.g., year 1850) greenhouse gas conditions and a second experiment under 4xCO2 conditions, similar to the typical CMIP protocol. The GEOS model is demonstrated to have a climate sensitivity of a 2.6 K increase in global mean surface temperature for an equivalent doubling of CO2 from pre-industrial conditions, in line with current CMIP models. This configuration of the GEOS model is suitable for application to multidecadal to century-long simulations of climate system response to changing greenhouse gas levels. We report here our evaluation of the climate diagnostics of this configuration and discuss future directions for work with this model configuration, including an ongoing twenty-first century projection experiment based on the Chemistry-Climate Model Intercomparison project protocol.

models↗
DOE OSTI2021

Too big, too small, or just right? A benchmark assessment of density functional theory for predicting the spatial extent of the electron density of small chemical systems

Multipole moments are the first-order responses of the energy to spatial derivatives of the electric field strength. The quality of density functional theory prediction of molecular multipole moments thus characterizes errors in modeling the electron density itself, as well as the performance in describing molecules interacting with external electric fields. However, only the lowest non-zero moment is translationally invariant, making the higher-order moments origin-dependent. Therefore, instead of using the 3 × 3 quadrupole moment matrix, we utilize the translationally invariant 3 × 3 matrix of second cumulants (or spatial variances) of the electron density as the quantity of interest (denoted by K). The principal components of K are the square of the spatial extent of the electron density along each axis. A benchmark dataset of the principal components of K for 100 small molecules at the coupled cluster singles and doubles with perturbative triples at the complete basis set limit is developed, resulting in 213 independent K components. The performance of 47 popular and recent density functionals is assessed against this Var213 dataset. Several functionals, especially double hybrids, and also SCAN and SCAN0 predict reliable second cumulants, although some modern, empirically parameterized functionals yield more disappointing performance. The H, Li, and Be atoms, in particular, are challenging for nearly all methods, indicating that future functional development could benefit from the inclusion of their density information in training or testing protocols.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

Electromagnetic scattering by coated convex surfaces and wedges simulated by approximate boundary conditions

Asymptotic/high-frequency solutions are developed for analyzing the non-specular scattering mechanisms associated with coated convex surfaces and edges simulated by approximate boundary conditions. In particular, the standard impedance boundary conditions (SIBC's) and the second order generalized impedance boundary conditions (GIBC's) are employed for a characterization of the edge diffraction, creeping wave, and surface diffracted wave contributions. To study the creeping wave and surface diffracted wave mechanisms, rigorous UTD (uniform geometrical theory of diffraction) diffraction coefficients are developed for a convex coated cylinder simulated with SIBC's and GIBC's. The ray solutions obtained remain valid in the transition region and reduce uniformly to those in the deep lit and shadow regions. A uniform asymptotic solution is also presented for observations in the close vicinity of the cylinder. The diffraction coefficient for a convex cylinder are obtained via a generalization of the corresponding ones of the circular cylinder. To validate the asymptotic/high-frequency solution, integral equations are derived for both E and H-polarization and solved numerically using the method of moments. Results are presented for a single and three layered coated convex cylinder. Some insights are also provided on the accuracy of the employed GIBC's versus SIBC's for application to curved surfaces. To characterize the scattering by impedance wedges illuminated at skew incidence, diffraction coefficients are derived from an approximate solution of the governing functional difference equations. This solution exactly recovers the known ones for an impedance half plane or an arbitrary wedge at normal incidence, and to validate it for other wedge angles, a moment method code was used. Finally, to test the usefulness of the approximate skew incidence impedance wedge diffraction coefficient for three dimensional structures, equivalent currents are derived in the context of the physical theory of diffraction (PTD) for a finite length impedance wedge of arbitrary internal angles. These are incorporated in a standard general purpose PTD code and results are presented for a number of different impedance structures.

Syed, H. H.↗
DOE OSTI2022

An Energy-Integrated Analysis for Measuring the Anomalous Precession Frequenzc for the Muon g - 2 Experiment at Fermilab

In the search for physics beyond the Standard Model, the Muon g - 2 Experiment at Fermilab (E989) will make the most precise measurement of the anomalous magnetic moment of the muon, a μ . Improvements in precision come from both increased statistics and new techniques to significantly reduce previous systematic uncertainties. The muon a μ is determined by extracting both the anomalous spin precession frequency, ω a , and the average magnetic field sampled by the muons, B. Traditionally an energy threshold analysis method which requires reconstruction of decay positrons from the muon decay, μ + → e + ν e ν μ has been used to extract ω a . This thesis will describe a novel, energy-integrated analysis method for ωa extraction that does not require full reconstruction of positron hits and lowers the energy threshold therefore sampling more of the stored muons. This new technique has different sensitivities to effects like gain fluctuations and positron pulse pile-up. The analysis results will focus on data taken in 2019, the second year of running the experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1938

Compressible flow about symmetrical Joukowski profiles

The method of Poggi is employed for the determination of the effects of compressibility upon the flow past an obstacle. A general expression for the velocity increment due to compressibility is obtained. The general result holds whatever the shape of the obstacle; but, in order to obtain the complete solution, it is necessary to know a certain Fourier expansion of the square of the velocity of flow past the obstacle. An application is made to the case flow of a symmetrical Joukowski profile with a sharp trailing edge, fixed in a stream of an arbitrary angle of attack and with the circulation determined by the Kutta condition. The results are obtained in a closed form and are exact insofar as the second approximation to the compressible flow is concerned, the first approximation being the result for the corresponding incompressible flow. Formulas for lift and moment analogous to the Blasius formulas in incompressible flow are developed and are applied to thin symmetrical Joukowski profiles for small angles of attack.

Kaplan, Carl↗
DOE OSTI2022

Effect of Initial Water Saturation on The Performance of Fracturing Fluids With and Without Polyallylamine under Simulated EGS Conditions

Objectives/Scope: StimuFrac (US Patents 9,873,828 B2 and 9,447,315 B2), a CO2-reactive polymer aqueous solution [polyallylamine (PAA) 1wt% in water] combined with CO2, can be used as a less water-intensive fracturing fluid for enhanced geothermal systems (EGS). Our previous results show that in hot dry rock (HDR), PAA/CO2 fracturing fluids outperformed other fluids such as water, CO2, and CO2/water in generating large fractures with less fluid consumed. The objective of this work is to study the effect of initial water saturation of rock on the performance of StimuFrac fluid in ½ foot cubic rock samples and under representative EGS pressure/temperature conditions using cyclic and constant flow rate injection strategies. The fracturing results are compared with results using different fracturing fluids in terms of controlling fracture propagation rates, fracture hydraulic conductivity, breakdown pressures and volumes of fluids required. Methods/Procedures/Process: In all tests, water was initially injected into the rock to increase the water saturation before the fracturing processes to simulate actual geothermal reservoir conditions. For the cyclic injection, one complete cycle consisted of (1) a PAA slug (or water slug) injection followed by (2) CO2 injection to initiate the fracture. In the second step of the first cycle, the pressure of CO2 is increased until a maximum pressure is reached (fracture is initiated at this moment), and then the injection of CO2 is allowed for another 30 seconds to propagate the fracture. Then, the two-step cycle of PAA followed by CO2 injection (up to 2-4 mL/min) was continued. For the constant flow rate injection strategies, the initial water saturation was increased by injecting water at 1000 psi and 200°C for three days. After that, an initial slug of water, CO2, or PAA was injected and then fracturing was initiated and propagated by injecting CO2 at a constant flow rate. Applications/Significance/Novelty: The results of this study suggest that water saturation, especially near the wellbore region, will significantly affect the fracturing fluid transmission into the rock porous media by changing the relative permeability of CO2 or water, thus affecting the fracture initiation and propagation. In this study, fracturing with cyclic injection or constant flow rate injection methods were performed using three different kinds of fluids systems. These fluids are water, CO2, or CO2 with PAA. Splitting the rock samples in half after fracturing reveals that the fracture propagation is significantly limited under these high water saturated conditions compared to dry initial conditions: The fractures propagate less than 1/3 length of the distance from the wellbore to rock surface, and in some cases no fracture is generated. This may be caused by the fact that leak-off is dominating the fracturing process and the injected fluid flow rate is not high enough to overcome the leak-off even under high flow rate injection conditions. Additionally, CO2 could be leaking off into the wellbore annulus and this may be making it more difficult to generate pressure gradients away from the near-wellbore region.

Jian, Guoqing↗
DOE OSTI2022

Multifunctional Cu 2 TSiS 4 (T = Mn and Fe): Polar Semiconducting Antiferromagnets with Nonlinear Optical Properties

In this work, Cu 2 TSiS 4 (T = Mn and Fe) polycrystalline and single-crystal materials were prepared with high-temperature solid-state and chemical vapor transport methods, respectively. The polar crystal structure (space group Pmn2 1 ) consists of chains of corner-sharing and distorted CuS 4 , Mn/FeS 4 , and SiS 4 tetrahedra, which is confirmed by Rietveld refinement using neutron powder diffraction data, X-ray single-crystal refinement, electron diffraction, energy-dispersive X-ray spectroscopy, and second harmonic generation (SHG) techniques. Magnetic measurements indicate that both compounds order antiferromagnetically at 8 and 14 K, respectively, which is supported by the temperature-dependent (100–2 K) neutron powder diffraction data. Additional magnetic reflections observed at 2 K can be modeled by magnetic propagation vectors k = (1/2,0,1/2) and k = (1/2,1/2,1/2) for Cu 2 MnSiS 4 and Cu 2 FeSiS 4 , respectively. The refined antiferromagnetic structure reveals that the Mn/Fe spins are canted away from the ac plane by about 14°, with the total magnetic moments of Mn and Fe being 4.1(1) and 2.9(1) μ B , respectively. Both compounds exhibit an SHG response with relatively modest second-order nonlinear susceptibilities. Density functional theory calculations are used to describe the electronic band structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1956

The Interference Effects of a Body on the Spanwise Load Distributions of Two 45 Degree Sweptback Wings of Aspect Ratio 8.02 from Low-Speed Tests

Tests of two wing-body combinations have been conducted in the Langley 19-foot pressure tunnel at a Reynolds number of 4 x 10(exp 6) and a Mach number of 0.19 to determine the effects of the bodies on the wing span load distributions. The wings had 45 degrees sweepback of the quarter-chord line, aspect ratio 8.02, taper ratio 0.45, and incorporated 12-percent-thick airfoil sections streamwise. One wing was untwisted and uncambered whereas the second wing incorporated both twist and camber. Identical bodies of revolution, of 10:1 fineness ratio, having diameter-to-span ratios of 0.10, were mounted in mid-high-wing arrangements. The effects of wind incidence, wing fences, and flap deflection were determined for the plane uncambered wing. The addition of the body to the plane wing increased the exposed wing loading at a given lift coefficient as much as 10 percent with the body at 0 degrees incidence and 4 percent at 4 degrees incidence. The bending-moment coefficients at the wing-body juncture were increased about 2 percent with the body at 0 degrees incidence, whereas the increases were as much as 10 percent with the body at 4 degrees incidence. The spanwise load distributions due to the body on the plane wing as calculated by using a swept-wing method employing 19 spanwise lifting elements and control points generally showed satisfactory agreement with experiment. The spanwise load distributions due to body on the flapped plane wing and on the twisted and cambered wing were dissimilar to those obtained on the plane wing. Neither of the methods of calculation which were employed yielded distributions that agreed consistently with experiment for either the flapped plane wing or the twisted and cambered wing.

Martina, Albert P.↗
NASA NTRS2013

Calibration Variable Selection and Natural Zero Determination for Semispan and Canard Balances

Independent calibration variables for the characterization of semispan and canard wind tunnel balances are discussed. It is shown that the variable selection for a semispan balance is determined by the location of the resultant normal and axial forces that act on the balance. These two forces are the first and second calibration variable. The pitching moment becomes the third calibration variable after the normal and axial forces are shifted to the pitch axis of the balance. Two geometric distances, i.e., the rolling and yawing moment arms, are the fourth and fifth calibration variable. They are traditionally substituted by corresponding moments to simplify the use of calibration data during a wind tunnel test. A canard balance is related to a semispan balance. It also only measures loads on one half of a lifting surface. However, the axial force and yawing moment are of no interest to users of a canard balance. Therefore, its calibration variable set is reduced to the normal force, pitching moment, and rolling moment. The combined load diagrams of the rolling and yawing moment for a semispan balance are discussed. They may be used to illustrate connections between the wind tunnel model geometry, the test section size, and the calibration load schedule. Then, methods are reviewed that may be used to obtain the natural zeros of a semispan or canard balance. In addition, characteristics of three semispan balance calibration rigs are discussed. Finally, basic requirements for a full characterization of a semispan balance are reviewed.

Ulbrich, Norbert M.↗
NASA NTRS1988

Radiation and scattering from loaded microstrip antennas over a wide bandwidth

The integral equation and moment method solution is developed for two different antennas in the presence of an infinite grounded dielectric substrate. The first antenna is a rectangular microstrip patch antenna. This antenna is analyzed for excitation by an incident plane wave in free space and a vertical filament of uniform current in the dielectric. This antenna can be loaded by a lumped impedance in a vertical filament of uniform current extending from the patch through the dielectric to the ground plane. The radar cross section of the microstrip antenna is found from the plane wave excitation and shows good agreement to measurement for both an unloaded and loaded antenna. The input impedance is found from the current filament excitation. This is compared to the measured input impedance of a coaxially fed microstrip antenna and shows good agreement for both unloaded and loaded antennas when the dielectric substrate is much less than a wavelength. The second antenna is a vertical thin wire extending from the ground plane into or through the dielectric substrate. The mutual impedance between two imbedded monopoles is compared to a previous calculation.

Forrai, D. P.↗
DOE OSTI2023

Adjoint Waveform Tomography for Crustal and Upper Mantle Structure the Middle East and Southwest Asia for Improved Waveform Simulations Using Openly Available Broadband Data

We present a new model of radially anisotropic seismic wavespeeds for the crust and upper mantle of a broad region of the Middle East and Southwest Asia (MESWA) derived from adjoint waveform tomography. We inverted waveforms from 192 Global Centroid Moment Tensor earthquakes (MW 5.5-7.0) recorded by over 1000 openly available broadband seismic stations from permanent and temporary networks in the region. Spatial coverage of the available data is highly uneven due to earthquakes clustered along plate boundaries and sparse coverage of open seismic networks in the region. We considered three possible starting models: the SPiRaL global model (Simmons et al., 2021); MEC-1 (Kaviani et al., 2020); and CSEM2.0 (Noe et al., 2023). Because the SPiRaL model provides good fits to the observed waveforms measured by the time-bandwidth product of selected windows in several period bands, provides all the necessary parameters and covers the entire domain we used it for the starting model with the period band 50-100 seconds. Inversion iterations proceeded using time-frequency phase misfits in six stages and 54 total iterations reducing the minimum period to 30 seconds. Our final model, MESWA, provides improved waveform fits compared to the starting model for both the data used in the inversion and an independent validation data set of 66 events. Two metrics of waveform fit (the time-frequency phase misfit used in the optimization and normalized L2 misfit) were both reduced by nearly 60% for both data sets and MESWA provides significantly larger misfit reductions relative to the SPiRaL model than the MEC-1 or CSEM models. We also find that MESWA provides a larger time-bandwidth product of selected windows indicating that more information content of the observed waveforms is explained by MESWA than the other models. Our new model reveals tectonic features imaged by other studies and methods but in a new holistic model of shear and compressional wavespeeds (v S and v P , respectively) with anisotropy covering the crust and uppermost mantle of a larger domain. MESWA has smaller scale-length features and tends to sharpen some features relative to the SPiRaL starting model. Examples include: low crustal v S in the TurkishIranian Plateau, Zagros Mountains, Afghan Central Blocks and Sulaiman Fold Belt; low mantle vSfollowing divergent (Gulf of Aden, Red Sea) and transform (Dead Sea Fault) margins of the Arabian Plate; low and high v S in the mantle beneath the Arabian Shield and Platform, respectively. Low vS is imaged below Cenozoic volcanic centers of the Arabian Peninsula, the so-called Mecca-Madina-Nafud (MMN) Line. Positive anisotropy (v SH > v SV ) is inferred for asthenospheric depths across the region except where up/downwelling may influence fabric alignment (e.g. Afar, Red Sea, Arabian Shield). Elevated vS tracks Makran subduction under southeast Iran. MESWA resembles the SPiRaL model in its long-wavelength structure, but enhances shorter wavelengths features on the order of 200 km and smaller. The resulting model could be used for as a starting model for further improvements, say using waveforms from in-country seismic networks that are not openly available or smaller-scale studies targeting shorter period waveforms. The model also could be used for source characterization and moment tensor inversion to improve earthquake hazard studies and nuclear explosion monitoring.

58 GEOSCIENCES↗