Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Physics 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,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
DOE OSTI2010

Use of Forward Sensitivity Analysis Method to Improve Code Scaling, Applicability, and Uncertainty (CSAU) Methodology

Code Scaling, Applicability, and Uncertainty (CSAU) methodology was developed in late 1980s by US NRC to systematically quantify reactor simulation uncertainty. Basing on CSAU methodology, Best Estimate Plus Uncertainty (BEPU) methods have been developed and widely used for new reactor designs and existing LWRs power uprate. In spite of these successes, several aspects of CSAU have been criticized for further improvement: i.e., (1) subjective judgement in PIRT process; (2) high cost due to heavily relying large experimental database, needing many experts man-years work, and very high computational overhead; (3) mixing numerical errors with other uncertainties; (4) grid dependence and same numerical grids for both scaled experiments and real plants applications; (5) user effects; Although large amount of efforts have been used to improve CSAU methodology, the above issues still exist. With the effort to develop next generation safety analysis codes, new opportunities appear to take advantage of new numerical methods, better physical models, and modern uncertainty qualification methods. Forward sensitivity analysis (FSA) directly solves the PDEs for parameter sensitivities (defined as the differential of physical solution with respective to any constant parameter). When the parameter sensitivities are available in a new advanced system analysis code, CSAU could be significantly improved: (1) Quantifying numerical errors: New codes which are totally implicit and with higher order accuracy can run much faster with numerical errors quantified by FSA. (2) Quantitative PIRT (Q-PIRT) to reduce subjective judgement and improving efficiency: treat numerical errors as special sensitivities against other physical uncertainties; only parameters having large uncertainty effects on design criterions are considered. (3) Greatly reducing computational costs for uncertainty qualification by (a) choosing optimized time steps and spatial sizes; (b) using gradient information (sensitivity result) to reduce sampling number. (4) Allowing grid independence for scaled integral effect test (IET) simulation and real plant applications: (a) eliminate numerical uncertainty on scaling; (b) reduce experimental cost by allowing smaller scaled IET; (c) eliminate user effects. This paper will review the issues related to the current CSAU, introduce FSA, discuss a potential Q-PIRT process, and show simple examples to perform FSA. Finally, the general research direction and requirements to use FSA in a system analysis code will be discussed.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

Probing parity violation with the four-point correlation function of BOSS galaxies

Parity-violating physics in the early universe can leave detectable traces in late-time observables. While vector- and tensor-type parity violation can be observed in the B -modes of the cosmic microwave background, scalar-type signatures are visible only in the four-point correlation function (4PCF) and beyond. This work presents a blind test for parity violation in the 4PCF of the baryon oscillation spectroscopic survey (BOSS) CMASS sample, considering galaxy separations in the range [20; 160]h -1 Mpc The parity-odd 4PCF contains no contributions from standard Λ CDM physics and can be efficiently measured using recently developed estimators. Data are analyzed using both a nonparametric rank test (comparing the BOSS 4PCFs to those of realistic simulations) and a compressed χ 2 analysis, with the former avoiding the assumption of a Gaussian likelihood. These find similar results, with the rank test giving a detection probability of 99.6% (2.9 σ). This provides significant evidence for parity violation, from either cosmological sources or systematics. We perform a number of systematic tests: although these do not reveal any observational artifacts, we cannot exclude the possibility that our detection is caused by the simulations not faithfully representing the statistical properties of the BOSS data. Our measurements can be used to constrain physical models of parity violation. As an example, we consider a coupling between the inflaton and a U (1) gauge field and place bounds on the latter’s energy density, which are several orders of magnitude stronger than those previously reported. Upcoming probes such as DESI and Euclid will reveal whether our detection of parity violation is due to new physics, and strengthen the bounds on a variety of models.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRSDate not supplied

Transfer function models for using empirical and physics-based simulation signal response data

