Search NASA⌕ Search

SEARCH · Search NASA

Results for “force 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 883 records · Page 49

FARGO3D: A NEW GPU-ORIENTED MHD CODE

We present the FARGO3D code, recently publicly released. It is a magnetohydrodynamics code developed with special emphasis on the physics of protoplanetary disks and planet–disk interactions, and parallelized with MPI. The hydrodynamics algorithms are based on finite-difference upwind, dimensionally split methods. The magnetohydrodynamics algorithms consist of the constrained transport method to preserve the divergence-free property of the magnetic field to machine accuracy, coupled to a method of characteristics for the evaluation of electromotive forces and Lorentz forces. Orbital advection is implemented, and an N-body solver is included to simulate planets or stars interacting with the gas. We present our implementation in detail and present a number of widely known tests for comparison purposes. One strength of FARGO3D is that it can run on either graphical processing units (GPUs) or central processing units (CPUs), achieving large speed-up with respect to CPU cores. We describe our implementation choices, which allow a user with no prior knowledge of GPU programming to develop new routines for CPUs, and have them translated automatically for GPUs.

79 ASTRONOMY AND ASTROPHYSICS↗

Higher-order phonon scattering: advancing the quantum theory of phonon linewidth, thermal conductivity and thermal radiative properties

Phonon scattering plays a central role in the quantum theory of phonon linewidth, which in turn governs important properties including infrared spectra, Raman spectra, lattice thermal conductivity, thermal radiative properties, and also signifi- cantly affects other important processes such as hot electron relaxation. Since Maradudin and Fein’s classic work in 1962, three-phonon scattering had been considered as the dominant intrinsic phonon scattering mechanism and has seen tremendous advances. However, the role of the higher-order four-phonon scattering had been persistently unclear and so was ignored. The tremendous complexity of the formalism and computational challenges stood in the way, prohibiting the direct and quantitative treatment of four-phonon scattering. In 2016, a rigorous four-phonon scattering formalism was developed, and the prediction was realized using empirical potentials. In 2017, the method was extended using first-principles calculated force constants, and the thermal conductivities of boron arsenides (BAs), Si and diamond were predicted. The predictions for BAs were later confirmed by several independent experiments. Four-phonon scattering has since been investigated in a range of materials and established as an important intrinsic scattering mechanism for thermal transport and radiative properties. Specifically, four-phonon scattering is important when the fourth-order scattering potential or phase space becomes relatively large. The former scenario includes: (i) nearly all materials when the temperature is high; (ii) strongly anharmonic (low thermal conductivity) materials, including most rocksalt compounds, halides, hydrides, chalcogenides and oxides. The latter scenario includes: (iii) materials with large acoustic–optical phonon band gaps, such as XY compounds with a large atomic mass ratio between X and Y; (iv) two- dimensional materials with reflection symmetry, such as single-layer graphene, single-layer boron nitride and carbon nanotubes; and (v) phonons with a large density of states, such as optical phonons, which are important for Raman, infrared and thermal radiative properties. Four-phonon scattering is expected to gain broad interest in various technologically important materials for thermoelectrics, thermal barrier coatings, thermal energy storage, phase change, nuclear power, ultra-high temperature ceramics, infrared spectra, Raman spectra, radiative transport, hot electron relaxation and radiative cooling. Four-phonon scattering has been, and will continue to be, established as an important intrinsic phonon scattering mechanism beyond three-phonon scattering. The prediction of four-phonon scattering will transition from a breakthrough to a new routine in the next decade.

Feng, Tianli↗

Toward Linear Stability Analysis of Upper Plenum Flows

Linear stability analysis is a tool used to help identify the sensitivity of a particular domain to perturbations and evaluate the system dynamics as a transition, if any, occurs. The development of a time invariant base flow for future use in a linear stability analysis is the focus of this investigation. A test case, with a Reynolds number of 100, was brought to a steady state using the spectral element code Nek5000 and the method of selective frequency damping. The convergence criteria for this method is the volume averaged value of the forcing function which, after approximately 80,000 time-steps, decreased to 3.6x10{sup -4}. Analysis shows that the case has yet to have reached the tolerance desired, 1x10{sup -6}. With the simulation in the process of converging a preliminary perturbation analysis was conducted. Due to the low Reynolds number of this case, 100, it was hypothesized that the flow would be linearly stable. This preliminary stability analysis only focused on the largest fundamental mode and the energy growth of that perturbation over time. Upon evaluation of the perturbation energy, which was shown to be monotonically decreasing in time, it was concluded that, despite not yet being converged, the system is in fact linearly stable. Future work includes evaluating the domain at higher Reynolds numbers, solving the eigenvalue problem created when the linearized Navier-Stokes are written in operator form and calculating the energy evolution of the optimal perturbation. (authors)

