SEARCH · Engineering Papers

Results for “numerical 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,657 records · Page 92

All sources NASA NTRS DOE OSTI

NASA NTRS1998

Structural Similitude and Scaling Laws

Aircraft and spacecraft comprise the class of aerospace structures that require efficiency and wisdom in design, sophistication and accuracy in analysis and numerous and careful experimental evaluations of components and prototype, in order to achieve the necessary system reliability, performance and safety. Preliminary and/or concept design entails the assemblage of system mission requirements, system expected performance and identification of components and their connections as well as of manufacturing and system assembly techniques. This is accomplished through experience based on previous similar designs, and through the possible use of models to simulate the entire system characteristics. Detail design is heavily dependent on information and concepts derived from the previous steps. This information identifies critical design areas which need sophisticated analyses, and design and redesign procedures to achieve the expected component performance. This step may require several independent analysis models, which, in many instances, require component testing. The last step in the design process, before going to production, is the verification of the design. This step necessitates the production of large components and prototypes in order to test component and system analytical predictions and verify strength and performance requirements under the worst loading conditions that the system is expected to encounter in service. Clearly then, full-scale testing is in many cases necessary and always very expensive. In the aircraft industry, in addition to full-scale tests, certification and safety necessitate large component static and dynamic testing. Such tests are extremely difficult, time consuming and definitely absolutely necessary. Clearly, one should not expect that prototype testing will be totally eliminated in the aircraft industry. It is hoped, though, that we can reduce full-scale testing to a minimum. Full-scale large component testing is necessary in other industries as well, Ship building, automobile and railway car construction all rely heavily on testing. Regardless of the application, a scaled-down (by a large factor) model (scale model) which closely represents the structural behavior of the full-scale system (prototype) can prove to be an extremely beneficial tool. This possible development must be based on the existence of certain structural parameters that control the behavior of a structural system when acted upon by static and/or dynamic loads. If such structural parameters exist, a scaled-down replica can be built, which will duplicate the response of the full-scale system. The two systems are then said to be structurally similar. The term, then, that best describes this similarity is structural similitude. Similarity of systems requires that the relevant system parameters be identical and these systems be governed by a unique set of characteristic equations. Thus, if a relation or equation of variables is written for a system, it is valid for all systems which are similar to it. Each variable in a model is proportional to the corresponding variable of the prototype. This ratio, which plays an essential role in predicting the relationship between the model and its prototype, is called the scale factor.

Simitses, George J.↗

DOE OSTI2019

Estimation of conductivity for nanoporous materials

Methods and systems for estimating conductivity of clay mineral systems, and for applying the estimates in larger-scale analysis. Conductivity of the clay may be estimated by constructing a molecular model of an anhydrous charge-neutral clay, and then assigning a charge density by substitution of ions in the model of the clay structure. Counterions are inserted for charge neutrality, and water molecules are added to the model to reflect a selected level of hydration. Following assignment of force-field coefficients, molecular dynamics simulation provides data from which diffusion coefficients can be estimated. Application of the Nernst-Einstein relationship to the diffusion coefficients of the counterions provides the ion conductivity of the clay system. This conductivity can be used to derive a formation factor, and can be applied in direct numerical simulation analysis.

58 GEOSCIENCES↗

DOE OSTI2025

Remote quantification of Cm(III) and HNO 3 by fluorescence spectroscopy and chemometrics

A unique approach to remotely quantify Cm(III) (0–100 µg mL −1 ) in HNO 3 (1–12 M) using steady-state laser fluorescence spectroscopy and multivariate regression models was developed. Photoluminescence is amenable to remote measurements using fiber-optic cables and is sensitive to numerous lanthanide and actinide species. In-line measurements can provide feedback to support complex processing in harsh environments (e.g., hot cells) to help guide and optimize radiochemical separations. In this work, Cm(III) spectra were acquired remotely in a glove box as a function of HNO 3 concentration to better understand spectral characteristics and evaluate the utility of multivariate regression models in this system. Furthermore, the Cm(III) fluorescence peak shape, width, position, and intensity changed significantly as a function of HNO 3 concentration, likely because of the displacement of emission quenching inner-sphere water molecules and complexation with nitrate ions. Despite significant covariance and nonlinearity in the data, a D-optimal design strategy successfully minimized training set sample size and was used to build effective partial least squares regression models for Cm(III) and HNO 3 concentrations without a priori knowledge of solution conditions. Chemometrics for modeling complex fluorescence spectra are promising and may find widespread applicability for online analysis in numerous chemical systems found in the nuclear field.

