Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “civil engineering”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Conceptual Design Report for the MATHUSLA Long-Lived Particle Detector near CMS

We present the Conceptual Design Report (CDR) for the MATHUSLA (MAssive Timing Hodoscope for Ultra-Stable neutraL pArticles) long-lived particle detector at the HL-LHC, covering the design, fabrication and installation at CERN Point 5. MATHUSLA is a 40 m-scale detector with an air-filled decay volume that is instrumented with scintillator tracking detectors, to be located near CMS. Its large size, close proximity to the CMS interaction point and about 100 m of rock shielding from HL-LHC backgrounds allows it to detect LLP production rates and lifetimes that are one to two orders of magnitude beyond the ultimate sensitivity of the HL-LHC main detectors for many highly motivated LLP signals. Data taking is projected to commence with the start of HL-LHC operations. We present a new 40m design for the detector: its individual scintillator bars and wavelength-shifting fibers, their organization into tracking layers, tracking modules, tower modules and the veto detector; define a high-level design for the supporting electronics, DAQ and trigger system, including supplying a hardware trigger signal to CMS to record the LLP production event; outline computing systems, civil engineering and safety considerations; and present preliminary cost estimates and timelines for the project. We also conduct detailed simulation studies of the important cosmic ray and HL-LHC muon backgrounds, implementing full track/vertex reconstruction and background rejection, to ultimately demonstrate high signal efficiency and $\ll 1$ background event in realistic LLP searches for the main physics targets at MATHUSLA. This sensitivity is robust with respect to detector design or background simulation details. Appendices provide various supplemental information.

Aitken, Branden [Victoria U.]↗
DOE OSTI2024

Dynamic Wind Loading on CSP Collectors Caused by Turbulent Wind Fluctuations: Insights from a 2-Year Field Campaign

Concentrating Solar Power (CSP) is a promising solar technology for electricity generation with thermal energy storage and with the additional benefit of industrial heat production. Wind loading on CSP collector structures, such as parabolic troughs or heliostats, is one of the primary drivers of their structural design costs. In particular, dynamic wind loading is a major source of uncertainty in the collector design process, which heavily relies on wind tunnel testing. In the field, the turbulent nature of the incoming wind creates fluctuating loads (support structure loads and resulting mirror deflections) on the collectors, with impacts on fatigue lifetime and optical performance. As is well known, wind tunnel tests cannot entirely reproduce the complex turbulent wind conditions typically observed at full-scale plants. To shed light on this topic, NREL initiated a field campaign at the operational Nevada Solar One (NSO) powerplant that uses parabolic troughs as solar collectors. The aim of the project is a detailed characterization of prevailing wind and turbulence conditions and resulting operational loads on parabolic troughs. We use the published 2-year dataset of high-resolution combined wind and structural loads measurements [1] to characterize the dynamic structural wind response. For quantifying dynamic wind loading, we apply the concept of admittance functions, which are spectral transfer functions that couple the turbulent wind to resulting structural loads (aerodynamic admittance), and to the structural response (mechanical admittance). In practice, aerodynamic admittance describes which turbulent eddy sizes are effective in creating structural loads. The mechanical admittance describes in which frequency ranges these loads are reinforced or dampened by the structure. While these functions are an established concept in civil engineering, their recent application to a single full-scale heliostat [2] proved their broader applicability to CSP collectors. Here, we present a characterization of admittance functions for full-scale parabolic trough collectors and show how wind characteristics (mean wind speed and direction, turbulent kinetic energy, turbulent length scales), the sun-tracking trough angle, and row position alter the admittance functions. Further, we study to which extent the admittance functions are universal for a specific trough geometry and how our findings compare to reported heliostat results. References [1] https://data.openei.org/submissions/5938. [2] Blume, K., Roger, M., and Pitz-Paal, R. 2023b. "Simplified analytical model to describe wind loads and wind-induced tracking deviations of heliostats." Solar Energy, 256, 96-109. https://doi.org/10.1016/j.solener.2023.03.055.

admittance functions↗
DOE OSTI2025

Introduction to Digital Image Correlation (DIC) with annotated bibliography

Digital Image Correlation (DIC) is “a non-contact means of measuring motion and deformation using digital images of the object of interest”, capable of full-field measurements over large areas. Getting started in DIC can be difficult, since there are many different techniques, and the body of literature on DIC is extensive. IEEE Xplore alone lists nearly 11,000 references, including magazines, conference proceedings, journal articles, books, etc. A Google search returns around 145 million results. The intent of this annotated bibliography is to provide an entry point for beginners and to collect references that may be useful to more advanced practitioners.

