SEARCH · Engineering Papers

Results for “convergence analysis”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,423 records · Page 79

All sources NASA NTRS DOE OSTI

NASA NTRS1997

Multistage Schemes with Multigrid for Euler and Navier-Strokes Equations: Components and Analysis

A class of explicit multistage time-stepping schemes with centered spatial differencing and multigrids are considered for the compressible Euler and Navier-Stokes equations. These schemes are the basis for a family of computer programs (flow codes with multigrid (FLOMG) series) currently used to solve a wide range of fluid dynamics problems, including internal and external flows. In this paper, the components of these multistage time-stepping schemes are defined, discussed, and in many cases analyzed to provide additional insight into their behavior. Special emphasis is given to numerical dissipation, stability of Runge-Kutta schemes, and the convergence acceleration techniques of multigrid and implicit residual smoothing. Both the Baldwin and Lomax algebraic equilibrium model and the Johnson and King one-half equation nonequilibrium model are used to establish turbulence closure. Implementation of these models is described.

Swanson, R. C.↗

NASA NTRS1979

Viscid/inviscid interaction analysis of thrust augmenting ejectors

A method was developed for calculating the static performance of thrust augmenting ejectors by matching a viscous solution for the flow through the ejector to an inviscid solution for the flow outside the ejector. A two dimensional analysis utilizing a turbulence kinetic energy model is used to calculate the rate of entrainment by the jets. Vortex panel methods are then used with the requirement that the ejector shroud must be a streamline of the flow induced by the jets to determine the strength of circulation generated around the shroud. In effect, the ejector shroud is considered to be flying in the velocity field of the jets. The solution is converged by iterating between the rate of entrainment and the strength of the circulation. This approach offers the advantage of including external influences on the flow through the ejector. Comparisons with data are presented for an ejector having a single central nozzle and Coanda jet on the walls. The accuracy of the matched solution is found to be especially sensitive to the jet flap effect of the flow just downstream of the ejector exit.

Bevilacqua, P. M.↗

NASA NTRS1996

Sonic Boom Propagation Codes Validated by Flight Test

The sonic boom propagation codes reviewed in this study, SHOCKN and ZEPHYRUS, implement current theory on air absorption using different computational concepts. Review of the codes with a realistic atmosphere model confirm the agreement of propagation results reported by others for idealized propagation conditions. ZEPHYRUS offers greater flexibility in propagation conditions and is thus preferred for practical aircraft analysis. The ZEPHYRUS code was used to propagate sonic boom waveforms measured approximately 1000 feet away from an SR-71 aircraft flying at Mach 1.25 to 5000 feet away. These extrapolated signatures were compared to measurements at 5000 feet. Pressure values of the significant shocks (bow, canopy, inlet and tail) in the waveforms are consistent between extrapolation and measurement. Of particular interest is that four (independent) measurements taken under the aircraft centerline converge to the same extrapolated result despite differences in measurement conditions. Agreement between extrapolated and measured signature duration is prevented by measured duration of the 5000 foot signatures either much longer or shorter than would be expected. The duration anomalies may be due to signature probing not sufficiently parallel to the aircraft flight direction.

Poling, Hugh W.↗

NASA NTRS1996

Application of the kappa-omega Turbulence Model to Quasi-Three-Dimensional Turbomachinery Flows

Many computational fluid dynamics codes for turbomachinery use the Baldwin-Lomax (B-L) turbulence model. It is easy to implement in two dimensions and works well for predicting overall turbomachinery performance. However, it is awkward to implement in three dimensions, often has difficulty finding the length scale, has a crude transition model, and neglects freestream turbulence, surface roughness, and mass injection. The kappa-omega model developed by Wilcox is an appealing alternative for several reasons. First, it is the only two-equation model that can be integrated to the wall without requiring damping functions or the distance to the wall, and hence, should behave well numerically. Second, the effects of freestream turbulence, surface roughness, and mass injection are easily included in the model. Finally, transition can be simulated using the low Reynolds number version of the model. Menter applied the kappa-model to external flows and showed very good results for flows with adverse pressure gradients. Liu and Zheng described their implementation of the kappa-model in a cascade code that included an area change term to account for endwall convergence. They validated the model for a flat plate, and compared the B-L and kappa-models to measured surface pressures for a low-pressure turbine cascade. Since they did not use the low Reynolds number version of the model, their results showed problems resulting from early transition. In this Note the low Reynolds number kappa-model was incorporated in the author's quasi-three-dimensional turbomachinery analysis code. The code includes the effects of rotation, radius change, and stream-surface thickness variation, and also includes the B-L turbulence model. The kappa-omega model was implemented using many of Menter's recommendations and an implicit approximate-factorization scheme described by Baldwin and Barth. The model was tested for a transonic compressor with rotation and variable stream-surface radius and height, and for a transonic turbine vane with transition and heat transfer. Results were compared to the B-L model and to experimental data.

