Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “velocity model building”

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

Designing for Supply and Return Air System Interaction in Residential Buildings

Standard practice for HVAC return design has evolved from running a dedicated return to each room with a supply, to systems with returns in more centrally located areas of the home with jump ducts, transfer grilles, or simply door undercuts used as return air pathways from isolated rooms . In these latter systems, hallways and stairwells act as large open ducts for conveying air back to a single (or sometimes multiple) central return. When partition doors to an isolated room are closed, the return airflow resistance goes up, significantly impacting airflow balance if an adequate relief pathway is not installed. Airflow imbalances can lead to comfort and building durability issues and increased envelope leakage. To combat this, some jurisdictions have requirements for return air pathways. The supply system topology – or layout - impacts the airflow balance stability in response to adjustments of return pathway resistances. Branching supply topologies typically have reduced static pressure after each split. The static pressure at the final split will be lower than the primary supply plenum. Because of this, if there is a restriction in a room’s return path, supply airflow will tend to redistribute to adjacent ducts at the end of the branch. A properly designed trunk and branch supply system can effectively equalize static pressure by reducing the cross sectional area after each takeoff . Maintaining static pressure within the trunk will reduce the system’s sensitivity to changes in return paths. However, in practice, it is difficult to design and time consuming to install a complex supply plenum. Velocity effects and poor takeoff placement also impact airflow balance. This fact sheet considers return systems in three main categories: distributed, with a return duct to each room; multiple central, with one return grille on each floor of a home; and single central, with a single return grille located near the air handling unit. Three supply categories are also considered: radial splitter box, trunk and branch, and home-run with all ducts connecting directly to a central manifold. A complete description of the modeling work and results can be found in the companion technical report.

buildings↗
DOE OSTI2020

Can assembly bias explain the lensing amplitude of the BOSS CMASS sample in a Planck cosmology?

ABSTRACT In this paper, we investigate whether galaxy assembly bias can reconcile the 20–40 ${{\ \rm per\ cent}}$ disagreement between the observed galaxy projected clustering signal and the galaxy–galaxy lensing signal in the Baryon Oscillation Spectroscopic Survey CMASS galaxy sample. We use the suite of abacuscosmos lambda cold dark matter simulations at Planck best-fitting cosmology and two flexible implementations of extended halo occupation distribution (HOD) models that incorporate galaxy assembly bias to build forward models and produce joint fits of the observed galaxy clustering signal and the galaxy–galaxy lensing signal. We find that our models using the standard HODs without any assembly bias generalizations continue to show a 20–40 ${{\ \rm per\ cent}}$ overprediction of the observed galaxy–galaxy lensing signal. We find that our implementations of galaxy assembly bias do not reconcile the two measurements at Planck best-fitting cosmology. In fact, despite incorporating galaxy assembly bias, the satellite distribution parameter, and the satellite velocity bias parameter into our extended HOD model, our fits still strongly suggest a $\sim \! 34{{\ \rm per\ cent}}$ discrepancy between the observed projected clustering and galaxy–galaxy lensing measurements. It remains to be seen whether a combination of other galaxy assembly bias models, alternative cosmological parameters, or baryonic effects can explain the amplitude difference between the two signals.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A non‐intrusive domain‐decomposition model reduction method for linear steady‐state partial differential equations with random coefficients

Abstract Domain decomposition methods have been proved to be an effective strategy to reduce the dimension of parametric partial differential equations (PDEs). However, existing domain decomposition methods for parametric PDEs are usually intrusive, which means domain decomposition based solvers need to be implemented from scratch for each target parametric PDE. To address this issue, we develop a new non‐intrusive domain‐decomposition model reduction method for linear steady‐state PDEs with random‐field coefficients. As a variant of our previous work by Mu and Zhang, the new method only needs access to the final linear system, that is, the global stiffness matrix and the right hand side, of a deterministic PDE solver, in order to build a domain‐decomposition‐based reduced model without intrusive implementation from scratch. The key idea is to remove the interface condition between sub‐domains and rely on the correlation between columns of the linear system to couple the sub‐domains. The non‐intrusive feature enables the applicability of the proposed method to a broader class of uncertainty quantification problems, where many legacy codes/solvers can be fully reused by our method. Two numerical examples including diffusion equations with random diffusivity and convection‐dominated transport with random velocity, are provided to demonstrate the effectiveness and efficiency of our method.

Zhang, Guannan↗
DOE OSTI2021

