Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Smooth approximation”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

The COBE normalization for standard cold dark matter

The Cosmic Background Explorer Satellite (COBE) detection of microwave anisotropies provides the best way of fixing the amplitude of cosmological fluctuations on the largest scales. This normalization is usually given for an n = 1 spectrum, including only the anisotropy caused by the Sachs-Wolfe effect. This is certainly not a good approximation for a model containing any reasonable amount of baryonic matter. In fact, even tilted Sachs-Wolfe spectra are not a good fit to models like cold dark matter (CDM). Here, we normalize standard CDM (sCDM) to the two-year COBE data and quote the best amplitude in terms of the conventionally used measures of power. We also give normalizations for some specific variants of this standard model, and we indicate how the normalization depends on the assumed values on n, Omega(sub B) and H(sub 0). For sCDM we find the mean value of Q = 19.9 +/- 1.5 micro-K, corresponding to sigma(sub 8) = 1.34 +/- 0.10, with the normalization at large scales being B = (8.16 +/- 1.04) x 10(exp 5)(Mpc/h)(exp 4), and other numbers given in the table. The measured rms temperature fluctuation smoothed on 10 deg is a little low relative to this normalization. This is mainly due to the low quadrupole in the data: when the quadrupole is removed, the measured value of sigma(10 deg) is quite consistent with the best-fitting the mean value of Q. The use of the mean value of Q should be preferred over sigma(10 deg), when its value can be determined for a particular theory, since it makes full use of the data.

Bunn, Emory F.↗
NASA NTRS2006

Electric Fields, Cloud Microphysics, and Reflectivity in Anvils of Florida Thunderstorms

A coordinated aircraft - radar project that investigated the electric fields, cloud microphysics and radar reflectivity of thunderstorm anvils near Kennedy Space Center is described. Measurements from two cases illustrate the extensive nature of the microphysics and electric field observations. As the aircraft flew from the edges of anvils into the interior, electric fields very frequently increased abruptly from approx.1 to >10 kV/m even though the particle concentrations and radar reflectivity increased smoothly. The abrupt increase in field usually occurred when the aircraft entered regions with a reflectivity of 10 to 15 dBZ. It is suggested that the abrupt increase in electric field may be because the charge advection from the storm core did not occur across the entire breadth of the anvil and was not constant in time. Screening layers were not detected near the edges of the anvils. Some long-lived anvils showed subsequent enhancement of electric field and reflectivity and growth of particles, which if localized, might be a factor in explaining the abrupt change of field in some cases. Comparisons of electric field magnitude with particle concentration or reflectivity for a combined data set that included all anvil measurements showed a threshold behavior. When the average reflectivity, such as in a 3-km cube, was less than approximately 5 dBZ, the electric field magnitude was <3 kV/m. Based on these findings, the Volume Averaged Height Integrated Radar Reflectivity (VAHIRR) is now being used by NASA, the Air Force and Federal Aviation Administration in new Lightning Launch Commit Criteria as a diagnostic for high electric fields in anvils.

Dye, J. E.↗
DOE OSTI2023

Extreme sensitivity of higher-order interatomic force constants and thermal conductivity to the energy surface roughness of exchange-correlation functionals

In this Letter, we report that the fourth-order interatomic force constants (4th-IFCs) are significantly sensitive to the energy surface roughness of exchange-correlation (XC) functionals in density functional theory calculations. This sensitivity, which is insignificant for the second- (2nd-) and third-order (3rd-) IFCs, varies for different functionals in different materials and can cause misprediction of thermal conductivity by several times of magnitude. As a result, when calculating the 4th-IFCs using the finite difference method, the atomic displacement needs to be taken large enough to overcome the energy surface roughness, in order to accurately predict phonon lifetime and thermal conductivity. We demonstrate this phenomenon on a benchmark material (Si), a high-thermal conductivity material (BAs), and a low thermal conductivity material (NaCl). For Si, we find that the LDA, PBE, and PBEsol XC functionals are all smooth to the 2nd- and 3rd-IFCs but all rough to the 4th-IFCs. This roughness can lead to a prediction of nearly one order of magnitude lower thermal conductivity. For BAs, all three functionals are smooth to the 2nd- and 3rd-IFCs, and only the PBEsol XC functional is rough for the 4th-IFCs, which leads to a 40% underestimation of thermal conductivity. For NaCl, all functionals are smooth to the 2nd- and 3rd-IFCs but rough to the 4th-IFCs, leading to a 70% underprediction of thermal conductivity at room temperature. In conclusion, with these observations, we provide general guidance on the calculation of 4th-IFCs for an accurate thermal conductivity prediction.