Chima, Rodrick V.↗

NASA NTRS2004

Rotor Airloads Prediction Using Loose Aerodynamic Structural Coupling

This work couples a computational fluid dynamics (CFD) code and rotorcraft computational structural dynamics (CSD) code to calculate helicopter rotor airloads across a range of flight conditions. An iterative loose (weak) coupling methodology is used to couple the CFD and CSD codes on a per revolution, periodic basis. The CFD uses a high fidelity, Navier-Stokes, overset grid methodology with first principles-based wake capturing. Modifications are made to the CFD code for aeroelastic analysis. For a UH-60A Blackhawk helicopter, four challenging level flight conditions are computed: 1) low speed (u = 0.15) with blade-vortex interaction, 2) high speed (u = 0.37) with advancing blade negative lift, 3) high thrust with dynamic stall (u = 0.24), and 4) hover. Results are compared with UH-60A Airloads Program fight test data. Most importantly, for all cases the loose coupling methodology is shown to be stable, convergent, and robust with full coupling of normal force, pitching moment, and chord force. In comparison with flight test data, normal force and pitching moment magnitudes are in good agreement. For the high speed and dynamic stall cases a phase lag in comparison with the data is seen, nonetheless, the shapes of the curves are very good. Overall, the results are noteworthy improvement over lifting line aerodynamics used in rotorcraft comprehensive codes.

Potsdam, Mark↗

NASA NTRS2011

Geodetic Strain Analysis Tool

A geodetic software analysis tool enables the user to analyze 2D crustal strain from geodetic ground motion, and create models of crustal deformation using a graphical interface. Users can use any geodetic measurements of ground motion and derive the 2D crustal strain interactively. This software also provides a forward-modeling tool that calculates a geodetic velocity and strain field for a given fault model, and lets the user compare the modeled strain field with the strain field obtained from the user s data. Users may change parameters on-the-fly and obtain a real-time recalculation of the resulting strain field. Four data products are computed: maximum shear, dilatation, shear angle, and principal components. The current view and data dependencies are processed first. The remaining data products and views are then computed in a round-robin fashion to anticipate view changes. When an analysis or display parameter is changed, the affected data products and views are invalidated and progressively re-displayed as available. This software is designed to facilitate the derivation of the strain fields from the GPS and strain meter data that sample it to facilitate the understanding of the strengths and weaknesses of the strain field derivation from continuous GPS (CGPS) and other geodetic data from a variety of tectonic settings, to converge on the "best practices" strain derivation strategy for the Solid Earth Science ESDR System (SESES) project given the CGPS station distribution in the western U.S., and to provide SESES users with a scientific and educational tool to explore the strain field on their own with user-defined parameters.

Kedar, Sharon↗

NASA NTRS2022

The Impacts of Freshwater Input and Surface Wind Velocity on the Strength and Extent of a Large High Latitude River Plume

Arctic Ocean physical and biogeochemical properties are strongly influenced by freshwater input from land and through the Bering Strait, where the mean currents transport water northward from the Bering Sea. The Yukon River is one of the largest rivers in North America and the Arctic, contributing large quantities of freshwater and terrigenous material to the coastal ocean in the northern Bering Sea. However, a detailed analysis of the coastal hydrodynamics at the outflow of the river has not been conducted in this remote but regionally important river. A three-dimensional hydrodynamic model was built to represent the lower Yukon River and coastal ocean for the ice-free months in 7 years. On average, a large anticyclonic eddy persisted at the main outflow of the Yukon that recirculates water back toward the coast where the currents converge to form a mean northward transport along the delta. Interannual spatial variance in salinity was relatively small, while there was substantial variance in u and v current velocity. u velocity spatial variance was correlated to the volume of freshwater discharge across years, while v velocity spatial variance was correlated to the N–S wind velocity. During strong wind events, plume structure was substantially altered: southerly winds deepened the plume and enhanced northward transport, while northerly winds shoaled and strengthened the pycnocline, and reversed the flow toward the south. The variability in plume dispersion on short time scales due to wind forcing has implications for where terrigenous material is processed in and settles out of the water column.

