Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generation equation”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

OpenSesame tutorial

OpenSesame is a program for generating tabular equations of state (EOS), with capabilities for multiphase EOS construction. In this tutorial, we provide an overview of how to run OpenSesame to construct a multiphase EOS. We discuss some general features of OpenSesame, followed by a description of sample input files required for multiphase EOS construction. We also discuss how to extract data from EOS tables in order to compare to experimental data, with an example using the OpenSesame GUI. Lastly, we provide a description of how to generate ASCII-formatted EOS tables most often used by hydro code users.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Confining but chirally symmetric dense and cold matter

The possibility for existence of cold, dense chirally symmetric matter with confinement is reviewed. The answer to this question crucially depends on the mechanism of mass generation in QCD and interconnection of confinement and chiral symmetry breaking. This question can be clarified from spectroscopy of hadrons and their axial properties. Almost systematical parity doubling of highly excited hadrons suggests that their mass is not related to chiral symmetry breaking in the vacuum and is approximately chirally symmetric. Then there is a possibility for existence of confining but chirally symmetric matter. We clarify a possible mechanism underlying such a phase at low temperatures and large density. Namely, at large density the Pauli blocking prevents the gap equation to generate a solution with broken chiral symmetry. However, the chirally symmetric part of the quark Green function as well as all color non-singlet quantities are still infrared divergent, meaning that the system is with confinement. A possible phase transition to such a matter is most probably of the first order. This is because there are no chiral partners to the lowest lying hadrons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1990

Trees, bialgebras and intrinsic numerical algorithms

Preliminary work about intrinsic numerical integrators evolving on groups is described. Fix a finite dimensional Lie group G; let g denote its Lie algebra, and let Y(sub 1),...,Y(sub N) denote a basis of g. A class of numerical algorithms is presented that approximate solutions to differential equations evolving on G of the form: dot-x(t) = F(x(t)), x(0) = p is an element of G. The algorithms depend upon constants c(sub i) and c(sub ij), for i = 1,...,k and j is less than i. The algorithms have the property that if the algorithm starts on the group, then it remains on the group. In addition, they also have the property that if G is the abelian group R(N), then the algorithm becomes the classical Runge-Kutta algorithm. The Cayley algebra generated by labeled, ordered trees is used to generate the equations that the coefficients c(sub i) and c(sub ij) must satisfy in order for the algorithm to yield an rth order numerical integrator and to analyze the resulting algorithms.

Crouch, Peter↗
NASA NTRS2006

Intra-seasonal Oscillations Inferred from SABER (TIMED) and MLS (UARS) Temperature Measurements

In the zonal mean meridional winds of the upper mesosphere, intra-seasonal oscillations with periods between 1 and 4 months have been inferred from UARS measurements and independently predicted with the Numerical Spectral Model WSM). The wind oscillations tend to be confined to low latitudes and appear to be driven, at least in part, by small-scale gravity waves propagating in the meridional direction. Winds across the equator should generate, due to dynamical heating and cooling, temperature oscillations with opposite phase in the two hemispheres. Investigating this phenomenon, we have analyzed SABER temperatures from TIMED in the altitude range between 55 and 95 km to delineate with an empirical model, the year-long variability of the migrating tides and zonal mean components. The inferred seasonal variations of the diurnal tide, characterized by amplitude maxima near equinox, are in substantial agreement with UARS observations and results from the NSM. For the zonal mean, the dominant seasonal variations in the SABER temperatures, with annual (12 months) and semiannual (6 months) periodicities, agree well with those derived from UARS measurements. The intra-seasonal variations with periods between 2 and 4 months have amplitudes close to 2 K, almost half as large as those for the dominant seasonal variations. Their amplitudes are in qualitative agreement with the corresponding values inferred from UARS during different years. The SABER and UARS temperature variations reveal pronounced hemispherical asymmetries, consistent with meridional wind oscillations across the equator. The phase of the semi-annual temperature oscillations from the NSM agrees with the observations from UARS and SABER. But the amplitudes are systematically smaller, which may indicate that planetary waves are more important than is allowed for in the model. For the shorter-period intra-seasonal variations, which can be generated by gravity wave drag, the model results are generally in better agreement with the observations.

Huang, F. T.↗
DOE OSTI2014

WE-E-17A-07: Patient-Specific Mathematical Neuro-Oncology: Biologically-Informed Radiation Therapy and Imaging Physics

