Numerical validation of the volume penalization method in three-dimensional pseudo-spectral simulations

Periodic boundary condition is commonly used in the FE homogenization method for predicting the thermo-mechanical properties of composites. However, its existing numerical algorithms are generally inefficient. Linear displacement boundary condition is much easier and more efficient to implement for the Representative Volume Elements (RVEs) with conformal and nonconformal surface meshes. Aiming at replacing the periodic boundary condition by the linear displacement boundary condition to easily and efficiently predict the thermo-mechanical properties of composites, this paper compares the implementation efficiency and prediction accuracy of the periodic and linear displacement boundary conditions. Regarding the prediction of the elastic properties of composites, the balance of micro-energy and macro-energy of the RVEs imposed the periodic and linear displacement boundary conditions is checked, and the results demonstrate that these two boundary conditions guarantee the Hill lemma (i.e. the balance of micro-energy and macro-energy of the RVEs). Moreover, comparing with the experimental test and the multi-step micro-mechanical method, the effectiveness of the periodic and linear displacement boundary conditions to predict the elastic properties of composites is validated.

Abstract. Atmospheric convection over the subtropical western North Pacific (SWNP) during boreal summer varies on timescales of around 10 d, with significant effects on both local and remote circulation. One of less understood effects of this variability is its coupling with equatorial wave dynamics. This paper quantifies equatorial wave perturbations and their evolution throughout the lifecycle of SWNP convection using wave-space regression between outgoing longwave radiation over the SWNP region and spectral expansion coefficients of global tropospheric circulation from ERA5 reanalyses. The regression distinguishes between convection-coupled linear Rossby and Kelvin waves, and mixed Rossby-gravity (MRG) and inertia-gravity (IG) waves. The former two correspond to the Gill-type tropical circulation response to asymmetric heating. The results show that MRG and IG waves exhibit amplitudes comparable to those of the Gill-response component in the upper troposphere. In particular, MRG and IG waves dominate the cross-equatorial northerly flow over the Maritime Continent, with MRG waves playing the larger role. These findings suggest caution when applying the Gill solution to interpret circulation responses to asymmetric heating in model simulations. As SWNP convection intensifies, the MRG-wave northerly winds across the equator strengthen, while IG waves represent enhanced upper-tropospheric outflow over the SWNP region. By contrast, the combined effect of Kelvin and Rossby waves reinforces circulation on the equatorward flank of the anticyclone over the SWNP region, with Rossby wave easterlies being about three times stronger than those associated with Kelvin waves. The Rossby wave signal resembles the n=1 Rossby wave, with its Southern Hemisphere (SH) subtropical anticyclonic gyre forming over the southern Indian ocean during the decay phase of the SWNP convection. This gyre, together with the SH IG meridional flow provides a dynamical bridge linking SWNP convection and extratropical circulation during austral winter.

The planar flow of a Newtonian incompressible fluid in a T-shaped channel is investigated. Three fluid interaction models with solid walls are considered: no slip boundary condition, Navier slip boundary condition and slip boundary condition with slip yield stress. The fluid flow is provided by uniform pressure profiles at the boundary sections of the channel. The problem is numerically solved using a finite difference method based on the SIMPLE procedure. Characteristic flow regimes have been found for the described models of liquid interaction with solid walls. The estimation of the influence of the Reynolds number, pressure applied to the boundary sections and the parameters of these models on the flow pattern was performed. The criterial dependences describing main characteristics of the flow under conditions of the present work have been demonstrated.

A formal framework is established for the way in which cloud-resolving numerical models are used to investigate the role of precipitating cloud systems in climate and weather forecasting models. Emphasis is on models with periodic lateral boundary conditions that eliminate unrealistic numerically generated circulations caused by open boundary conditions in long-term simulations. Defined in this formalism is the concept of large-scale forcing and the cloud-environment interactions that are consistent with the periodic boundary conditions. Two-dimensional numerical simulations of the evolution of cloud systems during 1–7 September 1974 in Phase III of the Global Atmospheric Research Program Atlantic Tropical Experiment (GATE) are conducted. Based on the above formalism, a simple technique is used to force an anelastic cloud-resolving model with evolving large-scale horizontal wind field and large-scale forcing for the temperature and moisture obtained from the GATE data. The 7-day period selected i...