river plume↗

NASA NTRSDate not supplied

OVERFLOW Analysis of Supersonic Retropropulsion Testing on a Blunt Mars Entry Vehicle Concept

Supersonic retropropulsion (SRP) flow over a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) blunt-body vehicle was simulated using the Overflow Computational Fluid Dynamics (CFD) solver. Two nozzle configurations were tested (1E and 1F) as a subset of the seven total configurations in the Descent System Study (DSS) testing campaign. A generalized, Adaptive Mesh Refinement (AMR) shock capturing and plume refinement technique was developed and calibrated for producing automatic and unique grid systems optimized for any given testing condition. Solution independence from grid resolution was determined by successively increasing grid refinement until asymptotic convergence of the mean aerodynamic loads was observed. Overflow Reynolds-averaged Navier-Stokes (RANS) solutions produced realistic SRP flow phenomena, including Mach disk normal shocks contained in high-thrust plumes and bow shock triple-points. Solutions for the 1E nozzle configuration were steady. Conversely, a subset of the 1F test conditions were unsteady and exhibited periodic, non-sinusoidal expansions and contractions of the streamwise-oriented plume, resulting in unsteady aerodynamic loads on the vehicle. Initial investigations into the effect of turbulence modeling fidelity demonstrated significant differences between the aerodynamic loads simulated with RANS and Detached Eddy Simulation (DES). Multiple flow mechanisms were identified as root causes of these differences, including bow shock shape augmentation and reduced entrainment due to reduced turbulence in the DES simulations. Initial comparisons of pitching moment control authority and aerodynamic drag performance between the tested nozzle configurations demonstrated that the 1E configuration may be more advantageous for the Mars Entry, Descent, and Landing (EDL) task.

Supersonic Retropropulsion↗

NASA NTRSDate not supplied

OVERFLOW Analysis of Supersonic Retropropulsion Testing on a Blunt Mars Entry Vehicle Concept

Supersonic Retropropulsion↗

NASA NTRS1973

An investigation of a two-dimensional propulsive lifting system

Several aspects of the nonhomogeneous flow associated with a system combining lifting and propulsive requirements of an aircraft are considered by analytical and experimental methods. The basic geometry of the problem is that of two lifting surfaces with an actuator disk located between them. The principles governing flow with energy addition are examined. Basic equations and boundary conditions are developed for the complete inviscid and incompressible analysis for the two-dimensional case. The corresponding flow singularities are discussed and the integral equations which completely specify the system are derived. The two special cases of small and large energy addition are considered in detail including solutions. A numerical procedure is developed to solve the full problem including allowance for the wake deflection. Appropriate vorticity forms are used to represent the entire system. An iterative scheme is presented which rapidly converges to a solution for the magnitude and location of the system vorticity distributions. Forces and moments are evaluated on the propulsive lift system.

Shollenberger, C. A.↗

NASA NTRS1980

Stratospheric pollutant transports by planetary waves and tropical circulations

Radiosonde and U-2 aircraft data from the Intertropical Convergence Zone experiment of July, 1977 were analyzed. The radiosonde data above 15 km was plotted in the form of time-height sections so that wave motions in the meteorological variables (temperature, zonal wind, and meridional wind) could be easily delineated. The results showed the presence of a planetary 16 day Kelvin wave. More significantly, large amplitude gravity waves with vertical wavelengths on the order of 3-6 km were found. Spectral analysis of frequencies showed peaks in meridional wind variance at 2-3 days and 5 days, with the shorter periods stronger at higher altitudes. In the time domain, two instances of organized wave groups propagating upward into the stratosphere were isolated. Associated with one of these wave groups was a region of very large vertical shear in which either turbulence or small scale (10 km wavelength or less) gravity wave activity had been found by the U-2 aircraft.

Riegel, C. A.↗

NASA NTRS1988

Formation of the lamellar structure in Group IA and IIID iron meteorites

