Direct numerical simulation of autoignition in a non-premixed, turbulent medium

In this paper, direct numerical simulation of autoignition in an initially non-premixed medium under isotropic, homogeneous, and decaying turbulence is presented. The pressure-based method developed herein is a spectral implementation of the sequential steps followed in the predictor-corrector type of algorithms; it includes the effects of density fluctuations caused by spatial inhomogeneities ill temperature and species. The velocity and pressure field are solved in the spectral space while the scalars and density field are solved in the physical space. The presented results reveal that the autoignition spots originate and evolve at locations where (1) the composition corresponds to a small range around a specific mixture fraction, and (2) the conditional scaler dissipation rate is low. A careful examination of the data obtained indicates that the autoignition spots originate in the vortex cores, and the hot gases travel outward as combustion progresses. Hence, the applicability of the transient laminar flamelet model for this problem is questioned. The dependence of autoignition characteristics on parameters such as (1) die initial eddy-turnover time and (2) the initial ratio of length scale of scalars to that of velocities are investigated. Certain implications of new results on the conditional moment closure modeling are discussed.

A new parallel overlapped domain decomposition method for nonlinear dynamic finite element analysis

In this paper a new parallel algorithm for nonlinear transient dynamic analysis of large structures has been presented. An unconditionally stable Newmark-beta method (constant average acceleration technique) has been employed for time integration. The proposed parallel algorithm has been devised within the broad framework of domain decomposition techniques. However, unlike most of the existing parallel algorithms (devised for structural dynamic applications) which are basically derived using nonoverlapped domains, the proposed algorithm uses overlapped domains. The parallel overlapped domain decomposition algorithm proposed in this paper has been formulated by splitting the mass, damping and stiffness matrices arises out of finite element discretisation of a given structure. A predictor-corrector scheme has been formulated for iteratively improving the solution in each step. A computer program based on the proposed algorithm has been developed and implemented with message passing interface as software development environment. PARAM-10000 MIMD parallel computer has been used to evaluate the performances. Numerical experiments have been conducted to validate as well as to evaluate the performance of the proposed parallel algorithm. Comparisons have been made with the conventional nonoverlapped domain decomposition algorithms. Numerical studies indicate that the proposed algorithm is superior in performance to the conventional domain decomposition algorithms. (C) 2003 Elsevier Ltd. All rights reserved.

任意阶隐式指数时程差分多步法及其在非线性系统中的应用

对非线性系统提出了任意阶隐式指数时程差分多步法 ,实现了任意阶次指数时程差分预测 校正算法 .发展完善了指数时程差分法 .将新算法应用于非线性系统 ,取得了较好的效果 .数值结果表明隐式指数时程差分多步法很好地修正了显式指数时程差分多步法 ,隐式指数时程差分多步法是一种高精度、高效率的方法

Propagation of an arbitrary elliptically polarized few-cycle ultrashort laser pulse in resonant two-level quantum systems

We investigate the propagation of an arbitrary elliptically polarized few-cycle ultrashort laser pulse in resonant two-level quantum systems using an iterative predictor-corrector finite-difference time-domain method. It is shown that when the initial effective area is equal to 2 pi, the effective area will remain invariant during the course of propagation, and a complete Rabi oscillation can be achieved. However, for an elliptically polarized few-cycle ultrashort laser pulse, polarization conversion can occur. Eventually, the laser pulse will evolve into two separate circularly polarized laser pulses with opposite helicities.

Polarization evolution of a few-cycle ultrashort laser pulse propagating in a degenerate three-level medium

The propagation of an arbitrary polarized few-cycle ultrashort laser pulse in a degenerate three-level medium is investigated by using an iterative predictor-corrector finite-difference time-domain method. It is found that the polarization evolution of the ultrashort laser pulse is dependent not only on the initial atomic coherence of the medium but also on the polarization condition of the incident laser pulse. When the initial effective area is equal to 2 pi, complete linear-to-circular and circular-to-linear polarization conversion of few-cycle ultrashort laser pulses can be achieved due to the quantum interference effects between the two different transition paths.

Lau cavity and phase locking of laser arrays

The Lau cavity is the self-imaging cavity with a phase corrector under the Lau reimaging condition. The author proposes the use of the Lau cavity to utilize both the Talbot and the Lau effects for phase locking one-dimensional and two-dimensional diode-laser arrays into a single-lobe coherent beam. Analyses on the self-reproducing of a coherent lasing field and the reimaging of initial incoherent radiation are given.

Uso do espaço por Hypselartemon contusulus (Férussac, 1827) (Gastropoda, Streptaxidae) da Trilha da Parnaioca, Ilha Grande, Rio de Janeiro, Brasil