42 ENGINEERING↗
NASA NTRS1985

Lubrication of nonconformal contacts

Minimum film thickness results for piezoviscous-rigid regime of lubrication are developed for a compressible Newtonian fluid with Roelands viscosity. The results provide a basis for the analysis and design of a wide range of machine elements operating in the piezoviscous-rigid regime of lubrication. A new numerical method of calculating elastic deformation in contact stresses is developed using a biquadratic polynomial to approximate the pressure distribution on the whole domain analyzed. The deformation of every node is expressed as a linear combination of the nodal pressures whose coefficients can be combined into an influence coefficient matrix. This approach has the advantages of improved numerical accuracy, less computing time and smaller storage size required for influence matrix. The ideal elastohydrodynamic lubrication is extended to real bearing systems in order to gain an understanding of failure mechanisms in machine elements. The improved elastic deformation calculation is successfully incorporated into the EHL numerical scheme. Using this revised numerical technique and the flow factor model developed by Patir and Cheng (1978) the surface roughness effects on the elastohydrodynamic lubrication of point contact is considered. Conditions typical of an EHL contact in the piezoviscous-elastic regime entrained in pure rolling are investigated. Results are compared with the smooth surface solutions. Experiments are conducted to study the transient EHL effects in instrument ball bearings.

Jeng, Y. R.↗
DOE OSTI2024

FLOWERS AEP: An Analytical Model for Wind Farm Layout Optimization

Annual energy production (AEP) is commonly used in objective functions for wind farm layout optimization. AEP is proportional to wind farm power production integrated over an annual distribution of free-stream wind conditions. Physics-based estimates of wind farm power production typically rely on low-fidelity engineering wake models that approximate the steady-state wind farm flow field. AEP estimates are then obtained by performing independent simulations for discrete wind conditions and using rectangular quadrature to account for each condition's expected frequency of occurrence. Depending on the number of simulated discrete wind conditions, this numerical integral could be hampered by poor accuracy or high computational costs. The FLOWERS AEP model instead poses an analytical integral of the engineering wake model over the variable wind conditions, yielding a closed-form, analytical function for wind farm AEP. This paper derives the analytical functions for FLOWERS AEP and its derivatives with respect to turbine position, which are useful for gradient-based wind farm layout optimization, in nondimensional form. We then analyze the benefits of the FLOWERS AEP model over conventional reference models, focusing on its low cost, adequate wake loss predictions, and smooth design space. Although the FLOWERS approach is found to predict the exact value of AEP with some error relative to the reference model (within 14% on average), it dramatically reduces computation time by an order of magnitude, produces a qualitatively similar design space at relatively low resolution, and yields comparable optimal layouts. This significant speed improvement is critical in layout optimization applications, where determining an optimal layout in an efficient manner is more important than precise AEP prediction.

17 WIND ENERGY↗
NASA NTRS2014

Predicting Solar Cycle 24 Using a Geomagnetic Precursor Pair

We describe using Ap and F(10.7) as a geomagnetic-precursor pair to predict the amplitude of Solar Cycle 24. The precursor is created by using F(10.7) to remove the direct solar-activity component of Ap. Four peaks are seen in the precursor function during the decline of Solar Cycle 23. A recurrence index that is generated by a local correlation of Ap is then used to determine which peak is the correct precursor. The earliest peak is the most prominent but coincides with high levels of non-recurrent solar activity associated with the intense solar activity of October and November 2003. The second and third peaks coincide with some recurrent activity on the Sun and show that a weak cycle precursor closely following a period of strong solar activity may be difficult to resolve. A fourth peak, which appears in early 2008 and has recurrent activity similar to precursors of earlier solar cycles, appears to be the "true" precursor peak for Solar Cycle 24 and predicts the smallest amplitude for Solar Cycle 24. To determine the timing of peak activity it is noted that the average time between the precursor peak and the following maximum is approximately equal to 6.4 years. Hence, Solar Cycle 24 would peak during 2014. Several effects contribute to the smaller prediction when compared with other geomagnetic-precursor predictions. During Solar Cycle 23 the correlation between sunspot number and F(10.7) shows that F(10.7) is higher than the equivalent sunspot number over most of the cycle, implying that the sunspot number underestimates the solar-activity component described by F(10.7). During 2003 the correlation between aa and Ap shows that aa is 10 % higher than the value predicted from Ap, leading to an overestimate of the aa precursor for that year. However, the most important difference is the lack of recurrent activity in the first three peaks and the presence of significant recurrent activity in the fourth. While the prediction is for an amplitude of Solar Cycle 24 of 65 +/- 20 in smoothed sunspot number, a below-average amplitude for Solar Cycle 24, with maximum at 2014.5+/-0.5, we conclude that Solar Cycle 24 will be no stronger than average and could be much weaker than average.