In many situations, real or induced flaws such as tight cracks with known morphology cannot be manufactured in part geometry specimens or in real parts. Typically, surface fatigue cracks are manufactured in simple geometry specimens such as flat plates, dog-bone shaped flat or cylindrical specimens. If a nondestructive evaluation (NDE)technique is required to provide a reliably detectable flaw size, denoted asa90/95, for induced flaws in a part, then a direct method for qualifying the NDE procedure is to use appropriate induced flaw specimens and perform NDE procedure demonstration on the specimens. Probability of detection (POD)analysis of the empirical data may provide estimation of a90/95. This approach is described as direct POD demonstration testing, which may follow guidelines of MIL-HDBK-1823. This paper considers a case, where embedded tight cracklike induced flaws are to be detected reliably using a signal response based NDE procedure. Here, it is assumed that it is not practical to make surface or embedded induced flaw specimens in part geometry or configuration. Therefore, a direct POD demonstration testing cannot be undertaken. It is also assumed that simulation of signal response is possible for both surface and embedded induced flaws in part geometry specimens using a physics-based model. The proposed approach for NDE procedure qualification uses artificial flaws in simple geometry and part geometry specimens, and induced flaws in the same type of simple geometry specimens. Signal response data is taken on all sets of artificial and induced flaws in simple geometry and part geometry specimens. Moreover, simulated signal response data is generated for surface and embedded flaws. Thus, a case of five signal response versus flaw size datasets is considered. Three of the datasets are empirical and two datasets are physics model-based simulation datasets. A method of devising and using transfer function calculation dataset blocks to estimate the either the reliably detectable flaw size or the demonstration flaw size is provided.

Nondestructive evaluation↗
NASA NTRS2004

How Consistent are Recent Variations in the Tropical Energy and Water Cycle Resolved by Satellite Measurements?

One notable aspect of Earth's climate is that although the planet appears to be very close to radiative balance at top-of-atmosphere (TOA), the atmosphere itself and underlying surface are not. Profound exchanges of energy between the atmosphere and oceans, land and cryosphere occur over a range of time scales. Recent evidence from broadband satellite measurements suggests that even these TOA fluxes contain some detectable variations. Our ability to measure and reconstruct radiative fluxes at the surface and at the top of atmosphere is improving rapidly. One question is 'How consistent, physically, are these diverse remotely-sensed data sets'? The answer is of crucial importance to understanding climate processes, improving physical models, and improving remote sensing algorithms. In this work we will evaluate two recently released estimates of radiative fluxes, focusing primarily on surface estimates. The International Satellite Cloud Climatology Project 'FD' radiative flux profiles are available from mid-1983 to near present and have been constructed by driving the radiative transfer physics from the Goddard Institute for Space Studies (GISS) global model with ISCCP clouds and TOVS (TIROS Operational Vertical Sounder)thermodynamic profiles. Full and clear sky SW and LW fluxes are produced. A similar product from the NASA/GEWEX Surface Radiation Budget Project using different radiative flux codes and thermodynamics from the NASA/Goddard Earth Observing System (GEOS-1) assimilation model makes a similar calculation of surface fluxes. However this data set currently extends only through 1995. We also employ precipitation measurements from the Global Precipitation Climatology Project (GPCP) and the Tropical Rainfall Measuring Mission (TRMM). Finally, ocean evaporation estimates from the Special Sensor Microwave Imager (SSM/I) are considered as well as derived evaporation from the NCAR/NCEP Reanalysis. Additional information is included in the original extended abstract.

Robertson, F. R.↗
DOE OSTI2015

Leapfrogging of multiple coaxial viscous vortex rings

A recent theoretical study [Borisov, Kilin, and Mamaev, “The dynamics of vortex rings: Leapfrogging, choreographies and the stability problem,” Regular Chaotic Dyn. 18, 33 (2013); Borisov et al., “The dynamics of vortex rings: Leapfrogging in an ideal and viscous fluid,” Fluid Dyn. Res. 46, 031415 (2014)] shows that when three coaxial vortex rings travel in the same direction in an incompressible ideal fluid, each of the vortex rings alternately slips through (or leapfrogs) the other two ahead. Here, we use a lattice Boltzmann method to simulate viscous vortex rings with an identical initial circulation, radius, and separation distance with the aim of studying how viscous effect influences the outcomes of the leapfrogging process. For the case of two identical vortex rings, our computation shows that leapfrogging can be achieved only under certain favorable conditions, which depend on Reynolds number, vortex core size, and initial separation distance between the two rings. For the case of three coaxial vortex rings, the result differs from the inviscid model and shows that the second vortex ring always slips through the leading ring first, followed by the third ring slipping through the other two ahead. A simple physical model is proposed to explain the observed behavior.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Coupled Raman-Raman modes in the ionic Raman scattering process

