Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical benchmark”

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,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Optimal Financing Decision in a Contract Food Supply Chain with Capital Constraint

To solve the financing problem of the food producers, we consider a two‐echelon contract food supply chain composed of a family farm with capital constraints and a food processing enterprise. With no capital constraints as the benchmark model, we analyze optimal decisions of the family farm and the food processing enterprise in the case of bank financing with bank participation only and bank financing with “government, bank, and insurance” coparticipation. Then, we discuss how the risk of yield uncertainty influences the optimal decisions and profits of the family farm and the food processing enterprise under different financing situation. Meanwhile, the reason why the government subsidizes agriculture is explored, and the policy of minimum purchase price of the food is initiated when the market price is too low. Finally, the numerical examples and sensitivity analysis are presented. The results show that the bank financing with “government, bank, and insurance” coparticipation improves the welfare of supply chain members more obviously than the bank financing with bank participation only; when the rice price is too low, the policy of minimum purchase price of food is initiated, which increases the revenue and the growing enthusiasm of the family farm; the profits of the family farm and the food processing enterprise will decrease as the risk of yield uncertainty increases in the case of bank financing, and the risk of yield uncertainty will be reduced for the family farm when bank financing with “government, bank, and insurance” coparticipation.

Luo, Ying (ORCID:0000000158101443)↗
DOE OSTI2023

SCALE 6.2 Lattice Physics Performance Assessment

The US Nuclear Regulatory Commission relies on the lattice physics analysis capabilities of the SCALE code system to perform confirmatory licensing analyses. Either SCALE lattice physics code—TRITON/NEWT or Polaris—can be used to generate cross section data used by the PARCS nodal core simulator for full-core neutronics calculations. This report presents an assessment of the accuracy of SCALE lattice physics codes for preparation of lattice physics data that are used to support simulator codes such as the NRC’s PARCS, for UO 2 -mixed oxide (MOX)/Zr fueled light water reactor (LWR) analyses. Due to the nature of lattice physics calculations, critical reactor experiment benchmarks cannot be modeled in explicit detail in a lattice physics code. However, this limitation does not mean that these measurement data are not usable for lattice physics studies. Therefore, either geometry approximations or axial buckling must be implemented to determine the critical water height. These modeling limitations have led to development of a three-phase assessment strategy. In the first phase, selected critical experiment benchmarks are modeled using the SCALE 3D continuous-energy (CE) Monte Carlo (MC) code KENO, which is the most rigorous neutron transport method available in SCALE, with no approximations in the spatial, angle, or energy treatments. Biases and statistical uncertainties in quantities of interest such as k eff and pin power distributions are determined by comparing CE KENO results to experimental data. This first phase in the accuracy assessment is to demonstrate and establish the use of CE KENO as a reference solution for the second phase. In the second phase, 14 numerical test suites are used to compare SCALE lattice physics calculations with CE KENO as a reference solution. In the third phase, SCALE lattice physics depletion calculations are performed, and the spent fuel isotopic results are compared with available radiochemical assay measurements. Isotopic measurement comparisons provide quantitative assessment of isotopic density distribution predictions with the depletion models in SCALE lattice physics codes. This report documents results for all test suites. The assessment was performed using standard production techniques unless otherwise noted. Both TRITON/NEWT and Polaris exhibited acceptable accuracy for most test cases. For the few test cases in which acceptable accuracy criteria were not met, further code and data development are planned. The computer codes used in this assessment are as follows: SCALE 6.2 rev19189 (pre-release of 6.2.1) was used for CE KENO, TRITON/NEWT and Polaris (PWR only) calculations; SCALE 6.2.2 was used for Polaris calculations to address several updates in support of support boiling water reactor (BWR) lattice geometries. For this work, there is no difference between SCALE 6.2.1 and SCALE 6.2.2 beyond the inclusion of the Polaris BWR analysis capability; The ENDF/B-VII.1 continuous energy and 252 group libraries deployed in SCALE 6.2 were generated with AMPX 6.2 that is distributed with SCALE.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1988

Primary instabilities and bicriticality in flow between counter-rotating cylinders

