Modeling neutronic transients with Galerkin projection onto a greedy-sampled, POD subspace

2021 ◽  
Vol 162 ◽  
pp. 108487
Author(s):  
Rabab Elzohery ◽  
Jeremy Roberts
Keyword(s):  
Author(s):  
J. Marconi ◽  
P. Tiso ◽  
D. E. Quadrelli ◽  
F. Braghin

AbstractWe present an enhanced version of the parametric nonlinear reduced-order model for shape imperfections in structural dynamics we studied in a previous work. In this model, the total displacement is split between the one due to the presence of a shape defect and the one due to the motion of the structure. This allows to expand the two fields independently using different bases. The defected geometry is described by some user-defined displacement fields which can be embedded in the strain formulation. This way, a polynomial function of both the defect field and actual displacement field provides the nonlinear internal elastic forces. The latter can be thus expressed using tensors, and owning the reduction in size of the model given by a Galerkin projection, high simulation speedups can be achieved. We show that the adopted deformation framework, exploiting Neumann expansion in the definition of the strains, leads to better accuracy as compared to the previous work. Two numerical examples of a clamped beam and a MEMS gyroscope finally demonstrate the benefits of the method in terms of speed and increased accuracy.


2014 ◽  
Vol 50 (2) ◽  
pp. 613-616 ◽  
Author(s):  
Zifu Wang ◽  
Zuqi Tang ◽  
Thomas Henneron ◽  
Francis Piriou ◽  
Jean-Claude Mipo

2001 ◽  
Author(s):  
H. Gunes ◽  
K. Gocmen ◽  
L. Kavurmacioglu

Abstract The two-dimensional incompressible non-isothermal confined twin-jet flow has been numerically studied in the transitional flow regime by a finite volume technique. Results have been obtained for the velocity and temperature distributions close to the onset of temporal oscillations. Next, the proper orthogonal decomposition (POD) is applied to the instantaneous flow and temperature data to obtain POD-based basis functions for both velocity and temperature fields. These basis functions are capable to identify the coherent structures in the velocity and temperature fields. The low-dimensional Galerkin models of the full Navier-Stokes and energy equations are constructed by the Galerkin projection onto basis functions. Since the low-dimensional Galerkin models are much easier to analyze than the full governing equations, basic insights into important mechanisms of dynamically complex flow and heat transfer (e.g. flow instabilities) can be easily studied by these models. The numerical implications, the validity of the models and their performance characteristics are discussed.


AIAA Journal ◽  
2021 ◽  
pp. 1-17
Author(s):  
Patrick J. Blonigan ◽  
Francesco Rizzi ◽  
Micah Howard ◽  
Jeffrey A. Fike ◽  
Kevin T. Carlberg

Author(s):  
Nariman Ashrafi ◽  
Habib Karimi Haghighi

Stress analysis of Pseudo-Plastic flow between rotating cylinders is studied in the narrow gap limit. The Galerkin projection method is used to derive dynamical system from the conservation of mass and momentum equations. Flow parameters were obtained using IMSL and also verified by Mathematica Software. Stresses are computed in a wide range of the Pseudo-Plastic effects. Azimuthal stress was found to be far greater than other stress components. All stress components increased as Pseudo-Plasticity decreased. Furthermore, complete stress and viscosity maps are presented for different scenarios in the flow regime.


Author(s):  
Nariman Ashrafi

The nonlinear stability and bifurcation of the one-dimensional channel (Poiseuille) flow is examined for a Johnson-Segalman fluid. The velocity and stress are represented by orthonormal functions in the transverse direction to the flow. The flow field is obtained from the conservation and constitutive equations using the Galerkin projection method. Both inertia and normal stress effects are included. The stability picture is dramatically influenced by the viscosity ratio. The range of shear rate or Weissenberg number for which the base flow is unstable increases from zero as the fluid deviates from the Newtonian limit as decreases. Typically, two turning points are observed near the critical Weissenberg numbers. The transient response is heavily influenced by the level of inertia. It is found that the flow responds oscillatorily. When the Reynolds number is small, and monotonically at large Reynolds number when elastic effects are dominated by inertia.


2016 ◽  
Vol 800 ◽  
pp. 72-110 ◽  
Author(s):  
Richard Semaan ◽  
Pradeep Kumar ◽  
Marco Burnazzi ◽  
Gilles Tissot ◽  
Laurent Cordier ◽  
...  

We propose a hierarchy of low-dimensional proper orthogonal decomposition (POD) models for the transient and post-transient flow around a high-lift airfoil with unsteady Coanda blowing over the trailing edge. The modal expansion comprises actuation modes as a lifting method for wall actuation following Graham et al. (Intl J. Numer. Meth. Engng, vol. 44 (7), 1999, pp. 945–972) and Kasnakoğlu et al. (Intl J. Control, vol. 81 (9), 2008, pp. 1475–1492). A novel element is separate actuation modes for different frequencies. The structure of the dynamic model rests on a Galerkin projection using the Navier–Stokes equations, simplifying mean-field considerations, and a stochastic term representing the background turbulence. The model parameters are identified with a data assimilation (4D-Var) method. We propose a model hierarchy from a linear oscillator explaining the suppression of vortex shedding by blowing to a fully nonlinear model resolving unactuated and actuated transients with steady and high-frequency modulation of blowing. The models’ accuracy is assessed through the mode amplitudes and an estimator for the lift coefficient. The robustness of the model is physically justified, and then observed for the training and the validation dataset.


2017 ◽  
Vol 829 ◽  
pp. 345-363 ◽  
Author(s):  
Jason R. Picardo ◽  
R. Narayanan

In this study, we revisit Rayleigh’s visionary hypothesis (Rayleigh, Proc. R. Soc. Lond., vol. 29, 1879a, pp. 71–97), that patterns resulting from interfacial instabilities are dominated by the fastest-growing linear mode, as we study nonlinear pattern selection in the context of a linear growth (dispersion) curve that has two peaks of equal height. Such a system is obtained in a physical situation consisting of two liquid layers suspended from a heated ceiling, and exposed to a passive gas. Both interfaces are then susceptible to thermocapillary and Rayleigh–Taylor instabilities, which lead to rupture/pinch off via a subcritical bifurcation. The corresponding mathematical model consists of long-wavelength evolution equations which are amenable to extensive numerical exploration. We find that, despite having equal linear growth rates, either one of the peak-modes can completely dominate the other as a result of nonlinear interactions. Importantly, the dominant peak continues to dictate the pattern even when its growth rate is made slightly smaller, thereby providing a definite counter-example to Rayleigh’s conjecture. Although quite complex, the qualitative features of the peak-mode interaction are successfully captured by a low-order three-mode ordinary differential equation model based on truncated Galerkin projection. Far from being governed by simple linear theory, the final pattern is sensitive even to the phase difference between peak-mode perturbations. For sufficiently long domains, this phase effect is shown to result in the emergence of coexisting patterns, wherein each peak-mode dominates in a different region of the domain.


Sign in / Sign up

Export Citation Format

Share Document