Understanding Type Ia Supernova Distance Biases by Simulating Spectral Variations

In the next decade, transient searches from the Vera C. Rubin Observatory and the Nancy Grace Roman Space Telescope will increase the sample of known Type Ia supernovae (SNe Ia) from ~10 3 to 10 5 . With this reduction of statistical uncertainties on cosmological measurements, new methods are needed to reduce systematic uncertainties. Here, characterizing the underlying spectroscopic evolution of SN Ia remains a major systematic uncertainty in current cosmological analyses, motivating a new simulation tool for the next era of SN Ia cosmology: Build Your Own Spectral Energy Distribution (byosed). byosed is used within the SNANA framework to simulate light curves by applying spectral variations to model SEDs, enabling flexible testing of possible systematic shifts in SN Ia distance measurements. We test the framework by comparing a nominal Roman SN Ia survey simulation using a baseline SED model to simulations using SEDs perturbed with byosed, and investigating the impact of ignoring specific SED features in the analysis. These features include semiempirical models of two possible, predicted relationships: between SN ejecta velocity and light-curve observables, and a redshift-dependent relationship between SN Hubble residuals and host-galaxy mass. We analyze each byosed simulation using the SALT2 and BEAMS with Bias Corrections framework, and estimate changes in the measured value of the dark-energy equation-of-state parameter, w. We find a difference of Δw = –0.023 for SN velocity and Δw = 0.021 for redshift-evolving host mass when compared to simulations without these features. By using byosed for SN Ia cosmology simulations, future analyses (e.g., the Rubin and Roman SN Ia samples) will have greater flexibility to constrain or reduce such SN Ia modeling uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A Numerical Model for the Dynamics of Pickup Ions Outside the Heliopause and IBEX “Ribbon” Observation

The leading mechanism for the origin of the Interstellar Boundary Explorer (IBEX) “ribbon” of enhanced energetic neutral atoms (ENAs) from the outer heliosphere is the so-called secondary ENA process. In this mechanism, neutralized ions escape the heliosphere, and after several times of charge exchange collision,s some propagate back toward Earth in directions nearly perpendicular to the local interstellar magnetic field. However, the physical processes governing the distribution of the pickup ions (PUIs) outside the heliopause are still unclear. In this study, we build a new global model of the IBEX ribbon, where a key component is to calculate the dynamics of the PUIs outside the heliopause by solving the gyrophase-averaged focused transport equations on top of the background magnetohydrodynamic-kinetic model for the heliosphere-interstellar medium interaction. We discuss how the properties of the simulated ribbon change with different scattering parameters and show simulation results of the ENA sky maps, ribbon centers, ribbon profiles at several azimuthal slices, the ribbon ENA source region, and the velocity distribution of PUIs in the ENA source region for different scattering frequencies. Our results show that the model can provide reasonable ribbon sky maps comparable to the observed ribbon when the PUI mean free path is on the order of several thousand astronomical units.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1977

Neutrally stable atmospheric flow over a two-dimensional rectangular block

The phenomena of atmospheric flow over a two dimensional surface obstruction such as a building modeled as a rectangular block are analyzed by an approach using the Navier-Stokes equations with a two equation model of turbulence. The partial differential equations for the vorticity, stream function, turbulence kinetic energy, and turbulence length scale are solved by a finite difference technique. The predicted results are in agreement with the limited experimental data available. Current computed results show that the separation bubble originates from the upper front corner of the block and extends approximately 11.5 block heights behind the block. The decay of the mean velocity along the wake center line coincides almost perfectly with the experimental data. The vertical profiles of the mean velocity defect are also in reasonable agreement with wind tunnel results. Velocity profiles in the mixing region are shown to agree with the error function profile typically found in the shear layer. Details of the behavior of the turbulence kinetic energy and the turbulence length scale are also discussed.

Shieh, C. F.↗
NASA NTRS1987

Preflare energy build-up in a filament circuit

The May 16, 1981 flare observed at Debrecen is studied by extending to a fully three-dimensional model the two-dimensional Van Tend and Kuperus (1978) scenario for preflare energy build-up. It is shown that there are 10 to the 33rd ergs of free energy available to explain the subsequent large two-ribbon flare. As a result of the three-dimensional character of the present model, this estimate is an order of magnitude larger than that made by Van Tend. It is confirmed that the global form of the preflare circuit is highly important in determining the amount of energy stored in the preflare configuration. The present model gives correct predictions for the independently observed photospheric flow velocity and current strength in filaments.