The primary instabilities and bicritical curves for flow between counter-rotating cylinders have been computed numerically from the Navier-Stokes equations assuming axial periodicity. The computations provide values of the Reynolds numbers, wavenumbers, and wave speeds at the primary transition from Couette flow for radius ratios from 0.40-0.98. Particular attention has been focused on the bicritical curves that separate (as the magnitude of counter-rotation is increased) the transitions from Couette flow to flows with different azimuthal wavenumbers m and m + 1. This lays the foundation for further analysis of nonlinear mode interactions and pattern formation occurring along the bicritical curves and serves as a benchmark for experimental studies. Preliminary experimental measurements of transition Reynolds numbers and wave speeds presented here agree well with the computations from the mathematical model.

Langford, W. F.↗
NASA NTRS2013

Regional Climate and Variability of the Summertime Continental United States in Reanalyses

Understanding climate variability at regional scales is an important for research and societal needs. Atmospheric retrospective-analyses (or reanalyses) integrate multitudes of observing systems with numerical models to produce continuous data that include variables not easily observed, if at all. The breadth of variables as well as observational influence included in reanalyses make them ideal for investigating climate variability. In this paper, we assess NASA s Modern Era Retrospective-analysis for Research and Application (MERRA) regional variability in North America, specifically the United States, in conjunction with current satellite data reanalyses. Emphasis is placed on summertime precipitation because 1) it is a difficult parameter to capture in the most difficult season, 2) significant observational resources exist to benchmark comparisons, and 3) accurate assessment of precipitation variability is crucial to a multitude of sectors and applications. Likewise, we have also begun to evaluate surface air temperature. While precipitation biases are identified, year to year variability of the precipitation variations, in many cases, are quite reasonable. However, some spurious long term trends and sudden shifts in the time series are also identified. In surface air temperature, analysis of station observations provides ERA Interim a clear overall advantage. However, in a number of regions, all the reanalyses are quite comparable in variability and trend. In other regions, significant precipitation biases may occur, which has implications for the ancillary process data in a reanalysis, such as surface fluxes. We also characterize the reanalyses ability to capture variability related to ENSO. In general, the summertime variations of precipitation in the reanalyses are more highly correlated (positively) to ENSO (using ENSO34) than are the observations. The Northwestern US shows the largest positive correlations to ENSO34, and reanalyses agree with that, and correlate highly with the gauge observations there. Additional evaluation of the data assimilation using the MERRA Gridded Innovations and Observations (GIO) data, which consists of the assimilated observations as well as forecast an analysis error fields. This work represents an initial evaluation of using the MERRA reanalysis to study regional climate variations of the United States, and identifying its applied limitations

Bosilovich, M G.↗
DOE OSTI2024

A Self Consistent 2D Simulation of Coherent Synchrotron Radiation Effects on Beam Dynamics

An increasing interest in high quality and high current electron beams necessitates a thorough understanding and prediction of coherent synchrotron radiation effects. The self-interaction of charged particles in a beam undergoing synchrotron motion is a physically significant process that is all too often computationally intensive with very little analytical results to rely on for the general case. The coherent spectrum of this interaction is of utmost importance to the design of free electron lasers (FELs) and an accurate assessment is imperative for their design. This work presents a novel implementation to the numerical simulation of charged particle beams. The simulation is a self-consistent approach including the self-fields generated by the beam of which coherent synchrotron radiation effects are of primary interest. A particle-in-cell model is used where a planar beam sampled by point particles is deposited on an encompassing grid at each timestep. The electromagnetic fields are calculated on the grid using the retarded potentials according to causality. The electromagnetic forces from the fields are interpolated on each particle which in turn advance in time. The simulation is benchmarked against well-established results for coherent synchrotron radiation effects. In addition, studies are provided that show the convergence of simulation results for increasing resolution. A study into the transverse beam size effects on beam dynamics is performed as well as a proof of concept where the simulation is used by a genetic algorithm to optimize the design parameters of a beam lattice. The results of these studies in tandem verify the efficacy of the simulation for its practical use in accelerator design or the study of synchrotron radiation effects