42 ENGINEERING↗

Understanding syringeability and injectability of high molecular weight PEO solution through time-dependent force-distance profiles

Opioid medications play a vital role in treating moderate to severe pain. Unfortunately, many drug misusers and abusers attempt to alter the formulations or properties of these drugs by manipulation, (e.g., crushing, chewing, smoking, snorting, injecting). The intravenous (IV) route is most dangerous to abusers, as the drugs directly enter the circulatory system and produce intense euphoria. To obtain a full understanding of the impact of syringe factors (e.g., needle gauge size, needle length, syringe barrel size), on the ease of injection, we undertook a comprehensive assessment of syringeability and injectability of manipulated abuse-deterrent formulations (ADFs). A texture analyzer-based testing method was developed for the measurement of the resistance force of pulling, holding, and pushing phases of injections. Results showed that the finer needle gauge sizes required higher injection force to withdraw drug solutions. In addition, the syringed liquid volume was highly dependent on needle gauge size, holding time, and sample viscosity. In most cases, a lower needle gauge number and a longer holding time increased the syringed volume. Needle length was highly correlated to injection force (R 2 = 0.99). Using longer needles to inject drug solution requires greater force. Furthermore, large barrel size was correlated to pushing force (R 2 = 0.99); thus, increasing the difficulty of pushing the plunger of a large syringe with one hand. Lastly, relationships between injection force, sample viscosity, and testing conditions were elucidated using a mathematical model, which could be used in the future to assess and predict injection force of solution samples.

60 APPLIED LIFE SCIENCES↗

An analytic model for buoyancy resonances in protoplanetary disks

Zhu et al. found in three-dimensional shearing box simulations a new form of planet-disk interaction that they attributed to a vertical buoyancy resonance in the disk. We describe an analytic linear model for this interaction. We adopt a simplified model involving azimuthal forcing that produces the resonance and permits an analytic description of its structure. We derive an analytic expression for the buoyancy torque and show that the vertical torque distribution agrees well with the results of the Athena simulations and a Fourier method for linear numerical calculations carried out with the same forcing. The buoyancy resonance differs from the classic Lindblad and corotation resonances in that the resonance lies along tilted planes. Its width depends on damping effects and is independent of the gas sound speed. The resonance does not excite propagating waves. At a given large azimuthal wavenumber k{sub y} > h {sup –1} (for disk thickness h), the buoyancy resonance exerts a torque over a region that lies radially closer to the corotation radius than the Lindblad resonance. Because the torque is localized to the region of excitation, it is potentially subject to the effects of nonlinear saturation. In addition, the torque can be reduced by the effects of radiative heat transfer between the resonant region and its surroundings. For each azimuthal wavenumber, the resonance establishes a large scale density wave pattern in a plane within the disk.

79 ASTRONOMY AND ASTROPHYSICS↗

FORCE-DISPATCHES Integration - Initial Demonstration