In this artile, we employ a lattice Hamiltonian to examine the dynamics of coupled normal modes excited through a nonlinear phononic process, extending the model beyond infrared and first-order Raman active mode couplings to include interactions between symmetry-allowed Raman active modes. We examine the strength of the interactions between two Raman active modes on the targeted driven-mode dynamics to confirm that all symmetry allowed Raman modes interact with each other. We also apply the model to the correlated insulating ferromagnet YTiO 3 and present the resulting renormalization effects on the driven-mode dynamics from the anharmonic interactions. Owing to the dependence of the displacive amplitude of the Raman active mode on such interactions, especially for those modes with small amplitude, we suggest that models of anharmonic phononic coupling in materials with electronic, ferroic, or superconducting properties derived from competing Raman-like distortions should include these low-order terms in the equations of motion describing the excited phonons to obtain accurate physical models.

36 MATERIALS SCIENCE↗
NASA NTRS2008

U.S. Eastern Continental Shelf Carbon Cycling (USECoS): Modeling, Data Assimilation, and Analysis

Although the oceans play a major role in the uptake of fossil fuel CO2 from the atmosphere, there is much debate about the contribution from continental shelves, since many key shelf fluxes are not yet well quantified: the exchange of carbon across the land-ocean and shelf-slope interfaces, air-sea exchange of CO2, burial, and biological processes including productivity. Our goal is to quantify these carbon fluxes along the eastern U.S. coast using models quantitatively verified by comparison to observations, and to establish a framework for predicting how these fluxes may be modified as a result of climate and land use change. Our research questions build on those addressed with previous NASA funding for the USECoS (U.S. Eastern Continental Shelf Carbon Cycling) project. We have developed a coupled biogeochemical ocean circulation model configured for this study region and have extensively evaluated this model with both in situ and remotely-sensed data. Results indicate that to further reduce uncertainties in the shelf component of the global carbon cycle, future efforts must be directed towards 1) increasing the resolution of the physical model via nesting and 2) making refinements to the biogeochemical model and quantitatively evaluating these via the assimilation of biogeochemical data (in situ and remotely-sensed). These model improvements are essential for better understanding and reducing estimates of uncertainties in current and future carbon transformations and cycling in continental shelf systems. Our approach and science questions are particularly germane to the carbon cycle science goals of the NASA Earth Science Research Program as well as the U.S. Climate Change Research Program and the North American Carbon Program. Our interdisciplinary research team consists of scientists who have expertise in the physics and biogeochemistry of the U.S. eastern continental shelf, remote-sensing data analysis and data assimilative numerical models.

Mannino, Antonio↗
DOE OSTI2021

CFD Simulation of Helium Flow Loop Test Section

A helium flow loop is being assembled at Oak Ridge National Laboratory to analyze heat transfer enhancement for systems such as blanket and divertor components. To efficiently identify optimum geometries for heat transfer enhancement in these applications, simulation work is performed to optimize test section designs that are built and tested in the helium flow loop that operates at 4 MPa and a mass flow rate of 100 g/s. Different ribbed geometries that examine rib shape, rib height, rib orientation, rib spacing, and three dimensional orientation are modeled and simulated in STAR-CCM+ to compare their ability to remove heat and mitigate pressure drop. Following the simulations, models are selected and manufactured for the helium flow loop tests. Simulations initially focus on a hydrodynamic study to determine the appropriate mesh and physics models and then add a heat flux to analyze the heat transfer abilities of the models. The simulations are run in steady state and use a Reynolds-averaged Navier-Stokes k-ε turbulence model. The helium is modeled as an ideal gas. The simulation explores models of geometries that enhance the heat transfer and decrease pressure drop with an overall goal of increasing fluid collision with the wall. Enhanced geometries are simulated to select appropriate designs for manufacturing, and preliminary experimental results are used to validate the simulations. Furthermore, the factors that are being analyzed in the comparison between the experimental and the simulated results include matching thermocouple temperatures, pressure drop, roughness, and fluid velocity.

