基于CFD技术的鼻腔输药的流动特性研究

鼻腔输药作为一种新颖的输药方式,具有方便可靠的特点,不仅可以适用于局部而且适用于全身用药,已经受到广泛的重视。但由于人类的呼吸道复杂的立体结构使得鼻腔输药成为生物医药研究的一个长期的挑战。利用CFD(计算流体动力学,Computational Fluid Dynamics)技术对鼻腔输药的流动特性进行研究,对一个基于真实的人类鼻腔解剖的3D模型的输药的流动过程进行了分析,着重研究了不同的药物颗粒的流动和沉淀特性。利用CFD的商用软件FLUENT在层流情况下对不同的密度、直径和流动速率的药物颗粒在鼻腔中流动进行了数值模拟,并对计算结果进行了比较分析。CFD研究的结果为鼻腔输药这种新颖的药物输送方式的优化设计提供了依据。     

Probing the era of galaxy formation via TeV gamma ray absorption by the near infrared extragalactic background

We present models of the extragalactic background light (EBL) based on several scenarios of galaxy formation and evolution. We have treated galaxy formation with the Press-Schecter approximation for both cold dark matter (CDM) and cold+hot dark matter (CHDM) models, representing a moderate (z _f 3) and a late (z _f 1) era of galaxy formation respectively. Galaxy evolution has been treated by considering a variety of stellar types, different initial mass functions and star formation histories, and with an accounting of dust absorption and emission. We find that the dominant factor influencing the EBL is the epoch of galaxy formation. A recently proposed method for observing the EBL utilizing the absorption of 0.1 to 10 TeV gamma-rays from active galactic nuclei (AGN) is shown to be capable of discriminating between different galaxy formation epochs. The one AGN viewed in TeV light, Mrk 421, does show some evidence for a cutoff above 3 TeV; based on the EBL models presented here, we suggest that this is due to extinction in the source. The large absorption predicted at energies > 200 GeV for sources at z > 0.5 indicates that observations of TeV gamma-ray bursts (GRB) would constrain or eliminate models in which the GRB sources lie at cosmological distances.Now at University of Chicago, Dept. of Astronomy & Astrophysics.     

Discrete-time observability and estimability analysis for bearings-only target motion analysis

Observability in the context of bearings-only tracking (BOT) is still the subject of important literature. Different from previous approaches, where continuous-time analysis was considered, our approach relies on discrete-time analysis. It is then shown that this allows us to use directly and efficiently the simple formalisms of linear algebra. Using the direct approach, observability analysis is essentially reduced to basic considerations about subspace dimensions. Even if this approach is conceptually quite direct, it becomes more and more complex as the source-encounter scenario complexity increases. For complex scenarios, the dual approach may present some advantages essentially due to the direct use of multilinear algebra. New results about BOT observability for maneuvering sources are thus obtained. Observability analysis is then extended to unknown instants of source velocity changes. Even if observability analysis provides thorough insights about the algebraic structure of the BOT problem, the optimization of the observer maneuvers is essentially a control problem. Basic algebraic considerations prove that a relevant cost functional for this control problem is the determinant of the Fisher information matrix (FIM). So, a large part of this work is devoted to the analysis of this cost functional. Using multilinear algebra, general approximations of this functional are given. In order to involve only directly estimable parameters, the source bearing-rates are examined. Using these approximations, a general framework for optimizing the observer trajectory is derived which allow us to approximate the optimal sequence of controls. It is worth stressing that our approach does not require the knowledge of the source trajectory parameters and is still valid for a maneuvering source.     

A New Mathematical Formulation for Strapdown Inertial Navigation

A differential equation is developed for the orientation vector relating the body frame to a chosen reference frame. The time derivative of this vector is the sum of the inertially measurable angular velocity vector and of the inertially nonmeasurable noncommutativity rate vector. It is precisely this noncommutativity rate vector that causes the computational problems when numerically integrating the direction cosine matrix. The orientation vector formulation allows the noncommutativity contribution to be isolated and, therefore, treated separately and advantageously. An orientation vector mechanization is presented for a strap down inertial system. Further, an example is given of the applica tion of this formulation to a typical rigid body rotation problem.     

Medium PRF set selection using evolutionary algorithms

This paper presents a new and novel method of selecting multiple pulse repetition frequency (PRF) sets for use in medium PRF pulsed-Doppler radars. Evolutionary algorithms are used to minimise the blind areas in the range/Doppler space. The evolutionary algorithm allows optimal solutions to be generated quickly, far faster than with exhaustive searches, and is fully automatic, unlike existing techniques. The evolved solutions compare very favorably against the results of both an exhaustive search and existing published PRF set selection methods. This evolutionary approach to generation of PRF sets is a major advance in medium PRF radar design.     

Matched Filter Criteria and Range Weighting for Weather Radar