Martens, P. C. H.↗
NASA NTRS2018

NASA ACC High Energy Dynamic Impact Methodology and Outcomes

For High Energy Dynamic Impact (HEDI) events, testing to evaluate the structural response of primary aircraft structure for design and certification is both expensive and time consuming. This paper discusses current work seeking to assess, develop, and validate appropriate analytical models that accurately predict physical response, damage, and failure modes for large scale composite structures in HEDI events. Four state-of-the-art Progressive Damage Analysis (PDA) methods were employed for this phased project: LS-DYNA MAT162, LS-DYNA MAT261, Smoothed Particle Galerkin (SPG), and EMU Peridynamics. Multiple material systems were considered, namely T700/5208 textile-infusion triaxial braid, T800/AMD-825 textile-infusion triaxial braid, IM7/8552 uni-directional tape, and SPG 196-PW/8552 plain-weave fabric. Extensive ballistic impact testing was performed to support this activity and measured results were compared to predictive models for assessment using panel delamination, panel displacement, force at the load cells, and threshold velocity as measures. Ultimately, the work under this activity provided significant progress in advancing the state-of-the-art in the use of PDA for HEDI events. Each material model had favorable performance comparing to test in some parameters and needed improvement in others. With the lessons learned from this activity, significant progress was made in the ability to predict panel behavior for a more general case beyond the flat panel in a ballistic impact event. Subsequent Phase II of the NASA ACC HEDI effort will continue to build on the coupon testing, flat panel ballistic impact testing, and analysis performed to-date with application of the PDA methods for intended material selections to test articles with greater complexity of configuration, curvature, and scale. It is not the intention of this paper to present a full set of data, but rather to give an overview of the NASA HEDI effort and show a small representative subset of the test and analysis results.

Hunziker, Kenneth J.↗
DOE OSTI2021

Three dimensional low-mode areal-density non-uniformities in indirect-drive implosions at the National Ignition Facility

To achieve hotspot ignition, an inertial confinement fusion implosion must achieve high hotspot pressure that is inertially confined by a dense shell of DT fuel. This requires a symmetric implosion having high in-flight shell velocity and high areal density at stagnation. The size of the driver and scale of the capsule required can be minimized by maintaining a high efficiency of energy coupling from the imploding shell to the hotspot. Significant 3D low mode asymmetries, however, are commonly observed in indirect-drive implosions and reduce the coupling of shell kinetic energy to the hotspot. To better quantify the magnitudes and impacts of shell density asymmetries, we have developed new analysis techniques and analytic models [Hurricane et al., Phys. Plasmas 27(6), 062704 (2020)]. To build confidence in the underlying data, we have also developed an analytic neutron transport model to cross-compare two independent measurements of asymmetry, which shows excellent agreement across shots for mode-1 (ℓ = 1). This work also demonstrates that asymmetry can introduce potential sampling bias into down-scattered ratio measurements causing the solid-angle-average and uncertainty-weighted-average down-scattered ratios to differ significantly. Diagnosing asymmetries beyond mode-1 (ℓ > 1) presents significant challenges. Using new diagnostic instruments and analysis techniques, however, evidence of significant Legendre mode P2 (ℓ = 2, m = 0) and additional 3D asymmetries (ℓ > 1, m ≠ 0) are beginning to emerge from the high precision activation diagnostic data (real-time nuclear activation detectors) and down-scattered neutron imaging data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Estimating Lossy Compressibility of Scientific Data Using Deep Neural Networks

Simulation based scientific applications generate increasingly large amounts of data on high-performance computing (HPC) systems. To allow data to be stored and analyzed efficiently, data compression is often utilized to reduce the volume and velocity of data. However, a question often raised by domain scientists is the level of compression that can be expected so that they can make more informed decisions, balancing between accuracy and performance. In this letter, we propose a deep neural network based approach for estimating the compressibility of scientific data. To train the neural network, we build both general features as well as compressor-specific features so that the characteristics of both data and lossy compressors are captured in training. Our approach is demonstrated to outperform a prior analytical model as well as a sampling based approach in the case of a biased estimation, i.e., for SZ. However, for the unbiased estimation (i.e., ZFP), the sampling based approach yields the best accuracy, despite the high overhead involved in sampling the target dataset.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Physics-constrained, low-dimensional models for magnetohydrodynamics: First-principles and data-driven approaches