Integrated energy systems (IES) combine, in mutually beneficial ways, power from variable renewable energy sources and nuclear power plants (NPP) to improve economic viability under uncertain market and weather conditions. The open-source Framework for Optimization of Resources and Economics (FORCE) tool suite, developed at Idaho National Laboratory (INL), has enabled comprehensive modeling and simulation of IES. The capabilities within FORCE include grid portfolio optimization through the Holistic Energy Resource Optimization Network (HERON) and the transient process model analysis library HYBRID, among others. Continuous efforts and investments from the IES programs have been made to expand and improve the versatility of the FORCE toolset in fiscal year 2022. Code-coupling and cross-tool communication have been important methods for improving this versatility. This report focuses on an additional workflow in the HERON tool for capacity and dispatch stochastic optimization through integration with the external tool Design Integration and Synthesis Platform to Advance Tightly Coupled Hybrid Energy Systems (DISPATCHES). DISPATCHES was primarily developed by the National Energy Technology Laboratory, in collaboration with other national laboratories, which included INL, universities, and industry partners. It is coupled to a library of algebraic models for specific plant components, and to a framework for stochastic optimization different from that provided in the current Risk Analysis Virtual Environment (RAVEN)-running-RAVEN algorithm in HERON. HERON currently conducts stochastic optimization via an outer-inner loop: it optimizes over variable capacity on the outer loop, and at each step within the capacity parameter space, conducts an inner optimization over scenarios (of market signals, demand, and/or weather patterns) and hourly dispatch throughout a user-specified number of years. On the other hand, DISPATCHES conducts stochastic optimization via an “all-at-once” strategy in which capacity variables are optimized at the same level as dispatch variables, as all scenarios are considered at once. The latter method works especially well for projects of limited size and project length, as the necessary computational power and memory increases with the number of variables and scenarios. The new capability to use the DISPATCHES workflow in HERON enhances standalone simulations by leveraging FORCE tools—namely, the economic metrics from the Tool for Economic Analysis (TEAL) and reduced-order model (ROM) sampling from RAVEN. The initial demonstration of the DISPATCHES workflow simulates an existing nuclear-case flowsheet within the DISPATCHES repository—this models a NPP with a secondary revenue stream for hydrogen production. Electrical output from the plant is converted to hydrogen via a proton-exchange membrane (PEM) electrolyzer, hydrogen tanks are used for storage, and an additional turbine is added for hydrogen combustion. Continued work regarding this FORCE-DISPATCHES integration will include automatic generation of DISPATCHES models from HERON inputs, offering analysts the option of using either the RAVEN-runsRAVEN or DISPATCHES workflow to solve technoeconomic optimization problems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Reactive molecular dynamics models from ab initio molecular dynamics data using relative entropy minimization

Here, we present two new multiscale molecular dynamics (MS-RMD) models for the hydrated excess proton in water developed directly from ab initio molecular dynamics (AIMD) simulation data of the same system. The potential of mean force along the proton transfer reaction coordinate and radial distribution functions for the MS-RMD models are shown faithfully reproduce those of AIMD. The models are developed using an algorithm based on relative entropy minimization, thus demonstrating the ability of the method to rapidly generate accurate and highly efficient reactive MD force fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Growth of strontium ruthenate films by hybrid molecular beam epitaxy

We report on the growth of epitaxial Sr 2 RuO 4 films using a hybrid molecular beam epitaxy approach in which a volatile precursor containing RuO 4 is used to supply ruthenium and oxygen. The use of the precursor overcomes a number of issues encountered in traditional molecular beam epitaxy that uses elemental metal sources. Phase-pure, epitaxial thin films of Sr 2 RuO 4 are obtained. At high substrate temperatures, growth proceeds in a layer-by-layer mode with intensity oscillations observed in reflection high-energy electron diffraction. Films are of high structural quality, as documented by x-ray diffraction, atomic force microscopy, and transmission electron microscopy. In conclusion, the method should be suitable for the growth of other complex oxides containing ruthenium, opening up opportunities to investigate thin films that host rich exotic ground states.

36 MATERIALS SCIENCE↗

Robust global detection of forced changes in mean and extreme precipitation despite observational disagreement on the magnitude of change