Purpose: To demonstrate a method of generating patient-specific, biologically-guided radiation therapy (RT) plans and to quantify and predict response to RT in glioblastoma. We investigate the biological correlates and imaging physics driving T2-MRI based response to radiation therapy using an MRI simulator. Methods: We have integrated a patient-specific biomathematical model of glioblastoma proliferation, invasion and radiotherapy with a multiobjective evolutionary algorithm for intensity-modulated RT optimization to construct individualized, biologically-guided plans. Patient-individualized simulations of the standard-of-care and optimized plans are compared in terms of several biological metrics quantified on MRI. An extension of the PI model is used to investigate the role of angiogenesis and its correlates in glioma response to therapy with the Proliferation-Invasion-Hypoxia- Necrosis-Angiogenesis model (PIHNA). The PIHNA model is used with a brain tissue phantom to predict tumor-induced vasogenic edema, tumor and tissue density that is used in a multi-compartmental MRI signal equation for generation of simulated T2- weighted MRIs. Results: Applying a novel metric of treatment response (Days Gained) to the patient-individualized simulation results predicted that the optimized RT plans would have a significant impact on delaying tumor progression, with Days Gained increases from 21% to 105%. For the T2- MRI simulations, initial validation tests compared average simulated T2 values for white matter, tumor, and peripheral edema to values cited in the literature. Simulated results closely match the characteristic T2 value for each tissue. Conclusion: Patient-individualized simulations using the combination of a biomathematical model with an optimization algorithm for RT generated biologically-guided doses that decreased normal tissue dose and increased therapeutic ratio with the potential to improve survival outcomes for treatment of glioblastoma. Simulated T2-MRI is shown to be consistent with known physics of MRI and can be used to further investigate biological drivers of imaging-based response to RT.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Characterization of scatter in digital mammography from physical measurements

Purpose: That scattered radiation negatively impacts the quality of medical radiographic imaging is well known. In mammography, even slight amounts of scatter reduce the high contrast required for subtle soft-tissue imaging. In current clinical mammography, image contrast is partially improved by use of an antiscatter grid. This form of scatter rejection comes with a sizeable dose penalty related to the concomitant elimination of valuable primary radiation. Digital mammography allows the use of image processing as a method of scatter correction that might avoid effects that negatively impact primary radiation, while potentially providing more contrast improvement than is currently possible with a grid. For this approach to be feasible, a detailed characterization of the scatter is needed. Previous research has modeled scatter as a constant background that serves as a DC bias across the imaging surface. The goal of this study was to provide a more substantive data set for characterizing the spatially-variant features of scatter radiation at the image detector of modern mammography units. Methods: This data set was acquired from a model of the radiation beam as a matrix of very narrow rays or pencil beams. As each pencil beam penetrates tissue, the pencil widens in a predictable manner due to the production of scatter. The resultant spreading of the pencil beam at the detector surface can be characterized by two parameters: mean radial extent (MRE) and scatter fraction (SF). The SF and MRE were calculated from measurements obtained using the beam stop method. Two digital mammography units were utilized, and the SF and MRE were found as functions of target, filter, tube potential, phantom thickness, and presence or absence of a grid. These values were then used to generate general equations allowing the SF and MRE to be calculated for any combination of the above parameters. Results: With a grid, the SF ranged from a minimum of about 0.05 to a maximum of about 0.16, and the MRE ranged from about 3 to 13 mm. Without a grid, the SF ranged from a minimum of 0.25 to a maximum of 0.52, and the MRE ranged from about 20 to 45 mm. The SF with a grid demonstrated a mild dependence on target/filter combination and kV, whereas the SF without a grid was independent of these factors. The MRE demonstrated a complex relationship as a function of kV, with notable difference among target/filter combinations. The primary source of change in both the SF and MRE was phantom thickness. Conclusions: Because breast tissue varies spatially in physical density and elemental content, the effective thickness of breast tissue varies spatially across the imaging field, resulting in a spatially-variant scatter distribution in the imaging field. The data generated in this study can be used to characterize the scatter contribution on a point-by-point basis, for a variety of different techniques.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS1981

On numerical nonlinear analysis of highly flexible spinning cantilevers