We have developed a rigorous self-consistent approach for the quantization of electromagnetic field in inhomogeneous structures. The approach is based on utilization of the scattering matrix of the system. Instead of the use of standard periodic Born-Karman boundary conditions, we use the quantization condition implying equating eigenvalues of the scattering matrix (S-matrix) of the system to unity (S-quantization). In the trivial case of uniform medium boundary condition for S-quantization is nothing but periodic boundary condition. S-quantization allows calculating modification of the spontaneous emission rate for arbitrary inhomogeneous structure and direction of the emitted radiation. S-quantization solves the long-standing problem coupled to normalization of the quasi-stationary electromagnetic modes. Examples of application of S-quantization for the calculation of spontaneous emission rate for the cases of Bragg reflector and microcavity are demonstrated.

We study finite-size corrections to the free energy of free-fermion models on a torus with periodic, twisted, and fixed boundary conditions. Inside the critical (striped-incommensurate) phase, the free energy densityf(N, M) on anN×M square lattice with periodic (or twisted) boundary conditions scales asf(N, M)=f ∞−A(s)/(NM)+.... We derive exactly the finite-size-scaling (FSS) amplitudesA(s) as a function of the aspect ratios=M/N. These amplitudes are universal because they do not depend on details of the free-fermion Hamiltonian. We establish an equivalence between the FSS amplitudes of the free-fermion model and the Coulomb gas system with electric and magnetic defect lines. The twist angle generates magnetic defect lines, while electric defect lines are generated by competition between domain wall separation and system size. The FSS behavior of the free-fermion model is consistent with predictions of the theory of conformal invariance with the conformal chargec=l. For instance, the FSS amplitude on an infinite cylinder with fixed boundary conditions is found to be one-quarter of that with periodic boundary conditions. Finally, we conjecture the exact form of the FSS amplitudes for an interacting-fermion model on a torus. Numerical calculations employing the Bethe Ansatz confirm our conjecture in the infinite-cylinder limit.

We illustrate the effect of boundary conditions on the evolution of structure in fuzzy dark matter. Scenarios explored include the evolution of single, ground-state equilibrium solutions of the Schr\"odinger-Poisson system, relaxation of a Gaussian density fluctuation, mergers of two equilibrium configurations, and the random merger of many solitons. For comparison, each scenario is evolved twice, with isolation vs periodic boundary conditions, the two commonly used to simulate isolated systems and structure formation, respectively. Replacing isolation boundary conditions (implemented by an absorbing sponge at large radius) by periodic boundary conditions changes the domain topology and dynamics of each scenario by affecting the outcome of gravitational cooling. With periodic boundary conditions, the initial ground-state equilibrium solution and Gaussian fluctuation each evolve toward the single equilibrium solitonic core of the isolated case, but surrounded by an envelope, or tail, in which additional mass is distributed nearly uniformly, unlike the isolated versions. The case of head-on, binary merger introduces additional effects caused by the pull on the system due to the infinite network of periodic images along each axis of the domain. Adding angular momentum to this binary merger results in a tail with polynomial profile when using a periodic domain. Finally, the 3D merger of many, randomly placed solitonic cores of different mass makes a solitonic core surrounded by a tail with power-law-like profile, for periodic boundary conditions, while producing a core with a much sharper falloff in the isolated case. This suggests the conclusion of earlier work, that the ground-state equilibrium solution is an attractor for the asymptotic state, is true even in 3D and for more general circumstances than previously considered, but only if gravitational cooling is able to carry mass and energy off to infinity, which isolation boundary conditions allow, but periodic ones do not.

<p style='text-indent:20px;'>In this paper, the chaotic oscillations of the initial-boundary value problem of linear hyperbolic partial differential equation (PDE) with variable coefficients are investigated, where both ends of boundary conditions are nonlinear implicit boundary conditions (IBCs). It separately considers that IBCs can be expressed by general nonlinear boundary conditions (NBCs) and cannot be expressed by explicit boundary conditions (EBCs). Finally, numerical examples verify the effectiveness of theoretical prediction.</p>

