Strong and weak constraint variational assimilations for reduced order fluid flow modeling

In this work we propose and evaluate two variational data assimilation techniques for the estimation of low order surrogate experimental dynamical models for fluid flows. Both methods are built from optimal control recipes and rely on proper orthogonal decomposition and a Galerkin projection of the Navier Stokes equation. The techniques proposed differ in the control variables they involve. The first one introduces a weak dynamical model defined only up to an additional uncertainty time-dependent function whereas the second one, handles a strong dynamical constraint in which the dynamical system’s coefficients constitute the control variables. Both choices correspond to different approximations of the relation between the reduced basis on which is expressed the motion field and the basis components that have been neglected in the reduced order model construction. The techniques have been assessed on numerical data and for real experimental conditions with noisy particle image velocimetry data.     

Evolution of multilevel order in supramolecular assemblies

The process of self-assembly at multiple length scales of bis-urea substituted toluene on a Au(111) surface was studied by low temperature scanning tunneling microscopy. Pattern formation is controlled by specific hydrogen bonds between these molecules but also by significantly weaker lateral coupling between the resulting supramolecular polymers and a quasiepitaxial interlocking with the substrate. The ordered assemblies exhibit a tunnel transparency. Our experiments indicate the necessity of multiple interactions of different strengths for obtaining ordered structures with hierarchical levels of organization.     

Orientational order of molecular assemblies on inorganic crystals

Surfactant micelles form oriented arrays on crystalline substrates although registration is unexpected since the template unit cell is small compared to the size of a rodlike micelle. Interaction energy calculations based on molecular simulations reveal that orientational energy differences on a molecular scale are too small to explain matters. With atomic force microscopy, we show that orientational ordering is a dynamic, multimolecule process. Treating the cooperative processes as a balance between van der Waals torque on a large, rodlike micellar assembly and Brownian motion shows that orientation is favored.     

Prediction of displacement and stress fields of a notched panel with geometric nonlinearity by reduced order modeling

The focus of this investigation is on a first assessment of the predictive capabilities of nonlinear geometric reduced order models for the prediction of the large displacement and stress fields of panels with localized geometric defects, the case of a notch serving to exemplify the analysis. It is first demonstrated that the reduced order model of the notched panel does indeed provide a close match of the displacement and stress fields obtained from full finite element analyses for moderately large static and dynamic responses (peak displacement of 2 and 4 thicknesses). As might be expected, the reduced order model of the virgin panel would also yield a close approximation of the displacement field but not of the stress one. These observations then lead to two “enrichment” techniques seeking to superpose the notch effects on the virgin panel stress field so that a reduced order model of the latter can be used. A very good prediction of the full finite element stresses, for both static and dynamic analyses, is achieved with both enrichments.     

Uncertainty-based experimental validation of nonlinear reduced order models

This paper presents the first results of a combined experimental–computational investigation focused on the validation of reduced order models of geometrically nonlinear structures in the presence of uncertainty. The validation approach considered here is based on the premise that the model is valid if the experimental results can be considered as random sample responses of the stochastic system of which the reduced order model is the mean. For the situation considered here, the power spectra of the experiments should lie within the 2nd and 98th percentiles of the response (forming the uncertainty band) of the stochastic model. Nominally clamped–clamped beams are considered to demonstrate the entire process. The construction of two mean reduced order models and their stochastic counterparts are first performed. Then, the validation effort is carried out by comparing experimentally obtained power spectra and their corresponding computational uncertainty bands. This process leads, for both reduced order models, to a very good representation of the important upper envelope (98th percentile) of the experimental data but a less good fit of the lower envelope (2nd percentile).     

Forced response statistics of mistuned bladed disks: a stochastic reduced basis approach

This paper presents a stochastic reduced basis approach for predicting the forced response statistics of mistuned bladed-disk assemblies. In this approach, the system response in the frequency domain is represented using a linear combination of complex stochastic basis vectors with undermined coefficients. The terms of the preconditioned stochastic Krylov subspace are used here as basis vectors. Two variants of the stochastic Bubnov-Galerkin scheme are employed for computing the undetermined terms in the reduced basis representation, which arise from how the condition for orthogonality between two random vectors is interpreted. Explicit expressions for the response quantities can then be derived in terms of the random system parameters, which allow for the possibility of efficiently computing the response statistics in the post-processing stage. Numerical studies are presented for mistuned cyclic assemblies of mono-coupled single-mode components. It is demonstrated that the accuracy of the response statistical moments computed using stochastic reduced basis methods can be orders of magnitude better than classical perturbation methods.     

