Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure modeling”

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 1,117 records · Page 62

Vadose Zone Model for B-3 Pond for Composite Analysis

The objectives of the vadose modeling for the updated Hanford Site composite analysis (CA) are to simulate the flow and transport of water and radionuclide releases from the surface to the water table and to provide radionuclide transfer rates for the plateau to river (P2R) model, version 8.3. Water additions include natural recharge and water discharged to the ground as a result of industrial processes associated with Hanford Site operations. Contaminant sources include radionuclides in water discharged to the ground during operations and radionuclides disposed “dry” in solid waste burial grounds or other means. The simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The scope of this ECF is to document the development and results of the B-3 Pond vadose zone model. The B-3 Pond model covers the main lobe of the 216-B-3 Pond System located in a topographic low east of the 200 East Area. The pond system operated from 1945 to 1997 and received effluents from the Plutonium Uranium Extraction Plant, B Plant, 284-E Powerhouse, and other 200 East Area facilities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Bidirectional reflectance of oceanic waters: A comparison of modeled and measured upward radiance fields

The bidirectional reflectance of oceanic waters is conveniently described in a normalized way by forming the ratio of the upwelling irradiance E(sub u) to any upwelling radiance L(sub u)(theta prime, phi). This ratio, Q (theta prime, theta(sub 0), (phi(sub 0) - phi), where theta prime, phi are the nadir and azimuth angles for the upward radiance and theta(sub 0), phi(sub 0) are the zenith and azimuth angles of the Sun, has been determined from measurements at sea and computed via Monte Carlo simulations using the inherent optical properties measured in the field and appropriate boundary conditions (clear sky, no wind, varying Sun angle). Experimental ad computed Q values are in excellent agreement. This successful comparison confirms the importance of the bidirectional character of ocean reflectance, already pointed out from a purely numerical approach without field validation, and corroborates the extended range of the Q variations. The later point is of importance when interpreting the marine signals detected by an ocean color satellite-borne sensor. The validation is extended by considering the historical data for the radiance distributions in Lake Pend Oreille determined at various depths. The closure issue in ocean optics is examined by solving the direct problem of radiative transfer and through a model-data comparison in terms of radiance field.

Morel, Andre↗

Reconstruction of Electron and Ion Distribution Functions in a Magnetotail Reconnection Diffusion Region

In the diffusion region of magnetotail reconnection, particle distributions are highly structured, exhibiting triangular shapes and multiple striations that deviate dramatically from the Maxwellian distribution. Fully kinetic simulations have been demonstrated to be capable of producing the essential structures of the observed distribution functions, yet are computationally not feasible for 3D global simulations. The fluid models used for large-scale simulations, on the other hand, do not have the kinetic physics necessary for describing reconnection accurately. Our study aims to bridge fully kinetic and fluid simulations by quantifying the information required to capture the non-Maxwellian features in the distributions underlying the closures used in the fluid code. Here, we compare the results of fully kinetic simulations with observed electron velocity distributions in a magnetotail reconnection diffusion region and use the maximum entropy model to reconstruct electron and ion distributions using various numbers of moments obtained from the simulation. Our results indicate that using only local moments, the maximum entropy model can reproduce many of the features of the distributions: (1) the electron outflow distribution with a tilted triangular structure is reproduced with 21 or more moments in agreement with Ng et al. and (2) counterstreaming distributions can be captured with the 35-moment model when the separation in velocity space between the populations is large.

58 GEOSCIENCES↗

An experimental documentation of a separated trailing-edge flow at a transonic Mach number

A detailed experiment on the separated flow field at a sharp trailing edge is described and documented. The separated flow is a result of sustained adverse pressure gradients. The experiment was conducted using an elongated airfoil-like model at a transonic Mach number and at a high Reynolds number of practical interest. Measurements made include surface pressures and detailed mean and turbulence flow quantities in the region just upstream of separation to downstream into the near-wake, following wake closure. The data obtained are presented mostly in tabular form. These data are of sufficient quality and detail to be useful as a test case for evaluating turbulence models and calculation methods.

Viswanath, P. R.↗

SOLPS-ITER modeling of a dedicated divertor for negative triangularity operation on DIII-D

The design of a new dedicated divertor for negative triangularity (NT) operation on DIII-D with neutral baffles and pumping is informed by SOLPS-ITER transport modeling. This dedicated NT divertor is the latest step in a progression of NT shapes with various divertor characteristics explored on DIII-D, including NT shapes at reduced triangularity and a campaign with stronger shaping that included new armored components on the outboard side. SOLPS simulations played a key role in these divertor designs. Interpretive simulations, using cross-field diffusivities constrained by experimental data in the NT Shelf shape were used to inform the design of the 2023 armor campaign components. A similar procedure used armor campaign data to predict conditions for the dedicated NT divertor. The predictive simulations were used to assess the divertor fluxes, detachment threshold, pumped flux, and neutral leakage. For the dedicated NT divertor, SOLPS simulations and two-point-modeling were used to show the relative impact of magnetic topology (mainly longer connection length) and divertor closure on the divertor conditions relative to the armor campaign. It is predicted that the dedicated NT divertor reaches detachment (measured by target ion flux rollover) at a lower upstream density (≈(1.75−−2.4)×1⁢019m−3) as compared to the armor campaign shape. For the preliminary design geometry, divertor closure reduces the neutral leakage by ≈10%. Parametric optimization indicating further ≈20%–60% improvement in the leakage flux and recycled flux crossing the pump entrance is possible for relatively minor changes to the divertor and baffle layout.

Lore, Jeremy [ORNL] (ORCID:000000029192465X)↗

Nucleonic shells and nuclear masses

The binding energy of an isotope is a sensitive indicator of the underlying shell structure as it reflects the net energy content of a nucleus. Since magic nuclei are significantly lighter, or more bound, compared to their neighbors, the presence of nucleonic shell structure makes an imprint on nuclear masses. In this study, using a carefully designed binding-energy indicator, we catalog the appearance of spherical and deformed shell and subshell closures throughout the nuclear landscape. After presenting experimental evidence for shell and subshell closures as seen through the lens of nuclear masses, we study the ability of global nuclear mass models to predict local binding-energy variations related to shell effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Literature Review of Documented Persistent Secondary Uranium Sources at DOE LM Sites

The U.S. Department of Energy Office of Legacy Management (LM) has the responsibility to protect human health and the environment at 101 legacy uranium sites across the nation by implementing long-term surveillance and maintenance requirements. Included in these responsibilities is to provide for the safe and environmentally sound disposal, long-term stabilization, and control of uranium mill tailings in a manner that minimizes or eliminates radiation health hazards to the public at 21 Uranium Mill Tailings Radiation Control Act (UMTRCA) Title I sites; 6 UMTRCA Title II sites; and 9 Comprehensive Environmental Response, Compensation, and Liability Act/Resource Conservation and Recovery Act sites. These sites are former uranium mills (processing sites), or sites with wastes from processing sites encapsulated in engineered disposal cells (disposal sites), or sites where components of atomic weapons were manufactured. During site characterization, groundwater contamination was often discovered and evaluated or modeled for future movement with the inclusion of uranium sorption in the groundwater plume. Much of this original work did identify uranium processing wastes as the original contaminant source. However, more recent work has identified the existence and influence of persistent secondary uranium sources that may prevent uranium in groundwater from attenuating as quickly as initial evaluations or modeling indicated. By understanding these processes, better decisions can be made to choose the appropriate remediation and closure strategies, identify long-term surveillance and maintenance activities, and more accurately estimate life-cycle costs. The purpose of this report is to summarize the literature available on documented persistence secondary uranium sources at LM sites as part of the Applied Studies and Technology Persistent Secondary Contaminant Sources Project.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Uncertainty quantification for equations of state: copper as an example

Equations of state are essential for providing a fundamental description of materials properties in thermodynamic equilibrium and are used to provide closure relations for hydrodynamics simulations. Generally, equations of state rely on simple physics-based parameterized materials models to inform on the free energy of a material through out a given thermodynamic state space. Historically the parameters of these models have been tuned by hand to fit various experimental data. However, modern optimization and uncertainty quantification techniques allow us to quickly test thousands of parameter combinations and obtain meaningful uncertainty estimates on the parameters, opening opportunities for assessing systematic uncertainties in experiments, assessing model adequacy, and more. In this report, we use Bayesian inference to fit the solid (fcc) equation of state of copper. We focus on fitting five different experimental datasets, including the isobaric density, isobaric heat capacity, room temperature isotherm, principal isentrope, and principal Hugoniot. We fit all five data types simultaneously, and then explore the extent to which combinations of 2 subsets of the 5 datasets can constrain the EOS parameters, as compared to the fit to all 5. This information is useful for investigating the extent to which different datasets can con strain EOS models and thereby help guide experimental investigations in order to best constrain the EOS. We also discuss ways that the methodologies can be used to investigate systematic discrepancies between experiments, as well as how the methods can be used to assess model uncertainty. The framework we develop is general, in that it can be used with a variety of optimization or uncertainty quantification techniques and with a variety of data sources, including both experimental and ab-inito data.

97 MATHEMATICS AND COMPUTING↗

Crack growth in Ti-8Al-1Mo-1V with real-time and accelerated flight-by-flight loading

Crack growth in Ti-8Al-1Mo-1V was measured and calculated for real-time and accelerated simulations of supersonic airplane loading and heating. Crack-growth rates calculated on the assumption that an entire flight could be represented by a single cycle predicted the experimental rates poorly. Calculated crack-growth rates were slower than the experimental rates for all tests with flight-by-flight loading. For room-temperature accelerated tests, the calculated rates agreed well with the experimental rates; but the calculations became progressively less accurate for progressively more complex test conditions (tests that included elevated temperature). Calculations of crack growth using the crack-closure concept can probably be improved through study of crack-opening stresses using finite-element models that account for variable-amplitude loading, residual stresses, and temperature effects. The calculations of crack growth could also be improved through detailed studies of material properties and interactions among stress, temperature, and time as appropriate for the real-time operating conditions of a supersonic transport airplane.

Imig, L. A.↗

Closure of the Africa-Eurasia-North America plate motion circuit and tectonics of the Gloria fault

The current motions of the African, Eurasian, and North American plates are examined. The problems addressed include whether there is resolvable motion of a Spitsbergen microplate, the direction of motion between the African and North American plates, whether the Gloria fault is an active transform fault, and the implications of plate circuit closures for rates of intraplate deformation. Marine geophysical data and magnetic profiles are used to construct a model which predicts about 4 mm/yr slip across the Azores-Gibraltar Ridge, and west-northwest convergence near Gibraltar. The analyzed data are consistent with a rigid plate model with the Gloria fault being a transform fault.

Argus, Donald F.↗

Modeling of turbulent shear flows

The current progress is documented in the research and development of modeling techniques for turbulent shear flows. These include a two-scale model for compressible turbulent flows and a new energy transfer model. The former represents the status of the efforts to identify compressibility effects in turbulence. The energy transfer model refines a weakly nonlinear wave model developed earlier, which models directly the turbulent large structures. The objective of these activities is to develop second-order closures for compressible turbulent flows.

Liou, William W.↗

The Pathway to a Safe and Effective Medication Formulary for Exploration Spaceflight

PURPOSE: Exploration space missions pose several challenges to providing a comprehensive medication formulary designed to accommodate the size and space limitations of the spacecraft; while addressing the individual medications needs and preferences of the Crew; the negative outcome of a degrading inventory over time, the inability to resupply before expiration dates; and the need to properly forecast the best possible medication candidates to treat conditions that will occur in the future. METHODS: The Pharmacotherapeutics Discipline has partnered with the Exploration Medical Capabilities (ExMC) Element to develop and propose a research pathway that is comprehensively focused on evidence-based models and theories, as well as on new diagnostic tools and treatments or preventive measures aimed at closure of the Med02 “Pharmacy” Gap; defined in the Human Research Program’s (HRP) risk-based research strategy. The Med02 Gap promotes the challenge to identify a strategy to ensure that medications used to treat medical conditions during exploration space missions are available, safe, and effective. It is abundantly clear that pharmaceutical intervention is an essential component of risk management planning for astronaut healthcare during exploration space. However, the quandary still remains of how to assemble a formulary that is comprehensive enough to prevent or treat anticipated medical events; and is also chemically stable, safe, and robust enough to have sufficient potency to last for the duration of an exploration space mission. In cases where that is not possible, addressing this Gap requires exploration of novel drug development techniques, dosage forms, and dosage delivery platforms that enhance chemical stability as well as therapeutic effectiveness. RESULTS: The proposed research pathway outlines the steps, processes, procedures, and a research portfolio aimed at identifying a capability that will provide a safe and effective pharmacy for any specific exploration Design Reference Mission (DRM). The proposed approach to building this research portfolio is to seek research projects that concentrate on four major focus areas; (1) Formulary selection, (2) Formulary potency and shelf life, (3) Formulary safety and toxicity, and (4) Novel technology and innovation such as portable real-time chemical analysis innovative drug therapies and dosage and delivery platforms. CONCLUSION: The research pathway has been completed and presented to the HRP. In spring 2017, it is scheduled to be reviewed by a panel of pharmaceutical and clinical experts that will evaluate the scientific merit and operational feasibility of the research pathway, as well as make suggestions for any warranted additions or improvements. Once finalized, the ExMC Element will proceed with the execution of this research pathway with the goal of gathering as much data, and learning as much as possible, to provide a safe and effective pharmaceutical formulary for use during exploration missions.

Daniels, V. R.↗

Simulating the Clamped Tapered Beam Specimen Under Quasi-Static and Fatigue Loading Using Floating Node Method

As part of the NASA Advanced Composites Project (ACP), a sub-element has been designed to provide validation data for progressive damage analysis models. The clamped tapered beam is a cross-ply laminated composite specimen designed to validate the simulation of the onset of matrix cracks and their interaction with delaminations, including delamination migration. A tapered geometry was used to localize the first damage occurrence in the tapered region, without prescribing an initial crack. The boundary and loading conditions were chosen to favor delamination growth and subsequent migration after the first damage occurrence. The typical sequence of events consists of a matrix crack located at the tapered region, leading to delamination onset, followed by delamination growth and subsequent delamination migration to a different interface via a dominant matrix crack. The Clamped Tapered Beam (CTB) was tested in both quasi-static and fatigue regimes. The results obtained are used in this study to assess and validate a methodology based on the Floating Node Method (FNM) implemented as an Extended Interface Element. In this methodology, quasi-static and fatigue damage formation and development are modeled by combining FNM to represent crack networks, with Directional Cohesive Zone Elements (DCZE) and Virtual Crack Closure Technique (VCCT), respectively. Qualitatively, the methodology is capable of predicting the sequence of events and overall failure morphology. Quantitatively, the simulation results generally bound the experimental data, based on the range of the characterization data used. In this paper, the results from quasi static and fatigue simulations are compared and correlated with experimental data.

Seshadri, B. R.↗

Vadose Zone Model for B-3C Pond for Composite Analysis

The objectives of the vadose modeling for the updated Hanford Site composite analysis (CA) are to simulate the flow and transport of water and radionuclide releases from the surface to the water table and to provide radionuclide transfer rates for the plateau to river (P2R) model, version 8.3. Water additions include natural recharge and water discharged to the ground as a result of industrial processes associated with Hanford Site operations. Contaminant sources include radionuclides in water discharged to the ground during operations and radionuclides disposed “dry” in solid waste burial grounds or other means. The simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The scope of this ECF is to document the development and results of the B-3C Pond vadose zone model. The B-3C Pond model covers the 216-B-3C Expansion Pond of the 216-B-3 Pond System located in a topographic low east of the 200 East Area. The pond system operated from 1945 to 1997 and received effluents from the Plutonium Uranium Extraction Plant, B Plant, 284-E Powerhouse, and other 200 East Area facilities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Scaling Behavior of a Turbulent Kinetic Energy Closure Scheme for the Stably Stratified Atmosphere: A Steady-State Analysis

Abstract We present a turbulent kinetic energy (TKE) closure scheme for the stably stratified atmosphere in which the mixing lengths for momentum and heat are not parameterized in the same manner. The key difference is that, while the mixing length for heat tends toward the stability independent mixing length for momentum in neutrally stratified conditions, it tends toward one based on the Brunt–Väisälä time scale and square root of the TKE in the limit of large stability. This enables a unique steady-state solution for TKE to be obtained, which we demonstrate would otherwise be impossible if the mixing lengths were the same. Despite the model’s relative simplicity, it is shown to perform reasonably well with observational data from the 1999 Cooperative Atmosphere–Surface Exchange Study (CASES-99) using commonly employed model constants. Analyzing the scaling behavior of the nondimensional velocity and potential temperature gradients, or of the stability (correction) functions, reveals that for large stability the present model scales in the same manner as the first-order operational scheme of Viterbo et al. Alternatively, it appears as a blend of two cases of the TKE closure scheme of Baas et al. Critically, because a unique steady-state TKE can be obtained, the present model avoids the nonphysical behavior identified in one of the cases of Baas et al.

Meteorology & Atmospheric Sciences↗

Modeling Progressive Damage Using Local Displacement Discontinuities Within the FEAMAC Multiscale Modeling Framework

A method for performing progressive damage modeling in composite materials and structures based on continuum level interfacial displacement discontinuities is presented. The proposed method enables the exponential evolution of the interfacial compliance, resulting in unloading of the tractions at the interface after delamination or failure occurs. In this paper, the proposed continuum displacement discontinuity model has been used to simulate failure within both isotropic and orthotropic materials efficiently and to explore the possibility of predicting the crack path, therein. Simulation results obtained from Mode-I and Mode-II fracture compare the proposed approach with the cohesive element approach and Virtual Crack Closure Techniques (VCCT) available within the ABAQUS (ABAQUS, Inc.) finite element software. Furthermore, an eccentrically loaded 3-point bend test has been simulated with the displacement discontinuity model, and the resulting crack path prediction has been compared with a prediction based on the extended finite element model (XFEM) approach.

Ranatunga, Vipul↗

Carrying the mass flux terms exactly in the first and second moment equations of compressible turbulence

In compressible turbulence models, it is assumed that the Favre-mean velocities are suitable approximations to the Reynolds-mean velocities in order to close unknown terms. This neglects, in the mean momentum and energy equations, the contribution to the stress and work terms by the mean of the fluctuating Favre velocity, a quantity proportional to the turbulent mass flux. As the stress and work terms do not introduce any new unknown correlations requiring closure in either k-epsilon or Reynolds stress closures and because the exact form of the terms can, with little additional work, be carried there is no need to make any modeling assumptions. In the Reynolds stress equations the viscous terms appear naturally in Reynolds variables while the problem is posed in Favre variables. In the process of splitting the viscous terms into the viscous transport terms, carried in Favre variables, and the dissipation terms, carried in Reynolds variables, important contributions from the mass flux appear. The accurate accounting of these terms is important for any consistent near wall modeling and the retention of the mass flux terms is important in complex compressible turbulent flows.

Ristorcelli, J. R.↗