predictions↗
NASA NTRSDate not supplied

Visualization and 3D Mesh Generation for ASSEMBLERS Range of Motion

In an effort to create a base and foundation for ASSEMBLERS with an informed layout, data representing the range of motion of varying heights of Stewart platforms were analyzed. 3D shapes were generated first to be approximately representative of these data in the form of rotationally symmetric shells, making use of 2D alpha shapes to create a border that did not account for asymmetry inherent in the true range of motion. Further work was performed to generate mesh files more fully representative of reach, creating 3D alpha shapes based directly on point clouds generated from uniformly random platform positions. These shapes were initially not representative of the entire range of motion of a platform stack due to the low likelihood of the most extreme positions being added to a given point cloud but were more representative of the asymmetry inherent in the possible positions of a Stewart platform as they no longer relied on a single 2D contour. The same algorithm for generating 3D alpha shapes was then applied to point clouds designed to represent extreme positions, generating larger smooth shells with six clear planes of symmetry. All shells were converted to STLs to aid in further modeling of range of motion for base layouts and in visualizing reachable area for varying stack heights, successfully creating closed shells representative of typical and extreme positions for platform stacks of varying heights.

Jonah DeGuire Vanke↗
NASA NTRS1993

The three-dimensional structure of cumulus clouds over the ocean. 1: Structural analysis

Thermal channel (channel 6, 10.4-12.5 micrometers) images of five Landsat thematic mapper cumulus scenes over the ocean are examined. These images are thresholded using the standard International Satellite Cloud Climatology Project (ISCCP) thermal threshold algorithm. The individual clouds in the cloud fields are segmented to obtain their structural statistics which include size distribution, orientation angle, horizontal aspect ratio, and perimeter-to-area (PtA) relationship. The cloud size distributions exhibit a double power law with the smaller clouds having a smaller absolute exponent. The cloud orientation angles, horizontal aspect ratios, and PtA exponents are found in good agreement with earlier studies. A technique also is developed to recognize individual cells within a cloud so that statistics of cloud cellular structure can be obtained. Cell structural statistics are computed for each cloud. Unicellular clouds are generally smaller (less than or equal to 1 km) and have smaller PtA exponents, while multicellular clouds are larger (greater than or equal to 1 km) and have larger PtA exponents. Cell structural statistics are similar to those of the smaller clouds. When each cell is approximated as a quadric surface using a linear least squares fit, most cells have the shape of a hyperboloid of one sheet, but about 15% of the cells are best modeled by a hyperboloid of two sheets. Less than 1% of the clouds are ellipsoidal. The number of cells in a cloud increases slightly faster than linearly with increasing cloud size. The mean nearest neighbor distance between cells in a cloud, however, appears to increase linearly with increasing cloud size and to reach a maximum when the cloud effective diameter is about 10 km; then it decreases with increasing cloud size. Sensitivity studies of threshold and lapse rate show that neither has a significant impact upon the results. A goodness-of-fit ratio is used to provide a quantitative measure of the individual cloud results. Significantly improved results are obtained after applying a smoothing operator, suggesting the eliminating subresolution scale variations with higher spatial resolution may yield even better shape analyses.

Kuo, Kwo-Sen↗
DOE OSTI2020

Matrix-free subcell residual distribution for Bernstein finite elements: Monolithic limiting

This paper is focused on the aspects of limiting in residual distribution (RD) schemes for high-order finite element approximations to advection problems. Both continuous and discontinuous Galerkin methods are considered in this work. Discrete maximum principles are enforced using algebraic manipulations of element contributions to the global nonlinear system. The required modifications can be carried out without calculating the element matrices and assembling their global counterparts. The components of element vectors associated with the standard Galerkin discretization are manipulated directly using localized subcell weights to achieve optimal accuracy. Low-order nonlinear RD schemes of this kind were originally developed to calculate local extremum diminishing predictors for flux-corrected transport (FCT) algorithms. In the present paper, we incorporate limiters directly into the residual distribution procedure, which makes it applicable to stationary problems and leads to well-posed nonlinear discrete problems. To circumvent the second-order accuracy barrier, the correction factors of monolithic limiting approaches and FCT schemes are adjusted using smoothness sensors based on second derivatives. Finally, the convergence behavior of presented methods is illustrated by numerical studies for two-dimensional test problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Chiral spin symmetry and the QCD phase diagram