Plasmas are highly nonlinear and multiscale, motivating a hierarchy of models to understand and describe their behavior. However, there is a scarcity of plasma models of lower fidelity than magnetohydrodynamics (MHD), although these reduced models hold promise for understanding key physical mechanisms, efficient computation, and real-time optimization and control. Galerkin models, obtained by projection of the MHD equations onto a truncated modal basis, and data-driven models, obtained by modern machine learning and system identification, can furnish this gap in the lower levels of the model hierarchy. This work develops a reduced-order modeling framework for compressible plasmas, leveraging decades of progress in projection-based and data-driven modeling of fluids. We begin by formalizing projection-based model reduction for nonlinear MHD systems. To avoid separate modal decompositions for the magnetic, velocity, and pressure fields, we introduce an energy inner product to synthesize all of the fields into a dimensionally consistent, reduced-order basis. Next, we obtain an analytic model by Galerkin projection of the Hall-MHD equations onto these modes. We illustrate how global conservation laws constrain the model parameters, revealing symmetries that can be enforced in data-driven models, directly connecting these models to the underlying physics. We demonstrate the effectiveness of this approach on data from high-fidelity numerical simulations of a three-dimensional spheromak experiment. Finally, this manuscript builds a bridge to the extensive Galerkin literature in fluid mechanics and facilitates future principled development of projection-based and data-driven models for plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

JHTDB-wind: a web-accessible large-eddy simulation database of a wind farm with virtual sensor querying