Actinide↗

DOE OSTI2017

Interface tracking simulations of bubbly flows in PWR relevant geometries

[en] Highlights: • Simulations were performed for turbulent bubbly flows in PWR subchannel geometry. • Liquid turbulence is fully resolved by direct numerical simulation approach. • Bubble behavior is captured using level-set interface tracking method. • Time-averaged single- and two-phase turbulent flow statistical quantities are obtained. - Abstract: The advances in high performance computing (HPC) have allowed direct numerical simulation (DNS) approach coupled with interface tracking methods (ITM) to perform high fidelity simulations of turbulent bubbly flows in various complex geometries. In this work, we have chosen the geometry of the pressurized water reactor (PWR) core subchannel to perform a set of interface tracking simulations (ITS) with fully resolved liquid turbulence. The presented research utilizes a massively parallel finite-element based code, PHASTA, for the subchannel geometry simulations of bubbly flow turbulence. The main objective for this research is to demonstrate the ITS capabilities in gaining new insight into bubble/turbulence interactions and assisting the development of improved closure laws for multiphase computational fluid dynamics (M-CFD). Both single- and two-phase turbulent flows were studied within a single PWR subchannel. The analysis of numerical results includes the mean gas and liquid velocity profiles, void fraction distribution and turbulent kinetic energy profiles. Two sets of flow rates and bubble sizes were used in the simulations. The chosen flow rates corresponded to the Reynolds numbers of 29,079 and 80,775 based on channel hydraulic diameter (D_h) and mean velocity. The finite element unstructured grids utilized for these simulations include 53.8 million and 1.11 billion elements, respectively. This has allowed to fully resolve all the turbulence scales and the deformable interfaces of individual bubbles. For the two-phase flow simulations, a 1% bubble volume fraction was used which resulted in 17 bubbles in the smaller case and 262 bubbles in the larger case. In the larger simulation case the size of the resolved bubbles is 0.65 mm in diameter, and the bulk mesh cell size is about 30 microns. Those large-scale simulations provide new level of details previously unavailable and were enabled by the excellent scaling performance of our two-phase flow solver and access to the state-of-the-art supercomputing resources. The presented simulations used up to 256 thousand processing threads on the IBM BG/Q supercomputer “Mira” (Argonne National Laboratory).

Nuclear Science & Technology↗

DOE OSTI2015

Rubble masonry response under cyclic actions: The experience of L’Aquila city (Italy)

Several methods of analysis are available in engineering practice to study old masonry constructions. Two commonly used approaches in the field of seismic engineering are global and local analyses. Despite several years of research in this field, the various methodologies suffer from a lack of comprehensive experimental validation. This is mainly due to the difficulty in simulating the many different kinds of masonry and, accordingly, the non-linear response under horizontal actions. This issue can be addressed by examining the local response of isolated panels under monotonic and/or alternate actions. Different testing methodologies are commonly used to identify the local response of old masonry. These range from simplified pull-out tests to sophisticated in-plane monotonic tests. However, there is a lack of both knowledge and critical comparison between experimental validations and numerical simulations. This is mainly due to the difficulties in implementing irregular settings within both simplified and advanced numerical analyses. Similarly, the simulation of degradation effects within laboratory tests is difficult with respect to old masonry in-situ boundary conditions. Numerical models, particularly on rubble masonry, are commonly simplified. They are mainly based on a kinematic chain of rigid blocks able to perform different “modes of damage” of structures subjected to horizontal actions. This paper presents an innovative methodology for testing; its aim is to identify a simplified model for out-of-plane response of rubbleworks with respect to the experimental evidence. The case study of L’Aquila district is discussed.