The effective weighting function for weather radar is defined. This weighting function considers the effects of both the transmitted signal and the receiver filter. It is used to assign effective ranges to samples taken at prescribed times. For uniformly distributed targets it is shown that "signal"-to-noise ratio depends on the receiver filter, transmitted signal envelope, and receiver noise spectral density. Maximization of this signal-to-noise ratio when range resolution constraints are imposed is discussed, and a receiver design approach specifically adapted for Doppler weather radars is developed.     

A Correlation Polarimeter for Noise-Like Signals

Optimum estimation (tracking) of the polarization plane of a linearly polarized electromagnetic wave is determined when the signal is a narrow-band Gaussian random process with a polarization plane angle which is also a Gaussian random process. This model is Compared to previous work and is applicable to space communication. The estimator performs a correlation operation similar to an amplitude -comparison monopulse angle tracker, giving the name correlation polarimeter. Under large signal-to-noise ratio (SNR), the estimator is causal. Performance of the causal correlation polarimeter is evaluated for arbitrary SNR. Optimum precorrelation filtering is determined. With low SNR, the performance of this system is far better than that of previously developed systems. Practical implementation is discussed. A scheme is given to reduce the effect of linearly polarized noise.     

Magnetosphere Imaging Instrument (MIMI) on the Cassini Mission to Saturn/Titan

The magnetospheric imaging instrument (MIMI) is a neutral and charged particle detection system on the Cassini orbiter spacecraft designed to perform both global imaging and in-situ measurements to study the overall configuration and dynamics of Saturn’s magnetosphere and its interactions with the solar wind, Saturn’s atmosphere, Titan, and the icy satellites. The processes responsible for Saturn’s aurora will be investigated; a search will be performed for substorms at Saturn; and the origins of magnetospheric hot plasmas will be determined. Further, the Jovian magnetosphere and Io torus will be imaged during Jupiter flyby. The investigative approach is twofold. (1) Perform remote sensing of the magnetospheric energetic (E > 7 keV) ion plasmas by detecting and imaging charge-exchange neutrals, created when magnetospheric ions capture electrons from ambient neutral gas. Such escaping neutrals were detected by the Voyager l spacecraft outside Saturn’s magnetosphere and can be used like photons to form images of the emitting regions, as has been demonstrated at Earth. (2) Determine through in-situ measurements the 3-D particle distribution functions including ion composition and charge states (E > 3 keV/e). The combination of in-situ measurements with global images, together with analysis and interpretation techniques that include direct “forward modeling’’ and deconvolution by tomography, is expected to yield a global assessment of magnetospheric structure and dynamics, including (a) magnetospheric ring currents and hot plasma populations, (b) magnetic field distortions, (c) electric field configuration, (d) particle injection boundaries associated with magnetic storms and substorms, and (e) the connection of the magnetosphere to ionospheric altitudes. Titan and its torus will stand out in energetic neutral images throughout the Cassini orbit, and thus serve as a continuous remote probe of ion flux variations near 20R _S (e.g., magnetopause crossings and substorm plasma injections). The Titan exosphere and its cometary interaction with magnetospheric plasmas will be imaged in detail on each flyby. The three principal sensors of MIMI consists of an ion and neutral camera (INCA), a charge–energy–mass-spectrometer (CHEMS) essentially identical to our instrument flown on the ISTP/Geotail spacecraft, and the low energy magnetospheric measurements system (LEMMS), an advanced design of one of our sensors flown on the Galileo spacecraft. The INCA head is a large geometry factor (G ∼ 2.4 cm² sr) foil time-of-flight (TOF) camera that separately registers the incident direction of either energetic neutral atoms (ENA) or ion species (≥5^∘ full width half maximum) over the range 7 keV/nuc < E < 3 MeV/nuc. CHEMS uses electrostatic deflection, TOF, and energy measurement to determine ion energy, charge state, mass, and 3-D anisotropy in the range 3 ≤ E ≤ 220 keV/e with good (∼0.05 cm² sr) sensitivity. LEMMS is a two-ended telescope that measures ions in the range 0.03 ≤ E ≤ 18 MeV and electrons 0.015 ≤ E≤ 0.884 MeV in the forward direction (G ∼ 0.02 cm² sr), while high energy electrons (0.1–5 MeV) and ions (1.6–160 MeV) are measured from the back direction (G ∼ 0.4 cm² sr). The latter are relevant to inner magnetosphere studies of diffusion processes and satellite microsignatures as well as cosmic ray albedo neutron decay (CRAND). Our analyses of Voyager energetic neutral particle and Lyman-α measurements show that INCA will provide statistically significant global magnetospheric images from a distance of ∼60 R _S every 2–3 h (every ∼10 min from ∼20 R _S). Moreover, during Titan flybys, INCA will provide images of the interaction of the Titan exosphere with the Saturn magnetosphere every 1.5 min. Time resolution for charged particle measurements can be < 0.1 s, which is more than adequate for microsignature studies. Data obtained during Venus-2 flyby and Earth swingby in June and August 1999, respectively, and Jupiter flyby in December 2000 to January 2001 show that the instrument is performing well, has made important and heretofore unobtainable measurements in interplanetary space at Jupiter, and will likely obtain high-quality data throughout each orbit of the Cassini mission at Saturn. Sample data from each of the three sensors during the August 18 Earth swingby are shown, including the first ENA image of part of the ring current obtained by an instrument specifically designed for this purpose. Similarily, measurements in cis-Jovian space include the first detailed charge state determination of Iogenic ions and several ENA images of that planet’s magnetosphere.This revised version was published online in July 2005 with a corrected cover date.     