42 ENGINEERING↗
DOE OSTI2021

Understanding protein-complex assembly through grand canonical maximum entropy modeling

Inside a cell, heterotypic proteins assemble in inhomogeneous, crowded systems where the abundance of these proteins vary with cell types. While some protein complexes form putative structures that can be visualized with imaging, there are far more protein complexes that are yet to be solved because of their dynamic associations with one another. Nevertheless, it is possible to infer these protein complexes through a physical model. However, it is often not clear to physicists what kind of data from biology is necessary for such a modeling endeavor. Here, we aim to model these clusters of coarse-grained protein assemblies from multiple subunits through the constraints of interactions among the subunits and the chemical potential of each subunit. We obtained the constraints on the interactions among subunits from the known protein structures. We inferred the chemical potential that dictates the particle number distribution of each protein subunit from the knowledge of protein abundance from experimental data. Guided by the maximum entropy principle, we formulated an inverse statistical mechanical method to infer the distribution of particle numbers from the data of protein abundance as chemical potentials for a grand canonical multicomponent mixture. Using grand canonical Monte Carlo simulations, we captured a distribution of high-order clusters in a protein complex of succinate dehydrogenase with four known subunits. The complexity of hierarchical clusters varies with the relative protein abundance of each subunit in distinctive cell types such as lung, heart, and brain. When the crowding content increases, we observed that crowding stabilizes emergent clusters that do not exist in dilute conditions. We, therefore, proposed a testable hypothesis that the hierarchical complexity of protein clusters on a molecular scale is a plausible biomarker of predicting the phenotypes of a cell.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2010

Hadronic Physics: an Outlook

A brief outlook, in two senses, is presented for hadronic physics. The likely near term future for experiment and lattice effort is sketched and I speculate on future directions in theory. I also look out at other fields, presenting a short review of QCD ideas in ''Beyond the Standard Model'' physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Wave Energy

After reflecting on the history of ocean wave energy research and technology development, the most prominent classification schemes for wave energy converter (WEC) technology concepts are produced, highlighting the diversity of WEC technologies that were or are under development at different levels of maturity and intensity. Three advanced WEC technology examples are provided in this article. The ocean wave energy resource is addressed with respect to its global energy potential, local variability, and issues associated with the established forms of spectral power presentation. The spectrum of theoretical modeling approaches is described, along with the range of simulation tasks and related tools. This article also covers the specific aspects of the physical modeling and scaled model testing of WES systems and subsystems. Following a comparison of the engineering and design challenges of wave energy and wind energy technology development, the WEC technology development status and process is described and analyzed. Finally, the requirements for successful WEC technology development for achieving economic WEC technology and energy production under commercial conditions are discussed, including a number of possible development and innovation pathways.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Modeling rarefied gas chemistry with QuiPS, a novel quasi-particle method

The goal of this work is to build up the capability of quasi-particle simulation (QuiPS), a novel flow solver, such that it can adequately model the rarefied portion of an atmospheric reentry trajectory. Direct simulation Monte Carlo (DSMC) is the conventional solver for such conditions, but struggles to resolve transient flows, trace species, and high-level internal energy states due to stochastic noise. Quasi-particle simulation (QuiPS) is a novel Boltzmann solver that describes a system with a discretized, truncated velocity distribution function. The resulting fixed-velocity, variable weight quasi-particles enable smooth variation of macroscopic properties. The distribution function description enables the use of a variance-reduced collision model, greatly minimizing expense near equilibrium. This work presents the addition of a neutral air chemistry model to QuiPS and some demonstrative 0D simulations. The explicit representation of internal distributions in QuiPS reveals some of the flaws in existing physics models. Overall, variance reduction, a key feature of QuiPS, can greatly reduce expense of multi-dimensional calculations, but is only cheaper when the gas composition is near chemical equilibrium.

42 ENGINEERING↗
DOE OSTI2024