Detection and attribution (D&A) of forced precipitation change are challenging due to internal variability, limited spatial, and temporal coverage of observational records and model uncertainty. These factors result in a low signal-to-noise ratio of potential regional and even global trends. Here, we use a statistical method – ridge regression – to create physically interpretable fingerprints for the detection of forced changes in mean and extreme precipitation with a high signal-to-noise ratio. The fingerprints are constructed using Coupled Model Intercomparison Project phase 6 (CMIP6) multi-model output masked to match coverage of three gridded precipitation observational datasets – GHCNDEX, HadEX3, and GPCC – and are then applied to these observational datasets to assess the degree of forced change detectable in the real-world climate in the period 1951–2020. We show that the signature of forced change is detected in all three observational datasets for global metrics of mean and extreme precipitation. Forced changes are still detectable from changes in the spatial patterns of precipitation even if the global mean trend is removed from the data. This shows the detection of forced change in mean and extreme precipitation beyond a global mean trend is robust and increases confidence in the detection method's power as well as in climate models' ability to capture the relevant processes that contribute to large-scale patterns of change. We also find, however, that detectability depends on the observational dataset used. Not only coverage differences but also observational uncertainty contribute to dataset disagreement, exemplified by the times of emergence of forced change from internal variability ranging from 1998 to 2004 among datasets. Furthermore, different choices for the period over which the forced trend is computed result in different levels of agreement between observations and model projections. These sensitivities may explain apparent contradictions in recent studies on whether models under- or overestimate the observed forced increase in mean and extreme precipitation. Lastly, the detection fingerprints are found to rely primarily on the signal in the extratropical Northern Hemisphere, which is at least partly due to observational coverage but potentially also due to the presence of a more robust signal in the Northern Hemisphere in general.

54 ENVIRONMENTAL SCIENCES↗

Quantifying the Sources of Intermodel Spread in Equilibrium Climate Sensitivity

This paper clarifies the causes of intermodel differences in the global-average temperature response to doubled CO 2 , commonly known as equilibrium climate sensitivity (ECS). The authors begin by noting several issues with the standard approach for decomposing ECS into a sum of forcing and feedback terms. This leads to a derivation of an alternative method based on linearizing the effect of the net feedback. Consistent with previous studies, the new method identifies shortwave cloud feedback as the dominant source of intermodel spread in ECS. This new approach also reveals that covariances between cloud feedback and forcing, between lapse rate and longwave cloud feedbacks, and between albedo and shortwave cloud feedbacks play an important and previously underappreciated role in determining model differences in ECS. Finally, defining feedbacks based on fixed relative rather than specific humidity (as suggested by Held and Shell) reduces the covariances between processes and leads to more straightforward interpretations of results.

54 ENVIRONMENTAL SCIENCES↗

Radiation forces and torque on a rigid elliptical cylinder in acoustical plane progressive and (quasi)standing waves with arbitrary incidence

This paper presents two key contributions; the first concerns the development of analytical expressions for the axial and transverse acoustic radiation forces exerted on a 2D rigid elliptical cylinder placed in the field of plane progressive, quasi-standing, or standing waves with arbitrary incidence. The second emphasis is on the acoustic radiation torque per length. The rigid elliptical cylinder case is important to be considered as a first-order approximation of the behavior of a cylindrical fluid column trapped in air because of the significant acoustic impedance mismatch at the particle boundary. Based on the rigorous partial-wave series expansion method in cylindrical coordinates, non-dimensional acoustic radiation force and torque functions are derived and defined in terms of the scattering coefficients of the elliptic cylinder. A coupled system of linear equations is obtained after applying the Neumann boundary condition for an immovable surface in a non-viscous fluid and solved numerically by matrix inversion after performing a single numerical integration procedure. Computational results for the non-dimensional force components and torque, showing the transition from the progressive to the (equi-amplitude) standing wave behavior, are performed with particular emphasis on the aspect ratio a/b, where a and b are the semi-axes of the ellipse, the dimensionless size parameter, as well as the angle of incidence ranging from end-on to broadside incidence. The results show that the elliptical geometry has a direct influence on the radiation force and torque, so that the standard theory for circular cylinders (at normal incidence) leads to significant miscalculations when the cylinder cross section becomes non-circular. Moreover, the elliptical cylinder experiences, in addition to the acoustic radiation force, a radiation torque that vanishes for the circular cylinder case. The application of the formalism presented here may be extended to other 2D surfaces of arbitrary shape, such as Chebyshev cylindrical particles with a small deformation, stadiums (with oval shape), or other non-circular geometries.

42 ENGINEERING↗

Methods for Estimating 2D Cloud Size Distributions from 1D Observations