58 GEOSCIENCES↗

DOE OSTI2015

The effects of phase decorrelation on the dynamics of the bump-on-tail instability

The dynamics of the bump-on-tail instability has been studied. The novelty of the work is the analysis of how the bump-on-tail dynamics is affected by an extrinsic stochastisation of the phase of the wave-particle interaction; here referred to as phase decorrelation. For this purpose, a nonlinear Monte Carlo model has been developed. When the characteristic time scale for macroscopic phase decorrelation becomes shorter than time scales of nonlinear wave-particle dynamics, the system may be described quasilinearly, with the phase decorrelation being replaced by a quasilinear diffusion coefficient in particle energy. A purely quasilinear Monte Carlo model, which is typically less computationally demanding than the fully nonlinear description due to the reduced dimensionality of phase space, has been developed for comparison. In this paper, parameter regimes, where the nonlinear and the quasilinear descriptions quantitatively agree on a macroscopic level, have been investigated, using combined theoretical and numerical analyses. Qualitative effects on the macroscopic dynamics by the presence of phase decorrelation and/or by structures of the energy distribution function in the proximity of the wave-particle resonance are also studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2016

Survey of Galileo Plasma Observations in Jupiter's Plasma Sheet

The plasma science (PLS) Instrument on the Galileo spacecraft (orbiting Jupiter from December 1995 to September 2003) measured properties of the ions that were trapped in the magnetic field. The PLS data provide a survey of the plasma properties between approx. 5 and 30 Jupiter radii [R(sub J)] in the equatorial region. We present plasma properties derived via two analysis methods: numerical moments and forward modeling. We find that the density decreases with radial distance by nearly 5 orders of magnitude from approx. 2 to 3000 cm(exp.-3) at 6R(sub j) to approx. 0.05cm(sub -3) at 30 R(sub j). The density profile did not show major changes from orbit to orbit, suggesting that the plasma production and transport remained constant within about a factor of 2. The radial profile of ion temperature increased with distance which implied that contrary to the concept of adiabatic cooling on expansion, the plasma heats up as it expands out from Io's orbit (where TI is approx.60-80 eV) at approx. 6R(sub j) to a few keV at 30R(sub j).There does not seem to be a long-term, systematic variation in ion temperature with either local time or longitude. This latter finding differs from earlier analysis of Galileo PLS data from a selection of orbits. Further examination of all data from all Galileo orbits suggests that System Ill variations are transitory on timescales of weeks, consistent with the modeling of Cassini Ultraviolet Imaging Spectrograph observations. The plasma flow is dominated by azimuthal flow that is between 80% and 100% of corotation out to 25 R(sub j).

Bagenal, Fran↗

DOE OSTI2014

Quasi-static three-dimensional magnetic field evolution in solar active region NOAA 11166 associated with an X1.5 flare

We study the quasi-static evolution of coronal magnetic fields constructed from the non-linear force-free field (NLFFF) approximation aiming to understand the relation between the magnetic field topology and ribbon emission during an X1.5 flare in active region (AR) NOAA 11166. The flare with a quasi-elliptical and two remote ribbons occurred on 2011 March 9 at 23:13 UT over a positive flux region surrounded by negative flux at the center of the bipolar AR. Our analysis of the coronal magnetic structure with potential and NLFFF solutions unveiled the existence of a single magnetic null point associated with a fan-spine topology and is co-spatial with the hard X-ray source. The footpoints of the fan separatrix surface agree with the inner edge of the quasi-elliptical ribbon and the outer spine is linked to one of the remote ribbons. During the evolution, the slow footpoint motions stressed the field lines along the polarity inversion line and caused electric current layers in the corona around the fan separatrix surface. These current layers trigger magnetic reconnection as a consequence of dissipating currents, which are visible as cusp-shaped structures at lower heights. The reconnection process reorganized the magnetic field topology whose signatures are observed at the separatrices/quasi-separatrix layer structure in both the photosphere and the corona during the pre-to-post flare evolution. In agreement with previous numerical studies, our results suggest that the line-tied footpoint motions perturb the fan-spine system and cause null point reconnection, which eventually causes the flare emission at the footpoints of the field lines.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2002