The general nonlinear discretized equations of motion of spinning elastic solids and structures are derived as a set of nonlinear ordinary differential equations for the case when the strain-displacement and velocity-displacement relations are nonlinear up to the second order. It is shown that the cost of generation of such equations is proportional to the fourth power of the number of degrees of freedom. A computer program is written to automatically generate the equations for the case of spinning cantilevers with initial imperfections. The types and the number of the coordinate functions used in the trial solution are parameters of the program.

Utku, S.↗
NASA NTRS1988

Grid generation for general 2-D regions using hyperbolic equations

A method for applying a hyperbolic grid generation scheme to the construction of meshes in general 2-D regions has been developed. This approach, which follows the theory developed by Steger and Chaussee (1980) and the algorithm outlined by Kinsey and Barth (1984), is based on improving local grid control. This is accomplished by adding an angle control source term to the equations and using a new algorithm for computing the volume source term. These modifications lead to superior methods for fixing the 'local' problems of hyperbolic grid generation, namely, propagation of initial discontinuities and formation of grid shocks (crossing grid lines). More importantly, a method for solving the global problem of constraining the grid with more than one boundary (internal grid generation) has been developed. These algorithms have been implemented in an interactive grid generation program and the results for several geometries are presented and discussed.

Cordova, Jeffrey Q.↗
NASA NTRS2012

Evaluation of Delamination Growth Characterization Methods Under Mode I Fatigue Loading

Reliable delamination characterization data for laminated composites are needed for input to analytical models of structures to predict delamination. The double-cantilevered beam (DCB) specimen is used with laminated composites to measure fracture toughness, G(sub Ic), delamination onset strain energy release rate, and growth rate data under cyclic loading. In the current study, DCB specimens of IM7/8552 graphite/epoxy supplied by two different manufacturers were tested in static and fatigue to compare the measured characterization data from the two sources, and to evaluate a proposed ASTM standard for generating Paris Law equations. Static results were used to generate compliance calibration constants for the fatigue data, and a delamination resistance curve, G(sub IR), which was used to determine the effects of fiber-bridging on delamination growth. Specimens were tested in fatigue at a cyclic G(sub Imax) level equal to 50, 40 or 30% of G(sub Ic), to determine a delamination onset curve and delamination growth rate. The delamination onset curve equations had similar exponents and the same trends. Delamination growth rate was calculated by fitting a Paris Law to the da/dN versus G(sub Imax) data. Both a 2-point and a 7-point data reduction method were used and the Paris Law equations were compared. To determine the effects of fiber-bridging, growth rate results were normalized by the delamination resistance curve for each material and compared to the non-normalized results. Paris Law exponents were found to decrease by 31% to 37% due to normalizing the growth data. Normalizing the data also greatly reduced the amount of scatter between the different specimens. Visual data records from the fatigue testing were used to calculate individual compliance calibration constants from the fatigue data for some of the specimens. The resulting da/dN versus G(sub Imax) plots showed much improved repeatability between specimens. Gretchen

Murri, Gretchen B.↗
DOE OSTI2014

The Milky Way tomography with Sloan Digital Sky Survey. V. Mapping the dark matter halo

We present robust constraints from the Sloan Digital Sky Survey (SDSS) on the shape and distribution of the dark matter halo within the Milky Way (MW). Using the number density distribution and kinematics of SDSS halo stars, we probe the dark matter distribution to heliocentric distances exceeding ∼10 kpc and galactocentric distances exceeding ∼20 kpc. Our analysis utilizes Jeans equations to generate two-dimensional acceleration maps throughout the volume; this approach is thoroughly tested on a cosmologically derived N-body+SPH simulation of a MW-like galaxy. We show that the known accelerations (gradients of the gravitational potential) can be successfully recovered in such a realistic system. Leveraging the baryonic gravitational potential derived by Bovy and Rix, we show that the gravitational potential implied by the SDSS observations cannot be explained, assuming Newtonian gravity, by visible matter alone: the gravitational force experienced by stars at galactocentric distances of ∼20 kpc is as much as three times stronger than what can be attributed to purely visible matter. We also show that the SDSS data provide a strong constraint on the shape of the dark matter halo potential. Within galactocentric distances of ∼20 kpc, the dark matter halo potential is well described as an oblate halo with axis ratio q{sub DM}{sup Φ}=0.7±0.1; this corresponds to an axis ratio q{sub DM}{sup ρ}∼0.4±0.1 for the dark matter density distribution. Because of our precise two-dimensional measurements of the acceleration of the halo stars, we can reject several MOND models as an explanation of the observed behavior.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Probabilistic learning and updating of a digital twin for composite material systems

