Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupling method”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25

Assessing North Atlantic Tropical Cyclone Rainfall Hazard Using Engineered-Synthetic Storms and a Physics-Based Tropical Cyclone Rainfall Model

In this study, we design a statistical method to couple observations with a physics-based tropical cyclone (TC) rainfall model (TCR) and engineered-synthetic storms for assessing TC rainfall hazard. We first propose a bias-correction method to minimize the errors induced by TCR via matching the probability distribution of TCR-simulated historical TC rainfall with gauge observations. Then we assign occurrence probabilities to engineered-synthetic storms to reflect local climatology, through a resampling method that matches the probability distribution of a newly proposed storm parameter named rainfall potential (POT) in the synthetic dataset with that in the observation. POT is constructed to include several important storm parameters for TC rainfall such as TC intensity, duration, and distance and environmental humidity near landfall, and it is shown to be correlated with TCR-simulated rainfall. The proposed method has a satisfactory performance in reproducing the rainfall hazard curve in various locations in the continental United States; it is an improvement over the traditional joint probability method (JPM) for TC rainfall hazard assessment.

54 ENVIRONMENTAL SCIENCES↗

Some experiences with the viscous-inviscid interaction approach

Methods for simulating compressible viscous flow using the viscid-inviscid interaction approach are described. The formulations presented range from the more familiar full-potential/boundary-layer interaction schemes to a method for coupling Euler/Navier-Stokes and boundary-layer algorithms. An effort is made to describe the advantages and disadvantages of each formulation. Sample results are presented which illustrate the applicability of the methods.

Vandalsem, W. R.↗

High accuracy computation of fluid-structure interaction in transonic cascades

A coupling strategy for simulating fluid-structure interaction phenomena is formulated and applied to the prediction of flutter in transonic cascades. The flow is governed by the Euler equations and discretized using a finite volume flux-splitting scheme. The structure is modeled using an isoparametric finite element formulation. The coupling strategy successfully reconciles these two formulations at the fluid-structure interface by enforcing both kinematic and kinetic boundary conditions. In particular, the conservation laws applicable to the combined fluid-structure system are preserved across the interface. Since the primary mechanism driving aeroelastic phenomena involves energy exchange occurring at the interface, this highly accurate coupling mechanism is believed to improve the predictive capability of the scheme. The coupled equations are advanced simultaneously in time using an implicit time integration method. Results obtained using the coupling method are presented for cascade geometries operating in transonic flow.

Boschitsch, Alexander↗

Fundamentals of wildlife dosimetry and lessons learned from a decade of measuring external dose rates in the field

Methods for determining the radiation dose received by exposed biota require major improvements to reduce uncertainties and increase precision. We share our experiences in attempting to quantify external dose rates to free-ranging wildlife using GPS-coupled dosimetry methods. The manuscript is a primer on fundamental concepts in wildlife dosimetry in which the complexities of quantifying dose rates are highlighted, and lessons learned are presented based on research with wild boar and snakes at Fukushima, wolves at Chornobyl, and reindeer in Norway. GPS-coupled dosimeters produced empirical data to which numerical simulations of external dose using computer software were compared. Our data did not support a standing paradigm in risk analyses: Using averaged soil contaminant levels to model external dose rates conservatively overestimate the dose to individuals within a population. Following this paradigm will likely lead to misguided recommendations for risk management. The GPS-dosimetry data also demonstrated the critical importance of how modeled external dose rates are impacted by the scale at which contaminants are mapped. When contaminant mapping scales are coarse even detailed knowledge about each animal’s home range was inadequate to accurately predict external dose rates. Importantly, modeled external dose rates based on a single measurement at a trap site did not correlate to actual dose rates measured on free ranging animals. These findings provide empirical data to support published concerns about inadequate dosimetry in much of the published Chernobyl and Fukushima dose-effects research. Furthermore, our data indicate that a huge portion of that literature should be challenged, and that improper dosimetry remains a significant source of controversy in radiation dose-effect research.

61 RADIATION PROTECTION AND DOSIMETRY↗

The Impact of Ocean Data Assimilation on Seasonal-to-Interannual Forecasts: A Case Study of the 2006 El Nino Event