Duffin, Dallan [Old Dominion Univ., Norfolk, VA (U↗
DOE OSTI2025

LLM Benchmarking with LLaMA2: Evaluating Code Development Performance Across Multiple Programming Languages

The rapid evolution of large language models (LLMs) has opened new possibilities for automating various tasks in software development. This paper evaluates the capabilities of the LLaMA 2-70B model in automating these tasks for scientific applications written in commonly used programming languages. Using representative test problems, we assess the model's capacity to generate code, documentation, and unit tests, as well as its ability to translate existing code between commonly used programming languages. Our comprehensive analysis evaluates the compilation, runtime behavior, and correctness of the generated and translated code. Additionally, we assess the quality of automatically generated code, documentation, and unit tests. Here, our results indicate that while LLaMA 2-70B frequently generates syntactically correct and functional code for simpler numerical tasks, it encounters substantial difficulties with more complex, parallelized, or distributed computations, requiring considerable manual corrections. We identify key limitations and suggest areas for future improvements to better leverage AI-driven automation in scientific computing workflows.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2011

Experimental and Numerical Investigation of Reduced Gravity Fluid Slosh Dynamics for the Characterization of Cryogenic Launch and Space Vehicle Propellants

As space programs increasingly investigate various options for long duration space missions the accurate prediction of propellant behavior over long periods of time in microgravity environment has become increasingly imperative. This has driven the development of a detailed, physics-based understanding of slosh behavior of cryogenic propellants over a range of conditions and environments that are relevant for rocket and space storage applications. Recent advancements in computational fluid dynamics (CFD) models and hardware capabilities have enabled the modeling of complex fluid behavior in microgravity environment. Historically, launch vehicles with moderate duration upper stage coast periods have contained very limited instrumentation to quantify propellant stratification and boil-off in these environments, thus the ability to benchmark these complex computational models is of great consequence. To benchmark enhanced CFD models, recent work focuses on establishing an extensive experimental database of liquid slosh under a wide range of relevant conditions. In addition, a mass gauging system specifically designed to provide high fidelity measurements for both liquid stratification and liquid/ullage position in a micro-gravity environment has been developed. This pUblication will summarize the various experimental programs established to produce this comprehensive database and unique flight measurement techniques.

Walls, Laurie K.↗
DOE OSTI2024

Efficient sensitivity analysis of the thermal profile in powder bed fusion of metals using hypercomplex automatic differentiation finite element method

Rapid cyclic temperature fluctuation occurring in powder bed fusion of metals using a laser beam (PBF-LB/M) influences the formation of flaws in printed parts. Consequently, there is a pressing need to enhance the quality of printed parts by developing innovative methodologies that can predict thermal histories and help uncover the intricate relationships between process parameters and thermal profiles. Sensitivity Analysis (SA) emerges as an essential tool for this, offering the potential for process optimization and enhanced quality control. Nonetheless, conventional SA methodologies often incur in excessive computational costs and potential numerical approximation errors. Here, to address this technical challenge, we present a novel method for SA that integrates the HYPercomplex-based Automatic Differentiation (HYPAD) technique with transient thermal simulations conducted via the finite element method (FEM). Leveraging this methodology, we efficiently and accurately perform SA for PBF-LB/M processes in a post-processing step. Compared to traditional methods like Finite Differences (FD), HYPAD-FEM required 96 % less computational time for obtaining sensitivities for 22 process parameters, under a comparative study conducted within the context of the 2018–02 AM benchmark of the National Institute of Standards and Technology. In summary, HYPAD-FEM offers superior efficiency and accuracy in SA over conventional methods, delivering the best sensitivity of a model without the need for step-size selection and problem or parameter-based implementations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling betatron radiation using particle-in-cell codes for plasma wakefield accelerator diagnostics

The analysis of plasma wakefield acceleration experimental measurements, particularly in the characterization of photons emitted through the betatron radiation mechanism, requires the development of accurate numerical models. These computational models are crucial for supporting modern instrumentation designed to measure the single-shot, double-differential angular-energy radiation spectra resulting from interactions between beams and plasmas. Motivated by the needs of such applications, this paper presents detailed numerical models of betatron radiation generated in beam-plasma acceleration experiments. These models are based on the integration of the Liénard-Wiechert (LW) potentials, applied to computed particle trajectories. The particle trajectories are generated using three distinct methods: first, by tracking particles through idealized fields in the blowout regime of PWFA; second, by obtaining trajectories using the fast quasistatic particle-in-cell (PIC) code quickpic; and third, obtaining trajectories from the fully self-consistent PIC code osiris. To ensure the accuracy and reliability of these models, the paper includes various benchmark tests using analytical expressions, as well as employing the PIC code epoch, which takes an alternative approach by using a Monte Carlo quantum electrodynamics (QED)-based radiation model. Additionally, the paper presents simulations of the expected experimental betatron radiation spectra, taking into account parameters relevant to PWFA and plasma photocathode experiments at the SLAC FACET-II facility.

Yadav, M. [University of California, Los Angeles, ↗
NASA NTRS1991

A computer code for multiphase all-speed transient flows in complex geometries. MAST version 1.0

The operation of the MAST code, which computes transient solutions to the multiphase flow equations applicable to all-speed flows, is described. Two-phase flows are formulated based on the Eulerian-Lagrange scheme in which the continuous phase is described by the Navier-Stokes equation (or Reynolds equations for turbulent flows). Dispersed phase is formulated by a Lagrangian tracking scheme. The numerical solution algorithms utilized for fluid flows is a newly developed pressure-implicit algorithm based on the operator-splitting technique in generalized nonorthogonal coordinates. This operator split allows separate operation on each of the variable fields to handle pressure-velocity coupling. The obtained pressure correction equation has the hyperbolic nature and is effective for Mach numbers ranging from the incompressible limit to supersonic flow regimes. The present code adopts a nonstaggered grid arrangement; thus, the velocity components and other dependent variables are collocated at the same grid. A sequence of benchmark-quality problems, including incompressible, subsonic, transonic, supersonic, gas-droplet two-phase flows, as well as spray-combustion problems, were performed to demonstrate the robustness and accuracy of the present code.

Chen, C. P.↗
DOE OSTI2011

Oscillator strengths and transition probabilities for allowed and forbidden transitions in Fe XIX

An extensive set of oscillator strengths, line strengths, and radiative decay rates for the allowed and forbidden transitions in Fe XIX is presented. They correspond to 1626 fine structure levels of total angular momenta 0{<=}J{<=}8 of even and odd parities with 2{<=}n{<=}10, 0{<=}l{<=}9, 0{<=}L{<=}10, and (2S+1)=1, 3, 5. In contrast, the compiled table of the National Institute for Standards and Technology (NIST) lists only 63 observed levels. A total of 289,291 electric dipole allowed transitions are presented. They were obtained in the close coupling approximation using the relativistic Breit-Pauli R-matrix method. The wavefunction expansion included 15 levels of the configurations 2s{sup 2}2p{sup 3}, 2s2p{sup 4}, and 2p{sup 5} of the Fe XX core. The calculated fine structure levels are assigned with spectroscopic identifications using quantum defect analysis. Comparison with the observed energies shows very good agreement, the largest difference being less than 4%. The transitions also compare well with the compiled data by NIST and recent calculations. The forbidden transitions of the electric quadrupole and octupole, and magnetic dipole and quadrupole, type are presented for the 379 levels of the configurations 2s{sup 2}2p{sup 4}, 2s2p{sup 5}, 2p{sup 6}, 2s{sup 2}2p{sup 3}3s, 2s{sup 2}2p{sup 3}3p, 2s{sup 2}2p{sup 3}3d, 2s{sup 2}2p{sup 3}4s, 2s{sup 2}2p{sup 3}4p, 2s{sup 2}2p{sup 3}4d, 2s{sup 2}2p{sup 3}4f, 2s2p{sup 4}3s, 2s2p{sup 4}3p, 2s2p{sup 4}3d, 2s2p{sup 4}4s, 2s2p{sup 4}4p, and 2s{sup 2}2p{sup 2}3s{sup 2} of Fe XIX. They correspond to a total of 66,619 transitions. These results have been obtained from relativistic Breit-Pauli atomic structure calculations using the program SUPERSTRUCTURE. The forbidden transition probabilities show very good agreement with those compiled by NIST. - Highlights: {yields} Presents the most complete (n up to 10) set of transitions for Fe XIX. {yields} Considers both allowed and forbidden transitions. {yields} Large number of new energy levels are found. {yields} Benchmarks existing accurate energies and transitions. {yields} Completeness and accuracy are crucial to solve its abundance in the Sun.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

In-plane crashworthiness of bio-inspired hierarchical honeycombs

Biological tissues like bone, wood, and sponge possess hierarchical cellular topologies, which are lightweight and feature an excellent energy absorption capability. Here we present a system of bio-inspired hierarchical honeycomb structures based on hexagonal, Kagome, and triangular tessellations. The hierarchical designs and a reference regular honeycomb configuration are subjected to simulated in-plane impact using the nonlinear finite element code LS-DYNA. The numerical simulation results show that the triangular hierarchical honeycomb provides the best performance compared to the other two hierarchical honeycombs, and features more than twice the energy absorbed by the regular honeycomb under similar loading conditions. We also propose a parametric study correlating the microstructure parameters (hierarchical length ratio r and the number of sub cells N) to the energy absorption capacity of these hierarchical honeycombs. The triangular hierarchical honeycomb with N = 2 and r = 1/8 shows the highest energy absorption capacity among all the investigated cases, and this configuration could be employed as a benchmark for the design of future safety protective systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Revenue-Maximizing Shared Parking and Electric Vehicle Charging Management in Multi-Unit Dwellings

In urban areas, searching for parking and electric vehicle (EV) charging can result in cruising, congestion, and environmental externalities. Recognizing the business opportunity of offering private parking and charging infrastructure access within multi-unit dwellings (MUDs) during daytime, we model a shared parking and EV charging management system. We maximize the revenue of MUD charging hubs in mixed land use, catering to public demand. Our approach accounts for the objectives of the two stakeholders involved: a demand model is fitted on the choices of EV charging users, and the supply model optimizes the allocation of parking and charging requests in an MUD parking lot. A binary integer linear programming model for the allocation of parking and charging spaces with a rolling horizon is integrated with matching rules that handle both parking and charging requests. In our numerical experiments in a neighborhood of Chicago, Illinois, we estimate the performance of the MUD parking and charging system with metrics that include revenue, number of matchings, and utilization rates. At any given time, MUDs with lower prices attract more charging requests, particularly those of longer duration, resulting in higher revenue and greater charging utilization. Dynamic pricing facilitates a more equitable distribution of requests; as MUD parking lots reach capacity and their fees increase, other MUDs become more competitive, attracting additional requests. Comparing our method against first-come-first-served and optimal-solution benchmarks, we demonstrate our model’s effectiveness in dynamically managing mixed parking and charging demand in MUD charging hubs.

electric vehicle, multi-unit dwelling, charging in↗
NASA NTRS1989

The effects of Reynolds number, rotor incidence angle, and surface roughness on the heat transfer distribution in a large-scale turbine rotor passage

A combined experimental and computational program was conducted to examine the heat transfer distribution in a turbine rotor passage geometrically similiar to the Space Shuttle Main Engine (SSME) High Pressure Fuel Turbopump (HPFTP). Heat transfer was measured and computed for both the full-span suction and pressure surfaces of the rotor airfoil as well as for the hub endwall surface. The primary objective of the program was to provide a benchmark-quality data base for the assessment of rotor passage heat transfer computational procedures. The experimental portion of the study was conducted in a large-scale, ambient temperature, rotating turbine model. Heat transfer data were obtained using thermocouple and liquid-crystal techniques to measure temperature distributions on the thin, electrically-heated skin of the rotor passage model. Test data were obtained for various combinations of Reynolds number, rotor incidence angle and model surface roughness. The data are reported in the form of contour maps of Stanton number. These heat distribution maps revealed numerous local effects produced by the three-dimensional flows within the rotor passage. Of particular importance were regions of local enhancement produced on the airfoil suction surface by the main-passage and tip-leakage vortices and on the hub endwall by the leading-edge horseshoe vortex system. The computational portion consisted of the application of a well-posed parabolized Navier-Stokes analysis to the calculation of the three-dimensional viscous flow through ducts simulating the a gas turbine passage. These cases include a 90 deg turning duct, a gas turbine cascade simulating a stator passage, and a gas turbine rotor passage including Coriolis forces. The calculated results were evaluated using experimental data of the three-dimensional velocity fields, wall static pressures, and wall heat transfer on the suction surface of the turbine airfoil and on the end wall. Particular attention was paid to an accurate modeling of the passage vortex and to the development of the wall boundary layers including crossflow.

Blair, Michael F.↗
DOE OSTI2015

Effects of interactions on the generalized Hong–Ou–Mandel effect

We numerically investigate the influence of interactions on the generalized Hong–Ou–Mandel (HOM) effect for bosonic particles in a (quasi-)one-dimensional set-up with weak harmonic confinement and show results for the cases of N = 2, N = 3 and N = 4 bosons interacting with a beam splitter, whose role is played by a δ-barrier. In particular, we focus on the effect of attractive interactions and compare the results with the repulsive case, as well as with the analytically available results for the non-interacting case (that we use as a benchmark). In addition, we observe a fermionization effect both for growing repulsive and attractive interactions, i.e., the dip in the HOM coincidence count is progressively smeared out, for increasing interaction strengths. The role of input asymmetries is also explored.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Fast and robust all-electron density functional theory calculations in solids using orthogonalized enriched finite elements

Here, we present a computationally efficient approach to perform systematically convergent real-space all-electron Kohn-Sham density functional theory calculations for solids using an enriched finite element (FE) basis. The enriched FE basis is constructed by augmenting the classical FE basis with atom-centered numerical basis functions, comprising of atomic solutions to the Kohn-Sham problem. Notably, to improve the conditioning, we orthogonalize the enrichment functions with respect to the classical FE basis, without sacrificing the locality of the resultant basis. In addition to improved conditioning, this orthogonalization procedure also renders the overlap matrix block diagonal, greatly simplifying its inversion. Subsequently, we use a Chebyshev polynomial based filtering technique to efficiently compute the occupied eigenspace in each self-consistent field iteration. We demonstrate the accuracy and efficiency of the proposed approach on periodic unit cells and supercells. The benchmark studies show a staggering 130× speedup of the orthogonalized enriched FE basis over the classical FE basis. We also present a comparison of the orthogonalized enriched FE basis with the linearized augmented plane-wave + local orbitals basis, both in terms of accuracy and efficiency. Notably, we demonstrate that the orthogonalized enriched FE basis can handle large system sizes of ~10 000 electrons. Finally, we observe good parallel scalability of our implementation with 92% efficiency at 22× speedup for a system with 620 electrons

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

An $\bar{F}$ Meshfree Treatment for Nearly Incompressible Materials

The computational modeling of nearly incompressible materials is a difficult task for many numerical methods, and even after several decades of investigation, it is still an active research area. This report seeks to address the treatment of incompressible materials in meshfree methods using a synergistic combination of two treatments. The first treatment is an $\bar{F}$ method, where the decomposed dilatational and deviatoric parts are calculated over different smoothing domains. The second treatment “activates” additional nodes throughout the domain to increase the flexibility of the model. We implement this synergistic combination in the context of the reproducing kernel particle method (RKPM) and present results for the Cook’s membrane benchmark problem. The results are compared with those using the composite tet10 finite element with a volume-averaged J formulation. We show that the combined treatment is an effective way to deal with nearly incompressible materials in a meshfree framework and compares well with other highly-effective treatments.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2001

Gust Acoustic Response of a Swept Rectilinear Cascade Using The Space-Time CE/SE Method

The benchmark problem 3 in Category 3 of the third Computational Aero-Acoustics (CAA) Workshop sponsored by NASA Glenn Research Center is solved using the space-time conservation element and solution element (CE/SE) method. This problem concerns the unsteady response of a rectilinear swept cascade to an incident gust. The acoustic field generated by the interaction of the gust with swept at plates in the cascade is computed by solving the 3D nonlinear Euler equations using the space-time CE/SE method. A parallel version of the 3D CE/SE Euler solver is employed to obtain numerical solutions for several sweep angles. Numerical solutions are presented and compared with the analytical solutions.

Wang, X. Y.↗