Systematic derivation of order parameters through reduced density matrices

A systematic method for determining order parameters for quantum many-body systems on lattices is developed by utilizing reduced density matrices. This method allows one to extract the order parameter directly from the wave functions of the degenerate ground states without the aid of empirical knowledge, and thus opens a way to explore unknown exotic orders. The applicability of this method is demonstrated numerically or rigorously in models that are considered to exhibit dimer, scalar chiral, and topological orders.     

Higher order parabolic approximations of the reduced wave equation

Asymptotic solutions of order k⁻ⁿ are developed for the reduced wave equation. Here k is a dimensionless wave number and n is the arbitrary order of the approximation. These approximations are an extension of geometric acoustics theory and provide corrections to that theory in the form of multiplicative functions which satisfy parabolic partial differential equations. These corrections account for the diffraction effects caused by variation of the field normal to the ray path and the interaction of these transverse variations with the variation of the field along the ray. The theory is applied to the example of radiation from a piston, and it is demonstrated that the higher order approximations are more accurate for decreasing values of k.     

Assessment of reduced models for the detection of modal interaction through rotor stator contacts

Interactions through direct contact between blade-tips and outer casings in modern turbomachines require complex formulations and subsequent expensive computational efforts when the classical finite element method is considered. The construction of reduced-order models through component mode synthesis techniques usually improves the computational efficiency and may be used for fast parameter studies yielding a better knowledge of the phenomena of interest.In this highly nonlinear framework, the present study is dedicated to the investigation of the capabilities of fixed- and free-interface reduction strategies to handle accurately such problems through a realistic 2D model and complements former results involving a direct modal projection with respective strong kinematic restrictions.The equations of motion are solved using an explicit time integration scheme together with the Lagrange multiplier method where friction is accounted for. The presented work discusses the notions of both displacement and motion convergences and the possibility to conduct fast parameter studies with the use of relevant reduction bases. It also shows that kinematic restrictions artificially enhance the detection of modal interactions.     

Development of a reduced biodiesel surrogate model for compression ignition engine modeling

Methylbutanoate (MB), a C₄ methyl ester, represents the simplest surrogate that captures the chemical effects of the ester moiety in biodiesel and biodiesel surrogates. An updated chemical kinetic model has been developed to characterize the ignition and flame characteristics of MB. The mechanistic elements within this model that relate to the MB and smaller ester/oxygenate sub-mechanisms are drawn from the prototypical Fisher et al. model and from more recent theory and modeling efforts. The MB model development which is based on an iterative procedure involving global sensitivity analyses to identify elementary reactions that govern ignition and subsequent high level ab initio based theoretical updates to these reaction rates are presented. The MB model makes reasonable predictions of ignition delays and laminar flame speeds.The C₅–C₇ submechanisms from the LLNL n-heptane (NH) model were merged with the present MB model to obtain a detailed chemical kinetics model for a surrogate blend representing biodiesel. The detailed MB-NH model (661 species) was reduced using graph based techniques. The robust reduction techniques employed result in a reduced model (145 species) that is in good agreement with the detailed model over a wide range of conditions. 3-D compression ignition (CI) engine simulations utilizing this reduced chemistry model for MB-NH blends as a surrogate for biodiesel show good agreement with the experimental data suggesting the utility of this model for predictions of combustion and emission characteristics of biodiesel in realistic CI engine simulations.     

Reduced order modeling of head related impulse responses for virtual acoustic displays

This study investigated the use of reduced order head related impulse response (HRIR) models to improve the computational efficiency in acoustic virtual displays. State space models of varying order were generated from zero-elevation HRIRs using a singular value decomposition technique. A source identification experiment was conducted under anechoic conditions in which three subjects were required to localize sounds in the front horizontal plane. The sounds were either (1) real sources (emitted by individual loudspeakers in a semi-circular array), (2) virtual sources generated from the original HRIRs, or (3) virtual sources generated using reduced order state space models. All virtual sources were created by simultaneous activation of two loudspeakers at +/- 30 degrees using a virtual source imaging technique based on either the measured or modeled HRIRs. The errors in the perceived direction of the virtual sources generated from the reduced order models were compared to errors in localization using the original HRIRs. The results demonstrate that a very significant reduction in model size can be achieved without significantly affecting the fidelity of the virtual display of horizontally placed sources.     