A Possible Transfer Mechanism for the 11-Year Solar Cycle to the Lower Stratosphere

Observational evidence of the 11-year solar cycle (SC) modulation of stratosphere temperatures and winds from the ERA-40 dataset is reviewed, with emphasis on the Northern winter hemisphere. A frequency modulation of sudden warming events is noted, with warmings occurring earlier in solar minimum periods than in solar maximum periods. The observed interaction between the influence of the SC and the quasi biennial oscillation (QBO) on the frequency of sudden warmings is noted as a possible clue for understanding their mechanism of influence. A possible transfer route for the 11-year solar cycle from the equatorial stratopause region to the lowest part of the stratosphere is proposed, via an influence on sudden warming events and the associated induced meridional circulation. SC and QBO composites of zonal wind anomalies show anomalous wind distributions in the subtropical upper stratosphere in early winter. Mechanistic model experiments are reviewed that demonstrate a sensitivity of sudden warmings to small wind anomalies in this region. Various diagnostics from these experiments are shown, including EP fluxes and their divergence and also the synoptic evolution of the polar vortex, in order to understand the mechanism of the influence. Some recent GCM experiments to investigate the SC/QBO interaction are also described. They simulate reasonably well the observed SC/QBO interaction of sudden warming events and appear to support the hypothesis that tropical/subtropical upper stratospheric wind anomalies are an important influence on the timing of sudden warmings.     

Deep Impact: Working Properties for the Target Nucleus – Comet 9P/Tempel 1

In 1998, Comet 9P/Tempel 1 was chosen as the target of the Deep Impact mission (A’Hearn, M. F., Belton, M. J. S., and Delamere, A., Space Sci. Rev., 2005) even though very little was known about its physical properties. Efforts were immediately begun to improve this situation by the Deep Impact Science Team leading to the founding of a worldwide observing campaign (Meech et al., Space Sci. Rev., 2005a). This campaign has already produced a great deal of information on the global properties of the comet’s nucleus (summarized in Table I) that is vital to the planning and the assessment of the chances of success at the impact and encounter. Since the mission was begun the successful encounters of the Deep Space 1 spacecraft at Comet 19P/Borrelly and the Stardust spacecraft at Comet 81P/Wild 2 have occurred yielding new information on the state of the nuclei of these two comets. This information, together with earlier results on the nucleus of comet 1P/Halley from the European Space Agency’s Giotto, the Soviet Vega mission, and various ground-based observational and theoretical studies, is used as a basis for conjectures on the morphological, geological, mechanical, and compositional properties of the surface and subsurface that Deep Impact may find at 9P/Tempel 1. We adopt the following working values (circa December 2004) for the nucleus parameters of prime importance to Deep Impact as follows: mean effective radius = 3.25± 0.2 km, shape – irregular triaxial ellipsoid with a/b = 3.2± 0.4 and overall dimensions of ∼14.4 × 4.4 × 4.4 km, principal axis rotation with period = 41.85± 0.1 hr, pole directions (RA, Dec, J2000) = 46± 10, 73± 10 deg (Pole 1) or 287± 14, 16.5± 10 deg (Pole 2) (the two poles are photometrically, but not geometrically, equivalent), Kron-Cousins (V-R) color = 0.56± 0.02, V-band geometric albedo = 0.04± 0.01, R-band geometric albedo = 0.05± 0.01, R-band H(1,1,0) = 14.441± 0.067, and mass ∼7×10¹³ kg assuming a bulk density of 500 kg m⁻³. As these are working values, {i.e.}, based on preliminary analyses, it is expected that adjustments to their values may be made before encounter as improved estimates become available through further analysis of the large database being made available by the Deep Impact observing campaign. Given the parameters listed above the impact will occur in an environment where the local gravity is estimated at 0.027–0.04 cm s⁻² and the escape velocity between 1.4 and 2 m s⁻¹. For both of the rotation poles found here, the Deep Impact spacecraft on approach to encounter will find the rotation axis close to the plane of the sky (aspect angles 82.2 and 69.7 deg. for pole 1 and 2, respectively). However, until the rotation period estimate is substantially improved, it will remain uncertain whether the impactor will collide with the broadside or the ends of the nucleus.