Abstract The two-dimensional (2D) size distribution of clouds in the horizontal plane plays a central role in the calculation of cloud cover, cloud radiative forcing, convective entrainment rates, and the likelihood of precipitation. Here, a simple method is proposed for calculating the area-weighted mean cloud size and for approximating the 2D size distribution from the 1D cloud-chord lengths measured by aircraft and vertically pointing lidar and radar. This simple method (which is exact for square clouds) compares favorably against the inverse Abel transform (which is exact for circular clouds) in the context of theoretical size distributions. Both methods also perform well when used to predict the size distribution of real clouds from a Landsat scene. When applied to a large number of Landsat scenes, the simple method is able to accurately estimate the mean cloud size. As a demonstration, the methods are applied to aircraft measurements of shallow cumuli during the Routine ARM Aerial Facility (AAF) Clouds with Low Optical Water Depths (CLOWD) Optical Radiative Observations (RACORO) campaign, which then allow for an estimate of the true area-weighted mean cloud size.

54 ENVIRONMENTAL SCIENCES↗

Impact assessment of land use planning driving forces on environment

Land use change may exert a negative impact on environmental quality. A state–impact–state (SIS) model describing a state transform under certain impacts has been integrated into land use planning (LUP) environmental impact assessment (LUPEA). This logical model is intuitive and easy to understand, but the exploration of impact is essential to establish the indicator system and to identify the scope of land use environmental impact when it is applied to a specific region. In this study, we investigated environmental driving forces from land use planning (LUPF), along with the conception, components, scope, and impact of LUPF. This method was illustrated by a case study in Zoucheng, China. Through the results, we concluded that (1) the LUPF on environment are impacts originated from the implementation of LUP on a regional environment, which are characterized by four aspects: magnitude, direction, action point, and its owner; (2) various scopes of LUPF on individual environmental elements based on different standards jointly define the final scope of LUPEA; (3) our case study in Zoucheng demonstrates the practicability of this proposed approach; (4) this method can be embedded into LUPEA with direction, magnitudes, and scopes of the LUPF on individual elements obtained, and the identified indicator system can be directly employed into LUPEA and (5) the assessment helps to identify key indicators and to set up a corresponding strategy to mitigate the negative impact of LUP on the environment, which are two important objectives of strategic environmental assessment (SEA) in LUP. - Highlights: • Environmental driving forces from land use planning (LUPF) are investigated and categorized. • Our method can obtains the direction, magnitudes and scopes of environmental driving forces. • The LUPEA scope is determined by the combination of various scopes of LUPF on individual elements. • LUPF assessment can be embedded into LUPEA. • The method can help to identify key indicators and set up a strategy to mitigate negative environmental impact.

54 ENVIRONMENTAL SCIENCES↗

Poster - 32: Atlas Selection for Automated Segmentation of Pelvic CT for Prostate Radiotherapy

Atlas based-segmentation has recently been evaluated for use in prostate radiotherapy. In a typical approach, the essential step is the selection of an atlas from a database that the best matches of the target image. This work proposes an atlas selection strategy and evaluate it impacts on final segmentation accuracy. Several anatomical parameters were measured to indicate the overall prostate and body shape, all of these measurements obtained on CT images. A brute force procedure was first performed for a training dataset of 20 patients using image registration to pair subject with similar contours; each subject was served as a target image to which all reaming 19 images were affinity registered. The overlap between the prostate and femoral heads was quantified for each pair using the Dice Similarity Coefficient (DSC). Finally, an atlas selection procedure was designed; relying on the computation of a similarity score defined as a weighted sum of differences between the target and atlas subject anatomical measurement. The algorithm ability to predict the most similar atlas was excellent, achieving mean DSCs of 0.78 ± 0.07 and 0.90 ± 0.02 for the CTV and either femoral head. The proposed atlas selection yielded 0.72 ± 0.11 and 0.87 ± 0.03 for CTV and either femoral head. The DSC obtained with the proposed selection method were slightly lower than the maximum established using brute force, but this does not include potential improvements expected with deformable registration. The proposed atlas selection method provides reasonable segmentation accuracy.

60 APPLIED LIFE SCIENCES↗

An augmented Lagrangian filter method

Here, we introduce a filter mechanism to enforce convergence for augmented Lagrangian methods for nonlinear programming. In contrast to traditional augmented Lagrangian methods, our approach does not require the use of forcing sequences that drive the first-order error to zero. Instead, we employ a filter to drive the optimality measures to zero. Our algorithm is flexible in the sense that it allows for equality-constrained quadratic programming steps to accelerate local convergence. We also include a feasibility restoration phase that allows fast detection of infeasible problems. We provide a convergence proof that shows that our algorithm converges to first-order stationary points. We provide preliminary numerical results that demonstrate the effectiveness of our proposed method.