Study of Fluid Flow Control in Protein Crystallization using Strong Magnetic Fields

An important component in biotechnology, particularly in the area of protein engineering and rational drug design is the knowledge of the precise three-dimensional molecular structure of proteins. The quality of structural information obtained from X-ray diffraction methods is directly dependent on the degree of perfection of the protein crystals. As a consequence, the growth of high quality macromolecular crystals for diffraction analyses has been the central focus for biochemists, biologists, and bioengineers. Macromolecular crystals are obtained from solutions that contain the crystallizing species in equilibrium with higher aggregates, ions, precipitants, other possible phases of the protein, foreign particles, the walls of the container, and a likely host of other impurities. By changing transport modes in general, i.e., reduction of convection and sedimentation, as is achieved in "microgravity", researchers have been able to dramatically affect the movement and distribution of macromolecules in the fluid, and thus their transport, formation of crystal nuclei, and adsorption to the crystal surface. While a limited number of high quality crystals from space flights have been obtained, as the recent National Research Council (NRC) review of the NASA microgravity crystallization program pointed out, the scientific approach and research in crystallization of proteins has been mainly empirical yielding inconclusive results. We postulate that we can reduce convection in ground-based experiments and we can understand the different aspects of convection control through the use of strong magnetic fields and field gradients. Whether this limited convection in a magnetic field will provide the environment for the growth of high quality crystals is still a matter of conjecture that our research will address. The approach exploits the variation of fluid magnetic susceptibility with concentration for this purpose and the convective damping is realized by appropriately positioning the crystal growth cell so that the magnetic susceptibility force counteracts terrestrial gravity. The general objective is to test the hypothesis of convective control using a strong magnetic field and magnetic field gradient and to understand the nature of the various forces that come into play. Specifically we aim to delineate causative factors and to quantify them through experiments, analysis and numerical modeling. Once the basic understanding is obtained, the study will focus on testing the hypothesis on proteins of pyruvate dehydrogenase complex (PDC), proteins E1 and E3. Obtaining high crystal quality of these proteins is of great importance to structural biologists since their structures need to be determined. Specific goals for the investigation are: 1. To develop an understanding of convection control in diamagnetic fluids with concentration gradients through experimentation and numerical modeling. Specifically solutal buoyancy driven convection due to crystal growth will be considered. 2. To develop predictive measures for successful crystallization in a magnetic field using analyses and numerical modeling for use in future protein crystal growth experiments. This will establish criteria that can be used to estimate the efficacy of magnetic field flow damping on crystallization of candidate proteins. 3. To demonstrate the understanding of convection damping by high magnetic fields to a class of proteins that is of interest and whose structure is as yet not determined. 4. To compare quantitatively, the quality of the grown crystals with and without a magnetic field. X-ray diffraction techniques will be used for the comparative studies. In a preliminary set of experiments, we studied crystal dissolution effects in a 5 Tesla magnet available at NASA Marshall Space Flight Center (MSFC). Using a Schlieren setup, a 1mm crystal of Alum (Aluminum-Potassium Sulfate) was introduced in a 75% saturated solution and the resulting dissolution plume was observed. The experiment was conducted both in the presence and absence of a magnetic field gradient. The magnet produces a gradient field of approx. 1 Tesla2/cm. Image analysis of the recorded images indicated an enhanced plume velocity that was of the order of the measurement limit. For this experiment, both the gradient and gravity fields are in the same direction resulting in an enhanced effective gravity that tends to accelerate the observed plume velocity. While the results are not conclusive, pending further tests, it clearly points out the inadequacy of the MSFC magnet for conducting protein crystallization experiments and the need for a stronger magnet. In spacebased experiments, however, where the gravitational effects are small, only a weak magnetic field will be required to control or mitigate the effects of convective contamination.