This study investigates the impact of four different ocean analyses on coupled forecasts of the 2006 El Nino event. Forecasts initialized in June 2006 using ocean analyses from an assimilation that uses flow-dependent background error covariances are compared with those using static error covariances that are not flow dependent. The flow-dependent error covariances reflect the error structures related to the background ENSO instability and are generated by the coupled breeding method. The ocean analyses used in this study result from the assimilation of temperature and salinity, with the salinity data available from Argo floats. Of the analyses, the one using information from the coupled bred vectors (BV) replicates the observed equatorial long wave propagation best and exhibits more warming features leading to the 2006 El Nino event. The forecasts initialized from the BV-based analysis agree best with the observations in terms of the growth of the warm anomaly through two warming phases. This better performance is related to the impact of the salinity analysis on the state evolution in the equatorial thermocline. The early warming is traced back to salinity differences in the upper ocean of the equatorial central Pacific, while the second warming, corresponding to the mature phase, is associated with the effect of the salinity assimilation on the depth of the thermocline in the western equatorial Pacific. The series of forecast experiments conducted here show that the structure of the salinity in the initial conditions is important to the forecasts of the extension of the warm pool and the evolution of the 2006 El Ni o event.

Yang, Shu-Chih↗

Returning CP-observables to the frames they belong

Optimal kinematic observables are often defined in specific frames and then approximated at the reconstruction level. We show how multi-dimensional unfolding methods allow us to reconstruct these observables in their proper rest frame and in a probabilistically faithful way. We illustrate our approach with a measurement of a CP-phase in the top Yukawa coupling. Our method makes use of key advantages of generative unfolding, but as a constructed observable it fits into standard LHC analysis frameworks.

Physics↗

Salt Sampling FY21 Technical Report

The goal of the salt sampling program at Argonne is to develop and deploy automated molten salt sampling approaches to enable high-precision in-process salt sample analysis to improve the timeliness of sampling-based accountancy measurements. Tools currently under development in support of this goal include (1) a modular vacuum sampler with an accompanying sample handling method for coupling vacuum sampling with high-precision at-line sample analysis, (2) a pneumatic sample generator that enables high-throughput sample analysis to improve the precision of existing analytical techniques, and (3) a windowless flow cell to enable on-line optical analysis of molten salt in a sampling loop. Compared to point sampling approaches (i.e., dip probes), vacuum sampling systems and on-line sampling loops facilitate access to a larger cross-section of a process fluid. This is known to improve the characterization of the process fluid by producing more representative samples and by enabling the analysis of a larger cross section of the fluid. A vacuum sampling approach for molten salts eliminates the risk of dross contamination of samples and avoids the use of moving parts in the salt. In FY21, two methods for integrating a vacuum sampler with a pneumatic sample generator were tested. These included direct fluidic coupling and coupling using a solid salt transfer mechanism. Solid salt transfer was ultimately selected over fluidic coupling, primarily to enable the transport of samples over longer distances to support automated at-line integration with high-precision techniques (such as microcalorimetry) that cannot withstand the extreme conditions near an electrorefining process. To facilitate rapid solid salt coupling, new mechanisms were developed for rapidly charging and discharging salt sample tubes at the vacuum sampler and pneumatic sample generator, respectively. While the charging mechanism will be deployed in FY22, the tube transfer method and discharge mechanism were tested in FY21. These were deployed at one of Argonne’s engineering-scale electrorefiners to implement at-line high-throughput pneumatic micro-sample generation capabilities. The method was used to generate precise uranium- and lanthanide-bearing electrorefiner micro-samples with the specific dimensions requested by researchers at Los Alamos National Laboratory for use in testing their novel microcalorimeter x-ray techniques. The solid salt transfer mechanism proved not only to be an effective means of integrating the precision sample generator with vacuum sampling, but also improved the performance of the sampler generator. Because the modular sampling approach described here eliminates the need for new high-radiation sample handling capabilities, salt-wetted seals, salt-wetted moving parts, and heated transfer lines outside the electrorefiner, it will address most of the remaining technical challenges for the at-line deployment of high-precision analytical techniques. This will enable significant reductions in the time delay for sampling-based accountancy measurements by eliminating the need for manual off-line sample processing and analysis. On-line optical analysis of molten salt in a sampling loop would provide complementary information to at-line and in-situ techniques. In FY21, an open-aperture molten salt gravity flow cell with windowless optical access to flowing salt was successfully demonstrated. Future work should include the refinement and performance testing of the on-line and at-line sampling tools, integration of additional analysis techniques, stakeholder outreach and collaboration, evaluation of the integrated methods, and analyses to determine how the various tools might fit into an integrated safeguards monitoring system of unattended near real time monitoring tools.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Imaging Inhomogeneities From Dry-Coupled Ultrasonic Scans