Lattice QCD simulations with chirally symmetric quarks have recently established approximate SU(2) CS and SU(2N F ) symmetries of the quantum effective action in a temperature range above the chiral crossover T ch , in which color-electric interactions between quarks dominate the dynamics. We show that such an intermediate temperature range between the chirally broken and plasma regimes is fully consistent with published screening mass spectra, which demonstrate the breakdown of thermal perturbation theory at the crossover between the partonic and the chiral spin symmetric regime at T s ~(2–3)T ch . From the known behavior of screening masses with baryon chemical potential, we deduce qualitatively how this chiral spin symmetric band extends into the QCD phase diagram. In the cold and dense region, we propose parity doubled baryons as possible candidates for chiral spin symmetric matter. This represents a special case of quarkyonic matter with confinement and restored chiral symmetry, and can smoothly transform to quark matter at sufficiently high densities. Finally, we discuss the potential of dilepton spectra to identify such matter forms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1986

Thermodynamic evaluation of transonic compressor rotors using the finite volume approach

A method was developed which calculates two-dimensional, transonic, viscous flow in ducts. The finite volume, time marching formulation is used to obtain steady flow solutions of the Reynolds-averaged form of the Navier Stokes equations. The entire calculation is performed in the physical domain. The method is currently limited to the calculation of attached flows. The features of the current method can be summarized as follows. Control volumes are chosen so that smoothing of flow properties, typically required for stability, is now needed. Different time steps are used in the different governing equations to improve the convergence speed of the viscous calculations. A new pressure interpolation scheme is introduced which improves the shock capturing ability of the method. A multi-volume method for pressure changes in the boundary layer allows calculations which use very long and thin control volumes. A special discretization technique is also used to stabilize these calculations. A special formulation of the energy equation is used to provide improved transient behavior of solutions which use the full energy equation. The method is then compared with a wide variety of test cases. The freestream Mach numbers range from 0.075 to 2.8 in the calculations. Transonic viscous flow in a converging diverging nozzle is calculated with the method; the Mach number upstream of the shock is approximately 1.25. The agreement between the calculated and measured shock strength and total pressure losses is good. Essentially incompressible turbulent boundary layer flow in a adverse pressure gradient is calculated and the computed distribution of mean velocity and shear stress are in good agreement with the measurements. At the other end of the Mach number range, a flat plate turbulent boundary layer with a freestream Mach number of 2.8 is calculated using the full energy equation; the computed total temperature distribution and recovery factor agree well with the measurements when a variable Prandtl number is used through the boundary layer.

Nicholson, S.↗
NASA NTRS1992

Adaptive filtering of Echelle spectra of distant Quasars

The study of the Ly alpha - forest of distant (approximately greater than 3) Quasars is an important tool in obtaining a more detailed picture of the distribution of matter along the line of sight and thus of the general distribution of matter in the Universe and is therefore of important cosmological significance. Obviously, this is one of the tasks where spectral resolution plays an important role. The spectra used were obtained with the EFOSC at the ESO 3.6m telescope. Applying for the data reduction the standard Echelle procedure, as it is implemented for instance in the MIDAS-package, one uses stationary filters (e.g. median) for noise and cosmic particle event reduction in the 2-dimensional Echelle image. These filters are useful if the spatial spectrum of the noise reaches essentially higher frequencies then the highest resolution features in the image. Otherwise the resolution in the data will be degraded and the spectral lines smoothed. However, in the Echelle spectra the highest resolution is already in the range of one or a few pixels and therefore stationary filtering means always a loss of resolution. An Echelle reduction procedure on the basis of a space variable filter described which recognizes the local resolution in the presence of noise and adapts to it is developed. It was shown that this technique leads to an improvement in resolution by a factor of 2 with respect to standard procedures.

Priebe, A.↗
NASA NTRS2006

The Hard X-ray 20-40 keV AGN Luminosity Function