Ramachandran, Narayanan↗

DOE OSTI2024

Liquid lithium divertor analysis using coupled plasma material interaction model

A liquid lithium divertor can improve performance of future fusion devices by creating efficient power exhaust and improving the energy confinement via pumping of the hydrogen isotopes. In addition, significantly higher heat fluxes can be handled if controlled vapor shielding is used to redistribute the divertor heat flux over a wider area. Design and optimization of such a system calls for an analysis model which includes a strong two-way coupling between the plasma and divertor material. The incoming plasma heat and particle flux will affect the divertor surface temperature, which is a defining factor of the lithium evaporative and sputtered flux going into the plasma. Results of the coupled model based on the plasma edge code SOLPS-ITER and the computational fluid dynamics (CFD) code ANSYS-CFX will be presented for different configurations. An analytical slab flow model is used as a heat transfer boundary condition for SOLPS, defining particle flux from the wall via calculation of the surface temperature. At the final step, results of the SOLPS analysis are verified using a 3D CFD magnetohydrodynamics (MHD) analysis which uses heat and particle flux from SOLPS as a boundary condition. In addition to plasma heat flux, both analytical and CFD temperature models include several plasma material interaction effects, such as lithium evaporation, condensation and sputtering based on deuterium target flux. New adatom sputtering model based on the available experimental data is presented. Analytical model is expanded to include free surface axisymmetric configurations. Results of parametric studies of the divertor configurations with different lithium inlet temperature and velocity will be presented leading to the optimal design resulting in the lowest possible lithium contamination in the core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Fuel Thermo-physical Characterization Project: Evaluation of Models to Calculate Thermal Diffusivity of Layered Composites

The Office of Material Management and Minimization Fuel Thermo-physical Characterization Project at Pacific Northwest National Laboratory (PNNL) is tasked with using PNNL facilities and processes to receive irradiated low enriched uranium-molybdenum fuel plate samples and perform analyses in support of the Office of Material Management and Minimization Reactor Conversion Program. This work is in support of the Fuel Development Pillar that is managed by Idaho National Laboratory. A key portion of the scope associated with this project was to measure the thermal properties of fuel segments harvested from plates that were irradiated in the Advanced Test Reactor. Thermal diffusivity of samples prepared from the fuel segments was measured using laser flash analysis. Two models, one developed by PNNL and the other developed by the Technische Universität München (TUM), were evaluated to extract the thermal diffusivity of the uranium-molybdenum alloy from measurements made on the irradiated, layered composites. The experimental data of the “TC” irradiated fuel segment was evaluated using both models considering a three-layer and five-layer system. Both models are in acceptable agreement with one another and indicate that the zirconium diffusion barrier has a minimal impact on the overall thermal diffusivity of the monolithic U-Mo fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1994

Zero potential vorticity envelopes for the zonal-mean velocity of the Venus/Titan atmospheres

The diagnostic analysis of numerical simulations of the Venus/Titan wind regime reveals an overlooked constraint upon the latitudinal structure of their zonal-mean angular momentum. The numerical experiments, as well as the limited planetary observations, are approximately consistent with the hypothesis that within the latitudes bounded by the wind maxima the total Ertel potential vorticity associated with the zonal-mean motion is approximately well mixed with respect to the neutral equatorial value for a stable circulation. The implied latitudinal profile of angular momentum is of the form M equal to or less than M(sub e)(cos lambda)(exp 2/Ri), where lambda is the latitude and Ri the local Richardson number, generally intermediate between the two extremes of uniform angular momentum (Ri approaches infinity) and uniform angular velocity (Ri = 1). The full range of angular momentum profile variation appears to be realized within the observed meridional - vertical structure of the Venus atmosphere, at least crudely approaching the implied relationship between stratification and zonal velocity there. While not itself indicative of a particular eddy mechanism or specific to atmospheric superrotation, the zero potential vorticity (ZPV) constraint represents a limiting bound for the eddy - mean flow adjustment of a neutrally stable baroclinic circulation and may be usefully applied to the diagnostic analysis of future remote sounding and in situ measurements from planetary spacecraft.