Method of imaging spatial distribution of selected physical properties and microstructure of specimen of material based on dry-coupled contact ultrasonic pulse/echo scanning. Ultrasonic transducer scanned across top surface of specimen. At each of many positions on two-dimensional grid on top surface, ultrasonic pulse/echo measurements taken and processed. Offers rapid, nondestructive alternative to destructive metallographic sectioning to obtain picture of inhomogeneity of specimen.

Roth, Don J.↗

An interactive three-dimensional boundary-layer method for transonic flow over swept wings

A three-dimensional laminar/turbulent boundary-layer method is developed for transonic flow over swept wings. The governing equations and curvature terms are written for a nonorthogonal body-oriented curvilinear coordinate system. The viscous method is coupled to the full-potential inviscid code, FLO-30, through a displacement-surface interaction. Typically, for transonic Mach numbers and moderate sweep angles, between 5 and 12 viscous updates are required for convergence. The method is applied to a variety of wing planforms, and the results are also compared to those obtained with the three-dimensional integral inviscid-viscous interaction method TAWFIVE (Streett and Melson, 1983).

Woodson, Shawn H.↗

A cell-centered AMR-ALE framework for 3D multi-material hydrodynamics. Part II: linesweep ALE rezoning for nonconformal block-structured AMR meshes

The simulation of flows presenting contact discontinuities, vorticity, and large variations in spatial scales can be performed in a framework coupling Arbitrary Lagrangian Eulerian (ALE) algorithms and Adaptive Mesh Refinement (AMR). This coupling requires adaptation of ALE rezoning techniques to meshes containing nonconformal nodes arising from both the AMR topology and the junction of mesh blocks. Here, in this paper, we present an ALE rezoning strategy that is compatible with such meshes, and that can also act as a disentangling algorithm. Emphasis is put on an algorithm that respects intrinsic Lagrangian mesh properties in order to preserve accuracy around discontinuities. To that end, we adapt the weighted linesweep algorithm to nonconformal block-structured AMR meshes. Then, we present control parameters introduced in the method for it to be applicable in practical situations. Notably, the method is coupled to a specific metric optimization in order to palliate some shortcomings of the linesweep method. Finally, numerical test cases are presented that feature the capabilities of the ALE-AMR algorithm for flows that present discontinuities, vorticity, and a variety of scales. Notably, we show that our ALE-AMR algorithm gives results at least similar to Euler-AMR, but provides better accuracy in cases where discontinuities are involved, thanks to a method that respects the Lagrangian features of the mesh. Additionally, it enables Euler-AMR-like computations on domains with temporally varying domain boundaries.

Adaptive mesh refinement↗

Ground State Energy Functional with Hartree–Fock Efficiency and Chemical Accuracy

We introduce the deep post Hartree–Fock (DeePHF) method, a machine learning-based scheme for constructing accurate and transferable models for the ground-state energy of electronic structure problems. DeePHF predicts the energy difference between results of highly accurate models such as the coupled cluster method and low accuracy models such as the Hartree–Fock (HF) method, using the ground-state electronic orbitals as the input. It preserves all the symmetries of the original high accuracy model. The added computational cost is less than that of the reference HF or DFT and scales linearly with respect to system size. We examine the performance of DeePHF on organic molecular systems using publicly available data sets and obtain the state-of-art performance, particularly on large data sets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

CpFe(CO) 2 Radical Generated from Dinuclear [CpFe(CO) 2 ] 2 and Mononuclear (Cp)(CO) 2 Fe(H): Density Functional Theory Is Accurate for One, But Not Both

Density functional theory (DFT) methods remain the most practical approach to calculating properties and reaction mechanisms of transition metal complexes. While the accuracy of DFT methods has been evaluated for some properties of mononuclear organometallic complexes there has been a general lack of evaluation for dinuclear organometallic complexes, in particular bonding changes related to reaction mechanisms. Here, this work evaluated DFT and coupled cluster methods for the accuracy of calculating the CpFe(CO) 2 radical (Fp•) generated from dinuclear [CpFe(CO) 2 ] 2 (Fp 2 ) and mononuclear [(Cp)(CO) 2 Fe(H)] (Fp-H). This transition metal radical fragment was evaluated because dinuclear complexes built with it have recently shown a variety of unique reactions but has proven challenging to accurately calculate with DFT methods. Here we show that DFT methods provide a surprising wide range of fragmentation energies for Fp 2 and lower and mid rung DFT methods as well as DLPNO–CCSD(T) perform well for this dissociation energy. The highest rung double-hybrid methods have a large range in the Fp 2 dissociation energy, and the energy greatly depends on the amount of MP2 correlation energy included. For generating Fp• from Fp-H the lower and mid rung methods that worked well for Fp 2 showed significant error. Double-hybrid methods unfortunately are only accurate for the Fe–H bond if they are very inaccurate for the Fp 2 dissociation energy. While DLPNO–CCSD(T) is not perfect, and not close to chemically accurate for the Fe–H bond, it does provide reasonable accuracy for both Fp 2 and Fp-H dissociation energies.