We have compiled a complete, significance limited extragalactic sample based on approximately 25,000 deg(sup 2) to a limiting flux of 3 x 10(exp -11) ergs per square centimeter per second. (approximately 7,000 deg(sup 2)) to a flux limit of 10(exp -11) ergs per square centimeter per second)) in the 20 - 40 keV band with INTEGRAL. We have constructed a detailed exposure map to compensate for effects of non-uniform exposure. The flux-number relation is best described by a power-law with a slope of alpha = 1.66 plus or minus 0.11. The integration of the cumulative flux per unit area leads to f(sub 20-40 keV) = 2.6 x 10(exp -10) ergs per square centimeter per second per sr(sup -1) which is about 1% of the known 20-40 keV X-ray background. We present the first luminosity function of AGN in the 20-40 keV energy range, based on 68 extragalactic objects detected by the imager IBIS/ISGRI on-board INTEGRAL. The luminosity function shows a smoothly connected two power-law form, with an index of gamma (sub 1) = 0.9 below, and gamma (sub 2) = 2.2 above the turn-over luminosity of L(sub *), = 4.6 x 10(sup 43) ergs per second. The emissivity of all INTEGRAL AGNs per unit volume is W(sub 20-40keV)(greater than 10(sup 41) ergs per second) = 2.8 x 10(sup 38) ergs per second h(sup 3)(sub 70) Mpc(sup -3). These results are consistent with those derived in the 2-20keV energy band and do not show a significant contribution by Compton-thick objects. Because the sample used in this study is truly local (z(raised bar) = 0.022)), only limited conclusions can be drawn for the evolution of AGNs in this energy band. But the objects explaining the peak in the cosmic X-ray background are likely to be either low luminosity AGN (L(sub x) less than 10(sup 41) ergs per second) or of other type, such as intermediate mass black holes, clusters, and star forming regions.

Beckmann, V.↗
NASA NTRS2007

Fracture Strength of Single-Crystal Silicon Carbide Microspecimens at Room and Elevated Temperature

Three shapes of tensile specimens were tested--curved with a very low stress concentration factor and straight with either a circular hole or an elliptical hole. The nominal thickness was 125 micron with a net section 100 micron wide; the overall length of these microspecimens was 3.1 mm. They were fabricated by an improved version of deep reactive ion etching, which produced specimens with smooth sidewalls and cross-sections having a slightly trapezoidal shape that was exaggerated inside the holes. The novel test setup used a vertical load train extending into a resistance furnace. The specimens had wedge-shaped ends which fit into ceramic grips. The fixed grip was mounted on a ceramic post, and the movable grip was connected to a load cell and actuator outside the furnace with a ceramic-encased nichrome wire. The same arrangement was used for tests at 24 and at 1000 C. The strengths of the curved specimens for two batches of material (made with slightly different processes) were 0.66+/-0.12 GPa and 0.45+/-0.20 GPa respectively at 24 C with identical values at 1000 C. The fracture strengths of the circular-hole and elliptical-hole specimens (computed from the stress concentration factors and measured loads at failure) were approximately 1.2 GPa with slight decreases at the higher temperature. Fractographic examinations showed failures initiating on the surface--primarily at corners. Weibull predictions of fracture strengths for the hole specimens based on the properties of the curved specimens were reasonably effective for the circular holes, but not for the elliptical holes.

Nemeth, Noel N.↗
NASA NTRS1994

Numerical study of the effects of icing on viscous flow over wings

An improved hybrid method for computing unsteady compressible viscous flows is presented. This method divides the computational domain into two zones. In the outer zone, the unsteady full-potential equation (FPE) is solved. In the inner zone, the Navier-Stokes equations are solved using a diagonal form of an alternating-direction implicit (ADI) approximate factorization procedure. The two zones are tightly coupled so that steady and unsteady flows may be efficiently solved. Characteristic-based viscous/inviscid interface boundary conditions are employed to avoid spurious reflections at that interface. The resulting CPU times are less than 60 percent of that required for a full-blown Navier-Stokes analysis for steady flow applications and about 60 percent of the Navier-Stokes CPU times for unsteady flows in non-vector processing machines. Applications of the method are presented for a rectangular NACA 0012 wing in low subsonic steady flow at moderate and high angles of attack, and for an F-5 wing in steady and unsteady subsonic and transonic flows. Steady surface pressures are in very good agreement with experimental data and are essentially identical to Navier-Stokes predictions. Density contours show that shocks cross the viscous/inviscid interface smoothly, so that the accuracy of full Navier-Stokes equations can be retained with a significant savings in computational time.

Sankar, L. N.↗
NASA NTRS2002

Thermal and Nonthermal X-ray Emission from the Forward Shock in Tycho's Supernova Remnant

