A computational fluid dynamics （CFD） approach is used to study the respiratory airflow dynamics within a human upper airway. The airway model which consists of the airway from nasal cavity, pharynx, larynx and trach...A computational fluid dynamics （CFD） approach is used to study the respiratory airflow dynamics within a human upper airway. The airway model which consists of the airway from nasal cavity, pharynx, larynx and trachea to triple bifurcation is built based on the CT images of a healthy volunteer and the Weibel model. The flow character- istics of the whole upper airway are quantitatively described at any time level of respiratory cycle. Simulation results of respiratory flow show good agreement with the clinical mea- sures, experimental and computational results in the litera- ture. The air mainly passes through the floor of the nasal cavity in the common, middle and inferior nasal meatus. The higher airway resistance and wall shear stresses are distrib- uted on the posterior nasal valve. Although the airways of pharynx, larynx and bronchi experience low shear stresses, it is notable that relatively high shear stresses are distrib- uted on the wall of epiglottis and bronchial bifurcations. Besides, two-dimensional fluid-structure interaction models of normal and abnormal airways are built to discuss the flow-induced deformation in various anatomy models. The result shows that the wall deformation in normal airway is relatively small.展开更多
The large finite element global stiffness matrix is an algebraic, discreet, even-order, differential operator of zero row sums. Direct application of the, practically convenient, readily applied, Gershgorin’s eigenva...The large finite element global stiffness matrix is an algebraic, discreet, even-order, differential operator of zero row sums. Direct application of the, practically convenient, readily applied, Gershgorin’s eigenvalue bounding theorem to this matrix inherently fails to foresee its positive definiteness, predictably, and routinely failing to produce a nontrivial lower bound on the least eigenvalue of this, theoretically assured to be positive definite, matrix. Considered here are practical methods for producing an optimal similarity transformation for the finite-elements global stiffness matrix, following which non trivial, realistic, lower bounds on the least eigenvalue can be located, then further improved. The technique is restricted here to the common case of a global stiffness matrix having only non-positive off-diagonal entries. For such a matrix application of the Gershgorin bounding method may be carried out by a mere matrix vector multiplication.展开更多
基金supported by the National Natural Science Foundation of China (10472025 10672036+1 种基金 10872043)Natural Science Foundation of Liaoning Province, China (20032109).
文摘A computational fluid dynamics （CFD） approach is used to study the respiratory airflow dynamics within a human upper airway. The airway model which consists of the airway from nasal cavity, pharynx, larynx and trachea to triple bifurcation is built based on the CT images of a healthy volunteer and the Weibel model. The flow character- istics of the whole upper airway are quantitatively described at any time level of respiratory cycle. Simulation results of respiratory flow show good agreement with the clinical mea- sures, experimental and computational results in the litera- ture. The air mainly passes through the floor of the nasal cavity in the common, middle and inferior nasal meatus. The higher airway resistance and wall shear stresses are distrib- uted on the posterior nasal valve. Although the airways of pharynx, larynx and bronchi experience low shear stresses, it is notable that relatively high shear stresses are distrib- uted on the wall of epiglottis and bronchial bifurcations. Besides, two-dimensional fluid-structure interaction models of normal and abnormal airways are built to discuss the flow-induced deformation in various anatomy models. The result shows that the wall deformation in normal airway is relatively small.
文摘The large finite element global stiffness matrix is an algebraic, discreet, even-order, differential operator of zero row sums. Direct application of the, practically convenient, readily applied, Gershgorin’s eigenvalue bounding theorem to this matrix inherently fails to foresee its positive definiteness, predictably, and routinely failing to produce a nontrivial lower bound on the least eigenvalue of this, theoretically assured to be positive definite, matrix. Considered here are practical methods for producing an optimal similarity transformation for the finite-elements global stiffness matrix, following which non trivial, realistic, lower bounds on the least eigenvalue can be located, then further improved. The technique is restricted here to the common case of a global stiffness matrix having only non-positive off-diagonal entries. For such a matrix application of the Gershgorin bounding method may be carried out by a mere matrix vector multiplication.