Analytical EM, light microscopy, and electron microprobe analysis are used to study the lamellar plessite structure of Group IA and IIID iron meteorites. The alpha lamellae in IIID structures contained a compositional gradient from 6.1 + or - 0.7 wt pct Ni at the center of the alpha lamellae to 3.6 + or - 0.5 wt pct at the alpha/gamma interface. For the Group IA irons, compositions of 4 wt pct Ni in alpha and about 48 wt pct Ni in gamma are found. Convergent beam electron diffraction was used to characterize the orientation relations at the alpha/gamma interface in the lamellar regions of both Group IA and IIID. The phase transformations responsible for the observed lamellar structure in the IA and IIID chemical groups were also investigated.

Kowalik, J. A.↗

NASA NTRS1993

A comparison of locally adaptive multigrid methods: LDC, FAC and FIC

This study is devoted to a comparative analysis of three 'Adaptive ZOOM' (ZOom Overlapping Multi-level) methods based on similar concepts of hierarchical multigrid local refinement: LDC (Local Defect Correction), FAC (Fast Adaptive Composite), and FIC (Flux Interface Correction)--which we proposed recently. These methods are tested on two examples of a bidimensional elliptic problem. We compare, for V-cycle procedures, the asymptotic evolution of the global error evaluated by discrete norms, the corresponding local errors, and the convergence rates of these algorithms.

Khadra, Khodor↗

NASA NTRS2012

SE-FIT

The mathematical theory of capillary surfaces has developed steadily over the centuries, but it was not until the last few decades that new technologies have put a more urgent demand on a substantially more qualitative and quantitative understanding of phenomena relating to capillarity in general. So far, the new theory development successfully predicts the behavior of capillary surfaces for special cases. However, an efficient quantitative mathematical prediction of capillary phenomena related to the shape and stability of geometrically complex equilibrium capillary surfaces remains a significant challenge. As one of many numerical tools, the open-source Surface Evolver (SE) algorithm has played an important role over the last two decades. The current effort was undertaken to provide a front-end to enhance the accessibility of SE for the purposes of design and analysis. Like SE, the new code is open-source and will remain under development for the foreseeable future. The ultimate goal of the current Surface Evolver Fluid Interface Tool (SEFIT) development is to build a fully integrated front-end with a set of graphical user interface (GUI) elements. Such a front-end enables the access to functionalities that are developed along with the GUIs to deal with pre-processing, convergence computation operation, and post-processing. In other words, SE-FIT is not just a GUI front-end, but an integrated environment that can perform sophisticated computational tasks, e.g. importing industry standard file formats and employing parameter sweep functions, which are both lacking in SE, and require minimal interaction by the user. These functions are created using a mixture of Visual Basic and the SE script language. These form the foundation for a high-performance front-end that substantially simplifies use without sacrificing the proven capabilities of SE. The real power of SE-FIT lies in its automated pre-processing, pre-defined geometries, convergence computation operation, computational diagnostic tools, and crash-handling capabilities to sustain extensive computations. SE-FIT performance is enabled by its so-called file-layer mechanism. During the early stages of SE-FIT development, it became necessary to modify the original SE code to enable capabilities required for an enhanced and synchronized communication. To this end, a file-layer was created that serves as a command buffer to ensure a continuous and sequential execution of commands sent from the front-end to SE. It also establishes a proper means for handling crashes. The file layer logs input commands and SE output; it also supports user interruption requests, back and forward operation (i.e. undo and redo), and others. It especially enables the batch mode computation of a series of equilibrium surfaces and the searching of critical parameter values in studying the stability of capillary surfaces. In this way, the modified SE significantly extends the capabilities of the original SE.

Chen, Yongkang↗

DOE OSTI2023

Optimization Algorithm for Criticality Experiment Design Using Whisper