This paper introduces JHTDB-wind (https://turbulence.idies.jhu.edu/datasets/windfarms, last access: 11 November 2025), a publicly accessible database containing large-eddy simulation (LES) data from wind farms. Building on the framework of the Johns Hopkins Turbulence Database (JHTDB), which hosts direct numerical simulation (DNS) and some LES datasets of canonical turbulent flows, JHTDB-wind stores the 4D space–time history of the flow and provides users the ability to access and query the data via a web-based virtual sensor interface. The initial dataset comprises LES results from a large wind farm with 10×6 turbines, modeled using a filtered actuator line method, under conventionally neutral atmospheric conditions. These data comprise 1 h (hour) of flow field data (velocity, pressure, potential temperature deviation, subgrid-scale (SGS) eddy viscosity, and turbine forces, approximately 15 TB (terabytes) and wind turbine data – including both turbine-level operational quantities and blade-level aerodynamic quantities (approximately 1.3 TB) – stored in Zarr and Parquet formats, respectively. Data retrieval is facilitated by the giverny Python package, allowing remote users to query the database in Python or MATLAB (C and Fortran support are available for flow field data). This paper details the simulation setup and demonstrates data access through examples that analyze wind farm flow structures and turbine performance. The framework is extensible to future datasets, including the JHTDB-wind diurnal cycle simulation analyzed in Xiao et al. (2025).

17 WIND ENERGY↗
NASA NTRS1995

The microscopic state of the solar wind: Links between composition, velocity distributions, waves and turbulence

An overview is given of the microscopic state of the solar wind with emphasis on recent Ulysses high-latitude observations and previous Helios in-ecliptic observations. The possible links between composition, ionization state. velocity distribution functions of electrons, protons and heavy ions. kinetic plasma waves and MHD-scale turbulence are elaborated. Emphasis is placed on a connection of interplanetary kinetic-scale phenomena with their generating microscopic processes in the corona. The fast streams seem to consist of mesoscale pressure-balanced plasma filaments and magnetic flux tubes, reminiscent of the supergranular-size structures building the open corona, from which copious Alfven waves emanate. The wind from the magnetically structured and active corona shows developed compressive turbulence and considerable abundance and ionization state variations. Some modelling attempts to explain the observed element fractionation are briefly reviewed. The causes of the nonthermal particle features. such as proton-ion differential streaming, ion beams. temperature anisotropies, and skewed distributions associated with collisionless heat conduction, are ultimately to be searched in the fact, that the corona is never quiescent but fundamentally variable in space and time. Consequently, the radial evolution of the internal state of the wind resembles at all latitudes a complicated relaxation process, in the course of which the free (in comparison with LTE conditions) particle kinetic energy is converted into plasma waves and MHD turbulence on a wide range of scales. This leads to intermittent wave-particle interactions and unsteady anomalous transport, mixed with the weak effects of the rare Coulomb collisions. Spherical expansion and large-scale inhomogeneity forces the wind to attain microscopically a complex internal state of dynamic statistical equilibrium.

Marsch, E.↗
NASA NTRS2001

Dendritic Alloy Solidification Experiment (DASE)

A space experiment, and supporting ground-based research, is proposed to study the microstructural evolution in free dendritic growth from a supercooled melt of the transparent model alloy succinonitrile-acetone (SCN-ACE). The research is relevant to equiaxed solidification of metal alloy castings. The microgravity experiment will establish a benchmark for testing of equiaxed dendritic growth theories, scaling laws, and models in the presence of purely diffusive, coupled heat and solute transport, without the complicating influences of melt convection. The specific objectives are to: determine the selection of the dendrite tip operating state, i.e. the growth velocity and tip radius, for free dendritic growth of succinonitrile-acetone alloys; determine the growth morphology and sidebranching behavior for freely grown alloy dendrites; determine the effects of the thermal/solutal interactions in the growth of an assemblage of equiaxed alloy crystals; determine the effects of melt convection on the free growth of alloy dendrites; measure the surface tension anisotropy strength of succinon itrile -acetone alloys establish a theoretical and modeling framework for the experiments. Microgravity experiments on equiaxed dendritic growth of alloy dendrites have not been performed in the past. The proposed experiment builds on the Isothermal Dendritic Growth Experiment (IDGE) of Glicksman and coworkers, which focused on the steady growth of a single crystal from pure supercooled melts (succinonitrile and pivalic acid). It also extends the Equiaxed Dendritic Solidification Experiment (EDSE) of the present investigators, which is concerned with the interactions and transients arising in the growth of an assemblage of equiaxed crystals (succinonitrile). However, these experiments with pure substances are not able to address the issues related to coupled heat and solute transport in growth of alloy dendrites.

Beckermann, C.↗
DOE OSTI2020

Simulations of Laser Experiments to Study the Origin of Cosmic Magnetic Fields

Magnetic fields are ubiquitous in the universe; however, their origin is not fully understood. Cosmologists and astrophysicists have proposed a variety of ways in which small seed magnetic fields could be created. It is widely thought that the much larger values of the cosmic magnetic fields that we observe are a result of the amplification of these seed fields by the nonlinear turbulent dynamo mechanism. Such a mechanism had not yet been demonstrated in a controlled laboratory environment. We conceived experiments designed to demonstrate and study the nonlinear turbulent dynamo mechanism in the laboratory. These experiments characterize the distribution of turbulent energy among the velocity, magnetic field, and density fluctuations, providing a comprehensive picture of the energy cascade in a magnetized, turbulent plasma. The experiments build on the pathfinder experiments we have conducted on the Vulcan laser at the Rutherford-Appleton Laboratory in the UK and the OMEGA laser at Laboratory for Laser Energetics at Rochester. They utilize the high-intensity lasers at the Omega Laser Facility, the National Ignition Facility (NIF) at LLNL and the Laser Megajoule (LMJ) Facility in France – the largest laser facilities in the world. The goal of this project was to design and model the highly demanding experiments through simulation campaigns using FLASH, a highly capable radiation-MHD code we have developed, and large-scale 3D simulations on the Mira supercomputer at ANL. The simulations were vital to ensuring the experiments achieve the strong turbulence and large magnetic Reynolds numbers required for the nonlinear turbulent dynamo mechanism to operate. The simulations were also critical to determining when to fire the diagnostics, since the experiments last tens of nanoseconds but the signals of the strongly amplified magnetic fields last only a few nanoseconds. Finally, high-fidelity, validated FLASH simulations were crucial to interpreting the results of the experiments. By combining theory, experiments, and simulations, we were able to succeed in demonstrating and characterizing the nonlinear turbulent dynamo mechanism. The effort leveraged (1) the emergence of the Flash Center for Computational Science as a leader in numerical modeling of academic High Energy Density Physics (HEDP) experiments; (2) the FLASH code, which has been developed by the Flash Center and is a leadership-class, multi- physics, publicly available community code that has been applied successfully in many laser-driven experiments; (3) the expertise acquired through pathfinder experiments we performed on the Vulcan laser at the Rutherford-Appleton Laboratory in the UK and on the Omega laser facility; (4) the experience we have gained in designing and interpreting these experiments using validated simulations done with FLASH; and (5) the close collaborations we have established with outstanding experimental groups in the US and abroad. The broader impacts of the project are two-fold: (1) the work furthered the transformation of the academic community’s ability to design and analyze HEDP experiments at large laser facilities that is happening through the availability of FLASH and its widespread adoption by the academic community – a community that previously had had limited access to open and validated hydro/MHD simulation tools for HEDP experiments; and (2) the work trained young scientists to design and interpret HEDP experiments using validated simulations – a critical national need.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2011

Kalman Orbit Optimized Loop Tracking

Under certain conditions of low signal power and/or high noise, there is insufficient signal to noise ratio (SNR) to close tracking loops with individual signals on orbiting Global Navigation Satellite System (GNSS) receivers. In addition, the processing power available from flight computers is not great enough to implement a conventional ultra-tight coupling tracking loop. This work provides a method to track GNSS signals at very low SNR without the penalty of requiring very high processor throughput to calculate the loop parameters. The Kalman Orbit-Optimized Loop (KOOL) tracking approach constitutes a filter with a dynamic model and using the aggregate of information from all tracked GNSS signals to close the tracking loop for each signal. For applications where there is not a good dynamic model, such as very low orbits where atmospheric drag models may not be adequate to achieve the required accuracy, aiding from an IMU (inertial measurement unit) or other sensor will be added. The KOOL approach is based on research JPL has done to allow signal recovery from weak and scintillating signals observed during the use of GPS signals for limb sounding of the Earth s atmosphere. That approach uses the onboard PVT (position, velocity, time) solution to generate predictions for the range, range rate, and acceleration of the low-SNR signal. The low- SNR signal data are captured by a directed open loop. KOOL builds on the previous open loop tracking by including feedback and observable generation from the weak-signal channels so that the MSR receiver will continue to track and provide PVT, range, and Doppler data, even when all channels have low SNR.

Young, Lawrence E.↗
NASA NTRS1976

Turbulent atmospheric flow over a backward-facing step

The phenomenon of atmospheric shear layer separation over a man-made structure such as a building (modeled as a backward-facing step) has been analyzed theoretically by (1) solving the two-dimensional equations of motion in the two variables, stream function and vorticity, and by (2) employing an approximate integral technique. Boundary conditions for the undisturbed flow are that of the turbulent atmospheric shear flow over a rough terrain. In the first approach a two-equation model of turbulence was used. In the second approach an approximate technique was utilized in an attempt to describe the details of the flow in the recirculation zone behind the step. The results predict velocity profiles in sufficient detail that the presence of the corner eddy in the region of negative surface pressure gradient is evident. The magnitude of the reversed flow velocity in the recirculation eddy has been found to agree with that found from experiments. Also, a surface eddy viscosity distribution has been an outgrowth of the method which realistically follows the magnitude of the surface pressure gradient distribution as found experimentally.

Kaul, U. K.↗
DOE OSTI2020

Strong Lensing Model of SPT-CL J0356–5337, a Major Merger Candidate at Redshift 1.0359

We present an analysis of the mass distribution inferred from strong lensing by SPT-CL J0356-5337, a cluster of galaxies at redshift $z=1.0359$ revealed in the follow-up of the SPT-SZ clusters. The cluster has an Einstein radius of ${\theta }_{{\rm{E}}}\simeq $ 14'' for a source at z = 3 and a mass within 500 kpc of ${M}_{500\mathrm{kpc}}=4.0\pm 0.8\times {10}^{14}\,$ ${M}_{\odot }$ . Our spectroscopic identification of three multiply imaged systems ($z=2.363$, $z=2.364$, and $z=3.048$), combined with HSTF606W-band imaging allows us to build a strong lensing model for this cluster with an rms of $\leqslant 0\buildrel{\prime\prime}\over{.} 3$. Our modeling reveals a two-component mass distribution in the cluster. One mass component is dominated by the brightest cluster Galaxy and the other component, separated by ~170 kpc, contains a group of eight red elliptical galaxies confined in a ~9'' (~70 kpc) diameter circle. We estimate the mass ratio between the two components to be between 1:1.25 and 1:1.58. In addition, spectroscopic data reveal that these two near-equal mass cores have only a small velocity difference of ~300 km s-1 between the two components. This small radial velocity difference suggests that most of the relative velocity takes place in the plane of the sky, and implies that SPT-CL J0356-5337 is a major merger with a small impact parameter seen face-on. We also assess the relative contributions of Galaxy-scale halos to the overall mass of the core of the cluster and find that within 800 kpc from the brightest cluster Galaxy about 27% of the total mass can be attributed to visible and dark matter associated with galaxies, whereas only 73% of the total mass in the core comes from cluster-scale dark matter halos.

79 ASTRONOMY AND ASTROPHYSICS↗