We present Chandra CCD images of Tycho's supernova remnant that delineate its outer shock, seen as a thin, smooth rim along the straight northeastern edge and most of the circular western half. The images also show that the Si and S ejecta are highly clumpy, and have reached the forward shock at numerous locations. Most of the X-ray spectra that we examine along the rim show line emission from Si and S, which in some cases must come from ejecta; the continuum is well represented by either thermal or nonthermal models. In the case that the continuum is assumed to be thermal, the temperatures at the rim are all similar at about 2 keV, and the ionization ages are very low because of the overall weakness of the line emission. Assuming shock velocities inferred from radio and X-ray expansion measurements, these temperatures are substantially below those expected for equilibration of the electron and ion temperatures; electron to mean temperature ratios of approximately less than 0.1 - 0.2 indicate at most modest collisionless heating of the electrons at the shock. The nonthermal contribution to these spectra may be important, however, and may account for as many as half of the counts in the 4-6 keV energy range, based on an extrapolation of the hard X-ray spectrum above 10 keV.

Hwang, Una↗
DOE OSTI2026

Back-to-back dijet production in DIS at arbitrary Bjorken x: TMD gluon distributions to twist-3 accuracy

We derive the gluon transverse-momentum-dependent (TMD) operator structure of back-to-back\\\\r\\\\nquark–antiquark dijet production in deep inelastic scattering at arbitrary Bjorken-x to twist-3 ac\\\\r\\\\ncuracy. Working at leading order in the strong coupling and in the kinematic regime where the\\\\r\\\\ntransverse momentum imbalance of the jets is much smaller than their individual transverse mo\\\\r\\\\nmenta, we perform a systematic gradient expansion of the quark propagator in a background gluon\\\\r\\\\nfield. This expansion organizes multiple interactions with the target in terms of longitudinal Wilson\\\\r\\\\nlines and gauge-invariant field-strength insertions, yielding a TMD description valid beyond the\\\\r\\\\nstrict high-energy eikonal (x → 0) approximation. We obtain explicit cross sections for longitudi\\\\r\\\\nnally and transversely polarized virtual photons, identifying all contributing gluon TMD operators\\\\r\\\\nup to twist-3, including structures involving F+−, Fij, and three-gluon correlators. The full lon\\\\r\\\\ngitudinal phase eixP+z− associated with Bjorken-x is retained throughout. In the small-x limit,\\\\r\\\\nour results reproduce the known sub-eikonal expressions obtained in the Color Glass Condensate\\\\r\\\\nframework, establishing a direct connection between the general-x TMD expansion and high-energy\\\\r\\\\nfactorization. We further reduce the operator basis using equations of motion, minimizing the num\\\\r\\\\nber of independent nonperturbative matrix elements entering the cross section. This work provides\\\\r\\\\na systematic foundation for extending TMD analyses of dijet production beyond leading twist, es\\\\r\\\\ntablishing a unified operator framework valid at arbitrary Bjorken-x that smoothly interpolates\\\\r\\\\nbetween moderate- and small-x descriptions of gluon TMDs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Insights into Slip-Rate Time Functions, Rupture Parameter Correlations, and Ground Motions from Validated Multicycle Earthquake Ruptures

Earthquake strong-motion predictions using kinematic source modeling require knowledge of the slip-rate functions (SRFs) along the rupture and their distinct characteristics in asperities, background (off-asperity) areas and near the surface. Here, in this study we analyzed SRFs from well-validated, self-consistent, and fully dynamic rupture models from earthquake cycles obeying a rate-and-state friction law, from our companion study (Galvez et al., 2021). The shapes of SRFs in asperities are well described by the regularized Yoffe function (RYF), which has only two parameters: rise time T r and smoothing time T s , which control the generation of long- and short-period ground motions, respectively. In background areas, we demonstrate that, in addition to the primary rupture, multiple secondary ruptures may also nucleate from rupture heterogeneities related to asperities, resulting in SRFs with multiple peaks. Because it is impossible to fit a multiple-peak SRF by the single-peak RYF, we describe SRFs in background areas in an effective way by fitting their amplitude spectra with the RYF spectra. Such spectrally effective RYFs capture salient aspects of seismic-wave generation and can be used in rupture generators for strong motion prediction. We found that small T s values correlate with small characteristic weakening distances, large peak slip rates (PSRs), and large rupture velocities. T r values are larger in background areas and smaller in asperities. Within the shallow aseismic zone, T s values approximately quadruple whereas T r values approximately double. Because of this dominant T s increase, PSR values decrease in the near-surface zone. These features indicate that the generation of strong motions by the near-surface portions of the rupture is negligible in the studied scenarios.

Geosciences↗