42 ENGINEERING↗
DOE OSTI2026

1016516285-AA - Handling Frame Assembly Table Structural and Seismic Analysis

Within the NIF (National Ignition Facility) sustainment project, LLNL is replacing all 1,728 blast shields with new ones. Before the new blast shield glass can be installed on NIF, it must pass through the AMOL (Amplifier Main Optics Loop). The AMOL process involves cleaning and coating the glass, assembling it into a Line Replaceable Unit (LRU), and performing testing and inspection. This work will be conducted in the OPFX cleanroom facility (Building 391, Room 1250) and requires several pieces of custom equipment. One key piece of equipment within the loop is the Handling Frame Assembly Table (see Figure 1.1). Its purpose is to assist operators in installing the blast shield glass into handling frames (1014442709) (see Figure 1.2). The table is also used to store handling frame assembly components. The primary components stored are the L-frames (1015208271), which are kept in a drawer underneath the table (see Figure 1.3). The purpose of this safety note is to ensure that the system does not pose a risk of injury to personnel during a seismic event. The system owner is OMST Engineering & Maintenance.

42 ENGINEERING↗
DOE OSTI2026

Design and Development of the Handling Frame Assembly Table for the National Ignition Facility Blast Shields

The Handling Frame Assembly Table is a custom designed piece of equipment developed to support the NIF Sustainment Project. The goal of the NIF Sustainment Project is to refurbish the National Ignition Facility (NIF), the world’s largest and most energetic laser, so it can continue advancing research in fusion energy, national security, and high energy density physics. A key aspect of this effort involves replacing all 1,728 blast shields on the NIF, which requires specialized equipment such as the Handling Frame Assembly Table to support the production of new blast shields. This project follows a systems engineering approach that guides the design sequence. This report documents the design process from initial concept through final design, including stakeholder analysis, requirements development, conceptual design, final design, and design validation. Future work will focus on building and commissioning the system.

42 ENGINEERING↗
DOE OSTI2025

Selective mass scaling for single-layer thick shell elements in DYNA3D

Hexahedral elements can be adapted to model thin and moderately thick structures by neglecting the coupling of through-thickness stress, resulting in a fully three-dimensional, but simplified, state of stress. These specialized elements, often referred to as “thick” or “solid” shells, are generally employed to model thin-walled structures using continuum mechanics-based material models. In explicit dynamics simulations, where computational speed is important, these elements are integrated with a single quadrature point and a set of anti-hourglassing (stabilizing) forces. Thick shells, by definition, have a thickness dimension smaller than their in-plane dimensions, and this small thickness often determines the stable time step size in simulations, despite the mechanics being approximated. To alleviate this limitation while retaining the relevant dynamics of thin-walled structures, selective mass scaling (SMS), or selective mass “augmentation,” has been proposed in the literature. In this technical report, we explore the application of SMS to single-layer thick shells in the simulation software DYNA3D.

42 ENGINEERING↗
DOE OSTI2026

Calculations Document

Figure 1 shows a simplified model schematic of the Tetracage target design. In this figure the green lines represent the Tetracage yarns, which consist of two parallel yarns on the top and two parallel yarns on the bottom with the pairs in a cross hatched pattern. For simplicity, as we are only considering vertical forces, we will perform calculations without consideration for the two pairs of yarns being perpendicular to each other. Also worth noting is the model shows the yarns extending directly to the capsule then closely hugging it – in reality, the yarns behave as you would expect where they hug the capsule only as much as the geometry of their deflection allows, similarly to if you held a tennis ball with a pair of parallel strings anchored at opposing ends. The dashed black lines show the path the yarns take straight across from slot glue point to slot glue point before the target is closed and the capsule deforms them. The variable r corresponds to the capsule outer radius. The variable L corresponds to the hohlraum inner radius. The variable a corresponds to the yarn offset from the center of the diagnostic band. The variable d corresponds to the wire deflection vertically.

42 ENGINEERING↗
DOE OSTI2026

Design and Development of the Handling Frame Assembly Table for the National Ignition Facility Blast Shields