A combinatorially complete epistatic fitness landscape in an enzyme active site

Protein engineering often targets binding pockets or active sites which are enriched in epistasis—nonadditive interactions between amino acid substitutions—and where the combined effects of multiple single substitutions are difficult to predict. Few existing sequence-fitness datasets capture epistasis at large scale, especially for enzyme catalysis, limiting the development and assessment of model-guided enzyme engineering approaches. We present here a combinatorially complete, 160,000-variant fitness landscape across four residues in the active site of an enzyme. Assaying the native reaction of a thermostable β-subunit of tryptophan synthase (TrpB) in a nonnative environment yielded a landscape characterized by significant epistasis and many local optima. These effects prevent simulated directed evolution approaches from efficiently reaching the global optimum. There is nonetheless wide variability in the effectiveness of different directed evolution approaches, which together provide experimental benchmarks for computational and machine learning workflows. The most-fit TrpB variants contain a substitution that is nearly absent in natural TrpB sequences—a result that conservation-based predictions would not capture. Thus, although fitness prediction using evolutionary data can enrich in more-active variants, these approaches struggle to identify and differentiate among the most-active variants, even for this near-native function. Overall, this work presents a large-scale testing ground for model-guided enzyme engineering and suggests that efficient navigation of epistatic fitness landscapes can be improved by advances in both machine learning and physical modeling.

biocatalysis↗
DOE OSTI2017

On the galaxy–halo connection in the EAGLE simulation

Empirical models of galaxy formation require assumptions about the correlations between galaxy and halo properties. These may be calibrated against observations or inferred from physical models such as hydrodynamical simulations. In this Letter, we use the EAGLE simulation to investigate the correlation of galaxy size with halo properties. We motivate this analysis by noting that the common assumption of angular momentum partition between baryons and dark matter in rotationally supported galaxies overpredicts both the spread in the stellar mass–size relation and the anticorrelation of size and velocity residuals, indicating a problem with the galaxy–halo connection it implies. We find the EAGLE galaxy population to perform significantly better on both statistics, and trace this success to the weakness of the correlations of galaxy size with halo mass, concentration and spin at fixed stellar mass. Here by, using these correlations in empirical models will enable fine-grained aspects of galaxy scalings to be matched.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1975

Development of a three-dimensional time-dependent flow field model

A three-dimensional, time-dependent mathematical model to represent Mobile Bay was developed. Computer programs were developed which numerically solve the appropriate conservation equations for predicting bay and estuary flow fields. The model is useful for analyzing the dispersion of sea water into fresh water and the transport of sediment, and for relating field and physical model data.

Farmer, R. C.↗
NASA NTRS1993

Calculations of combustion response profiles and oscillations

The theory and procedures for determining the characteristics of pressure oscillations in rocket engines with prescribed burning rate oscillations are presented. Pressure and velocity oscillations calculated using this procedure are presented for the Space Shuttle Main Engine (SSME) to show the influence of baffles and absorbers on the burning rate oscillations required to achieve neutral stability. Results of calculations to determine local combustion responses using detailed physical models for injection, atomization, and vaporization with gas phase oscillations in baffled and unbaffled SSME combustors are presented. The contributions of the various physical phenomena occurring in a combustor to oscillations in combustion response were determined.

Priem, Richard J.↗
NASA NTRS1995

Computational Modeling of the Ames 11-Ft Transonic Wind Tunnel in Conjunction with IofNEWT