O comportamento espacial dos indivíduos é um componente chave para se entender a dinâmica de população dos organismos e esclarecer o potencial de migração e dispersão das espécies. Vários fatores afetam a atividade de locomoção de moluscos terrestres, como temperatura, luz, umidade, época do ano, tamanho da concha, sexo, estratégia reprodutiva, idade, densidade de coespecíficos e disponibilidade de alimento. Um dos métodos usados para estudar deslocamento de gastrópodes terrestres é o de marcação-recaptura. Gastrópodes terrestres se prestam a este tipo de estudo por causa de (1) seu reduzido tamanho, (2) fácil manejo, (3) fácil captura e (4) pequenas distâncias de deslocamento e, consequentemente, reduzidas áreas de vida. Estes organismos servem como modelo para o estudo de ecologia espacial e dispersão. Estudos de população, investigando o uso do espaço, a distribuição espacial, a densidade populacional e a área de vida são escassos para moluscos terrestres e ainda mais raros em áreas naturais tropicais. Nosso objeto de estudo é Hypselartemon contusulus (Férussac, 1827), um molusco terrestre carnívoro, da família Streptaxidae, muito abundante na serrapilheira, em trechos planos de mata secundária na Trilha da Parnaioca, Ilha Grande, Rio de Janeiro. A espécie é endêmica para o estado do Rio de Janeiro. Seu tamanho é de até 7,2 mm de altura, apresentando 6 a 7 voltas. Neste trabalho estudamos as variáveis temperatura ambiente, temperatura do solo, umidade do ar, luminosidade, profundidade do folhiço, tamanho do animal, densidade de co-específicos e densidade de presas, relacionando estes dados ecológicos ao deslocamento observado em Hypselartemon contusulus. Uma das hipóteses de trabalho é que estas variáveis afetam seu deslocamento. O trabalho foi realizado na Ilha Grande, situada ao sul do Estado do Rio de Janeiro, no município de Angra dos Reis. Os animais foram capturados e marcados com um código individual pintado na concha com corretor ortográfico líquido e caneta nanquim. As distâncias de deslocamento, em cm, foram registradas medindo-se as distâncias entre marcadores subsequentes. Os resultados encontrados indicam que o método utilizado é eficaz para marcar individualmente Hypselartemon contusulus em estudos de médio prazo (até nove meses). Sugerimos o uso deste método de marcação para estudos com gastrópodes terrestres ameaçados de extinção, como algumas espécies das famílias Bulimulidae, Megalobulimidae, Streptaxidae e Strophocheilidae. Hypselartemon contusulus não mantém uma distância mínima de seus vizinhos, é ativo ao longo de todo o ano e ao longo do dia, demonstrando atividade de locomoção e predação. Não foram encontrados animais abrigados sob pedra ou madeira morta. Não foram observados locais de atividade em oposição a lugares de repouso/abrigo. Beckianum beckianum (Pfeiffer, 1846) foi a presa preferencial. A densidade populacional variou de 0,57 a 1,2 indivíduos/m2 entre as campanhas de coleta. A espécie desloca-se, em média, 26,57 17,07 cm/24h, na Trilha da Parnaioca, Ilha Grande. A área de vida de H. contusulus é pequena, sendo de, no máximo, 0,48 m2 em três dias e 3,64 m2 em 79 dias. O deslocamento da espécie variou ao longo do ano, mas esta variação não é afetada pelas variáveis ecológicas estudadas. Este é, portanto, um comportamento plástico em H. contusulus e, provavelmente, controlado por fatores endógenos.

Two-orbit switch-pitch structures

The Hoberman 'switch-pitch ' ball is a transformable structure with a single folding and unfolding path. The underlying cubic structure has a novel mechanism that retains tetrahedral symmetry during folding. Here, we propose a generalized class of structures of a similar type that retain their full symmetry during folding. The key idea is that we require two orbits of nodes for the structure: within each orbit, any node can be copied to any other node by a symmetry operation. Each member is connected to two nodes, which may be in different orbits, by revolute joints. We will describe the symmetry analysis that reveals the symmetry of the internal mechanism modes for a switch-pitch structure. To follow the complete folding path of the structure, a nonlinear iterative predictor-corrector algorithm based on the Newton method is adopted. First, a simple tetrahedral example of the class of two-orbit structures is presented. Typical configurations along the folding path are shown. Larger members of the class of structures are also presented, all with cubic symmetry. These switch-pitch structures could have useful applications as deployable structures.

Low-rank optimization with trace norm penalty

The paper addresses the problem of low-rank trace norm minimization. We propose an algorithm that alternates between fixed-rank optimization and rank-one updates. The fixed-rank optimization is characterized by an efficient factorization that makes the trace norm differentiable in the search space and the computation of duality gap numerically tractable. The search space is nonlinear but is equipped with a Riemannian structure that leads to efficient computations. We present a second-order trust-region algorithm with a guaranteed quadratic rate of convergence. Overall, the proposed optimization scheme converges superlinearly to the global solution while maintaining complexity that is linear in the number of rows and columns of the matrix. To compute a set of solutions efficiently for a grid of regularization parameters we propose a predictor-corrector approach that outperforms the naive warm-restart approach on the fixed-rank quotient manifold. The performance of the proposed algorithm is illustrated on problems of low-rank matrix completion and multivariate linear regression. © 2013 Society for Industrial and Applied Mathematics.