Many criticality experiments performed to aid in nuclear data evaluation are designed to maximize the sensitivity of the system’s effective neutron multiplication factor to a certain nuclide reaction pair over an energy region of interest. This is typically done by evaluating possible designs in a transport code such as MCNP and selecting the one with the highest desired sensitivity. A designer has many tools to try to maximize this sensitivity such as different moderators, reflectors, fuels, and geometries. This balancing act of identifying a critical and maximally sensitive system become very computationally expensive as more variables are added and higher precisions are desired. In order to identify these optimal configurations more efficiently a Particle Swarm Optimization (PSO) algorithm coupled with MCNP has been developed by Los Alamos National Laboratory (LANL). This algorithm has been used to design two upcoming criticality experiments that will be performed at the National Criticality Experiments Research Center (NCERC), located at the Nevada National Security Site, and operated by LANL, the only general-purpose critical experiments laboratory in the United States. PSO uses a population (swarm) of candidate solutions (particles) on a search space of dimensions such as moderator and reflector thicknesses or enrichments and concentrations. These particles move around the search space from generation to generation according to simple rules. Eventually, the swarm converges on the configuration that is both critical and maximally sensitive to a piece of nuclear data. PSO is well suited for criticality experiments as the algorithm is agnostic to the underlying physics, meaning it is effective on many different experimental setups. This algorithm has been modified to maximize the nuclear data similarity coefficient between an application case and an experiment aimed at replicating the application case using WHISPER, a nuclear criticality safety analysis tool. This allows for the efficient design of critical experiments informed by nuclear data sensitives.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1974

Finite-element analysis of vibration and flutter of cantilever anisotropic plates

The hybrid stress finite element method is used to study the effects of filament angle and orthotropicity parameter on the vibration and flutter characteristics of cantilevered anisotropic plates. The results indicate a generally strong lack of monotonic dependence on filament angle. Also, the natural frequency in certain cases involving the first few modes, can become relatively insensitive to both filament angle and orthotropicity parameter for a range of filament angle beyond 70 deg. Values of critical dynamic pressure are obtained by a modal approach, in which the mode shapes are obtained by the hybrid stress method. Convergence of the modal method is rather rapid for the configurations analyzed, and a comparison of the method with an exact solution for the case of an isotropic simply supported plate shows that use of six modes gives excellent agreement.

Rossettos, J. N.↗

DOE OSTI2021

Pricing and Energy Trading in Peer-to-Peer Zero Marginal-Cost Microgrids

Efforts to utilize 100% renewable energy in community microgrids require new approaches to energy markets and transactions to efficiently address periods of scarce energy supply. In this paper we contribute to the promising approach of peer-to-peer (P2P) energy trading in two main ways: analysis of a centralized, welfare-maximizing economic dispatch that characterizes optimal price and allocations, and a novel P2P system for negotiating energy trades that yields physically feasible and at least weakly Pareto-optimal outcomes. Our main results are 1) that optimal pricing is insufficient to induce agents with batteries to take optimal actions, 2) a novel P2P algorithm to addresses this while keeping private information, 3) a formal proof that this algorithm converges to the centralized solution in the case of two agents negotiating for a single period, and 4) numerical simulations of the P2P algorithm performance with up to 10 agents and 24 periods that show it converges on average to total welfare within 0.1% of the social optimum in on the order of 10s to 100s of iterations, increasing with the number of agents, time periods, and total storage capacity.

batteries↗

DOE OSTI2021

Code and Solution Verification Assessment of the CTF Thermal Hydraulic Subchannel Code

CTF is a thermal-hydraulics subchannel code jointly developed by Oak Ridge National Laboratory and North Carolina State University. Over the past seven years, the Consortium for Advanced Simulation of Light Water Reactors (CASL) has made a significant investment in developing CTF so it can be used to model light water reactors, including nominal operating conditions, departure from nucleate boiling analysis, and transients ranging from loss of flow to reactivity insertion accidents. In addition to implementing new modeling capabilities and developing the user input and output interface, extensive work has been performed to improve the code’s quality assurance program, resulting in a development process that conforms with NQA-1 requirements. The CASL program follows the Predictive Capability Maturity Model (PCMM) approach for assessing code quality, which emphasizes performing code verification(ensuring the code converges to the correct answer) and solution verification (ensuring the code converges for the intended application). Code and solution verification are used to identify uncertainty errors introduced by numerical approximations in the code and are important for demonstrating that the model is coded without error, which is an important aspect of the Best Estimate plus Uncertainty method. This paper presents a comprehensive overview of the code and solution verification testing that has been performed on CTF. A top-down approach is taken in which the intended CTF applications are presented, followed by the code features required for their modeling. These features are then linked to the applicable code and solution verification tests that demonstrate proper functioning. Past testing efforts are summarized, and new tests are added to help close gaps in the presented test matrix. Rather than performing “one-off” exercises, these tests are added to the automated CTF regression test suite to ensure continual code quality.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