Transient-free convection in a porous enclosure having heat-conducting solid walls of finite thickness under conditions of convective heat exchange with an environment was studied numerically. A heat source of constant temperature was located at the bottom of the cavity. The governing equations in porous volume formulated in dimensionless variables such as the temperature and vector potential functions within the Darcy–Boussinesq approach and the transient three-dimensional heat conduction equation based on the Fourier hypothesis for solid walls with corresponding initial and boundary conditions were solved using an iterative implicit finite-difference method. The main objective was to investigate the influence of the Rayleigh number 103 ≤ Ra ≤ 106, the Darcy number 10−5 ≤ Da ≤ 10−3, the thermal conductivity ratio 1 ≤ k1,2 ≤ 20, the solid wall thickness ratio 0.1 ≤ l/L ≤ 0.3, and the dimensionless time 0 ≤ τ ≤ 200 on the fluid flow and heat transfer. Comprehensive analysis of the effects of these key parameters on the average Nusselt number at the heat source surface was conducted.

Entropy generation due to conjugate natural convection in enclosures bounded by vertical solid walls with different thicknesses

Entropy analysis due to conjugate-buoyant flow in a right-angle trapezoidal enclosure filled with a porous medium bounded by a solid vertical wall

Cellular automaton model within the fundamental-diagram approach reproducing some findings of the three-phase theory

In this article, a new cellular automaton traffic flow model, in which randomisation effect is enhanced with the decrease of time gap and the long time stopped vehicle has large randomisation probability, is presented. Both periodic and open boundary conditions are considered. Under periodic boundary condition, the fundamental diagram and the properties of synchronised flow are studied in detail. Under open boundary condition, different congested patterns induced by an on-ramp are analysed. The results are consistent with the well-known results of Kerner's three-phase traffic theory published before.