Allison, Michael↗

DOE OSTI2021

Theory of Trotter Error with Commutator Scaling

The Lie-Trotter formula, together with its higher-order generalizations, provides a simple approach to decomposing the exponential of a sum of operators. Despite significant effort, the error scaling of such product formulas remains poorly understood. We develop a theory of Trotter error that overcomes the limitations of truncating the Baker-Campbell-Hausdorff expansion. Our analysis directly exploits the commutativity of operator summands, producing tighter error bounds for both real- and imaginary-time evolutions. Whereas previous work achieves similar goals for systems with geometric locality or Lie-algebraic structure, our approach holds in general. We give a host of improved algorithms for digital quantum simulation and quantum Monte Carlo methods, nearly matching or even outperforming the best previous results. Our applications include: (i) a simulation of second-quantized plane-wave electronic structure, nearly matching the interaction-picture algorithm of Low and Wiebe; (ii) a simulation of $k$-local Hamiltonians almost with induced one-norm scaling, faster than the qubitization algorithm of Low and Chuang; (iii) a simulation of rapidly decaying power-law interactions, outperforming the Lieb-Robinson-based approach of Tran et al.; (iv) a hybrid simulation of clustered Hamiltonians, dramatically improving the result of Peng, Harrow, Ozols, and Wu; and (v) quantum Monte Carlo simulations of the transverse field Ising model and quantum ferromagnets, tightening previous analyses of Bravyi and Gosset. We obtain further speedups using the fact that product formulas can preserve the locality of the simulated system. Specifically, we show that local observables can be simulated with complexity independent of the system size for power-law interacting systems, which implies a Lieb-Robinson bound nearly matching a recent result of Tran et al. Our analysis reproduces known tight bounds for first- and second-order formulas. We further investigate the tightness of our bounds for higher-order formulas. For quantum simulation of a one-dimensional Heisenberg model with an even-odd ordering of terms, our result overestimates the complexity by only a factor of $5$. Our bound is also close to tight for power-law interactions and other orderings of terms. This suggests that our theory can accurately characterize Trotter error in terms of both the asymptotic scaling and the constant prefactor.

quantum computing, numerical analysis↗

DOE OSTI2019

Development of a Novel Accelerator for Neutron Transport Solution Using the Galerkin Spectral Element Methods (Final Report)

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to develop and implement advanced acceleration schemes for the DOE NEAMS neutronics code PROTEUSSN. The project team included the University of Massachusetts Lowell, The Ohio State University, University of Michigan, and Argonne National Laboratory. The PROTEUS code is a high-fidelity capable deterministic neutron transport code based on unstructured finite element meshes, which solves the steady-state and transient neutron transport problem using the 2nd-order discrete ordinate method (SN2ND), the method of characteristics (MOC), and the advanced nodal transport method (NODAL). The DSA scheme has been implemented in PROTEUS to speed up SN transport calculations. The existing DSA scheme employs a consistent finite element formulation on the same fine mesh structure as the SN solution While the use of consistent discretizations makes the DSA method effective, it does not necessarily make it efficient because the numerical solution of a discretized elliptic diffusion problem itself can be costly, particularly when the problem size becomes large. We developed a new discontinuous Galerkin (DG) discretization of the diffusion equation, called DG-DSA, which can effectively and efficiently accelerate the SN transport iterations. As compared with the previous work, the novelty of our method is that the diffusion equation is solved on a coarse-mesh grid using the DG methods, and the DG diffusion discretization incorporates local hp adaptation, i.e., local adaptation of mesh size and/or polynomial degree, based on local total cross section (or optical thickness). Therefore, the resulting number of degrees of freedom (DOF) of the DG discretization is much less than the conventional consistent DSA discretizations, and thus DG-DSA can achieve significant improvement in computational efficiency. We implemented this scheme in PROTEUS-SN. In addition, we developed a new stabilization scheme named linear prolongation CMFD (lpCMFD). A novel feature of this scheme is that the conventional flat flux ratio–based scaling approach is replaced with a linear interpolation of the scalar flux differences at the coarse-mesh cell edges between the neutron transport and CMFD calculations. Fourier convergence analysis and numerical results show that lpCMFD is unconditionally stable and more effective than other CMFD based on acceleration schemes such as pCMFD, odCMFD. We also developed a new nonlinear diffusion acceleration scheme for solving neutron transport equations. This scheme, called LR-NDA, employs a local refinement approach on the framework of CMFD by solving a local boundary value problem of the scalar flux on the coarse-mesh structure to replace the piecewise constant scalar flux obtained by CMFD. The refined flux is then used to update the scalar flux in the neutron transport source iteration. It has been demonstrated that that LR-NDA is much more effective and stable than CMFD for a wide range of optical thicknesses. LR-NDA is a local adaptive method, which means LR-NDA does not necessarily require local refinement for all the coarse-mesh cells on the problem domain, i.e., it can be used only for relatively optically thick regions where the standard CMFD scheme would encounter the convergence problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1996