Hybrid numerical methods for exponential models of growth

A series of novel numerical methods for the exponential models of growth are proposed. Based on these methods, hybrid predictor-corrector methods are constructed. The hybrid numerical methods can increase the accuracy and the computing speed obviously, as well as enlarge the stability domain greatly. (c) 2005 Published by Elsevier Inc.

水下机器人神经网络自适应逆控制

水下环境的复杂性以及自身模型的不确定性,给水下机器人的控制带来很大困难。针对水下机器人的特点和控制方面所存在的问题,提出了基于预测 校正控制策略的水下机器人神经网络自适应逆控制结构及训练算法。通过在线辨识系统的前向模型,估计出系统的Jacobian矩阵,然后采用预报误差法实现控制器的自适应。同时,为了提高系统对于外扰的鲁棒性,在伪线性回归算法的基础上,在评价函数中引入微分项。理论分析和仿真结果表明,与原来的算法相比,微分项的引入改善了系统对于外扰的鲁棒性和动态性能。

Dynamic fluid-structure interaction using finite volume unstructured mesh procedures

A three dimensional finite volume, unstructured mesh method for dynamic fluid-structure interation is described. The broad approach is conventional in that the fluid and structure are solved sequentially. The pressure and viscous stresses from the flow algorithm provide load conditions for the solid algorithm, whilst at the fluid structure interface the deformed structure provides boundary condition from the structure to the fluid. The structure algorithm also provides the necessary mesh adaptation for the flow field, the effect of which is accounted for in the flow algorithm. The procedures described in this work have several novel features, namely: * a single mesh covering the entire domain. * a Navier Stokes flow. * a single FV-UM discretisation approach for both the flow and solid mechanics procedures. * an implicit predictor-corrector version of the Newmark algorithm. * a single code embedding the whole strategy. The procedure is illustrated for a three dimensional loaded cantilever in fluid flow.

Electron-impact excitation of CrII: A theoretical calculation of effective collision strengths for optically allowed transitions.

In this paper, we present electron-impact excitation collision strengths and Maxwellian averaged effective collision strengths for the complicated iron-peak ion Cr II. We consider specifically the allowed lines for transitions from the 3d(5) and 3d(4)4s even parity configuration states to the 3d(4)4p odd parity configuration levels. The parallel suite of R-Matrix packages, RMATRX II, which have recently been extended to allow for the inclusion of relativistic effects, were used to compute the collision cross sections. A total of 108 LS pi/280 J pi levels from the basis configurations 3d(5), 3d(4)4s, and 3d(4)4p were included in the wavefunction representation of the target including all doublet, quartet, and sextet terms. Configuration interaction and correlation effects were carefully considered by the inclusion of seven more configurations and a pseudo-corrector (4d) over bar type orbital. The 10 configurations incorporated into the Cr II model thus listed are 3d(5), 3d(4)4s, 3d(4)4p, 3d(3)4s(2), 3d(3)4p(2), 3d(3)4s4p, 3d(4)(4d) over bar, 3d(3)4s (4d) over bar, 3d(3)4p (4d) over bar, and 3d(3)(4d) over bar (2), constituting the largest Cr II target model considered to date in a scattering calculation. The Maxwellian averaged effective collision strengths are computed for a wide range of electron temperatures 2000-100,000 K which are astrophysically significant. Care has been taken to ensure that the partial wave contributions to the collision strengths for these allowed lines have converged with "top-up" from the Burgess-Tully sum rule incorporated. Comparisons are made with the results of Bautista et al. and significant differences are found for some of the optically allowed lines considered.

Use of arc-length method for capturing mode jumping in postbuckling aerostructures

The arc-length method has become a widely established solution technique for studying nonlinear structural behavior. By augmenting the set of nonlinear equilibrium equations with a constraint equation, which is a function of both the displacements and load increment, it is capable of traversing limit points. Numerous investigations have shown that highly nonlinear behavior such as sharp "snap-backs" can still lead to numerical difficulties. Two practical examples are presented to assess the effectiveness of this solution technique in capturing secondary instabilities in postbuckling structures, which present themselves as abrupt mode jumps. Although the first example poses no special difficulties, in the second case the nonlinear procedure fails to converge. An improvement to the method's formulation is suggested, which accounts for the residual forces that are usually neglected, when proceeding to the next increment once convergence is reached on the current increment. The choice of a correct load increment at the first iteration, within a predictor-corrector scheme, is central to the method's effectiveness. Current strategies for a choice of this load increment are discussed and are shown to be no longer consistent with the modified formulation; therefore, a new approach is proposed.

On various low-hardware-complexity LMS algorithms for adaptive I/Q correction in quadrature receivers

In this paper, the performance and convergence time comparisons of various low-complexity LMS algorithms used for the coefficient update of adaptive I/Q corrector for quadrature receivers are presented. We choose the optimum LMS algorithm suitable for low complexity, high performance and high order QAM and PSK constellations. What is more, influence of the finite bit precision on VLSI implementation of such algorithms is explored through extensive simulations and optimum wordlengths established.