The Handling Frame Assembly Table is a custom designed piece of equipment developed to support the NIF Sustainment Project. The goal of the NIF Sustainment Project is to refurbish the National Ignition Facility (NIF), the world’s largest and most energetic laser, so it can continue advancing research in fusion energy, national security, and high energy density physics. A key aspect of this effort involves replacing all 1,728 blast shields on the NIF, which requires specialized equipment such as the Handling Frame Assembly Table to support the production of new blast shields. This project follows a systems engineering approach that guides the design sequence. This report documents the design process from initial concept through final design, including stakeholder analysis, requirements development, conceptual design, final design, and design validation. Future work will focus on building and commissioning the system.

42 ENGINEERING↗
DOE OSTI2026

Proliferated Resilient Economical Half-Meter Aperture Space Telescopes (PREEMPT)

The PREEMPT LDRD was motivated by a major national need for better space-based imaging systems that are both high performing and affordable for the US Government. Current optical payloads for intelligence, surveillance and reconnaissance (ISR) and space domain awareness (SDA) missions can cost hundreds of millions of dollars and take years to develop, which makes it difficult to build the large constellations required for persistent, world-wide coverage. To address this, the team aimed to advance a different kind of large-aperture (>25 cm) telescope, called a monolithic telescope, in which key optical surfaces are built into a single piece of fused silica. This design greatly reduces payload and spacecraft complexity, the need for precision focus actuators, improves mechanical and thermal robustness, and lowers cost when compared with traditional Cassegrain telescopes that rely on many precisely aligned components. The project focused on five primary thrusts. The first thrust was to advance the concept of a V10, 25 cm, monolithic telescope forward from optical design to flight-ready stage. This was accomplished in partnership with Optimax, who delivered the first test unit in the early stages of the LDRD. The team developed several technologies necessary for this optic to be integrated into a flight demonstration. These include carbon fiber housings, highly detailed structural and thermal models and stress-reducing elastic averaging Hirth groove designs. These technologies resulted in a successful maturation of the optic, which is now slated to fly in late 2026/early 2027 for a demonstration mission. The second thrust was to advance the manufacturability of these optics. In collaboration with NIF’s optical manufacturing shop, we reduced polishing time from 480 hours to 65 hours through the implementation of optimized processes and new tools. The NIF team utilized a conceptual V8 (18 cm) optic to demonstrate this optimization, though it can be applied to the rest of the monolithic optic portfolio. Third, the team developed the first conceptual V20 (50 cm) payload, which is slated to be the next generation of LLNL optical payload systems. A set of structural, dynamic and thermal simulations were performed to identify potential challenges in the future development of this payload. Early-stage simulations suggest the payload is feasible, though thermal management will be the key focus area to maintain optimal performance. Fourth, a non-linear model of Viton was developed, to further enhance the reliability of our structural and dynamic models for future payloads. Viton acts as the primary interface material between the optic and its housing. Lastly, the team focused on successfully displaying the feasibility of using additively manufactured metal composites for optical space payloads. The team successfully demonstrated layer by layer deposition of Al-SiC composites, which have highly tunable structural and coefficient of thermal expansion (CTE) properties. These are crucial for optical payloads because CTE mismatch is one of the causes for degraded optical performance for telescopes in orbit. Overall, the work showed that monolithic telescopes could become a practical, lower-cost path to high-resolution space imaging for both national security and scientific missions.

42 ENGINEERING↗
DOE OSTI2025

Impact of Electromagnetic Fields on Gypsum and Silica Scaling in Reverse Osmosis

Electromagnetic field (EMF) is a cost-effective, simple, and energy-efficient method for scale control in reverse osmosis (RO) systems. However, its effects on gypsum and silica scaling, as well as the underlying mechanisms, remain poorly understood. This study systematically investigates the effects of EMF treatment on gypsum and silica scaling in RO systems, utilizing synthetic brackish water and natural RO concentrate (ROC) from a desalination facility. For gypsum, EMF changes the crystal morphology, resulting in the formation of a porous, less compact scaling layer. It is more readily removed through hydraulic flushing (HF), enhancing scaling reversibility and water recovery. In the case of silica scaling, EMF promotes homogeneous polymerization in the bulk solution, producing larger silica particles that inhibit the formation of a dense, cross-linked gel layer on the membrane surface, mitigating flux decline. This study thus demonstrates EMF’s effectiveness in controlling gypsum scaling in undersaturated feedwaters when combined with HF and in mitigating silica scaling under both HF and non-HF conditions for supersaturated feedwaters. These findings underscore EMF’s versatility as a nonchemical approach for scale control in RO desalination and show its substantial potential to enhance membrane performance and operational efficiency in real-world water treatment applications.