Technical advances in Computational Fluid Dynamics have now made it possible to simulate complex three-dimensional internal flows about models of various size placed in a Transonic Wind Tunnel. TWT wall interference effects have been a source of error in predicting flight data from actual wind tunnel measured data. An advantage of such internal CFD calculations is to directly compare numerical results with the actual tunnel data for code assessment and tunnel flow analysis. A CFD capability has recently been devised for flow analysis of the NASA/Ames 11-Ft TWT facility. The primary objectives of this work are to provide a CFD tool to study the NASA/Ames 11-Ft TWT flow characteristics, to understand the slotted wall interference effects, and to validate CFD codes. A secondary objective is to integrate the internal flowfield calculations with the Pressure Sensitive Paint data, a surface pressure distribution capability in Ames' production wind tunnels. The effort has been part of the Ames IofNEWT, Integration of Numerical and Experimental Wind Tunnels project, which is aimed at providing further analytical tools for industrial application. We used the NASA/Ames OVERFLOW code to solve the thin-layer Navier-Stokes equations. Viscosity effects near the model are captured by Baldwin-Lomax or Baldwin-Barth turbulence models. The solver was modified to model the flow behavior in the vicinity of the tunnel longitudinal slotted walls. A suitable porous type wall boundary condition was coded to account for the cross-flow through the test section. Viscous flow equations were solved in generalized coordinates with a three-factor implicit central difference scheme in conjunction with the Chimera grid procedure. The internal flow field about the model and the tunnel walls were descretized by the Chimera overset grid system. This approach allows the application of efficient grid generation codes about individual components of the configuration; separate minor grids were developed to resolve the model and overset onto a main grid which discretizes the interior of the tunnel test section. Individual grid components axe not required to have mesh boundaries joined in any special way to each other or to the main tunnel grid. Programs have been developed to rotate the model about the tunnel pivot point and rotation axis, similar to that of the tunnel turntable mechanism for adjusting the pitch of the physical model in the test section.

Djomehri, M. Jahed↗
NASA NTRS2006

Earth Satellite Population Instability: Underscoring the Need for Debris Mitigation

A recent study by NASA indicates that the implementation of international orbital debris mitigation measures alone will not prevent a significant increase in the artificial Earth satellite population, beginning in the second half of this century. Whereas the focus of the aerospace community for the past 25 years has been on the curtailment of the generation of long-lived orbital debris, active remediation of the current orbital debris population should now be reconsidered to help preserve near-Earth space for future generations. In particular, we show in this paper that even if launch operations were to cease today, the population of space debris would continue to grow. Further, proposed remediation techniques do not appear to offer a viable solution. We therefore recommend that, while the aerospace community maintains the current debris-limiting mission regulations and postmission disposal procedures, future emphasis should be placed on finding new remediation technologies for solving this growing problem. Since the launch of Sputnik 1, space activities have created an orbital debris environment that poses increasing impact risks to existing space systems, including human space flight and robotic missions (1, 2). Currently, more than 9,000 Earth orbiting man-made objects (including many breakup fragments), with a combined mass exceeding 5 million kilograms, are tracked by the US Space Surveillance Network and maintained in the US satellite catalog (3-5). Three accidental collisions between cataloged satellites during the period from late 1991 to early 2005 have already been documented (6), although fortunately none resulted in the creation of large, trackable debris clouds. Several studies conducted during 1991-2001 demonstrated, with assumed future launch rates, the unintended growth potential of the Earth satellite population, resulting from random, accidental collisions among resident space objects (7-13). In some low Earth orbit (LEO) altitude regimes where the number density of satellites is above a critical spatial density, the production rate of new satellites (i.e., debris) due to collisions exceeds the loss of objects due to orbital decay. NASA s evolutionary satellite population model LEGEND (LEO-to-GEO Environment Debris model), developed by the Orbital Debris Program Office at the NASA Lyndon B. Johnson Space Center, is a high fidelity three-dimensional physical model that is capable of simulating the historical satellite environment, as well as the evolution of future debris populations (14, 15). The subject study assumed no rocket bodies and spacecraft were launched after December 2004, and no future disposal maneuvers were allowed for existing spacecraft, few of which currently have such a capability. The rate of satellite explosions would naturally decrease to zero within a few decades as the current satellite population ages. The LEGEND future projection adopts a Monte Carlo approach to simulate future on-orbit explosions and collisions. Within a given projection time step, once the explosion probability is estimated for an intact object, a random number is drawn and compared with the probability to determine if an explosion would occur. A similar procedure is applied to collisions for each pair of target and projectile involved within the same time step. Due to the nature of the Monte Carlo process, multiple projection runs must be performed and analyzed before one can draw reliable and meaningful conclusions from the outcome. A total of fifty, 200-year future projection Monte Carlo simulations were executed and evaluated (16).

Liou, Jer-chyi↗