This paper presents an approach for characterizing and estimating statistical dependence between a large number of observables in a composite material system. Conditional regression is carried out using the estimated joint density function, permitting a systematic exploration of interdependence between fine scale and coarse observables that can be used for both prognosis and design of complex material systems. An example demonstrates the integration of experimental data with a computational database. The statistical approach is based on the probabilistic learning on manifolds recently developed by the authors. This approach leverages intrinsic structure detected through diffusion on graphs with projected stochastic differential equations to generate samples constrained to that structure.

42 ENGINEERING↗
DOE OSTI2026

Idealized Offshore Low‐Level Jets for Turbine Structural Impact Considerations

Low-level jets (LLJs) describe conditions in which the wind speed reaches a local maximum with respect to altitude near the surface; they have been observed intermittently in the US mid-Atlantic offshore environment. LLJs pose unique operating conditions for future wind turbines operating in the region, such as negative shear and locally strong veer, but they are not typically considered in existing turbine standards. This work builds upon recent research that explains the formation and evolution of the US mid-Atlantic LLJs through a simple analytical governing equation. We generate several LLJ inflow conditions with varying jet characteristics based on this analytical model and create monotonically sheared (MS) analogues with matched veer to assess the impacts of the LLJ on turbine performance and loading. Using aeroelastic simulations with these inflow conditions on the International Energy Agency 15-MW reference turbine, we find that the LLJ leads to a greater range of tower-top pitching and yawing moments, which could contribute to larger accumulated structural fatigue in components compared to MS inflow. These preliminary results demonstrate a path toward a unified set of test cases for low-level wind maxima that can inform the International Electrotechnical Commission standards related to offshore wind turbine design.

17 WIND ENERGY↗
DOE OSTI2017

Morse-Smale Analysis of Ion Diffusion in Ab Initio Battery Materials Simulations

Ab initio molecular dynamics (AIMD) simulations are increasingly useful in modeling, optimizing and synthesizing materials in energy sciences. In solving Schrödinger’s equation, they generate the electronic structure of the simulated atoms as a scalar field. However, methods for analyzing these volume data are not yet common in molecular visualization. The Morse-Smale complex is a proven, versatile tool for topological analysis of scalar fields. In this paper, we apply the discrete Morse-Smale complex to analysis of first-principles battery materials simulations. We consider a carbon nanosphere structure used in battery materials research, and employ Morse-Smale decomposition to determine the possible lithium ion diffusion paths within that structure. Our approach is novel in that it uses the wavefunction itself as opposed distance fields, and that we analyze the 1-skeleton of the Morse-Smale complex to reconstruct our diffusion paths. Furthermore, it is the first application where specific motifs in the graph structure of the complete 1-skeleton define features, namely carbon rings with specific valence. We compare our analysis of DFT data with that of a distance field approximation, and discuss implications on larger classical molecular dynamics simulations.

25 ENERGY STORAGE↗
DOE OSTI2012

Gauge invariance for a whole Abelian model

Light invariance is a fundamental principle for physics be done. It generates Maxwell equations, relativity, Lorentz group. However there is still space for a fourth picture be developed which is to include fields with same Lorentz nature. It brings a new room for field theory. It says that light invariance does not work just to connect space and time but it also associates different fields with same nature. Thus for the ((1/2),(1/2)) representation there is a fields family {l_brace}A{sub {mu}I}{r_brace} to be studied. This means that given such fields association one should derive its corresponding gauge theory. This is the effort at this work. Show that there is a whole gauge theory to cover these fields relationships. Considering the abelian case, prove its gauge invariance. It yields the kinetic, massive, trilinear and quadrilinear gauge invariant terms.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2011

TURBULENT CROSS-HELICITY IN THE MEAN-FIELD SOLAR DYNAMO PROBLEM