Du, Xuewei [New Mexico State University, Las Cruce↗
DOE OSTI2026

Data for Greenhouse Gas Accounting Procedures in Low Carbon Fuel Policies Overlook the Spatial Variability of Miscanthus-Derived Sustainable Aviation Fuel

Low carbon fuel policies such as the U.S. Renewable Fuel Standard (RFS), Canada Clean Fuel Regulations (CFR), and California Low Carbon Fuel Standard (LCFS) as well as the 45Z tax credit are intended to reduce greenhouse gas (GHG) emissions from transportation. Cellulosic feedstocks, optimized biorefineries, and favorable farming locations can significantly reduce biofuel carbon intensity (CI). Despite advances in field-to-fuel GHG monitoring and flexibility in resource allocation within biorefineries (e.g., governing net electricity production), rigid CI accounting procedures in current policies may limit CI responsiveness across candidate sites and processing facilities. This work examines a hypothetical biomass-to-sustainable aviation fuel (SAF) pathway using miscanthus and alcohol-to-jet (i) to demonstrate how GHG accounting requirements drive estimates of biofuel CIs and (ii) to explore potential CI and financial implications of scenario-specific life cycle assessment (LCA). Results demonstrate that GHG accounting using the CFR/LCFS can reasonably account for distinct levels of net electricity production by a biorefinery, but only the CFR yields similar CI sensitivity to spatially explicit factors (feedstock CI, grid electricity CI) as scenario-specific LCA: most GHG accounting frameworks do not capture CI variation across candidate sites in the United States. Ultimately, this work demonstrates the importance of LCA methodological specifications in low carbon fuel policies and tax credits.

Miscanthus↗
DOE OSTI2026

Characteristics of Fluid‐Solid Interaction Constitutive Models Within Poroelastodynamics at Higher Strain‐Rates and Large Deformations Implemented in 1D

The large deformation, mixed formulation, finite element (FE) modeling approach presented in Irwin et al. 2024 is extended herein to include improved constitutive models for representing dynamic solid-fluid interactions at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger overpressure magnitudes (𝒪⁡(1⁢0 2 )⁢kPa) within a biphasic soft porous material using Theory of Porous Media (TPM) at finite strain. Specifically, these constitutive modeling improvements are the following: (i) a more physically robust constitutive model for pore fluid seepage velocity via inclusion of pore fluid viscous stress, and (ii) a modified deformation-dependent-permeability model and updated hyperelastic constitutive model better suited for handling larger volumetric compressions and extensions. The novelty of the present work is mainly the contribution (i): inclusion of pore fluid viscous stress at higher strain-rate and large deformations, which requires 𝐶 1 continuity in the weak formulation, accomplished by employing Hermite cubic interpolation functions within a mixed nonlinear poromechanical finite element formulation. In (ii), the model is updated to weakly enforce solid phase incompressibility, such that this assumption is not violated numerically, which provides improved numerical stability for achieving larger overpressure magnitudes on 𝒪⁡(1⁢0 2 ) kPa, which were not achievable with the previous Kozeny–Carman model in Irwin et al. 2024. Also in (ii), the volumetric part of the solid skeleton free energy function is modified to ensure proper bounds on the solid skeleton Jacobian of deformation 𝐽 s related to incompressibility of the solid phase. Uniaxial strain, unidirectional flow examples at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger deformations (up to 0.2 (or 20%) nominal axial strain) demonstrate the improved physical representation—and numerical stability—of these constitutive model improvements.

42 ENGINEERING↗
DOE OSTI2026

Critical Review of Brazil Disk Techniques for Tensile Strength Characterization With an Emphasis on High Explosive Materials

Mechanical properties are a critical performance metric for many high explosive (HE) materials and tensile strength properties are particularly important. Direct tensile measurements using dogbone shaped samples are the gold standard but they have the disadvantage that they are fairly large and require samples machined from billets. Diametral compression, more commonly known as Brazil disk (BD) testing, is an indirect method for measuring tensile strength on smaller and more easily fabricated samples. A review of the BD literature is presented with an emphasis on tensile strength measurements in high explosive materials. BD literature is reviewed in three primary areas: (i) rocks and concrete, (ii) pharmaceutical materials, and (iii) high explosive materials. The literature for rocks/concrete is extensive and dates back over 80 years; despite this there is no consensus on the validity/accuracy of the BD technique or the optimal variant of the BD technique to employ. The pharmaceutical literature is the opposite, being limited in scope and quantity of studies. BD literature on high explosive materials falls in between, not as impressive as in the rocks/concrete community but more substantiative than in the pharmaceutical community. After the review of the literature practical parameters for HE BD testing and recommended future work is discussed.

Brazil disk↗
DOE OSTI2026

Nonlinear viscoelastic response of silicone additively manufactured direct ink write (DIW) foams under repetitive compression

To investigate dynamic fatigue behavior of foam in military protective applications, such as helmets, additively manufactured (AM) foams were compressively strained into the plateau region using a reduced design of experiments. A simple power law was found to govern the decline in dynamic stiffness (complex modulus) as the foams underwent the purchase order requirement of 10,000 cycles of small deformation in the plateau region. This rate of decline was newly found to correlate with the degree of nonlinearity in the material’s deformation, quantified using total harmonic distortion. Materials with low nonlinearity exhibited relatively stable stiffness across cycles, while those with high nonlinearity experienced greater losses. The observed nonlinearity depended on both applied stress and strain rate. A strong linear correlation (R 2 = 0.78) was identified between second-order nonlinearity and the time-dependent stiffness response. Two lattice structures were examined: face-centered tetragonal (FCT) and simple cubic (SC). The SC material exhibited higher total harmonic distortion (5%) and lower stiffness retention than the FCT (2%). In conclusion, these results suggest that for cyclic compression applications in a wide variety of industries such as packaging, personal protective equipment, or aerospace, selecting materials with lower stress and greater structural uniformity can enhance the stability of dynamic performance.

42 ENGINEERING↗
DOE OSTI2023

On the convergence of statistics in simulations of stationary incompressible turbulent flows

When reporting statistics from simulations of statistically stationary chaotic phenomenon, it is important to verify that the simulations are time-converged. This condition is connected with the statistical error or number of digits with which statistics can be reliably reported. In this work we consider numerical experiments of low Reynolds number incompressible homogeneous and isotropic turbulence as a model problem to investigate statistical convergence over finite simulation times. Specifically, we investigate the time integration requirements that allow meaningful reporting of the statistical error associated with finiteness of the temporal domain. We address two key questions: (1) How long should a simulation be performed in terms of large eddy time, and (2) How should the simulation time be divided among temporal windows over which a quantity of interest is estimated so that its statistical error could be reliably reported? We find that reliable reporting of statistical errors requires simulations on the order of 10 4 large eddy times, which is orders of magnitude longer than typically performed. Additionally, data post-processing should employ windows of at least ten times the large eddy time scale, with the most robust computation of statistical error of the mean requiring window sizes of an additional factor of ten. For practical simulations, we demonstrate that it is possible to estimate the statistical error within a factor of two under a less stringent condition in which a minimum of four windows with size at least ten large eddy times are used. In conclusion, our observations for homogeneous isotropic turbulence are also shown to hold in turbulent channel flow.

42 ENGINEERING↗
DOE OSTI2025

Adaptive Interface-PINNs (AdaI-PINNs) for transient diffusion: Applications to forward and inverse problems in heterogeneous media

We model transient diffusion in heterogeneous materials using a novel physics-informed neural networks framework (PINNs) termed Adaptive interface physics-informed neural networks or AdaI-PINNs (Roy et al. arXiv preprint arXiv:2406.04626, 2024). AdaI-PINNs utilize different activation functions with trainable slopes tailored to each material region within the computational domain, allowing for a fully automated and adaptive PINNs approach to model interface problems with strongly and weakly discontinuous solutions. To enhance its performance in highly heterogeneous transient diffusion systems, we prescribe a suite of robust practices, including appropriate non-dimensionalization of equations, a biased sampling method, Glorot initialization, and the hard enforcement of boundary and initial conditions. Here we evaluate the efficacy of the proposed method on several benchmark forward and inverse problems. Comparative studies on one-dimensional and two-dimensional benchmark problems reveal that the modified AdaI-PINNs outperform its unmodified counterpart, achieving root-mean-square errors that are at least two orders of magnitude better in forward problems. For inverse problems, the maximum errors in the approximated diffusion coefficients by modified AdaI-PINNs are four orders of magnitude better than those of the unmodified version. Additionally, modified AdaI-PINNs demonstrate improved stability in problems with large material mismatches.

42 ENGINEERING↗