97 MATHEMATICS AND COMPUTING↗

Impact of Short-Range Forces on Defect Production from High-Energy Collisions

Primary radiation damage formation in solid materials typically involves collisions between atoms that have up to a few hundred keV of kinetic energy. The distance between two colliding atoms can approach 0.05 nm during these collisions. At such small atomic separations, force fields fitted to equilibrium properties tend to significantly underestimate the potential energy of the colliding dimer. To enable molecular dynamics simulations of high-energy collisions, it is common practice to use a screened Coulomb force field to describe the interactions and to smoothly join this to the equilibrium force field at a suitable interatomic spacing. But, there is no accepted standard method for choosing the parameters used in the joining process, and our results prove that defect production is sensitive to how the force fields are linked. A new procedure is presented that involves the use of ab initio calculations to determine the magnitude and spatial dependence of the pair interactions at intermediate distances, along with systematic criteria for choosing the joining parameters. Results are presented for the case of nickel, which demonstrate the use and validity of the procedure.

74 ATOMIC AND MOLECULAR PHYSICS↗

Petaloid-shaped ZnO coated carbon felt as a controllable host to construct hierarchical Li composite anode

Prestoring Li metal into the 3D host via thermal infusion strategy has been considered as an effective approach on suppressing Li dendrites growth of Li-metal anode during cycling. However, the fast Li infusing process generally induces uncontrollable and undesirable structures during the formation of composite Li metal anode. It is widely accepted that both physical and chemical interactions between the molten Li and the hosts provide driving forces for Li infusing into the host, but which one plays the dominant role has not been explored clearly. In this work, a 3D carbon felt coated by petaloid zinc oxide (CFZO) is designed as a controllable host for prestoring Li metal via thermal infusion method. It is revealed that the morphology induced physical driving force plays a more important role than the chemical reaction between ZnO and molten Li in the thermal infusion process. The physical driving force originates from the capillarity effect, depending on the roughness of the host according to Wenzel's model. An adjustable Li infusing process can be controlled to fabricate the Li composite electrode with a hierarchical structure (with dense Li-rich layer and porous Li-poor layer). The unique structure of the CFZO-Li not only provides enough spaces for Li deposition, but also disperses local current density to induce a uniform Li deposition behavior, thus Li dendrites can be suppressed effectively. As a result, CFZO-Li anode exhibits promising cycle stability when coupled with both Li-containing (LiNi 0.8 Co 0.15 Al 0.05 O 2 ) and Li-free (sulfur) cathodes.

25 ENERGY STORAGE↗

Application of small-angle neutron scattering to the study of forces between magnetically chained monodisperse ferrofluid emulsion droplets

The optical magnetic chaining technique (MCT) developed by Leal-Calderon, Bibette and co-workers in the 1990 s allows precise measurements of force profiles between droplets in monodisperse ferrofluid emulsions. However, the method lacks an in-situ determination of droplet size and therefore requires the combination of separately acquired measurements of droplet chain periodicity versus an applied magnetic field from optical Bragg scattering and droplet diameter inferred from dynamic light scattering (DLS) to recover surface force-distance profiles between the colloidal particles. Compound refractive lens (CRL) focussed small-angle scattering (SANS) MCT should result in more consistent measurements of droplet size (form factor measurements in the absence of field) and droplet chaining period (from structure factor peaks when the magnetic field is applied); and, with access to shorter length scales, extend force measurements to closer approaches than possible by optical measurements. We report on CRL-SANS measurements of monodisperse ferrofluid emulsion droplets aligned in straight chains by an applied field perpendicular to the incident beam direction. Analysis of the scattering from the closely spaced droplets required algorithms that carefully treated resolution and its effect on mean scattering vector magnitudes in order to determine droplet size and chain periods to sufficient accuracy. At lower applied fields scattering patterns indicate structural correlations transverse to the magnetic field direction due to the formation of intermediate structures in early chain growth.

Colloid↗