density functional theory↗

Characterization of Tank 9H Annulus Sample in Support of Residual Material Inventory Determinations

The Savannah River National Laboratory (SRNL) was requested by Savannah River Mission Completion (SRMC) to provide sample preparation and characterization of the Tank 9H annulus sample in support of Residual Material Inventory Determinations. One Tank 9H sample in three vials [HTF-9-25-13, HTF-9-25-14 and HTF-9-25-15], with each vial containing approximately 200 mL of the Tank 9H annulus salt solution, were delivered to the SRNL Shielded Cells for sample preparation and characterizations in February 2025. The density of the “as-received” solution contained in each of the three Tank 9H annulus sample vials were determined followed by a solid-liquid separation on each one using 0.45-micron Nalgene® nylon filter membranes. The resulting filtrates were combined to form the Tank 9H annulus sample with a total volume of about 600 mL. The combined wet solid fractions, about a total of 4.8 grams of salt material, remaining on the filter membranes were air-dried in the Shielded Cells for 72 hours. The total weight of the air-dried solids was 2.1 grams. These air-dried solids were washed with deionized water (DI water) at a phase ratio of 60 mL DI water/gram of solids to recover insoluble solids, if any. No visible or measurable quantity of insoluble solids were recovered after DI water washing of the air-dried solids because the air-dried solids completely dissolved in the DI water. The solid fraction-wash water was not combined with the 600 mL of the filtrate solution, and the resulting solution was not screened or analyzed for radionuclides. Aliquot sample volumes of the undiluted Tank 9H annulus sample were sent to the SRNL analytical services groups for radionuclides, elementals, anions and total mercury analysis by various methods including radiochemical separations/counting methods, inductively coupled plasma-atomic emission spectroscopy (ICP-AES), and Inductively Coupled Plasma Mass Spectroscopy (ICP-MS) and special preparations. All sample analyses were performed in triplicate. This report presents the analytical characterization results for the Tank 9H annulus sample. The results are also reported where analytical methods yielded additional analytes, other than those requested by SRMC. In the characterization of the Tank 9H annulus sample, the detection limits for all the analytes, as specified in the Technical Task Request (TTR) and Task Technical and Quality Assurance Plan (TTQAP), were met.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A fully ro-vibrationally resolved corona model for the molecular hydrogen Fulcher-𝛼 system

A ro-vibrationally resolved corona model for the molecular hydrogen Fulcher-𝛼 transition based on the Yacora solver is introduced. The model couples 1365 ro-vibrational levels of the 𝑋 1 ⁢𝛴$^+_𝑔$, 𝑑 3⁢ 𝛱 𝑢 and 𝑎 3 ⁢𝛴$^+_𝑔$ states via electron impact excitation from the ground state into the 𝑑 3⁢ 𝛱 𝑢 state and subsequent spontaneous emission into the 𝑎 3 ⁢𝛴$+_𝑔$ state. For the process of electron impact excitation a set of 45260 fully ro-vibrationally resolved cross sections calculated with the molecular convergent close-coupling (MCCC) method in the adiabatic-nuclei formulation is applied. The MCCC cross sections are compared with results of a scaling method used in other works, demonstrating the need for dedicated ro-vibrationally resolved cross sections. By post-processing the model output, entire Fulcher spectra can be simulated. These spectra are benchmarked with measurements from an inductively coupled plasma discharge at a pressure between 1.1–10 Pa and a RF power of 700 W. The model results agree very well with the experiment both in relative shape and in absolute value. A non-ro-vibrationally resolved, purely electronically resolved collisional radiative model for molecular hydrogen (including further states and reaction channels) is applied to investigate the validity of the corona approximation for the benchmark plasmas and the relevance of possible process extensions of the corona model. Furthermore, the influence of collisional quenching, autoionization and predissociation is discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗

Enabling topography-resolving structural dynamic contact simulation