On the parameters identification of the Duffing's system by means of a reduced order observer

An on-line procedure for recovering the unknown parameters set of the Duffing's oscillator by means of a reduced order proportional integral observer is presented in this Letter. First, it is shown that the oscillator has the properties of being algebraically observable and algebraically identifiable with respect to a well-chosen output (which turns out to be the oscillator's position). Therefore, an extended differential parametrization of the output and its time derivatives can be obtained. This extended differential parametrization has the necessary information to estimate the output time derivatives and to recover the unknown parameters. The numerical implementation of this method is easily accomplished in a digital computer.     

A high order kinetic flux-vector splitting method for the reduced five-equation model of compressible two-fluid flows

We present a high order kinetic flux-vector splitting (KFVS) scheme for the numerical solution of a conservative interface-capturing five-equation model of compressible two-fluid flows. This model was initially introduced by Wackers and Koren (2004) [21]. The flow equations are the bulk equations, combined with mass and energy equations for one of the two fluids. The latter equation contains a source term in order to account for the energy exchange. We numerically investigate both one- and two-dimensional flow models. The proposed numerical scheme is based on the direct splitting of macroscopic flux functions of the system of equations. In two space dimensions the scheme is derived in a usual dimensionally split manner. The second order accuracy of the scheme is achieved by using MUSCL-type initial reconstruction and Runge–Kutta time stepping method. For validation, the results of our scheme are compared with those from the high resolution central scheme of Nessyahu and Tadmor [14]. The accuracy, efficiency and simplicity of the KFVS scheme demonstrate its potential for modeling two-phase flows.     

电感电流连续模式下Boost变换器的分数阶建模与仿真分析

基于分数阶微积分理论以及实际电容和实际电感在本质上是分数阶的事实,建立电感电流连续模式下Boost变换器的分数阶数学模型以及分数阶状态平均模型并进行理论分析,给出Boost变换器运行于电感电流连续模式下的参数条件.最后,基于Matlab/Simulink软件建立运行于电感电流连续模式下Boost变换器的分数阶仿真模型,并进行数值仿真分析以验证分数阶建模与理论分析的正确性. 关键词： Boost变换器 电感电流连续模式 分数阶微积分 数值仿真     

Finite element reduced order model for noise and vibration reduction of double sandwich panels using shunted piezoelectric patches

This work concerns the control of sound transmission through double laminated panels with viscoelastic core using semi-passive piezoelectric shunt technique. More specifically, the system consists of two laminated walls, each one composed of three layers and called sandwich panel with an air cavity in between. The external sandwich panel has a surface-mounted piezoelectric patches. The piezoelectric elements, connected with resonant shunt circuits, are used for the vibration damping of some specific resonance frequencies of the coupled system. Firstly, a finite element formulation of the fully coupled visco-electro-mechanical-acoustic system is presented. This formulation takes into account the frequency dependence of the viscoelastic material. A modal reduction approach is then proposed to solve the problem at a lower cost. In the proposed technique, the coupled system is solved by projecting the mechanical displacement unknown on a truncated basis composed by the first real short-circuit structural normal modes and the pressure unknown on a truncated basis composed by the first acoustic modes with rigid boundaries conditions. The few initial electrical unknowns are kept in the reduced system. A static correction is also introduced in order to take into account the effect of higher modes. Various results are presented in order to validate and illustrate the efficiency of the proposed finite element reduced order formulation.     

Conversion of high-order Laguerre-Gaussian beams into Bessel beams of increased, reduced or zeroth order by use of a helical axicon

We propose a new method for transformation of a Laguerre-Gaussian beam of azimuthal index l and radial index n = 0 (LG_l,0) into a vortex, diverging or nondiverging Bessel beam, which can have increased or decreased phase singularity order, or into a zeroth order Bessel beam, by means of a helical axicon. The Bessel beam divergence or nondivergence depends upon the waist position of the input Laguerre-Gaussian beam, regarding the plane where the helical axicon is situated.The expressions for the amplitude and the intensity distribution of the diffracted wave field, in the process of Fresnel diffraction, are deduced using the stationary phase method. The theoretical analysis for the vortex radius and the maximum propagation distance of the Bessel beams obtained is presented.