A three-dimensional unsteady compressible viscous flow solver has been developed, and applied to a low speed centrifugal compressor configuration. Good agreement with experiments have been obtained for a compressor tested at NASA Lewis Research Center. Comparisons with measured surface pressures, and performance data are given. For this compressor, stall was found to occur in the diffuser region. It is demonstrated that this stall may be eliminated, and stable operation may be restored by the use of bleed valves located on the diffuser walls. INTRODUCTION In many army applications (tanks and rotorcraft) small gas turbine engines are often used. The mass flow rate in these engines is small (below 10 Ibs/sec). The compression is provided by one or two centrifugal compressor stages, sometimes in combination with an axial compressor stage. Centrifugal compressors are well suited for handling such low mass flow rates. They are rugged and compact. However, the geometric complexity of the blades and the twisting blade passages give rise to a host of aerodynamic problems such as unsteadiness, establishment of cross flows, flow separation at the trailing edge of the compressor blade, separation over diffuser guide vanes, and unsteady meandering motion of a 'tip vortex' that originates in the vicinity of the gap between the rotor and the casing. In many instances, this unsteadiness leads to rotating stall and surge. Considerable interest exists in the engine community in understanding and eventually controlling these phenomena, particularly compressor stall and surge. Computational fluid dynamics methods provide an efficient way of studying the complex flow phenomena within centrifugal compressors. Thanks to the availability of post-processing and scientific visualization software, the computed flow fields can be extensively analyzed. It is also convenient to test within a numerical model various control strategies such as bleeding, energizing of incoming flow using jets, excitation /wiggling of guide vanes, etc. Passive mechanisms such as slots, slats and flaps can also be easily studied within the scope of a numerical method. These simulations are thus a useful and economical first step to detailed flow control experiments. Finally, CFD simulations provide a wealth of numerical and flow physics information which are useful in developing reduced order models. For these reasons, a 3-D CFD code capable of analyzing centrifugal compressor configurations has been developed by the present investigators. Rather than solving the flow equations in a rotating coordinate system, the present calculations are carried out in an inertial system, with grid velocity terms used to simulate a rotor. Inflow and outflow boundary conditions appropriate for internal flows in general and centrifugal compressors in particular are used. This paper is organized as follows. The flow solution algorithm is first described. Next, results in the form of velocity fields, surface pressure data and performance map are presented for a NASA low-speed centrifugal compressor. The performance of this compressor for off-design low mass flow ratio conditions is subsequently studied. Methods for improving the compressor performance during such conditions are discussed. NUMERICAL FORMULATION The 3-D unsteady compressible Reynolds averaged Navier-Stokes equations are solved numerically using a time marching scheme. This involves solving the governing equations at each time step by marching in time from an initial flow condition with appropriate boundary conditions. The governing equations are cast in a curvilinear body-fitted coordinate system (£,T|,£) and solved using an Alternating Direction Implicit (ADI) scheme 1. This implicit scheme has been validated for a number of fixed wing and rotary wing configurations both in steady and unsteady flight conditions^'*. The partial derivatives in the governing equations are calculated using finite difference approximations. The Graduate Research Assistant, Student Member AIAA. * Regents' professor, Associate Fellow AIAA. Copyright© 1997, American Institute of Aeronautics and Astronautics, Inc. inviscid flux terms are calculated using a third order upwind scheme and the viscous stress terms are calculated using standard central difference representations. The details of the upwind formulation are given by Roe. The time derivative is approximated using a two point first order backward difference. The finite difference representation of the governing Navier-Stokes equations at time level 'n' is components at the inflow boundary were set to zero, assuming that there was no swirl at the inlet. It is possible within the framework of the present boundary conditions to prescribe swirl, non-uniform total pressure and total temperature at the inlet. At every time step, the speed of sound a and the normal component of velocity un were found at the inflow boundary from the relations: AT 2a 2a — u = n — u n = I 8 f R + S S + S , T (1) Interior where Aq* is the change in the flow property vector q during adjacent time levels and At is the time step. The operators 8£, 5r| and 8^ are the standard central difference operators. The viscous terms are evaluated explicitly and the inviscid flux vectors F, G and H are calculated implicitly at time level n+1. A three factor approximate factorization is used to factorize the implicit coefficient matrix operator into block tri-diagonal matrices. The factorized system of equations is solved using the Thomas algorithm. Turbulence effects are modeled using the SpalartAllmaras model. INITIAL AND BOUNDARY CONDITIONS The centrifugal compressor is shown in figure 1. The boundary grids are also shown for a single blade passage. The volume grid is made of three segments: inlet region, the rotor, and the diffuser. The following boundary conditions were applied: Periodicity: In the present study, the flow through only a single blade passage has been studied, assuming blade-toblade periodicity. At the periodic boundaries, the flow properties are computed by averaging the properties on either side of the periodic boundary. Solid walls: At the solid walls, the no slip boundary condition was used. The velocity on the walls of the inlet and the diffuser were set to zero, while the velocity for grid points on the compressor blades and the shaft were set equal to Q x r. The pressure, density and temperature values at the solid surfaces were extrapolated from the interior using relationships such as 9p/3n=dT/9n=3p/3n=0, where n represents the direction normal to the solid surface. Inflow Boundary: At this boundary, the stagnation temperature T0 and the total pressure p0 were assumed known. In the simulation of a compressor drawing the air from the atmosphere, these quantities may be set based on ambient conditions. In a multistage compressor analysis, these quantities will be known at the downstream boundary of the previous stage. The tangential (vand w-) (2)

We propose a new boundary condition for the particle swarm optimization (PSO) method. The new boundary is duped periodic boundary condition (PBC) and is inspired by a molecular dynamics (MD) formulation of the PSO. Based on this formulation, the PBC is introduced formally by postulating that the swarm environment is similar to crystal-bulk materials studied using MD. A combination of soft and hard boundary modes is employed in the implementation of the hybrid version of the proposed PBC. The hybrid PBC constitutes a substantial alteration of the original scenario assumed in theoretical physics but is allowed here to suite the PSO environment. When applied to standard test functions and the practical problem of linear antenna array synthesis, the hybrid PBC proved to be superior to previously established boundary conditions.