Development and Application of Viscoplastic Models in Nonlinear Structural Analysis

A number of numerical solution technologies were developed for advanced analysis capabilities involving the finite element method in advanced constitutive modeling and structural analysis for improved and rational designs of aerospace propulsive systems. These technologies were incorporated into several advanced viscoplastic models and were applied to a wide spectrum of structural engineering problems involving extremes of thermal/mechanical loading. Results of the research performed under the grant were presented at a number of scientific meetings and conferences and have resulted in numerous publications.

Keith, Theo G., Jr.↗

DOE OSTI2021

Midfidelity Model Verification for a Point-Absorbing Wave Energy Converter with Linear Power Take-Off

In the preliminary design stage of a wave energy converter (WEC), fast and reliable simulation tools are required. High-fidelity numerical models are usually employed to study the wave-structure interaction, yet the computational cost is demanding. Instead, mid-fidelity models provide simulations in the order of real time. In this study, Uppsala University WEC operates in a relatively mild sea state and is modeled using WEC-Sim. The model is verified based on OpenFOAM simulations. To analyze the ability of mid-fidelity model to capture the WEC dynamics, the system is investigated separately in 1, 2 and 3 DoF. The contribution of viscous phenomena is examined. Moreover, linear and weakly non-linear solutions provided by WEC-Sim are studied. The results obtained indicate that the viscous effects in heave and surge motion can be neglected but not for the pitch. The weakly non-linear WEC-Sim solution successfully agrees with the CFD, whereas the linear solution could drive to misleading results.

CFD↗

DOE OSTI2021

Midfidelity Model Verification for a Point-Absorbing Wave Energy Converter with Linear Power Takeoff: Preprint

CFD↗

DOE OSTI2022

Midfidelity model verification for a point-absorbing wave energy converter with linear power take-off

In the preliminary design stage of a wave energy converter (WEC), researchers need fast and reliable simulation tools. High-fidelity numerical models are usually employed to study the wave-structure interaction, but the computational cost is demanding. As an alternative, midfidelity models can provide simulations in the order of real time. In this study, we operate Uppsala University’s WEC in a relatively mild sea state and model it using WEC-Sim. The model is verified based on OpenFOAM simulations. To analyze the ability of the midfidelity model to capture WEC dynamics, we investigate the system separately with 1, 2, and 3 degrees of freedom. We examine the contribution of viscous phenomena, and study both linear and weakly nonlinear solutions provided by WEC-Sim. Our results indicate that the viscous effects can be neglected in heave and surge motion, but not for pitch. We also find that the weakly nonlinear WEC-Sim solution successfully agrees with the computational fluid dynamics, whereas the linear solution could suggest misleading results.

16 TIDAL AND WAVE POWER↗