Damping of structures and systems is often dominated by frictional dissipation in connections, the prediction of which remains a longstanding scientific challenge. Previous studies have shown that the actual topography of contact interfaces may have a strong effect, especially in the partial slip/liftoff regime. We recently proposed a multi-scale method, which couples finite element and boundary element modeling. The primary benefit of this approach is that it permits to analyze the effect of the actual contact topography on the dynamics of jointed structures. While this multi-scale modeling method was initially developed for quasi-static analysis, we demonstrate herein how it can be used for time step integration and Harmonic Balance analysis. We cross-verify those fully dynamic analysis methods against each other and quasi-static results, for the S4 Beam benchmark. We compare the multi-scale method against state-of-the-art full-FE analysis, in terms of numerical damping and computational performance. Some discrepancy is found to be of physical origin. Depending on the load history, it is shown that the system settles to a slightly different equilibrium. Finally, transient multi-scale simulations enable the prediction of this interesting phenomenon, for the first time, for a structure with bolted joints.

Frictional-unilateral contact↗

Aeroelastic Analysis of Aircraft: Wing and Wing/Fuselage Configurations

A previously developed interface method for coupling aerodynamics and structures is used to evaluate the aeroelastic effects for an advanced transport wing at cruise and under-cruise conditions. The calculated results are compared with wind tunnel test data. The capability of the interface method is also investigated for an MD-90 wing/fuselage configuration. In addition, an aircraft trim analysis is described and applied to wing configurations. The accuracy of turbulence models based on the algebraic eddy viscosity formulation of Cebeci and Smith is studied for airfoil flows at low Mach numbers by using methods based on the solutions of the boundary-layer and Navier-Stokes equations.

Chen, H. H.↗

The Ozonide Anion: A Theoretical Study

Complete active space self-consistent field (CASSCF) and CASSCF second order perturbation theory (CASPT2) methods have been used for the geometry optimization and calculation of harmonic and fundamental frequencies of the ozonide ion O3(-) and the ozonide lithium complex Li(+)O3(-). For O3(-) harmonic frequencies have also been obtained using the coupled-cluster method, CCSD(T). Infrared intensities are computed from dipole moment derivatives at the CASSCF level. The predicted equilibrium geometry for O3(-) is: Roo = 1.361 Angstroms and angle ooo = 115.4 degrees, and the fundamental frequencies are: nu(sub 1) = 989 per centimeter, nu(sub 2) 556 per centimeter, nu(sub 3) = 870 per centimeter (experimental values are: Roo =1.36 plus or minus 0.02 Angstroms, angle (ooo) = 111.8 plus or minus 2.0 degrees, nu(sub 1) = 975(50) per centimeter, nu(sub 2) = 550(50) per centimeter, nu(sub 3) = 880(50) per centimeter). Corresponding data for the lithium ozonide complex have also been obtained. The presented data contradict the previous interpretation of the IR and Raman spectrum obtained after deposition of ozone in N2, argon, or neon matrices with atomic beams of alkali metals. The presence of the lithium cation raises the asymmetric stretch frequency to about 940 per centimeter, which is contradictory to assumptions made in the assignments of the matrix spectra. Calculations made in a dielectric medium strongly suggest that the effect of the matrix on the IR spectrum is small for O3(-) itself. The dissociation and atomization energies of O3(-) are found to be in agreement with experiment.

Borowski, Poitr↗

Realistic fuel spray modeling for gasoline direct injection engine applications

Fuel spray modeling plays a critical role during modern gasoline direct injection (GDI) engine development due to fuel injection’s dominant impact on engine performance and emissions as well as the complex physical processes involved. In engineering three-dimensional (3D) computational fluid dynamics (CFD) simulations, the liquid-phase fuel atomization, evaporation, and mixing are usually modeled with the discrete droplet model (DDM) adopting a Lagrangian approach for multiphase CFD simulations. To this end, general practices heavily depend on the reduced order characterization of the injector nozzle flow. However, such simplified injector modeling may lead to insufficient representations of the complex spray dynamics. To tackle this problem, this study proposes a novel workflow to numerically evaluate GDI sub-cooled and flash-boiling sprays under engine-relevant conditions using a side-mounted GDI injector together with real gasoline fuel properties. The workflow introduces a one-way coupling (OWC) method leveraging high-fidelity nozzle flow simulations to provide realistic boundary conditions to the Lagrangian injector model. The proposed workflow was first verified in a constant volume chamber (CVC) environment and then implemented in a practical GDI engine setup to study spray morphology, fuel-air mixing, and wall-wetting propensity. In addition, detailed comparison was performed between the OWC method and the conventional rate of injection (ROI) routine. In conclusion, quantitative analysis of spray characteristics was conducted to highlight possible source of discrepancies of the conventional ROI method.

33 ADVANCED PROPULSION SYSTEMS↗