We study the dynamical and statistical properties of turbulent cross-helicity (correlation of the aligned fluctuating velocity and magnetic field components). We derive an equation governing generation and evolution of the turbulent cross-helicity and discuss its meaning for the dynamo. Using the symmetry properties of the problem we suggest a general expression for the turbulent cross-helicity. Effects of the density stratification, large-scale magnetic fields, differential rotation, and turbulent convection are taken into account. We investigate the relative contribution of these effects to the cross-helicity evolution for two kinds of dynamo models of the solar cycle: a distributed mean-field model and a flux-transport dynamo model. We show that the contribution from the density stratification follows the evolution of the radial magnetic field, while large-scale electric currents produce a more complicated pattern of the cross-helicity of comparable magnitude. The pattern of the cross-helicity evolution strongly depends on details of the dynamo mechanism. Thus, we anticipate that direct observations of the cross-helicity on the Sun may serve for the diagnostic purpose of the solar dynamo process.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Nonlinear equation for Farley–Buneman waves in multispecies plasma

The nonlinear equation describing the Farley–Buneman (FB) waves in multispecies collisional plasmas is derived by employing the multiple-scale reduction analysis. It is shown that the presence of several ion species with different collisionalities and different ion masses removes the degeneracy of the nonlinear equation and generates the nonlinear terms resulting in wave steepening and wave breaking. This effect may be responsible for formation of one-dimensional coherent FB waves of a finite amplitude.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A Finite-Time Ensemble Method for Mixed Layer Model Comparison

Here, this work evaluates the fidelity of various upper-ocean turbulence parameterizations subject to realistic monsoon forcing and presents a finite-time ensemble vector (EV) method to better manage the design and numerical principles of these parameterizations. The EV method emphasizes the dynamics of a turbulence closure multimodel ensemble and is applied to evaluate 10 different ocean surface boundary layer (OSBL) parameterizations within a single-column (SC) model against two boundary layer large-eddy simulations (LES). Both LES include realistic surface forcing, but one includes wind-driven shear turbulence only, while the other includes additional Stokes forcing through the wave-average equations that generate Langmuir turbulence. The finite-time EV framework focuses on what constitutes the local behavior of the mixed layer dynamical system and isolates the forcing and ocean state conditions where turbulence parameterizations most disagree. Identifying disagreement provides the potential to evaluate SC models comparatively against the LES. Observations collected during the 2018 monsoon onset in the Bay of Bengal provide a case study to evaluate models under realistic and variable forcing conditions. The case study results highlight two regimes where models disagree 1) during wind-driven deepening of the mixed layer and 2) under strong diurnal forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

A finalized determinant for complete lignocellulose enzymatic saccharification potential to maximize bioethanol production in bioenergy Miscanthus

Miscanthus is a leading bioenergy crop with enormous lignocellulose production potential for biofuels and chemicals. However, lignocellulose recalcitrance leads to biomass process difficulty for an efficient bioethanol production. Hence, it becomes essential to identify the integrative impact of lignocellulose recalcitrant factors on cellulose accessibility for biomass enzymatic hydrolysis. In this study, we analyzed four typical pairs of Miscanthus accessions that showed distinct cell wall compositions and sorted out three major factors that affected biomass saccharification for maximum bioethanol production. Among the three optimal (i.e., liquid hot water, H 2 SO 4 and NaOH) pretreatments performed, mild alkali pretreatment (4% NaOH at 50 °C) led to almost complete biomass saccharification when 1% Tween-80 was co-supplied into enzymatic hydrolysis in the desirable Miscanthus accessions. Consequently, the highest bioethanol yields were obtained at 19% (% dry matter) from yeast fermentation, with much higher sugar–ethanol conversion rates by 94–98%, compared to the other Miscanthus species subjected to stronger pretreatments as reported in previous studies. By comparison, three optimized pretreatments distinctively extracted wall polymers and specifically altered polymer features and inter-linkage styles, but the alkali pretreatment caused much increased biomass porosity than that of the other pretreatments. Based on integrative analyses, excellent equations were generated to precisely estimate hexoses and ethanol yields under various pretreatments and a hypothetical model was proposed to outline an integrative impact on biomass saccharification and bioethanol production subjective to a predominate factor (CR stain) of biomass porosity and four additional minor factors (DY stain, cellulose DP, hemicellulose X/A, lignin G-monomer). Using four pairs of Miscanthus samples with distinct cell wall composition and varied biomass saccharification, this study has determined three main factors of lignocellulose recalcitrance that could be significantly reduced for much-increased biomass porosity upon optimal pretreatments. It has also established a novel standard that should be applicable to judge any types of biomass process technology for high biofuel production in distinct lignocellulose substrates. Hence, this study provides a potential strategy for precise genetic modification of lignocellulose in all bioenergy crops.

09 BIOMASS FUELS↗